精确解： p=[1,5,10]; x=0:1/19:1; for i=1:1:3 for j=1:1:20 y(i,j)=(exp(p(1,i)*19*x(1,j))-1)/(exp(p(1,i)*19)-1); end plot(x,y(i,:)); hold on ; end 由题对中心差分、一阶迎风、混合格式进行模块编程： 他们之间可以通用，只需更改ae 关于p 的函数即可： 程序如下: (1)中心差分 p=[1,5,10]; for i=1:1:3 ae=1-0.5*p(1,i); x/L (Φ-ΦL )/(Φ0-ΦL ) 精确解图像
aw=p(1,i)+ae; ap=ae+aw; for i=1:1:18 for j=1:1:20 a(i,j)=0; end end for i=1:1:18 j=i; a(i,j)=aw; a(i,j+1)=-ap; a(i,j+2)=ae; end for i=1:1:17 n=i+1; for m=i:-1:1 b(1,1)=a(m,n); a(m,n)=-a(i+1,n)/a(i+1,n)*b(1,1)+a(m,n); a(m,n+1)=-a(i+1,n+1)/a(i+1,n)*b(1,1)+a(m,n+1); a(m,n+2)=-a(i+1,n+2)/a(i+1,n)*b(1,1)+a(m,n+2); end end F(1)=0; F(20)=1; F(19)=(-a(1,20)*F(20)-a(1,1)*F(1))/a(1,19); for i=2:1:18 F(i)=(-a(i,20)*F(20)-a(i,19)*F(19))/a(i,i); end x=0:1/19:1; y(1,:)=F; plot(x,y); hold on end
取步长δx=0.02。已知x=0,Φ=0;x=1,Φ=1.令k=ρｕ/Γ计算结果图表： 程序及数据结果： 追赶法： #include a[i]=2+0.02*k; b[i]=4; c[i]=2-0.02*k; f[i]=0; } tdma(a,b,c,f,x); for(i=0;i 4-1 解：采用区域离散方法A 时；网格划分如右图。内点采用中心差分 23278.87769.9 T T T === 22d T T=0dx - 有 i+1i 12 2+T 0i i T T T x ---=? 将2点,3点带入 32122 2+T 0T T T x --=? 即321 209T T -+= 432322+T 0T T T x --=?4321322+T 0T T T x --=? 即4 321 209 T T T -+-= 边界点4 （1）一阶截差 由x=1 1dT dx =，得 431 3 T T -= （2）二阶截差 11B M M q x x x T T S δδλλ -=++ 所以 434111. 1. 36311 T T T =++ 即 431 22293 T T -= 采用区域离散方法B 22d T T=0dx - 由控制容积法 0w e dT dT T x dT dT ????--?= ? ????? 所以代入2点4点有 322121011336 T T T T T ----= 即 239 028T T -= 544431011363 T T T T T ----= 即 34599 02828T T T -+= 对3点采用中心差分有 432 32 2+T 013T T T --=?? ??? 即 23499 01919 T T T -+= 对于点5 由x=1 1dT dx =，得 541 6 T T -= （1）精确解求左端点的热流密度 由 ()2 1 x x e T e e e -= -+ 所以有 ()22 20.64806911x x x x dT e e q e e dx e e λ -====- +=-=++ （2）由A 的一阶截差公式 21 0.247730.743113 x T T dT q dx λ =-=-= =?= （3）由B 的一阶截差公式 0 0.21640 0.649213 x dT q dx λ =-=-= = （4）由区域离散方法B 中的一阶截差公式： 210.108460.6504()B B T T dT dx x δ-?? ==?= ? ?? 通过对上述计算结果进行比较可得：区域离散B 有控制容积平衡法建立的离散方程与区域离散方程A 中具有二阶精度的格式精确度相当！ 4-3 解：将平板沿厚度方向3等分，如图 参考文献 [1]金国淼等.除尘设备[M].北京:化学工业出版社,2002:1-300 [2]Louis E. Stein, Alex. C. Hoffmann.旋风分离器-原理、设计和工程应用 [M].北京,化学工业出版社,2004:1-78 [3]国家环保局标准处.中华人民共和国国家标准环境空气质量标准[J],油气田环境保护,1996(04 ) [4]姚玉英,黄凤廉,陈常贵等.化工原理[M].天津:天津大学出版社,1999:138 [5]舒帆.影响旋风除尘器除尘效率的因素分析[J],粮食加工.2008, 33 (3):73-75 [6]韩占忠,王敬,兰小平.FLUENT流体工程仿真计算实例与应用[M].北京：北京理工大学出版社,2004:20 [7]魏志军,张平.旋风分离器气相流场的数值模拟[J].北京理工大学学报.2000, 20 (5):19-21 [8]嵇鹰,张红波,田耀鹏等.进口位置对旋风分离器特性影响的数值模拟[J].金属矿山,2008, 387 (3):127-129 [9]岑可法,倪明江,骆仲泱等.循环流化床锅炉理论设计与运行[M].北京:中国电力出版社, 2002:511-540 [10]陈明绍，吴光兴，张大中等.除尘技术的基本原理与应用[M].北京：中国建筑工业出版社，1981:333-518 [11]钱付平,章名耀.基于边界层理论旋风分离器分离效率的改进模型[J],中国电机工程学报.2007, 27 (5):71-74 [12]Hoffmann A C, Stein L E. Gas cyclones and twirl tubes：principles，design and operation [M]. Springer-Verlag，Berlin，Heidelberg，2002，169． [13]Leith D, Licth W. The collection efficiency of cyclone type particle collector. A new theoretical approach[J]. AIChE Symp Series,1972，126 (68)：196-206． [14]Obermair S，Woisetschlager J，Staudinger G．Investigation of the flow pattern in different dust outlet geometries of a gas cyclone by laser Doppler anemometry[J].Powder Technology，2003，2-3 (138):239-251 [15]Zhao Bingtao．Development of a new method for evaluating cyclone 【2-1】一食品冷藏室由内层为19 mm 厚的松木，中层为软木层，外层为51 mm 厚的混凝土所组成。内壁面温度为-17.8 ℃,混凝土外壁面温度为29.4 ℃。松木、软木和混凝土的平均热导率分别为, 3, W/(m ·K)，要求该冷藏室的热损失为15W/m 2。求所需软木的厚度及松木和软木接触面处的温度。 解：三层平壁的导热。 1）所需软木的厚度2b 由 ∑=-=3141i i i b T T q λ 得 151 .0019.00433.0762.0051.08.174.29152+++=b 解得： m b 128.02= 2）松木和软木接触面处的温度3T 由 151 .0019 .08.17153+==T q 解得：9.153-=T ℃ 解题要点：多层平壁热传导的应用。 【2-2】为减少热损失，在外径为150 mm 的饱和蒸汽管道外加有保温层。已知保温材料的热导率λ=+ 198 T(式中T 为℃),蒸汽管外壁温度为180 ℃,要求保温层外壁温度不超过50 ℃,每米管道由于热损失而造成蒸汽冷凝的量控制在1×10-4 kg/(m ·s)以下,问保温层厚度应为多少(计算时可假定蒸汽在180 ℃下冷凝)。 解：保温层平均热导率为： )./(126.02 501801098.1103.04K m W =+??+=-λ 由于本题已知的是蒸汽管道外壁面温度，即保温层内壁面温度，故为一层导热。 由 )()(21 221r r Ln T T L Q -=λπ 得： )()(21 221r r Ln T T L Q -=πλ （1） 式中：m W L Wr L Q /9.2011 103.20191013 4=???==- 将其及其它已知数据代入式（1）得： )075 .0()50180(126.029.2012r Ln -??=π 解得：m r 125.02= mm m 5005.0075.0125.0==-=∴δ壁厚 解题要点：单层圆筒壁热传导的应用。 【2-8】烤炉内在烤一块面包。已知炉壁温度为175 ℃，面包表面的黑度为,表面温度为100 ℃,表面积为 5 m 2，炉壁表面积远远大于面包表面积。求烤炉向这块面包辐射 传递的热量。 解：两物体构成封闭空间，且21S S <<，由下式计算辐射传热量： W T T S Q 0.65)448373(0645.085.01067.5) (448424111012-=-????=-=-εσ 负号表示炉壁向面包传递热量。 解题要点：辐射传热的应用，两个灰体构成的封闭空间。 【2-10】在逆流换热器中，用初温为20 ℃的水将1.25 kg/s 的液体［比热容为 kJ/(kg ·K)、密度为850 kg/m 3 ］由80 ℃冷却到30 ℃。换热器的列管直径为Φ25 mm ×2.5 mm,水走管内。水侧和液体侧的对流传热系数分别为850 W/（m 2·K ）和1 700W/（m 2·K ）,污垢热阻可忽略。若水的出口温度不能高于50 ℃，求水的流量和换热器的传热面积。 第二章 稳态导热 本章重点：具备利用导热微分方程式建立不同边界条件下稳态导热问题的数学模型的能力 第一节 通过平壁的导热 1-1 第一类边界条件 研究的问题： （1）几何条件：设有一单层平壁，厚度为δ，其宽度、高度远大于其厚度（宽度、高度是厚度的10倍以上）。这时可认为沿高度与宽度两个方向的温度变化率很小，温度只沿厚度方向发生变化。（属一维导热问题） （2）物理条件：无内热源，材料的导热系数λ为常数。 （3） 边界条件：假设平壁两侧表面分别保持均匀稳定的温度 1w t 和2w t ，21w w t t >。（为第一类边界条件，同时说明过程是稳态的） 求：平壁的温度分布及通过平壁的热流密度值。 方法1 导热微分方程： 采用直角坐标系，这是一个常物性、无内热源、一维稳态导热问题（温度只在 x 方向变化）。 导热微分方程式为：022=dx t d （2-1） 边界条件为：10w x t t == ， 2w x t t ==δ （2-2） 对式（2-1）连续积分两次，得其通解： 21c x c t += （2-3） 这里1c 、2c 为常数，由边界条件确定 ，解得：?? ???=-= 11221w w w t c t t c δ （2-4） 最后得单层平壁内的温度分布为： x t t t t w w w δ 2 11-- = （2-5） 由于δ 、1w t 、2w t 均为定值。所以温度分布成线性关系，即温度分布曲线的斜率是常数（温度梯度）， const t t dx dt w w =-=δ 1 2 （2-6） 热流密度为：)(21w w t t dx dt q -=-=δ λ λ 2/m W （2-7） 若表面积为 A, 在此条件下 , 通过平壁的导热热流量则为 : t A qA ?==Φδ λ W （2-8） 考虑导热系数随温度变化的情况： 对于导热系数随温度线形变化，即)1(0bt +=λλ，此时导热微分方程为：0=?? ? ??dx dt dx d λ 解这个方程，最后得： ?? ? ???++-+?? ? ?? +=+)(211212121121 122w w w w w w t t b x t t bt t bt t δ 或 x t t t t b b t b t w w w w w δ 12211)(2112 2-??????+++??? ??+=??? ??+ 说明：壁内温度不再是直线规律，而是按曲线变化。 对上式求导得：??? ?????+??? ??-=)1/(222bt dx dt b dx t d 因为 01>+bt ，02 >?? ? ??dx dt 所以 0>b ? 02 2 第二章 3.导热系数为常数的无内热源的平壁稳态导热过程，试问，若平壁两侧给定边界条件Tw1和Tw2，为什么这一导热过程的温度分布与平壁的材料无关？相同的平壁厚度，不同的平壁材料，仍给定第一类边界条件，热流密度是否相同？ （1）温度分布为 12 1w w w t t t t x δ -=- (设12w w t t >) 其与平壁的材料无关的根本原因在 coust λ=（即常物性假设） ，否则t 与平壁的材料有关 （2）由 dt q dx λ =- 知，q 与平壁的材料即物性有关 6.同上题，若已知边界条件为第三类，即已知Tf1，h1，Tf2，h2.试倒通过空心球壁热量的计算公式和球壁的传热热阻。 9.某教室有一层厚度为240mm 的砖层和一厚度为20mm 的灰泥构层。现安装空调设备，并在内表面加贴一层硬泡某塑料，是导入室内的热量比原来少了80%。已知砖的导热系数λ=0.7W/(m*k)，灰泥为λ=0.58W/(m*k)，硬泡某塑料的导热系数为λ=0.06W/(m*k)，试求出硬泡某塑料厚度。 已 知 ： 12240,20mm mm δδ==， 120.7/(),0.58/()W m k W m k λλ=?=? 3210.06/(),0.2W m k q q λ=?= 求：3δ 解： 设两种情况下的内外面墙壁温度12w w t t 和保持不变， 且12w w t t > 由题意知：12 112 12 w w t t q δδλλ-= + 12 23 12123 w w t t q δδδλλλ-= ++ 再由： 210.2q q =，有 12 12 3 12 1212 123 0.2 w w w w t t t t δδδδδλλλλλ--=+++ 22 131 3 1 2 tw 1 q tw 2 1 1 λ1 2 λ2 tw 1 tw 2 q 1 1λ1 2λ 2 3λ 3 传热学杨世铭第四版第二章答案 第二章 思考题 1 试写出导热傅里叶定律的一般形式，并说明其中各个符号的意义。 答：傅立叶定律的一般形式为：n x t gradt q ??-=λλ＝－，其中：gradt 为空间某点的温度梯度；n 是通过该点的等温线上的法向单位矢量，指向温度升高的方向；q 为该处的热流密度矢量。 2 已知导热物体中某点在x,y,z 三个方向上的热流密度分别为y x q q ,及z q ，如何 获得该点的 热密度矢量？ 答：k q j q i q q z y x ?+?+?=，其中k j i ,,分别为三个方向的单位矢量量。 3 试说明得出导热微分方程所依据的基本定律。 答：导热微分方程式所依据的基本定律有：傅立叶定律和能量守恒定律。 4 试分别用数学语言将传热学术语说明导热问题三种类型的边界条件。 答：① 第一类边界条件：)(01ττf t w =>时， ② 第二类边界条件：)()(02τλτf x t w =??->时 ③ 第三类边界条件：)()(f w w t t h x t -=??-λ 5 试说明串联热阻叠加原则的内容及其使用条件。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各串联环节热阻的和。使用条件是对于各个传热环节的传热面积必须相等。 7.通过圆筒壁的导热量仅与内、外半径之比有关而与半径的绝对值无关，而通过球壳的导热量计算式却与半径的绝对值有关，怎样理解？ 答：因为通过圆筒壁的导热热阻仅和圆筒壁的内外半径比值有关，而通过球壳的导热热阻却和球壳的绝对直径有关，所以绝对半径不同时，导热量不一样。 6 发生在一个短圆柱中的导热问题，在下列哪些情形下可以按一维问题来处理？ 答：当采用圆柱坐标系，沿半径方向的导热就可以按一维问题来处理。 8 扩展表面中的导热问题可以按一维问题来处理的条件是什么？有人认为，只要扩展表面细长，就可按一维问题来处理，你同意这种观点吗？ 答：只要满足等截面的直肋，就可按一维问题来处理。不同意，因为当扩展表面的截面不均时，不同截面上的热流密度不均匀，不可看作一维问题。 9 肋片高度增加引起两种效果：肋效率下降及散热表面积增加。因而有人认为，随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热数流量反而会下降。试分析这一观点的正确性。 第二章 思考题 1 试写出导热傅里叶定律的一般形式，并说明其中各个符号的意义。 答：傅立叶定律的一般形式为： n x t gradt q ??-=λλ＝－，其中：gradt 为空间某点的温度梯度；n 是通过该点的等温线上的法向单位矢量，指向温度升高的方向；q 为该处的热流 密度矢量。 2 已知导热物体中某点在x,y,z 三个方向上的热流密度分别为y x q q ,及z q ，如何获得该点的 热密度矢量？ 答：k q j q i q q z y x ?+?+?=，其中k j i ,,分别为三个方向的单位矢量量。 3 试说明得出导热微分方程所依据的基本定律。 答：导热微分方程式所依据的基本定律有：傅立叶定律和能量守恒定律。 4 试分别用数学语言将传热学术语说明导热问题三种类型的边界条件。 答：① 第一类边界条件：)(01ττf t w =>时， ② 第二类边界条件： )()( 02τλτf x t w =??->时 ③ 第三类边界条件： )()( f w w t t h x t -=??-λ 5 试说明串联热阻叠加原则的内容及其使用条件。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各串联环节热阻的和。使用条件是对于各个传热环节的传热面积必须相等。 7.通过圆筒壁的导热量仅与内、外半径之比有关而与半径的绝对值无关，而通过球壳的导热量计算式却与半径的绝对值有关，怎样理解？ 答：因为通过圆筒壁的导热热阻仅和圆筒壁的内外半径比值有关，而通过球壳的导热热阻却和球壳的绝对直径有关，所以绝对半径不同时，导热量不一样。 6 发生在一个短圆柱中的导热问题，在下列哪些情形下可以按一维问题来处理？ 答：当采用圆柱坐标系，沿半径方向的导热就可以按一维问题来处理。 8 扩展表面中的导热问题可以按一维问题来处理的条件是什么？有人认为，只要扩展表面细长，就可按一维问题来处理，你同意这种观点吗？ 答：只要满足等截面的直肋，就可按一维问题来处理。不同意，因为当扩展表面的截面不均时，不同截面上的热流密度不均匀，不可看作一维问题。 9 肋片高度增加引起两种效果：肋效率下降及散热表面积增加。因而有人认为，随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热数流量反而会下降。试分析这一观点的正确性。 答：错误，因为当肋片高度达到一定值时，通过该处截面的热流密度为零。通过肋片的热流已达到最大值，不会因为高度的增加而发生变化。 10 在式（2-57）所给出的分析解中，不出现导热物体的导热系数，请你提供理论依据。 答：由于式（2-57）所描述的问题为稳态导热，且物体的导热系数沿x 方向和y 方向的数值相等并为常数。 11 有人对二维矩形物体中的稳态无内热源常物性的导热问题进行了数值计算。矩形的一个边绝热，其余三个边均与温度为f t 的流体发生对流换热。你能预测他所得的温度场的解吗？ 答：能，因为在一边绝热其余三边为相同边界条件时，矩形物体内部的温度分布应为关于绝热边的中心线对称分布。 习题 主讲陶文铨 西安交通大学能源与动力工程学院热流中心CFD-NHT-EHT CENTER 2009年9月7日,西安 数值传热学 第一章绪论 课程简介 1. 教材－《数值传热学》第二版，2001 2. 学时－45学时理论教学；10学时程序教学 3. 考核－平时作业/计算机大作业： 考试－40/60；考查－60/40 4. 方法－开放，参与，应用 5. 助手－郭东之，周文静，李兆辉 有关的主要国外期刊 1.Numerical Heat Transfer, Part A-Applications; Part B- Fundamentals 2.International Journal of Numerical Methods in Fluids. https://www.doczj.com/doc/082403146.html,puter & Fluids 4.Journal of Computational Physics 5.International Journal of Numerical Methods in Engineering 6.International Journal of Numerical Methods in Heat and Fluid Flow https://www.doczj.com/doc/082403146.html,puter Methods of Applied Mechanics and Engineering 8.Engineering Computations 9.Progress in Computational Fluid Dynamics 10. Computer Modeling in Engineering & Sciences (CMES) 11.ASME Journal of Heat Transfer 12.International Journal of Heat and Mass Transfer 13.ASME Journal of Fluids Engineering 14.International Journal of Heat and Fluid Flow 15.AIAA Journal 第一章: 1-1 对于附图所示的两种水平夹层，试分析冷、热表面 间热量交换的方式有何不同？如果要通过实验来测定夹层中流体的导热系数，应采用哪一种布置？ 解：（a ）中热量交换的方式主要有热传导和热辐射。 （b ）热量交换的方式主要有热传导，自然对流和热辐射。 所以如果要通过实验来测定夹层中流体的导热系数，应采用（ a ）布置。 1-2 一炉子的炉墙厚13cm ，总面积为20m 2 ，平均导热系数为 1.04w/m 〃k ，内外壁温分别是520 ℃及50 ℃。试计算通过炉墙的热损失。如果所燃用的煤的发热量是 2.09 ×10 4 kJ/kg ，问每天因热损失要用掉多少千克煤？ 解：根据傅利叶公式 每天用煤 1-3 在一次测定空气横向流过单根圆管的对流换热实验中，得到下列数据：管壁平均温度t w = 69 ℃，空气温度t f = 20 ℃，管子外径d= 14mm ，加热段长80mm ，输入加热段的功率8.5w ，如果全部热量通过对流换热传给空气，试问此时的对流换热表面传热系数多大？ 解：根据牛顿冷却公式 1-4宇宙空间可近似的看作0K 的真空空间。一航天器在太空中飞行，其外表面平均温度为250K ，表面发射率为0.7 ，试计算航天器单位表面上的换热量？ 解：航天器单位表面上的换热量 1-5附图所示的空腔由两个平行黑体表面组成，孔腔内抽成真空，且空腔的厚度远小于其高度与宽度。其余已知条件如图。表面 2 是厚δ= 0.1m 的平板的一侧面，其另一侧表面 3 被高温流体加热，平板的平均导热系数λ=17.5w/m ? K ，试问在稳态工况下表面3 的t w3 温度为多少？ 解： 表面1 到表面2 的辐射换热量= 表面2 到表面3 的导热量 第二章: 西安交通大学西安交通大学《《数值传热学数值传热学》》课程大作业 20140114 一. 题目 （1） 百叶窗翅片的二维模型如图1 所示。在流动与换热已经进入周期性充分发展的阶段，可以取 出一个翅片单元进行传热与流动阻力的分析计算。在稳态，层流，常物性，翅片温度恒定的条件下，对于表1给定的几何尺寸，进行Re ＝10－500 范围内的数值模拟，揭示每个计算单元的平均Nu 数与阻力系数f 与Re 的关系； Nu ，f 以及Re 定为：1 12()Re ;;0.5p m m m dp dx L u L h L f Nu u νρλ==?= 其中m u 为来流平均速度；m h 为每块条片的平均换热系数。 表1 几何参数 L1/mm Tp/mm Lp/mm Delta/mm /θ 30 18.6 30 1.5 25 图1 百叶窗翅片二维模型 图2 阶梯型逼近 二. 建议建议与要求与要求 1. 为便于处理流固耦合问题，计算可对图1中打阴影线的区域进行； 2. 可采用图2 所示的阶梯型网格处理倾斜的翅片； 3. 按照《西安交通大学学报》的论文格式撰写本报告； 4. 2014年4月30号前交课程论文到东三楼204房间。 三. 参考文献 [1] 陶文铨编著，数值传热学（第二版），2001， 西安交通大学出版社，节11.2 [2] Wang L B, Tao, W Q. Numerical analysis on heat transfer and fluid flow for arrays of non-uniform plate length aligned at angles to the flow direction. Int J Numerical Methods for Heat and Fluid Flow , 1997, 7(5,6):496 [3] Gong L. Li Z Y, He Y L, Tao W Q. Discussion on numerical treatment of periodic boundary condition for temperature. Numerical Heat Transfer, Part B , 2007, 52(5):429-448 热物理过程的数值模拟Numerical Simulation of Thermophysics Process 讲稿 主讲：李隆键 第一章概论 1．1流动与传热过程的予测方法及特点 流动、传热、燃烧问题是热工类各专业和机械类动力机械专业所研究和解决的主要问题之一，燃烧问题实际上是有化学反应的流动与传热问题，推而广之，在所有热物理过程中，几乎都涉及到流动、传热问题。 预测的重要性： ①在规定设计参数的相应的结构下，热物理过程是否满足要求，达到预定的指 标？要预测； ②优化设计，不同方案的比较，要预测； ③减少设计、生产、再设计和再生产的费用； ④减少设计更改； ⑤减少试验和测量次数。 问题的核心：速度场、温度场（传热量）、浓度场等。 一、热物理问题的予测方法：理论分析法、实验测定、数值模拟 1、理论分析 以数学分析为基础，求解描述热物理过程的定解问题，获得函数形式的解，表示求解区域内物理量连续分布的场（速度场、温度场、浓度场……）。 控制方程+单值条件（数学模型）→理论解（分析解，解析解） 根据解的准确程度，又可再分为： （1）精确分析解（严格解） 特点：函数形式的解；它在求解区域精确地满足定解问题。 具体解法：直接积分法、分离变量法、积分变换法、热源法、映射法。 （2）近似分析解法 特点：函数形式的解，在求解区域上近似地满足定解问题（但在总量上满足相应的守恒原理，动量守恒、动量守恒、能量守恒、质量守恒）。 具体解法：积分法（从积分方程出发） 变分近似解法 摄动法（从微分方程出发） 2、实验测定 （1）纯实验法 （2）相似理论实验法：同类相似，减少变量数目→减少工作量，得到规律性结 数值传热学（Numerical Heat Transfer,NHT）又称计算传热学(Computational Heat Transfer,CHT)，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点(称为节点)上的值的集合来代替，通过一定的原则建立起这些离散点变量值之间关系的代数方程（称为离散方程，discretization equation）,求解所建立起来的代数方程已获得求解变量的近似值。 一、数值传热学的研究作用与地位 数值传热学在最近20年中得到飞速的发展，除了计算机硬件工业的发展给它提供了坚实的物质基础外，还主要因为无论分析的方法或实验的方法都有较大的限制，例如由于问题的复杂性，既无法做分析解，也因费用的昂贵而无力进行实验测定，而数值计算的方法正具有成本较低和能模拟复杂或较理想的过程等优点。经过一定考核的数值计算软件可以拓宽实验研究的范围，减少成本昂贵的实验工作量。在给定的参数下用计算机对现象进行一次数值模拟相当于进行一次数值实验，历史上也曾有过首先由数值模拟发现新现象而后由实验予以证实的例子。在这里要指出对数值模拟结果准确度应持正确认识。计算机本身不能创造信息，发现规律，它只是把人们送入的信息按照计算者所选定的规律进行处理，加工而已。但一旦建立了实际问题合理的数学模型，数值模拟又能发挥很大的作用。由于它本身的一些固定优点，它以发展成为工业界进行CAD/CAM及过程控制的重要手段,在多种工程领域中得到广泛应用。例如：叶轮机器粘性三元流体的计算，电站锅炉炉堂内流场与温度场的模拟；大型初见凝固过程中温度场的预测； 热流问题的数值计算 物理问题数值解的基本思想 把原来在空间与时间坐标中连续的物理量的场 (如速度场,温度场,浓度场等),用一系列有限 个离散点(称为节点,node)上的值的集合来代替; 通过一定的原则建立起这些离散点上变量值之间关 系的代数方程 (称为离散方程,discretization 大规模科学计算的重要性 传热与流动问题数值计算是应用计算机求解热量传 递过程中的速度场,温度场等的分支学科,是大规模 科学计算的重要组成部分,其重要性不言而喻. 2005年美国总统顾问委员会向美国总统提出要大 力发展计算科学以确保美国在世界上的竞争能力. 波音公司实现了对航空发动机的网格数达10亿量 级的直接数值模拟,以研究所设计发动机的性能. 现代科学研究的三大基本方法及其关系 课程简介 1、二维非稳态导热微分方程：S Y T X T t T p +??+??=)2222c （λδδρ。对于时间步进（x 方向，y 方向）及空间而言，该方程为何种类型的方程？ 解： 将二维非稳态导热微分方程化为：0c 2222=+-??+??S t T Y T X T p δδρλλ （1）x 方向：0,0,a ===c b λ。则：04b 2=-=?ac ，所以该二维非稳态导热方程为抛物型方程。 （2）y 方向：0,0,a ===c b λ。则：04b 2=-=?ac ，所以该二维非稳态导热方程也为抛物型方程。 （3）对于空间而言，二维非稳态导热方程可知：,0,a b c λλ===则：2240b ac λ?=-=-<，所以该二维非稳态导热方程为椭圆型方程。 2、 （补充不可压、常物性的条件。写出守恒型和非守恒型控制方程，并推导二者关系。） 解：由题可知，该流体为不可压缩、常物性流体，而且是有内热源的二维问题。 守恒型控制方程： 质量守恒方程：0=??+??y v x u ； 由于流体自身条件，使得0==v u S S ，得 动量守恒方程：()()??? ? ????+??+??-=??+??22221y u x u v x p y vu x uu ρ ()()??? ? ????+??+??-=??+??22221y v x v v y p y vv x uv ρ ； 能量守恒方程：()()T S y T x T a y vT x uT +??? ? ????+??=??+??2222 . 非守恒型控制方程： 质量守恒方程：无非守恒型 动量守恒方程：??? ? ????+??+??-=??+??22221y u x u v x p y u v x u u ρ ??? ? ????+??+??-=??+??22221y v x v v y p y v v x v u ρ ； 考博详解与指导 考试科目代码考试科目名称参考书目： 1001英语《新世纪研究生英语教材--阅读B,C》戴炜栋，柴小平编，上海外语教育出版社 1002俄语①《基础俄语》(1-3册)北京外语学院编，外语教学与研究出版社②《大学俄语基础教程》(1-3册)张智罗，高等教育出版社 1003日语《新编日语》(1-3册)周平、陈小芬，上海外语教育出版社 1004德语①《大学德语》戴鸣钟，高等教育出版社②《新编大学德语》朱建华编，外语教学与研究出版社，2002年9月第一版 1005法语《法语》(1-3册)马晓宏，外语教育出版社 2001工程流体力学①《工程流体力学》，归柯庭汪军王秋颖，科学出版社，2004年②《工程流体力学》(第二版)，孔珑，中国电力出版社，2007年 2002传热学《传热学》杨世铭，高等教育出版社，2006年 2003计算方法《数值分析》李庆杨等编著，清华大学出版社，2008年 2004高等光学《近代光学》袁一方译，高等教育出版社，1987年 2005物理光学《物理光学》梁铨庭，机械工业出版社 2006传感器技术及应用①《传感器》强锡富主编，机械工业出版社，2004年7月第三版②《非电量电测技术》严钟豪等主编，机械工业出版社，2003年1月第二版 2007激光原理《激光原理及应用》(第1-4章，6章)清华大学出版社 2008普通物理(光学)《普通物理学》(光学部分)程守洙，人民教育出版社 2009仪器电路原理与应用①《仪器电路设计与应用》，郝晓剑等编著，电子工业出版社，2007年6月②《基于运算放大器和模拟集成电路的电路设计》，赛尔吉欧。佛朗哥著西安交通大学出版社，2004年8月第1版 传热学思考题参考答案 第一章： 1、用铝制水壶烧开水时，尽管炉火很旺，但水壶仍安然无恙。而一旦壶内的水烧干后水壶很快就被烧坏。试从传热学的观点分析这一现象。 答：当壶内有水时，可以对壶底进行很好的冷却（水对壶底的对流换热系数大），壶底的热量被很快传走而不至于温度升得很高；当没有水时，和壶底发生对流换热的是气体，因为气体发生对流换热的表面换热系数小，壶底的热量不能很快被传走，故此壶底升温很快，容易被烧坏。 2、什么是串联热阻叠加原则，它在什么前提下成立？以固体中的导热为例，试讨论有哪些情况可能使热量传递方向上不同截面的热流量不相等。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各 串联环节热阻的和。例如：三块无限大平板叠加构成的平壁。例如通过圆筒壁，对于各个传热环节的传 热面积不相等，可能造成热量传递方向上不同截面的热流量不相等。 第二章： 1、扩展表面中的导热问题可以按一维问题处理的条件是什么?有人认为,只要扩展表面细长,就可按一维问题处理,你同意这种观点吗? 答：条件：（1）材料的导热系数，表面传热系数以及沿肋高方向的横截面积均各自为常数（2）肋片温度在垂直纸面方向（即长度方向）不发生变化，因此可取一个截面（即单位长度）来分析（3）表面上的换热热阻远远大于肋片中的导热热阻，因而在任一截面上肋片温度可认为是均匀的（4）肋片顶端可视为绝热。并不是扩展表面细长就可以按一维问题处理，必须满足上述四个假设才可视为一维问题。 2、肋片高度增加引起两种效果:肋效率下降及散热表面积增加。因而有人认为随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热流量会下降，试分析该观点的正确性。 答：的确肋片高度增加会导致肋效率下降及散热表面积增加，但是总的导热量是增加的，只是增加的部分的效率有所减低，所以我们要选择经济的肋片高度。 第三章： 1、由导热微分方程可知，非稳态导热只与热扩散率有关,而与导热系数无关。你认为对吗？答：错，方程的边界条件有可能与λ有关，只有当方程为拉普拉斯方程和边界条件为第一边界条件时才与λ无关。 2、对二维非稳态导热问题，能否将表面的对流换热量转换成控制方程中的内热源产生的热量？ 答：不能，二维问题存在边界微元和内边界微元，内边界微元不一定与边界换热，所以不存在源项。 第四章： 1、在第一类边界条件下，稳态无内热源导热物体的温度分布与物体的导热系数是否有关？为什么？ 答：无关，因为方程为拉普拉斯方程，边界为第一边界条件均与λ无关。 2、非稳态导热采用显式格式计算时会出现不稳定性，试述不稳定性的物理含义。如何防止这种不稳定性？ 答：物理意义：显示格式计算温度时对时间步长和空间步长有一定的限制，否则会出现不合 数值传热学第二章作业 2—1： POWER=input('POWER=?'); L1=input('L1=?'); M1=input('M1=?'); XL=input('XL=?'); YL=input('YL=?'); for i=2:L1 XF(i)=XL*((i-2)/(L1-2))^POWER; end for j=2:M1 YF(j)=YL*((j-2)/(M1-2))^POWER; end X(1)=0; for i=2:L1-1 X(i)=(XF(i)+XF(i+1))/2; end X(L1)=XF(L1); Y(1)=0; for j=2:M1-1 Y(j)=(YF(j)+YF(j+1))/2; end Y(M1)=YF(M1); for j=2:M1-1 plot(X(1),Y(j),'b.'); plot(X(1),Y(2), 'b.'); plot(X(L1),Y(j),'b.'); hold on end for i=2:L1-1 for j=1:M1 plot(X(i),Y(j),'b.'); hold on end end for i=2:L1 m=[XF(i),XF(i)]; n=[0,M1]; plot(m,n,'b-.'); hold on end for j=2:M1 m=[YF(j),YF(j)]; n=[0,L1]; plot(n,m,'b-.'); hold on end xlabel('x'); ylabel('y'); title('POWER= ') 运行结果如下： 2—3： 解：由2 2 2 1()u 2u u u x x y η ???== =???得： 原方程的守恒形式为： 2 2 2()2u u x y η ??=?? 对方程两端在t ?时间间隔内对其控制容积积分，把可积的部分积出后得： 22()t t s n e w t u u dtdy +?-? ? = 2t t e w t n s u u dtdx y y η+?????????-?? ? ????????? ? ? 选定2 u 随y 而变化的型线，这里取为阶梯式，即在控制容积内沿y 方向不变，则 2222 ()=y ()t t t t s n e w e w t t u u dtdy u u dt +?+?-?-? ? ? 选定2 u 随t 而变化的规律，这里采用阶梯式显式，则 2 2 ()t t e w t y u u dt +??-? = ()()22t t e w u u t y ??-?????? 选定 u y ??随x 而变化的型线，这里取为阶梯式，即在控制容积内沿x 方向不变，则 22t t t t e w t t n s n s u u u u dtdx x dt y y y y ηη+?+?????????????????-=?-???? ? ? ? ?????????? ???????? ? ? 选定 u y ??随t 而变化的规律，这里采用阶梯显式，则 2t t t n s u u x dt y y η+??????????-?? ? ???? ?????? = 2t t n s u u t x y y η?? ??????-???? ? ??????????? 进一步选取u 随x,y 分段线性变化，则 22 2 2 E P e u u u += ， 22 2 w 2 W P u u u +=数值传热学陶文铨第四章作业
旋风分离器参考文献
第二章 传热习题答案
传热学讲义设计—第二章
传热学第二章答案
最新传热学杨世铭第四版第二章答案
传热学第二章答案
数值传热学chapter_1
传热学_杨茉_部分习题与解答
西安交通大学西安交通大学《 《《《数值传热学数值传热学
热物理过程的数值模拟-计算传热学1
关于数值传热学的调研报告..
数值传热学(课件)-1
Numerical Simulations of Thermal & Fluid Problems
第一章 绪论
主讲 陶文铨
西安交通大学能源与动力工程学院 热流中心 CFD-NHT-EHT CENTER 2007年10月16日, 西安
1/88
equation);求解所建立起来的代数方程以获得所求
解变量的近似解.
2/88
3/88
理论分析
Analytical
实验研究
Experimental
数值模拟
Numerical
4/88
1. 学时- 30学时理论教学;6学时计算机作业 2. 考核- 平时作业/计算机大作业/考试: 20/30/50 3. 方法- 理解,参与,应用 努力将与数学处理相对应的物理背景联系起来理解. 4. 助手- 于乐 5. 参考教材-《计算流体力学与传热学》,中国建筑 工业出版社,1991
5/88数值传热学作业-第一章
上海理工大学博士研究生入学考试参考书目
传热学思考题参考答案(陶文铨第四版)
数值传热学第二章作业
相关主题
文本预览