第4章_热传导问题数值解法.ppt

第4章 热传导问题的数值解法 4.1 导热问题数值求解的基本思想应用 4.1.1 基本思想 把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。 这些离散点上被求物理量的集合称为该物理量的数值解。 这一基本思想可用图4-1所示的求解过程的框图表示。 * * 第4章 热传导问题的数值解法 4.1 导热问题数值求解的基本思想应用 4.2 内节点离散方程的建立方法 4.3 边界节点离散方程的建立及代数方程的求解 4.4 非稳态导热问题的数值解法 4.1.2 导热问题数值解求解基本步骤 以图4-2a所示的二维矩形域内的稳态、无内热源、常物性的导热稳态为例。 （1）建立控制方程及定解条件。 控制方程－－描写物理问题的微分方程 边界条件 （2）区域离散化 如图4-2b所示，用一系列与坐标轴平行的网格线将求解区域划分成许多子区域，以网格线的交点作为需要确定温度值的空间位置，称为节点。节点的位置以该点在两个方向的坐标号m、n表示。 相邻节点间的距离称为步长，记为?x和?y。 每个节点都可以看成是以它为中心的一个小区域的代表，如图4-2b中的有阴影线的小区域即是节点（m，n）所代表的区域，它由相邻两节点连线的中垂线构成。 控制容积（元体）——节点所代表的区域。 （3）建立节点物理量的代数方程 离散方程——节点上物理量的代数方程，其建立是数值求解过程中的重要环节。当?x和?y时，对于节点（m，n），其代数方程为 （4）设立迭代初场 代数方程的求解方法：直接解法和迭代法。 在传热问题的有限差分解法中主要采用迭代法。 初场——对被求解的温度场预先假定的一个解，在求解过程中不断得到改进。 （5）求解代数方程组 借助于计算机。 线性问题－－常物性、无内热源 非线性问题－－变物性或有内热源 （6）解的分析 通过求得的温度场进一步计算热流密度等。 4.2 内节点离散方程的建立方法 建立内节点离散方程的方法有： 泰勒级数展开法 热平衡法 考察图4-3。 4.2.1 泰勒级数展开法 以节点（m，n）处的二阶偏导数为例，对节点（m+1，n)及（m-1，n）分别写出函数t对（m，n）点的泰勒级数展开式： 将式（a）、（b）相加得 将式（c）改写成如下形式： 截断误差 略去截断误差有 将上式代入到前述的控制方程，得节点（m，n）的离散方程为 表4-1 一阶、二阶导数的常用差分表达式 中心差分 4.2.2 热平衡法 采用热平衡法时，对每个节点所代表的元体用傅里叶定律直接写出其能量守恒的表达式。 如：节点（m-1，n）通过界面w传导到节点（m，n）的热流量为 对于所研究的问题，元体（m，n）的能量守恒方程为 节点（m，n）的离散方程为 4.3 边界节点离散方程的建立及代数方程的求解 4.3.1 边界节点离散方程的建立 利用热平衡法。 （1）位于平直边界上的节点 如图4-4所示。节点（m，n）的离散方程为 （2）外部角点 如图4-5所示。节点（m，n）的离散方程为 （3）内部角点 如图4-5所示。节点（m，n）的离散方程为 4.4 非稳态导热问题的数值解法 4.4.1 时间－空间区域的离散化 以一维问题为例，如图4-8所示。 x代表空间坐标 ?代表时间坐标 时间步长??——从一个时间层到下一个时间层的间隔。 空间网格线和时间网格线的交点，如（n,i），代表了时间－空 间区域中的一个节点的位置相应的温度记为 将函数t在节点（n,i+1）对节点（n,i）作泰勒展开，可得 4.4.2 一维平板非稳态导热的显式格式 扩散项取中心差分，非稳态项取向前差分，有 求解非稳态导热方程就是从已知的初始温度分布出发，根据边界条件以次求得以后各个时间层上的温度值，由上式可知，一旦i时层上各节点的温度已知，可立即求得i+1时层上各节点的温度，而不必联立方程，因而上式所代表的格式称为显式差分格式。