第四章热传导热问题的数值解法.pptVIP

▲数值求解的基本思想及常用的数值求解方法 ▲有限差分法 ▲节点离散方程的建立 ——泰勒级数展开法与热平衡法。 ▲节点离散方程(组)的求解 1、直接求解; 2、简接求解——高斯-赛德尔(Gauss- Seidel)迭代法* ▲非稳态导热问题数值求解的有关概念 常用的数值求解方法 4.2 内节点离散方程的建立方法 4.2.1 泰勒级数展开法 将函数t 的一阶导数及二阶导数在点(m,n)按泰勒 级数展开，并取其一阶近似，结果代入相应的导热微 分方程，即可得出该节点的离散方程。(详略) 4.3.3 代数方程的求解方法---迭代法 雅可比迭代法(简单迭代); 高斯—赛德尔(Gauss - Seidel ) 迭代法*; 超松弛迭代法; 块迭代；交替方向迭代等 1.高斯-赛德尔(Gauss- Seidel)迭代法 迭代原则 基本求解步骤 迭代过程已经收敛的判据 (即控制过程终止的方法) 迭代过程能否收敛的判据 (获得收敛解的条件) 其中，ε’ =10-3～10-6 迭代结果 （2）迭代结果 4-3 非稳态导热问题的数值解 非稳态导热与稳态导热问题的主要差别： ≠0 本节主要介绍： 非稳态项的离散（三种差分格式） ——向前差分、向后差分与中心差分 差分方程的两种格式 ——显式差分格式与隐式差分格式 显式格式的稳定性问题 一、非稳态项 离散的三种差分格式 ——向前差分、向后差分与中心差分 △τ n,i+1 n.,i n+1,i n-1,i n,i-1 △x n=1 n N i=1 i I X △τ：时间步长 △x ：几何步长 以一维非稳态导热问题为例： 1、向前差分： 2、向后差分: 3、中心差分: (4-11) (4-12) (4-13) (1)显式差分格式——扩散项用i 时层的值表示： ——内节点差分方程（显式） 优点：计算工作量小(不必解联立方程, 可直接求解) 缺点：对时间步长及空间步长有一定限制(稳定性条件) 二. 差分方程 (4-14a) 改写为： (4-14b) (2)隐式差分格式——扩散项用i+1时层的值表示： ——内节点差分方程（隐式） 优点：对步长无限制，所得的解稳定 缺点：计算工作量大—— 不能根据i时层各节点的温度直接由方程求 出(i+1)时层上各内点的温度值tn(i+1)，而必 须求解 (i+1)时层的一个联立方程。 (4-15) 三、步长△x、 △τ对显式格式稳定性的影响 ——稳定性限制条件 (1)稳定性限制的物理意义 i时刻点n的温度越高，则其相继时刻的 温度也较高；反之, i时刻点n的温度越低, 则 其相继时刻的温度也较低。 (2)稳定性限制条件 相应节点的显式差分方程中, 项前的 系数大于等于零 例题4-3：大平板对称冷却 * * 第四章 主要内容 重点：用热平衡法建立稳态导热问题的离散方程， 数值求解的高斯-赛德尔(Gauss- Seidel)迭代法 4.1.1 数值求解的基本思想(见P162)： 把原来在时间、空间坐标系中连续的 物理量的场，用有限个离散点上的值的集 合来代替，通过求解按一定方法建立起来 的关于这些值的代数方程（组），来获得 离散点上被求物理量的值（其集合称为该 物理量的数值解） 4-1 导热问题数值求解的基本思想 1、有限差分法* （finite-difference） ——对控制方程进行离散(近似) 2、有限元法（finite-element） ——对解进行离散(近似) 3、边界元法（boundary- element） 4、分子动力学模拟（MD） 4、解方程 , 并用节点的解的集合( 离散值 )来 代替原物体内的连续温度分布 有限差分法： 1、离散：将连续体用网格分割成有限单元体 2、取节点：以单元体的中心点代表该单元体 3、建立节点离散方程：对每一单元体按一定 方法，将针对微元体得出的导热微分方程 简转化成针对有限单元体的节点离散方程 ( 代数方程组 ) 基本概念： 单元体(控制容积)、网格线、节点、 界面线、步长 4.1.2 用有限差分法求解导热问题的基本步骤 是否收敛？ 是 否 建立控制方程及定解条件 确定节点（区