传热学-导热数值计算共46页
传热学导热问题的数值解法
导热问题的数值解法1 、重点内容：① 掌握导热问题数值解法的基本思路；② 利用热平衡法和泰勒级数展开法建立节点的离散方程。
2 、掌握内容：数值解法的实质。
3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。
由前述3 可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。
但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。
随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种：（１）有限差分法（ 2 ）有限元方法（ 3 ）边界元方法数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。
如：几何形状、边界条件复杂、物性不均、多维导热问题。
分析解法与数值解法的异同点：相同点：根本目的是相同的，即确定① t=f(x ，y ，z) ；②。
不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。
§4-1 导热问题数值求解的基本思想及内节点离散方程的建立实质对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。
该方法称为数值解法。
这些离散点上被求物理量值的集合称为该物理量的数值解。
2 、基本思路：数值解法的求解过程可用框图4-1 表示。
由此可见：1 ）物理模型简化成数学模型是基础；2 ）建立节点离散方程是关键；3 ）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。
一数值求解的步骤如图4-2 （a ），二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下：1 建立控制方程及定解条件控制方程：是指描写物理问题的微分方程针对图示的导热问题，它的控制方程（即导热微分方程）为：（a ）边界条件：x=0 时，x=H 时，当y=0 时，当y=W 时，区域离散化（确立节点）用一系列与坐标轴平行的网格线把求解区域划分成若干个子区域，用网格线的交点作为需要确定温度值的空间位置，称为节点( 结点) ，节点的位置用该节点在两个方向上的标号m ，n 表示。
《传热学》第4章-导热问题的数值解法
数值解法的基本思想
v 用导热问题所涉及的空间和时间区域内有限 个离散点(称为节点)的温度近似值来代替物体 内实际连续的温度分布，将连续温度分布函 数的求解问题转化为各节点温度值的求解问 题，将导热微分方程的求解问题转化为节点 温度代数方程的求解问题。因此，求解域的 离散化、节点温度代数方程组的建立与求解 是数值解法的主要内容。
= ti, j
−
∂t ∂x
i
,
j
∆x
+
∂2 ∂x
t
2
i, j
∆x 2 2!
−
∂ 3t ∂x 3
i, j
∆x 3 3!
+ ...
∂t ∂x
i,
j
=
ti, j
− ti−1, j ∆x
+ O(∆x)
一阶截差公式（向后差分）
ti+1, j
= ti, j
4适用于内节点和边界节点3二控制容积热平衡法0nsew根据导热付里叶定律对于垂直于纸面方向单位宽度而言01111??????????????????yttxyttxxttyxttyjijijijijijijijixttyjijiw?????1xttyjijie????1yttxjijis?????1yttxjijin????1二控制容积热平衡法如果选择步长??xy01111??????????????????yttxyttxxttyxttyjijijijijijijijitttttijijijijij???111140二维稳态导热均匀步长情况下的节点温度差分方程1上上式为内部节点温度差分方程二控制容积热平衡法2边界节点温度差分方程第一类边界条件边界节点温度已知
《传热学》第四章 导热数值解法基础
边界
2.第二类边界条件：
Байду номын сангаас
Δx＝Δy时简化为：
绝热边界：
3.第三类边界条件：
Δx＝Δy时简化为：
其他情况的节点方程 ——见教材表4－1
外拐角与内拐角节点
对流边界内部拐角节点热平衡：
节点方程式推导实例 ——对流边界外部拐角节点
Δx＝Δy时简化为：
数值导热离散方程组＝内节点离散方程＋边界节点离散方程
二、常用计算软件
1.MATLAB——矩阵计算软件
matlab软件主界面
2.FLUENT——流体流动通用数值计算软件
3. FLUENT AIRPAK ——人工环境系统分析软件，暖通空调专业和传热学领域必备软件
AIRPAK模拟温度场
第四章重点： 1.有限差分方程的建立 2.高斯－赛德尔迭代方法
谢谢观看
《传热学》
第四章 导热数值解法基础
本章研究的目的 ——利用计算机求解难以用 分析解求解的导热问题 基本思想 ——把原来在时间、空间坐 标系中连续的物理量的场， 用有限个离散点的值的集合 来代替，通过求解按一定方 法建立起来的关于这些值的 代数方程，来获得离散点 上被求物理量的值。 研究手段——有限差分法
数值导热离散方程组内节点离散方程边界节点离散方程三节点离散方程组的求解迭代法迭代法的原理离散方程组的求解方法消元法方程过多时计算机内存不足迭代法假定初值根据假定的初值求新值并重复此步骤若干次两次计算值足够接近认为达到真实值简单迭代法每次迭代时使用上次迭代的结果允许误差简单迭代法的缺点由于每次迭代中使用与真实值偏差较大的上次迭代的旧值使运算过程接近真实值的时间增加高斯赛德尔迭代法将本次迭代的最新结果立刻代入本次迭代过程计算其他未知值高斯赛德尔迭代法的优点由于每次迭代中使用与真实值偏差较小的本次迭代的新值使运算过程接近真实值的时间缩短第三节非稳态导热的数值计算一显式差分格式研究对象一维非稳态导热问题一维非稳态导热内节点差分方程
传热学：第四章 导热问题数值解法
t m,n
1 t m 1,n t m 1,n t m ,n 1 t m ,n 1 4
•二维导热问题；网格线；
沿x、y方向的间距为x、 y；网格单元。
每个节点温度就代表了它 所在网格单元的温度。 p(m,n)
•此方法求得的温度场
在空间上不连续。
•网格越细密、节点越多，结果越接近分析解 •网格越细密，计算所花时间越长
2) 数值计算法，把原来在时间和空间连续的物理量的
场，用有限个离散点上的值的集合来代替，通过求解
按一定方法建立起来的关于这些值的代数方程，从而
获得离散点上被求物理量的值；并称之为数值解；
3) 实验法 就是在传热学基本理论的指导下，采用实
验对所研究对象的传热过程进行测量的方法。 3 三种方法的特点 1) 分析法 a 能获得所研究问题的精确解，可以为实验和数值 计算提供比较依据；
t m,n 1 2t m,n t m,n 1 2t 同理： 2 y y 2 m,n
将以上两式代入导热微分方程得到节点（m,n）的温 度离散方程： t tm,n1 2tm,n tm,n1 m 1, n 2t m , n t m 1, n 0 2 2 x y
x y 上式可简化
第三类边界条件： y x
qw h(t f tm,n )
2hx 2hx x 2 tm1,n tm,n1 2 tf 0 tm,n 2
(3) 内部角点
y t m 1,n t m ,n y y qw 2 x x 2 t m ,n 1 t m ,n x x t m ,n 1 t m ,n x qw 2 y 2 y 3xy 0 4
传热学数值计算
Fe aE De 2
Thermal
Fw Fe aP Dw De aW aE Fe Fw 2 2
2、对方程的几点说明 由于连续性，Fe=Fw， aP aW aE（只是在流 场满足连续性条件时才具有这一性质）； 方程 aP P aE E aW W 隐含着分段线性分布的含
讨论只有对流项和扩散项存在时的一维稳态问题，控 制方程为：
u j ( ) S x j x j x j
d d d u ( ) dx dx dx
d 连续方程： u 0 dx
u const
任务：导出相应方程的离散化形式
义，也是熟知的中心差分格式（用左右节点值表示
界面上的值以及界面上的导数值）； 方程必须遵守四项基本法则，否则会产生灾难性的 结果。
2018-11-24
太 原 理 工 大 学
7 /70
Thermal
例如：设 De Dw 1， Fe Fw 4 若E、W给定，即可由离散方程求得P 。
即， F 2D时，有可能使 aE 或aW 为负 产生不切实际的结果
2018-11-24
太 原 理 工 大 学
8 /70
Thermal
这就是中心差分格式求解对流换热问题时仅限于低
Fw Fe Re(低的F/D)的原因 . aP Dw De aW aE Fe Fw 2 2
原通用方程可改写为
u j ( )S x j x j x j
对于已知的ρ、uj、Γ及S（常量）的分布，任何解及
+c 将同时满足方程，故系数和的法则仍然适用。
2018-11-24
太 原 理 工 大 学
传热学 -- 导热数值分析解法
•
else T(i,1)=(qw*0.5*b+0.5*h*b*tf+0.5*a*T(i-1,1)+0.5*a*T(i,j+1))/(a+0.5*h*b);
•
end
•
end
• end
• for j=2:5
•
for i=1:4
•
c(i,j)=T(i,j);
•
if i==1
•
T(i,j)=(a*T(i,j+1)+0.5*a*T(i,j-1)+a*T(i+1,j))/(2*a);•ຫໍສະໝຸດ end•end
•
end
• end
• for i=1:4
•
for j=1:5
•
z=max(T(i,j)-c(i,j));
•
eps=z;
•
end
•
end
• end
• surf(T);
下图为本题的温度分布
节点越多，温度场划分越细。 如下图以10mm为步长划分节点
谢谢观看！
根据题意进行边界条件分析及节点划分，如 下：
6 5 4
3 2
1
2
3
4
对每个节点按热平衡法列稳态方程:
1) 若 ti,j为平壁的内部节点，则有
ti-1,j-ti,j y+ ti+1,j -ti,j y ti,j-1-ti,j x+ ti,j+1-ti,j x 0
x
x
•
else if i>1&&i<4
•
T(i,j)=(a*T(i,j+1)+a*T(i+1,j)+a*T(i-1,j)+a*T(i,j-1))/(4*a);
热传导和导热系数的计算
热传导和导热系数的计算热传导是指热量在物体内部由高温区向低温区传递的过程，它是固体、液体和气体等物质的一种基本热传递方式。
热传导的计算通常涉及到导热系数这个物理量，它是一个材料特性，用来描述材料内部热量传递的能力。
一、热传导的基本公式1.一维稳态热传导：对于一维稳态热传导，热量在物体内部的传递可以用傅里叶定律来描述：[ q = -kA ]其中，( q ) 是单位面积的热流量（W/m^2），( k ) 是导热系数（W/m·K），( A ) 是物体的横截面积（m^2），( ) 是温度梯度（K/m）。
2.二维和三维稳态热传导：对于二维和三维稳态热传导，热量在物体内部的传递可以用傅里叶定律的微分形式来描述：[ = ]其中，( q ) 是单位体积的热流量（W/m^3），( t ) 是时间（s），( ) 是热扩散系数（m^2/s），( T ) 是温度（K或°C），( ) 是温度梯度的二阶导数。
二、导热系数的定义和影响因素导热系数（k）是描述材料内部热量传递能力的物理量，单位为W/m·K。
导热系数反映了材料在单位厚度、单位温差条件下，单位时间内通过单位面积的热量。
2.影响因素：a)材料的种类：不同材料的导热系数不同，金属的导热系数一般较大，而绝缘材料的导热系数较小。
b)温度：材料的导热系数随温度的变化而变化，一般情况下，随着温度的升高，导热系数增大。
c)湿度：对于多孔材料，湿度对导热系数有较大影响，湿度越大，导热系数越大。
d)孔隙率：对于多孔材料，孔隙率越大，导热系数越小。
三、常见材料的导热系数以下是一些常见材料的导热系数（单位：W/m·K）：1.金属：40-460（如铜：380，铝：237）2.木材：0.1-0.2（如松木：0.14，柚木：0.2）3.塑料：0.1-1.5（如聚乙烯：0.4，聚丙烯：1.0）4.玻璃：1-2（如普通玻璃：1.1，高强度玻璃：1.6）5.空气：0.026（在常温常压下）四、热传导和导热系数的应用1.建筑领域：热传导和导热系数的计算在建筑领域具有重要意义，可以用于设计保温层、隔热材料等，以提高建筑的能源效率。
传热学-第4章-导热数值解法基础
环境与能源工程学院 ( SEEE )
二阶导数的中心差分表达式 [1]-[2]相加
在表示温度对时间的一阶导数时只采用向前或向后差分 表达式——温度对时间的中心差分表达式求解非稳态导 表达式 中心差分 热问题将导致数值解的不稳定，参考：《工程传热学》 导致数值解的不稳定 以常物性、无内热源、二维稳态导热为例 P(I,j)
2ΔxqW 1 = (2ti −1, j + ti , j +1 + ti , j −1 + ) λ 4
第三类边界条件
i,j+1 i1,j i,j
Beijing University of Civil Engineering and Architecture ( BUCEA )
传热学 ( Heat transfer )
环境与能源工程学院：School of Environment and Energy Engineering -SEEE
第四章 导热数值解法基础 4-1 建立离散方程的方法 4-2 稳态导热的数值计算 4-3 非稳态导热的数值计算
ti +1, j ⎛ ∂ 2t ⎞ Δx 2 ⎛ ∂ 3t ⎞ Δx 3 ⎛ ∂t ⎞ = ti , j + ⎜ ⎟ Δx + ⎜ 2 ⎟ ⎜ ∂x ⎟ 2! + ⎜ ∂x 3 ⎟ 3! + LL ⎜ ⎟ ⎝ ∂x ⎠i , j ⎝ ⎠i , j ⎝ ⎠i , j
节点（i, j）一阶导数的向后差分表达式（一阶截差公式）；
P(i.j)
q ∂t ∂ 2t ∂ 2t = a( 2 + 2 ) + v ρc ∂τ ∂x ∂y
二维导热问题；网格；沿x、 y方向的间距为Δx、Δy； 网格单元 （Nodal point） 节点：网格线的交点；p(i,j) 节点
传热学 数值计算
数值计算
4、一无限大平壁厚度为 0.3m，其导热系数为λ =36.4 W/ (m·K)。平壁两侧表面均给定 为第三类边界条件，即 h 1 =60W/ (m 2 ·K)，t f 1 =25℃；h 2 =300W/ (m 2 ·K)，t f 2 =215℃。 当平壁中具有均匀的内热源 q v =2 10 W / m 时，试计算沿平壁厚度的稳态温度分布。
k
0
k k 1
11
k 2Fo (t 10 t f Bi ) (1 2Fo 2Bi Fo) t11 t11
0
四、
计算过程
⑴ 设定初值：
t (1~11) 35 ℃;
36W / m k ;
Bi h / ;
根据不同的 Fo 计算Δ τ ：
qv x 2
)/2
h2 x
h2 x TRB
qv x 2
) /(1
)
|T[i]-t[i]|<=EPS
NO
IT=IT+1
YES
打印“t[i]” ， “IT”
YES
IT>K
NO
停止
⑷ 程序与计算结果 #include"iostream.h" #include"iomanip.h" #include"math.h" #include"stdio.h" #define N 16 void main() { int M,i,IT,flag; //定义节点个数
计算结果：
各节点温度： 节点 1 2 411.24 10 422.65 3 420.33 11 414.23 4 427.22 12 403.62 5
传热学课件第四章 导热问题数值解法基础
i , j
t x
t i 1 , j t i , j x
0 x
2.一阶导级的向后差分表达式：舍去<2>式△x2后各项，则有：
i , j
t x
t i , j t i 1 , j x
0 x
第一节 建立离散方程的方法
二、泰勒级数展开法（有限差分法）
k 2 k 1
对 流 h t f t1 A
k k
显式
△x
C.内能增量△u：
u c
x 2
A t1
k
k 1
t1 /
k
△x/2
k hx
据热平衡A+B=C并整理得：
k f
t 2 t1
k
t
t1
k
1 2
c
x
2
t1
k 1
LP
△y
t i 1 , j t i , j x
t i , j 1 t i , j y
y 2
x 2
1
BP
1
x 2
y 2
EP h t f t i , j
△x
1
FP h t f t i , j
t x
t
2
2
x i , j 2!
2
t x
3
x i , j 3!
3
3.一阶导级的中心差分表达式：<1>-<2>式且忽略后项，则有：
i , j
t x
传热与流动的数值计算
1.2 传热与流动问题数值计算的基本思想及应用举例
1.2.1 数值解基本思想（基于连续介质假设）
把原来在空间与时间坐标中连续的物理量的场 （如速度场、温度场、浓度场等），用一系列有限 个离散点（称为节点，node)上的值的集合来代替； 通过一定的原则建立起这些离散点上变量值之间关 系的代数方程(称为离散方程，discretization equation)；求解所建立起来的代数方程以获得所求 解变量的近似值。
u v w 0 x y z
div( U ) 0 t
称为流动无散（度）条件 (Zero divergence)。
2. 动量守恒方程
对上图所示的微元体分别在三个坐标方向上应用 Newton第2定律（F=ma)在流体中的表现形式： [微元体内动量的增加率]=[作用在微元体上各种力之和] 假设流体中切应力与正应力满足Stokes假定：应 力与应变成线性关系，可得u-动量方程如下：
为流体的动力粘度 ， 称为流体的第2分子粘度。
上式右端部分可进一步转化：
v u p u u w (divU 2 ) [ ( )] [ ( )] Fx x x y x y z z x x
u u u u v w ( ) ( ) ( ) ( ) ( ) ( ) (divU ) x x y y z z x x y x z x x p Fx u u u x div( gradu ) Su grad (u ) i j k x y z
Elliptic
的函数。 椭圆型 （回流型） 抛物型 （边界层）
0,
b 4ac
2
0, 0,
Parabolic
导热量计算公式
导热量计算公式导热量是物质传热过程中的一个重要参数，通常用于评估材料的传热性能。
它能够指示材料的导热能力以及导热系数大小。
导热量的计算公式如下：导热量=导热系数×截面积×温差/导热长度其中，导热系数是材料的一个固有属性，截面积为传热过程中与热交换有关的面积，导热长度为传热过程中热量从一个介质传到另一个介质所需的长度。
下面详细介绍各个参数的意义和计算方法。
1. 导热系数导热系数（热导率）是物质的一种热物理性质，反映了物体的导热能力大小。
导热系数是由于物质分子内部热运动而引起的热量传递引起的，其单位是(W/(m·K))，即每秒每米温度差为1K时所传递的热量。
在实际中我们可以通过测定物质的热传导速率来获取它的导热系数。
2. 截面积截面积是指传热过程中与热交换有关的面积，单位通常是平方米（m²）。
它是传热过程中热量传递的一个重要参考参数。
物体传热过程中，截面积的大小与物体的形状和大小、传热介质和传热方法有关。
3. 温差温差是指传热过程中两个物体之间的温度差值，通常表示为△T，其单位是开尔文（K）或摄氏度（℃）。
在热传导过程中，高温区域的热量会向低温区域传递，导热过程的速率与温差成正比关系。
4. 导热长度导热长度是指热量从一个介质传到另一个介质所需的长度，通常表示为L，其单位通常是米（m）。
在热传导过程中，物体热量的传递距离越长，传热速率就会变慢。
通过以上四个参数的计算，可以得到物体导热量的大小。
导热量是热传导过程中热量传递的一个重要参考值，可以用于评估材料的导热能力和传热性能。
总结：公式：导热量=导热系数×截面积×温差/导热长度参数意义：1. 导热系数：反映了物体的导热能力大小；2. 截面积：传热过程中与热交换有关的面积；3. 温差：传热过程中两个物体之间的温度差值，与传热速率成正比；4. 导热长度：热量从一个介质传到另一个介质所需的长度，传热距离越长，传热速率越慢。
传热流体数值计算
1 傅立叶定律傅立叶定律是导热理论的基础。
其向量表达式为:q g r a d T λ=-⋅ (2-1)式中:q —热流密度，是向量，2/()Kcal m h ；gradT —温度梯度，是向量，℃/m ；λ—导热系数，又称热导率，/()Kcal mh C ； 式中的负号表示q 的方向始终与gradT 相反。
2 导热系数（thermal conductivity ）及其影响因素导热系数λ（/()Kcal mh C ）是一个比例常数，在数值上等于每小时每平方米面积上，当物体内温度梯度为1℃/m 时的导热量。
导热系数是指在稳定传热条件下，1m 厚的材料，两侧表面的温差为1度（K ，°C ）,在1秒内，通过1平方米面积传递的热量，用λ表示，单位为瓦/米·度，w/m·k （W/m·K,此处的K 可用℃代替）。
导热系数为温度梯度1℃/m ，单位时间通过每平方米等温面的热传导热流量。
单位是：W/(m·K)。
3．热传导微分方程推导 ♥ 在t 时刻w 界面的温度梯度为xT∂∂在t 时刻e 界面的温度梯度为dx x T x T dx x x Tx T 22∂∂+∂∂=∂∂∂∂+∂∂ 单位时间内六面体在x 方向流入的热流量为：dydz xT∂∂-λ； 单位时间内六面体在x 方向流出的热流量为：dydz dx x T x T ⎥⎦⎤⎢⎣⎡∂∂+∂∂-22λ；单位时间内六面体在x 方向流入的净热量为：dxdydz xT22∂∂λ 图3-1 微分单元体各面上进出流量示意图同理，单位时间内六面体在y 方向流入的净热量为：dxdydz yT22∂∂λ； 单位时间内六面体在y 方向流入的净热量为：dxdydz z T 22∂∂λ； 单位时间内流入六面体的总热量为：dxdydz z T y T xT ⎥⎦⎤⎢⎣⎡∂∂+∂∂+∂∂222222λ (3-1) 六面体内介质的质量为:dxdydz ρ。
热传导和导热系数的计算方法
热传导和导热系数的计算方法热传导是指热量在物体内部由高温区向低温区传递的过程，其本质是物体内部粒子（如电子、原子、分子）的振动和碰撞引起的能量传递。
热传导的计算方法主要包括傅里叶定律、导热系数的概念及其计算方法。
1.傅里叶定律傅里叶定律是热传导的基本定律，表述为：物体内部的热流密度q与温度梯度dT/dx之间存在以下关系：[ q = -k ]其中，q表示热流密度，单位为瓦特每平方米（W/m^2）；k表示导热系数，单位为瓦特每米·开尔文（W/m·K）；dT/dx表示温度梯度，单位为开尔文每米（K/m）。
2.导热系数导热系数是描述材料导热性能的一个物理量，定义为：在稳态热传导条件下，1米厚的物体，在两侧表面温差为1开尔文时，单位时间内通过单位面积的热量。
导热系数用符号k表示，其单位为瓦特每米·开尔文（W/m·K）。
导热系数的计算方法主要有：（1）实验测定：通过实验方法，如热线法、热板法等，测定材料的导热系数。
（2）理论计算：根据材料的微观结构和组成，运用热力学和物理学原理，计算导热系数。
例如，对于均匀多晶材料，导热系数可通过以下公式计算：[ k = ( k_1 + k_2 + k_3 ) ]其中，k1、k2、k3分别为材料三个方向上的导热系数。
3.热传导的计算方法热传导的计算方法主要包括以下步骤：（1）建立热传导模型：根据实际问题，假设物体为均匀、各向同性或各向异性，简化模型以便于计算。
（2）确定边界条件和初始条件：如物体表面的温度、热流密度等。
（3）选择合适的数学方法求解：如有限差分法、有限元法、解析法等。
（4）分析结果：根据计算得到的温度分布、热流密度等，分析问题的热传导特性。
总之，热传导和导热系数的计算方法是热力学和物理学中的重要知识点，掌握这些方法有助于我们更好地理解和解决实际中的热传导问题。
习题及方法：1.习题：一长方体铜块的尺寸为2m×1m×0.5m，左表面温度为100℃，右表面温度为0℃。
传热学-导热数值计算
对物理问题进行数值解法的基本思路可以概括 为：把原来在时间、空间坐标系中连续的物理量的 场，如导热物体的温度场等，用有限个离散点上的 值的集合来代替，通过求解按一定方法建立起来的 关于这些值的代数方程，来获得离散点上被求物理 量的值，该方法称为数值解法。 这些离散点上被求物理量值的集合称为该物理 量的数值解。
1 (tm 1,n tm 1,n tm ,n 1 tm,n 1 ) 4
一阶
4.2.2 控制容积平衡法(热平衡法)
基本思想：对每个有限大小的控制容积应用能量守
恒，从而获得温度场的代数方程组，它从基本物理
现象和基本定律出发，不必事先建立控制方程，依 据能量守恒和Fourier导热定律即可。 能量守恒： 流入控制体的总热流量＋控制体内热源生成热 ＝ 流出控制体的总热流量＋控制体内能的增量
x y
tm , n
y
2x 2 qw x 1 3x 2 (2tm 1,n 2tm ,n 1 tm ,n 1 tm 1,n ) 6 2
讨论关于边界热流密度的三种情况： （1）绝热边界
即令上式 qw 0 即可。
（2）qw 值不为零
流入元体，qw 取正，流出元体，qw 取负 （3）对流边界 此时 qw h(t f tm,n ) ，将此表达式代入上述方程，并 将此项中的
tm , n 1 (tm 1,n tm 1,n tm ,n 1 tm ,n 1 ) 4
（4） 设立迭代初场
代数方程组的求解方法有直接解法与迭
代解法，传热问题的有限差分法中主要采用
迭代法。采用迭代法求解时，需对被求的温
度场预先设定一个解，这个解称为初场，并
在求解过程中不断改进。
4.1.2 物理问题的数值求解过程
第四章-传热学数值计算方法资料
达到给定精度所需要的网格点数,以及这些网格点在计算域内应 采取的分布方式与所求问题的特性有关。
4. 采用仅几个网格点进行试探性计算，为弄清有关解 的情况提供了一个方便的途径。 也可来指导实验。
2020-11-2
太原理工大学
8 /72
§4.2-3 界面导热系数
1. 问题的提出
Thermal
通用离散方程式 aPTP aETE aWTW b
§4.2-1 基本方程
一维稳态导热问题的控制方程：
其中：S SC SPTP
d dx
k
dT dx
S
0
相应的离散化方程： aPTP aETE aWTW b
式中：
aE
ke
xe
aW
kw
xw
aP aE aW SPx b SCx
2020-11-2
太原理工大学
5 /72
Thermal
分布假设： dT 由T 对 x 的分段线性的变化算得；
2020-11-2
太原理工大学
3 /72
Thermal
④. 把热传导用作流体流动计算方案的基本组成部分 的做法有助于理解动量传递与热量传递之间的类 似性（用某种方法把速度与温度相比拟）。
本章内容将是流动与换热数值解的基础
2020-11-2
太原理工大学
4 /72
§4.2 一维稳态热传导
Thermal
式中：
aE
ke
xe
aW
kw
xw
aE、aW分别是节点E 与 P 和节点W 与P 间的热导，热导的大小
反映了周围节点对节点P的影响程度。系数aE、aW中分别含有交
界面导热系数ke与kw。当k 是x 的函数时，只知道kP 、kE、 kW，
导热系数的计算公式
导热系数的计算公式导热系数是物质传导热量的性质参数，它反映了物质对热量的传导能力。
在工程领域中，准确计算导热系数对于材料的选择和热传导问题的分析都具有重要的意义。
本文将介绍导热系数的计算公式及其相关内容。
导热系数的计算公式如下：λ = (Q × L) / (A × ΔT)其中，λ表示导热系数，Q表示热量，L表示传热长度，A表示传热面积，ΔT表示温差。
导热系数的计算公式是基于热传导定律推导出来的。
热传导定律指出，热量的传导速率正比于传热面积，与温差成正比，与传热长度成反比。
因此，当我们知道热量、传热面积、传热长度和温差时，就可以通过导热系数的计算公式计算出导热系数。
在实际应用中，我们常常需要根据材料的特性参数来计算导热系数。
各种材料的导热系数不同，这取决于材料的热传导性能。
一般来说，金属材料的导热系数较高，而绝缘材料的导热系数较低。
根据导热系数的大小，我们可以判断材料的热传导性能。
导热系数的计算公式还可以用于解决一些热传导问题。
例如，我们可以利用导热系数计算材料的热传导速率，从而了解材料的保温性能。
另外，通过导热系数的计算公式，我们还可以计算材料的热阻，用于评估材料的隔热性能。
需要注意的是，导热系数的计算公式是在一定条件下成立的。
在实际应用中，我们需要根据具体的情况来选择合适的计算公式。
例如，如果材料的导热系数随温度变化较大，我们需要考虑温度对导热系数的影响。
导热系数是物质传导热量的重要参数，它反映了物质对热量的传导能力。
通过导热系数的计算公式，我们可以计算出导热系数，并应用于热传导问题的分析和材料的选择。
在实际应用中，我们需要根据具体情况选择合适的计算公式，并注意计算条件对导热系数的影响。
通过合理应用导热系数的计算公式，我们可以更好地理解和应用导热系数的概念，为工程实践提供有益的参考。