五种离散格式Word文档格式.docx
《五种离散格式Word文档格式.docx》由会员分享,可在线阅读,更多相关《五种离散格式Word文档格式.docx(21页珍藏版)》请在冰豆网上搜索。
3.1区域离散 3
3.2方程离散 3
3.2.1中心差分格式 3
3.2.1迎风差分格式 4
3.2.3混合格式 4
3.2.4指数格式 5
3.2.5指数格式 6
3.2.6五种格式格式系数aEDe的表达式 6
4数值方法及程序流程 7
5计算结果验证及网格独立性考核 8
5.1计算结果验证 8
5.1网格独立性考核 11
6结果分析与讨论 14
7参考资料 15
附录 16
附录A:
计算环境及源程序 16
1计算题目
对于有源项的一维稳态空气对流-扩散传热方程:
ddxρuT=ddxλcdTdx+S
设源项S=0.5-100x,利用中心差分格式,一阶迎风格式,混合格式,指数格式,乘方格式求解在不同流速情况下,温度T的一维分布。
2数学物理模型
物理模型:
一维,稳态,有内热源,常物性
控制方程:
ddxρuT=ddxλcdTdx+Sλc=Γ
物理条件:
ρ,Γ=consts=0.5-100x
边界条件:
x=0T=300K;
x=1T=500K
由于在常物性均分网格的情况下网格贝克勒数与流速成正比,所以在流速不同的情况下求温度场即在网格贝克勒数不同的情况下求温度场。
3计算区域及方程离散
3.1区域离散
3.2方程离散
3.2.1中心差分格式
由控制方程:
ddxρuT=ddxΓdTdx+S
对一维模型方程在P控制容积内做积分,取分段线性型线,经整理可得:
(ρuT)e-ρuTw=ΓdTdxe-ΓdTdxw+SPΔxP
(ρu)eTE+TP2-ρuwTP+TW2=ΓeδxeTE-TP-ΓwδxwTP-TW+SPΔxP
整理可得:
aPTP=aETE+aWTW+SPΔxP
其中:
aE=De-12Fe=De(1-12P∆e)De=ΓeδxeFe=(ρu)e
aW=Dw+12Fw=Dw(1+12P∆w)Dw=ΓwδxwFw=(ρu)w
aP=aE+aW+Fe-Fw
为保证代数方程迭代求解的收敛性,我们要求计算中质量守恒一定要满足,于是:
Fe-Fw=0aP=aE+aW
3.2.1迎风差分格式
控制容积法定义:
界面上未知函数永远取上游节点之值
紧凑形式:
(ρuT)e=TP[Fe,0]-TE[-Fe,0]
(ρuT)w=TW[Fw,0]-TP[-Fw,0]
带入控制方程,化简得:
aE=De+-Fe,0=De(1+-P∆e,0)
aW=Dw+Fw,0=Dw(1+P∆w,0)
aP=aE+aW+Fe,0+-Fw,0--Fe,0-Fw,0=aE+aW+Fe-Fw
3.2.3混合格式
紧凑定义:
aEDe=[-P∆e,1-12P∆e,0]
aWDw=[0,1+12P∆w,P∆w]
3.2.4指数格式
含义:
根据一维模型的精确解相对应的离散方程构成的一种格式。
导出方法:
将精确解表示成三点变量间的代数方程。
T-T0TL-T0=expρuxΓ-1expρuLΓ-1=expPexL-1expPe-1
对流扩散总通量密度:
定义:
J=ρuT-ΓdTdxJ=const对控制容积P:
Je=Jw
将对流扩散方程的精确解带入J的定义式,化简得:
J=F(T0+T0-TLexpPe-1)
对于e,w界面写出总通量密度的解析表达式:
Je=FeT0+TP-TEexpPΔe-1Jw=Fw(T0+TW-TPexpPΔw-1)
带入Je=Jw,整理得:
aE=FeexpPΔe-1aW=FWexpPΔwexpPΔw-1
aP=aE+aW
3.2.5指数格式
3.2.6五种格式格式系数aEDe的表达式
4数值方法及程序流程
初始化变量:
x=0,T=300K;
x=1,T=500K
网格划分:
外界点法,取Δx=0.05
确定方程系数:
根据不同情况设网格贝克勒数,求出ae/D,上一点的ae与下一点aw间相差一个网格贝克勒数
求解方程:
高斯赛德尔迭代
5计算结果验证及网格独立性考核
5.1计算结果验证
图1pe=0.7时温度分布图
xCDFUDHBSEXPPL
0.00300.00300.00300.00300.00300.00
0.05348.14299.63299.82299.80299.80
0.10364.97299.37299.65299.61299.60
0.15376.30298.99299.36299.29299.29
0.20383.91298.44298.91298.82298.81
0.25388.99297.71298.27298.15298.14
0.30392.35296.75297.41297.26297.25
0.35394.54295.56296.30296.10296.09
0.40395.92294.13294.89294.66294.64
0.45396.75292.48293.18292.90292.89
0.50397.20290.66291.14290.84290.82
0.55397.36288.80288.80288.48288.47
0.60397.34287.12286.22285.92285.92
0.65397.16286.03283.57283.37283.39
0.70396.88286.27281.29281.29281.34
0.75396.50289.08280.25280.64280.75
0.80396.06296.64282.36283.41283.60
0.85395.56312.65291.61293.61293.93
0.90395.00343.43316.31319.38319.83
0.95426.89399.76373.78377.15377.60
1.00500.00500.00500.00500.00500.00
图1pe=0.7时的温度图线
表2pe=7.8时温度分布表
xFUDHBSEXPPL
0.00300.00300.00300.00300.00
0.05300.00300.00300.00300.00
0.10299.99299.99299.99299.99
0.15299.98299.98299.98299.98
0.20299.95299.96299.95299.95
0.25299.91299.92299.91299.91
0.30299.85299.86299.86299.86
0.35299.77299.78299.78299.78
0.40299.67299.68299.68299.68
0.45299.54299.56299.55299.55
0.50299.38299.40299.39299.39
0.55299.18299.21299.19299.19
0.60298.95298.98298.96298.96
0.65298.68298.72298.69298.69
0.70298.36298.41298.38298.38
0.75298.00298.05298.02298.02
0.80297.62297.65297.61297.61
0.85297.42297.19297.14297.14
0.90299.24296.68296.63296.63
0.95319.27296.10296.13296.06
1.00500.00500.00500.00500.00
表3pe=50时温度分布表
0.10300.00300.00300.00300.00
0.15300.00300.00300.00300.00
0.20299.99299.99299.99299.99
0.25299.99299.99299.99299.99
0.30299.98299.98299.98299.98
0.35299.97299.97299.97299.97
0.40299.95299.95299.95299.95
0.45299.93299.93299.93299.93
0.50299.91299.91299.90299.90
0.55299.88299.88299.87299.87
0.60299.84299.84299.84299.84
0.65299.80299.80299.80299.80
0.70299.75299.75299.75299.75
0.75299.69299.70299.69299.69
0.80299.63299.63299.63299.63
0.85299.56299.56299.55299.55
0.90299.56299.48299.47299.47
0.95