热能工程硕士答辩ppt

1、谱方法求解变物性介质内辐射换热谱方法求解变物性介质内辐射换热姓名：马健姓名：马健 指导老师：李本文指导老师：李本文 教授教授材料电磁过程研究教育部重点实验室材料电磁过程研究教育部重点实验室2022-5-2411.研究背景研究背景2.谱方法简介3.谱方法求解变物性介质内辐射换热u3.1二维变吸收系数和散射系数介质内辐射换热求解u3.2二维非稳态变导热系数导热与辐射耦合求解4.结论与展望2022-5-2421.研究背景近几十年，辐射换热在民用领域和国防高科技领域得到了广泛的应用和发展。2022-5-24钢坯加热导弹尾焰分析卫星散热玻璃制造热红外成像发动机燃烧31.研究背景在实际生活和工业生产中，介

2、质的物性参数往往是其组分、压强和温度的复杂函数，如熔融的玻璃、厚厚的大气层、污浊的水、锅炉汽包中的水蒸气、炉膛中的高温烟气、泡沫耐火材料等。2022-5-24近三十年来学者对变物性介质内辐射换热的数值研究取得了巨大的成就，发展了各种各样的数值计算方法。1.蒙特卡罗法2.有限元法3.离散坐标法4.球形谐波法5.格子玻尔兹曼方法6. 低精度方法低精度方法41.研究背景2.谱方法简介谱方法简介3.谱方法求解变物性介质内辐射换热u3.1二维变吸收系数和散射系数辐射换热求解u3.2二维非稳态变导热系数导热与辐射耦合求解4.结论与展望2022-5-2452.谱方法简介谱方法是一种高精度数值的计算方法，目前

3、已广泛应用于大气扩散、海洋洋流、天气预报等实际问题。2022-5-241i-2i-1ii+1i+2N有限差分法1i-2i-1ii+1i+2N谱方法谱方法的全局性谱方法的指数收敛特性62.谱方法简介谱方法与有限差分法和有限元法都不同，在谱方法中试探函数被取为无穷可微的整体函数，一般是Sturm-Liouville问题的特征函数。按照正交多项式的类型，可把谱方法分为Fourier谱方法、Chebyshev谱方法、Legendre谱方法和Hermit谱方法等。我们采用Chebyshev谱方法求解变物性介质辐射换热问题。2022-5-2471.研究背景2.谱方法简介3.谱方法求解变物性介质内辐射换热谱

4、方法求解变物性介质内辐射换热u3.1二维变吸收系数和散射系数介质内辐射换热求解u3.2二维非稳态变导热系数导热与辐射耦合求解4.结论与展望2022-5-2483.谱方法求解变物性介质内辐射换热3.1二维变吸收系数和散射系数介质内辐射换热求解2022-5-24物理模型Tb,bTr,rTl,lTt,txy(1,0)(0,0)(0,1)12, ,gasTf x y f x y9假设：(1) 四个壁面均为不透明灰壁面(2) 介质折射率均匀为1(3) 介质为灰介质3.谱方法求解变物性介质内辐射换热求解变吸收系数和散射系数介质内辐射换热就是求解辐射传递方程。辐射传递方程描述的是辐射能量在介质内部传播时与介

5、质相互作用的过程。2022-5-24数学模型辐射传递方程边界条件4,(),( ,)(,) 4sasab I I II d r r r 01,0 wwww bwwwwI IIdn r rr nn10在角度上采用SRAPN格式的离散坐标法离散，在空间上采用Chebyshev-Gauss-Lobatto配置点离散，得到一系列离散的代数方程，将代数方程写为如下的矩阵形式。2022-5-24控制方程离散11111,2,yxNNmmmmmmmilljjlilijijijllA IB IX IFmM式中各矩阵元素如下 1,maxmin1,maxmin,4,mmm,1224mijx ijmijy ijija

6、ijs ijMg ijs ij,mija ijijmADxxBDyyXTCkIw3.谱方法求解变物性介质内辐射换热11由于吸收系数和散射系数是空间位置的函数，我们很难对矩阵方程直接求解。故在引入边界条件后，我们运用Matlab中的点乘和张量积将二维矩阵方程一维化，然后直接求解。如下所示：矩阵方程求解1111211,111122112212,112111,111,211,11ExxyxxxyxyyyxyxyyxmmmmmmNNNNmmmmmmNNNNmmmmmmNNNNNNNNNBBBBBBBBBEEEIIAIEEEIAIIAEE111222111yyymmmmmmmmmNNN XIFXIFXI

7、F3.谱方法求解变物性介质内辐射换热2022-5-24122022-5-24133.谱方法求解变物性介质内辐射换热节点数和方向数独立性测试NxNyqr*880.174%16160.037%24240.013%32320.008%4040当M=112时，节点数独立性测试Mqr*721.487%1120.643%1600.291%2160.083%280当NxNy=2424时，方向数独立性测试四个壁面均为冷壁面（0K），壁面发射率为1.0，介质为各向同性散射介质，介质吸收系数和散射系数均为1，介质温度为1000K。当吸收系数a=0, 散射系数s=1.0时。热壁面温度为1000K，其余均为冷壁面（0

8、K），介质为冷介质（0K）。壁面发射率均为1。我们将采用谱方法计算的结果与采用相同参数和边界条件的文献相对比。最大相对误差仅为0.37%。2022-5-24结果验证3.谱方法求解变物性介质内辐射换热热壁面无量纲辐射热流与文献中的对比Ratzel A C, Howell J R. Two-dimensional radiation in absorbing-emitting-scattering media using the P-N approximation J, ASME Journal, 1982, 36(1): 1-52.14这里我们讨论两个算例，分析不同空间位置变化的吸收系数和散射系

9、数对于辐射传热的影响。2022-5-243.谱方法求解变物性介质内辐射换热算例分析Case 1我们首先研究的四个壁面均为冷壁面(0K)，壁面发射率为1.0。介质温度为1000K，介质吸收系数为空间位置的函数，介质散射系数为0。152022-5-243.谱方法求解变物性介质内辐射换热16不同空间分布的吸收系数对各壁面无量纲辐射热流分布的影响，（a）a=1.0a下壁面右壁面左壁面上壁面xyx=0.5y=0.5y=xy=-x2022-5-243.谱方法求解变物性介质内辐射换热17不同空间分布的吸收系数对各壁面无量纲辐射热流分布的影响，（b） a=1.0+yb下壁面右壁面左壁面上壁面xyx=0.520

10、22-5-243.谱方法求解变物性介质内辐射换热18不同空间分布的吸收系数对各壁面无量纲辐射热流分布的影响，（c）a=1.0+xc下壁面右壁面左壁面上壁面xyy=0.52022-5-243.谱方法求解变物性介质内辐射换热19不同空间分布的吸收系数对各壁面无量纲辐射热流分布的影响，（d） a=1.0+x+yd下壁面右壁面左壁面上壁面xyy=x2022-5-243.谱方法求解变物性介质内辐射换热20不同空间分布的吸收系数对各壁面无量纲辐射热流分布的影响，（e）a=1.0-x+0.5ye下壁面右壁面左壁面上壁面xy2022-5-243.谱方法求解变物性介质内辐射换热21不同空间分布的吸收系数对各壁面

11、无量纲辐射热流分布的影响，（f）a=1.0+x2+y2f下壁面右壁面左壁面上壁面xyy=x2022-5-243.谱方法求解变物性参数介质内辐射换热Case 2当介质处于辐射热平衡时介质吸收的能量等于介质发出的能量。这里，我们研究的是辐射热平衡时不同吸收系数和散射系数对辐射换热的影响。此时左壁面为热壁面，温度为1000K，其他三个壁面温度为500K，壁面发射率均为1。介质处于辐射热平衡中。222022-5-243.谱方法求解变物性介质内辐射换热辐射热平衡时不同空间分布吸收系数对介质温度的影响，（a）Y=0.5，（b）X=0.523ab2022-5-243.谱方法求解变物性介质内辐射换热辐射热平衡

12、时不同空间分布吸收系数对热壁面无量纲辐射热流的影响24y=0.52022-5-243.谱方法求解变物性介质内辐射换热3.2二维非稳态变导热系数导热与辐射耦合求解物理模型下壁面冷壁面热壁面上壁面xy(1,0)(0,0)(0,1)=0+(T-Tref)25假设：(1) 四个壁面均为不透明灰壁面(2) 介质折射率均匀为1(3) 介质为灰介质2022-5-243.谱方法求解变物性介质内辐射换热数学模型导热与辐射耦合问题的数学模型包括辐射传递方程和能量方程。辐射传递方程表述的是辐射能束在介质中传播并与介质相互作用的过程的。能量方程实质是带有辐射热源项的导热微分方程。*4*41,( , ) (1) ( )

13、( ,)( ,)4IIId *r r rr 辐射传递方程边界条件*4* 01( ,)( )( , )0wmwwwwwwwI I|dn r rr n n 262022-5-243.谱方法求解变物性介质内辐射换热能量方程22224*2*2*411111(1)( , )4crIdxyxyNr 边界条件*(0,)1(1,)0.5(,0,)0.5(,1,)0.5yy x x在论文中，我们求解了第一类边界条件和混合边界条件（上下绝热）下变导热系数导热与辐射耦合问题。在这里，我们只讲述谱方法求解第一类边界条件下变导热系数导热与辐射耦合问题。27控制方程中的无量纲参数：无量纲时间：=2t，导热辐射参数：Ncr

14、=0/4T3，光学厚度：= L散射反照率：=s/ ，变导热系数参数：= T/ 0。其中为衰减系数。2022-5-243.谱方法求解变物性介质内辐射换热控制方程的离散对辐射传递方程，在角度上采用SRAPN格式的离散坐标法离散，在空间上采用Chebyshev-Gauss-Lobatto配置点离散。而对于能量方程，时间上采用全隐格式的差分离散，在空间上依然采用Chebyshev-Gauss-Lobatto配置点离散。得到一系列的代数方程，我们把这些代数方程写成矩阵形式。282022-5-243.谱方法求解变物性介质内辐射换热矩阵形式的辐射传递方程11*111,2,yxNNmmmmmilljjlili

15、jllP IQ IRmM式中各矩阵元素如下：*(1)*,maxmin(1)*,maxmin(1)*,maxmin4, 12122(1)4*mxilmilmxilmmjlyjlMmoldm oldm mmijijijmDil xxPDil xxQDyyRIw292022-5-243.谱方法求解变物性介质内辐射换热矩阵形式的能量方程1111yxNNilljjlilijijijllA B WC式中各矩阵元素如下：*2(2)*,maxmin2(2)*,maxmin1,2141,* ,1(1)1,*,11maxmin221/ 111124-xilx iljly jlnoldijijNMnoldmnmno

16、ldijijkjx ikmkcrijADxxBDyyWIwDNxxC *2221(1)1,*,1maxmin1,211yNnoldnikijy jkknoldijDyy 302022-5-24矩阵方程求解3.谱方法求解变物性介质内辐射换热对于矩阵形式的辐射传递方程，我们采用二维Schur分解求解，能量方程我们采用张量积将二维矩阵方程一维化后直接求解，我们将能量方程中辐射热源项和温度的一阶导数平方项放在方程右边进行迭代求解。根据文献定义稳态标准是在每一空间点上都满足| | 10-4，在我们的计算中，选取无量纲时间步长数=10-4，当|n+1n | 10-8时，我们判断此时已达到稳态。我们所有的结

17、果都是达到稳态后的结果。312022-5-24323.谱方法求解变物性介质内辐射换热节点数和方向数独立性测试当M=160，=5.0时节点数独立性测试当NxNy=3232，=5.0时方向数独立性测试由于无量纲辐射热流对节点数和方向数的变化很敏感，所以选取了无量纲辐射热流最敏感的情形，变导热系数参数=0，导热辐射参数Ncr=0.01，壁面发射率=1.0，散射反照率=0。NxNyqR*881.373%5.226%16160.161%0.470%24240.072%0.132%32320.033%0.051%40400.004%0.014%4848MqR*1120.205%2.010%1600.087

18、%0.672%2160.048%0.401%2800.015%0.219%3520.006%0.092%4322022-5-243.谱方法求解变物性介质内辐射换热结果验证Ncrcenterline at x/Lpresent workIEM-MFDMCDM-FDMLBM-FVMLBM-DOM0.010.30.7910.7910.7890.7830.7770.50.7250.7250.7250.7250.7220.70.6650.6630.6660.6760.67600.7600.7590.7590.7610.50.6630.6630.6630.6630.6670.70.59

19、40.5900.5940.5960.5981.00.30.7360.7330.7370.7370.7380.50.6300.6300.6300.6300.6310.70.5640.5600.5640.5640.565当=0，=1.0，=0，=1.0时不同导热辐射参数下方腔中心线Y=0.5处三点的温度值与文献的对比332022-5-243.谱方法求解变物性介质内辐射换热 Mishra S C, Talukdar P, Trimis D, Durst F. Two-dimensional transient conduction and radiation heat transfer with t

20、emperature dependent thermal conductivity J, International Communication Heat and Mass Transfer, 2005, 32(3-4): 305-314.不同变导热系数参数下方腔中心线Y=0.5温度变化曲线与文献中的对比，（a）=1.0，（b）=-1.034ab2022-5-243.谱方法求解变物性介质内辐射换热变导热系数参数对传热的影响分析当变导热系数参数=0时，介质导热系数不随温度变化而变化，当变导热系数参数0时，介质导热系数随温度变化而变化。当越大时，介质导热系数越小。352022-5-243.谱方法求

21、解变物性介质内辐射换热当Ncr=0.1，=0，=1.0，=1.0时不同变导热系数参数下介质等温线图，（a）=-1.0，（b） =0，（c） =1.036abc2022-5-243.谱方法求解变物性介质内辐射换热不同变导热系数参数下介质中心线上温度分布图，（a）Y=0.5，（b） X=0.537ab2022-5-243.谱方法求解变物性介质内辐射换热不同变导热系数参数对热壁面无量纲热流分布的影响，（a）导热热流，（b） 辐射热流，（c）总热流38abc2022-5-243.谱方法求解变物性介质内辐射换热不同变导热系数参数对各壁面无量纲辐射热流所占比例的影响，（a）热壁面，（b） 冷壁面，（c）上

22、壁面39abc2022-5-24在本论文中，我们还详细讨论了导热辐射参数、散射反照率，壁面发射率、光学厚度、长宽比对介质温度和壁面热流的影响。并且我们还求解了混合边界条件下（上下壁面绝热）变导热系数导热与辐射耦合问题，并研究了变导热系数参数和散射相函数对介质温度和壁面热流的影响。3.谱方法求解变物性介质内辐射换热401.研究背景2.谱方法简介3.谱方法求解变物性介质内辐射换热u3.1二维变吸收系数和散射系数介质内辐射换热求解u3.2二维非稳态变导热系数导热与辐射耦合求解4.结论与展望结论与展望2022-5-24414.结论与展望1.采用Chebyshev谱方法求解了二维变吸收系数和散射系数介质辐射换热，并根据对称性验证了计算结果，表明计算结果是准确的。不同空间分布的吸收系数和散射系数对处于热平衡状态下介质温度影响不大，但改变了壁面热流的对称性，且随着吸收系数和散射系数的增加，热壁面辐射热流明显减少。2.采用Chebyshev谱方法求解了第一类边界条件和混合边界条件下二维变导热系数导热与辐射耦合问题。对常导热系数和变导热系数与文献中的结果进行了对比，发现即使采用很少的节点，谱方法也能达到很高的精度。随着变导热系数参数的增加，介质温度明显增加，热壁面辐射热流降低，导热热流也降低，各壁面辐射热流所占比例增加。20