1、4.1 导热问题数值求解的基本思想 4.2 计算区域离散化与节点方程的建立 4.3 温度场内节点离散方程的建立 4.4 边界节点离散方程的建立 4.5 节点方程的求解 4.6 *一维非稳态导热问题数值解法,第4章 导热问题的数值解,4.1.1 数值计算求解思想把原来在时间和空间连续的物理量的场，用有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起这些离散点上变量值之间的关系的代数方程组，通过求解代数方程组以获得场变量的近似值。,4.1.2 物理问题的数值求解过程,4.1.3 数值传热学中常用数值方法有限差分法（finite difference method，FDM ）有限容积法
2、(Finite volume method，FVM )有限单元法（finite element method，FEM ）边界元法（boundary element method，BEM ）,节点,节点,节点,有限单元,边界单元,有限差分法,有限单元法,边界元法,有限容积法,基本思想：将求解域划分为差分网格，用有限个离散节点代替连续的求解域。,有限差分法,数学基础：在每个节点上，将控制方程中各阶导数（微商）用相应的差分表达式（差商）来代替，从而在每个节点上形成一个代数方程，每个方程中包括了本节点及其附近一些节点上的未知值，求解这些代数方程即可获得近似数值解。,几何意义：是用函数在某区域内的平均变
3、化率代替函数的真实变化率。,基本思想：将求解区域划分成一系列控制容积，每个控制容积都有一个节点做代表。,有限容积法,数学基础：通过将守恒型的控制方程对控制容积做积分来导出离散方程。 特点：离散方程系数物理意义明确，导出的离散方程可以保证具有守恒特性。,4.2 计算区域离散化与节点方程的建立,区域离散化（domain discritization ）：对空间连续的计算区域进行剖分，将其划分为许多个子区域，并确定每个区域的节点，此过程又称为网格生成。 实质： 用一组有限个离散的点来代替原来连续的空间,4.2.1 区域离散化的实质与内容,区域离散化后几何要素：,界面,节点,元体,网格线,4.2.2
4、节点方程的建立,Taylor展开法各阶导数的差分表达式由Taylor级数展开而得； 多项式拟合法 控制容积平衡法将物理上的守恒方程直接应用于所研究的控制容积。,4.3 温度场内节点离散方程的建立,【例】二维矩形域内稳态，常物性的导热问题。,控制方程,采用均匀网格,区域离散,4.3.1 Taylor展开法,用节点(m,n)的温度tm,n来表示节点(m+1,n)的温度tm+1,n,用节点(m,n)的温度tm,n来表示节点(m-1,n)的温度tm-1,n,由(1)得：,一阶差分格式,若略去截断误差，则有,截断误差 （truncation error）,向前差分 （forward diffence）,
5、一阶精度,由(2)得：,向后差分 （backward diffence）,由(1)-(2)得：,中心差分 （central diffence）,同理：,若将式(1)和式(2)相加：,二阶差分格式,中心差分 （central diffence）,同理可得：,导热微分方程的离散,差分方程：,对于传热问题，也称热平衡法Energy balance method 。 思路：节点有限差分方程式通过对围绕这个节点的控制体应用能量守恒获得。,4.3.2控制容积平衡法,流入控制体的总热流量,稳态,二维内部节点,注：w,e,n,s表控制体界面,节点越多，假设的分段线性分布越接近真实的温度布。此时：,型线选择,内
6、热源：,离散方程,【引申思考】所求节点的温度前的系数与相邻节点温度前的系数有定量关系吗？,由离散方程,可知，当x=y且无内热源时，,【重要结论】所求节点的温度前的系数一定等于其他所有相邻节点温度前的系数之和。,4.4.1平直边界上的节点,4.4 边界节点离散方程的建立,4.4.2 外部角点,4.4.3 内部角点,4.4.4边界热流密度不同情况分析,【引申思考】当x=y时，试推导出不同边界热流密度下的最终离散表达式。,4.5 节点方程的求解,写出所有内节点和边界节点的温度离散方程。,n个未知节点温度，n个代数方程式方程组封闭,代数方程组的求解方法：直接解法、迭代解法,特点：每次迭代时总是使用节点
7、温度的最新值,将方程组(a)改写成,显性格式,过程：,4.5.1 高斯-赛德尔迭代法,根据第 k（k=0,1,）次迭代的数值,可以求得节点第k+1次迭代温度：,取初值,结论： 一般采用相对偏差小于规定数值的判据更合理； 当有接近于零的t 时，第三种判 据更合理,4.5.2 迭代收敛判据,第k次迭代得到的最大值.,4.6* 一维非稳态导热问题数值解法,4.6.1时间域的离散,时间域与空间域不耦合，故,【例】常物性无内热源一维导热问题:,采用控制容积平衡法。,控制方程,无内热源:,型线选择：分段线性型线,4.6.2节点方程的建立,内节点,离散方程,进一步改写：,显式格式,引入无量纲参数：,则离散方程描述为：,网格傅里叶数,无内热源:,边界节点,第三类边界条件,则离散方程描述为:,显式格式,网格毕渥数,引入无量纲参数：,41,特点：所列出的方程不需要联立迭代求解。,步进式推进：由已知的初始场分布，通过节点方程直接求解下一时间步长下的场分布。以此类推，可得所有时间步长下的场分布。此种方式称之为步进式推进。,4.6.3 节点方程的求解,对于显式格式，常采用步进式推进法。,
