下载文档原格式
导热问题的数值求解方法数值解法的基本思想是用空间和时间区域内有限个离散点（称为节点）上温度的近似值，代替物体内实际的连续温度分布，然后由导热方程和边界条件推导出各节点温度间的相互关系的代数方程组（称为离散方程），求解此方程组，得到节点上的温度值，此即物体中温度场的解。
只要节点分布的足够稠密，数值解就有足够的精度。
求解导热问题的数值方法有有限差分法及有限元法，近几年又发展了边界元法和有限分析法。
数值方法适用于求解各种导热问题，不管物体的几何形状有多复杂，不管线性或非线性问题，都能使用。
由于计算机的飞速发展，计算技术软件发展也很快，数值方法的的地位越来越重要。
1 数值求解的基本思路及稳态导热内节点离散方程的建立一、 解法的基本思路1、基本思路：数值解法的求解过程可用框图4-1表示。
由此可见：1）物理模型简化成数学模型是基础；2）建立节点离散方程是关键；3）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。
二、稳态导热中位于计算区域内部的节点离散方程的建立方法1、基本方法方法：①泰勒级数展开法；②热平衡法。
1）泰勒级数展开法如图4-3所示，以节点(m,n)处的二阶偏导数为例，对节点(m+1,n)及(m-1,n)分别写出函数t 对(m,n)点的泰勒级数展开式：对(m+1,n)：+∂∂∆+∂∂∆+∂∂∆+∂∂∆+=+444333,222,,,12462x t x x t x x t x x t x t t n m n m n m n m （a ）对（m-1,n ）：+∂∂∆+∂∂∆-∂∂∆+∂∂∆-=-444333,222,,,12462x t x x t x x t x x t xt t n m n m n m n m （b ）（a ）+（b ）得： +∂∂∆+∂∂∆+=+-+444,222,,1,1122x t x x t x t t t n m n m n m n m 变形为n m x t,22∂∂的表示式得：n m x t,22∂∂)(0222,1,,1x x t t t nm n m n m ∆+∆+-=-+ 上式是用三个离散点上的值计算二阶导数n m x t ,22∂∂的严格表达式，其中：)(02x ∆―― 称截断误差，误差量级为2x ∆在数值计算时，用三个相邻节点上的值近似表示二阶导数的表达式即可，则相应的略去)(02x ∆。
传热学知识点总结传热学知识点总结传热学，是研究热量传递规律的科学，是研究由温差引起的热能传递规律的科学。
大约在上世纪30年代，传热学形成了独立的学科。
以下是小编整理的传热学知识点总结，欢迎阅读！第一章§1-1 “三个W”§1-2 热量传递的三种基本方式§1-3 传热过程和传热系数要求：通过本章的学习，读者应对热量传递的三种基本方式、传热过程及热阻的概念有所了解，并能进行简单的计算，能对工程实际中简单的传热问题进行分析（有哪些热量传递方式和环节）。
作为绪论，本章对全书的主要内容作了初步概括但没有深化，具体更深入的讨论在随后的章节中体现。
本章重点：1.传热学研究的基本问题物体内部温度分布的计算方法热量的传递速率增强或削弱热传递速率的方法2.热量传递的三种基本方式(1).导热：依靠微观粒子的热运动而产生的热量传递。
传热学重点研究的是在宏观温差作用下所发生的热量传递。
傅立叶导热公式：(2).对流换热：当流体流过物体表面时所发生的热量传递过程。
牛顿冷却公式：(3).辐射换热：任何一个处于绝对零度以上的物体都具有发射热辐射和吸收热辐射的能力，辐射换热就是这两个过程共同作用的结果。
由于电磁波只能直线传播，所以只有两个物体相互看得见的部分才能发生辐射换热。
黑体热辐射公式：实际物体热辐射：3.传热过程及传热系数：热量从固壁一侧的流体通过固壁传向另一侧流体的过程。
最简单的传热过程由三个环节串联组成。
4.传热学研究的基础傅立叶定律能量守恒定律+ 牛顿冷却公式 + 质量动量守恒定律四次方定律本章难点1.对三种传热形式关系的理解各种方式热量传递的机理不同，但却可以（串联或并联）同时存在于一个传热现象中。
2.热阻概念的理解严格讲热阻只适用于一维热量传递过程，且在传递过程中热量不能有任何形式的损耗。
思考题：1.冬天经太阳晒过的棉被盖起来很暖和，经过拍打以后，效果更加明显。
为什么？2.试分析室内暖气片的散热过程。
第五章复习题1、试用简明的语言说明热边界层的概念。
答：在壁面附近的一个薄层内，流体温度在壁面的法线方向上发生剧烈变化，而在此薄层之外，流体的温度梯度几乎为零，固体表面附近流体温度发生剧烈变化的这一薄层称为温度边界层或热边界层。
2、与完全的能量方程相比，边界层能量方程最重要的特点是什么？答：与完全的能量方程相比，它忽略了主流方向温度的次变化率σα22x A ，因此仅适用于边界层内，不适用整个流体。
3、式（5—4）与导热问题的第三类边界条件式（2—17）有什么区别？答：=∂∆∂-=yyt th λ（5—4）)()(f w t t h h t-=∂∂-λ （2—11）式（5—4）中的h 是未知量，而式（2—17）中的h 是作为已知的边界条件给出，此外（2—17）中的λ为固体导热系数而此式为流体导热系数，式（5—4）将用来导出一个包括h 的无量纲数，只是局部表面传热系数，而整个换热表面的表面系数应该把牛顿冷却公式应用到整个表面而得出。
4、式（5—4）表面，在边界上垂直壁面的热量传递完全依靠导热，那么在对流换热中，流体的流动起什么作用？答：固体表面所形成的边界层的厚度除了与流体的粘性有关外还与主流区的速度有关，流动速度越大，边界层越薄，因此导热的热阻也就越小，因此起到影响传热大小5、对流换热问题完整的数字描述应包括什么内容？既然对大多数实际对流传热问题尚无法求得其精确解，那么建立对流换热问题的数字描述有什么意义？答：对流换热问题完整的数字描述应包括：对流换热微分方程组及定解条件，定解条件包括，（1）初始条件 （2）边界条件 （速度、压力及温度）建立对流换热问题的数字描述目的在于找出影响对流换热中各物理量之间的相互制约关系，每一种关系都必须满足动量，能量和质量守恒关系，避免在研究遗漏某种物理因素。
基本概念与定性分析5-1 、对于流体外标平板的流动，试用数量级分析的方法，从动量方程引出边界层厚度的如下变化关系式：x xRe 1~δ解：对于流体外标平板的流动，其动量方程为：221xy u v dx d y u v x y u ∂+-=∂∂+∂∂ρρ 根据数量级的关系，主流方的数量级为1，y 方线的数量级为δ则有2211111111δρδδv +⨯-=⨯+⨯ 从上式可以看出等式左侧的数量级为1级，那么，等式右侧也是数量级为1级， 为使等式是数量级为1，则v 必须是2δ量级。
第5章 非稳态问题的求解方法1.1 通用输运方程()()()()()t t f q Γv tφφρφρφφ,grad div div =++-=∂∂ ( 5-1 )5.1 显式Euler 方法考虑1D, 定速度，常物性，无源项的特例22xx u t ∂∂Γ+∂∂-=∂∂φρφφ ( 5-2 ) 时间向前，空间中心差分，得FD 与FV 相同形式代数方程()t x x u nin i n i n i n i nin i∆⎥⎦⎤⎢⎣⎡∆-+Γ+∆--+=-+-++21111122φφφρφφφφ( 5-3 ) 可写成()ni n i n i n i c d c d d 1112221-++⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-+-=φφφφ ( 5-4 ) 其中()xtu c and x t d ∆∆=∆Γ∆=2ρ ( 5-5 ) d 表示时间步长与特征扩散时间()Γ∆/2ξρ的比。
后者代表一个扰动由于扩散通过∆x 一段距离所需时间。
c 表示时间步长与特性对流传递时间x u ∆/的比。
后者代表一个扰动由于对流通过∆x 一段距离所需时间。
c 成为Courant number, 为CFD 中一个关键的参数。
此格式为时间为1阶精度，空间为2阶精度。
方程（4）内的系数在某些条件下，可能会是负值。
用矩阵表示：n n A φφ=+1 ( 5-6 )观察函数：()∑---=-=in i ni n n 211φφφφε( 5-7 )如果系数矩阵A 的本征值中有大于1，则ε随着n 的增加而增加。
如果本征值全部小于1，则ε是递减的。
一般本征值很难求得，对于本特例，它的解可用复数形式表示ji n n j e ασφ= ( 5-8 )其中，α为波数，可取任意值。
∙ 无条件发散：φn 无条件随n 增加→|σ|>1 ∙无条件稳定：φn 无条件随n 降低→|σ|<1代入差分方程，得到本征值为：()αασsin 2cos 21c i d +1-+= ( 5-9 )考虑特殊情况，∙ 无扩散：d=0, →σ >0, 无条件发散，充分条件∙无对流：c=0, →当cos α= -1时，σ最大，→d<1/2，无条件收敛，充分条件从另一个稳定条件考虑，要求系数矩阵A 的所有系数为正，可得到类似稳定性条件：（充分条件）d c d 2and 5.0<<( 5-10 )第一个条件要求()Γ∆<∆22x t ρ ( 5-11 )表示，每当∆x 减少一半，时间步长需减少到1/4. 第二个条件要求2Pe or2<<Γ∆cell xu ρ ( 5-12 )这同前述的用1D 稳态对流/扩散问题的CDS 要求是一致的。
