导热问题的数值解法.pptxVIP

1 、重点内容： ① 掌握导热问题数值解法的基本思路； ② 利用热平衡法和泰勒级数展开法建立节点的离散方程。 2 、掌握内容：数值解法的实质。 3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。 ; 求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，这些数值解法主要有以下几种： ;分析解法与数值解法的异同点： ?? 相同点：根本目的是相同的，即确定 ① t=f(x ， y ， z) ； ② 。 ?? 不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。 ;??数值解法的实质 ???? 对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。该方法称为数值解法。 ???? 这些离散点上被求物理量值的集合称为该物理量的数值解。 ;§4-1 导热问题数值求解的基本思想 及内部节点离散方程的建立;二维矩形域内稳态无内热源，常物性的导热问题;; 如图（a）所示二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下： ;（2）区域离散化（确立节点） ; （3）建立节点物理量的代数方程（离散方程） ;（4） 设立迭代初场 ???? 代数方程组的求解方法有直接解法与迭代解法，传热问题的有限差分法中主要采用迭代法。采用迭代法求解时，需对被求的温度场预先设定一个解，这个解称为初场，并在求解过程中不断改进。 ;（5） 求解代数方程组 ;2 ）非线性代数方程组：代数方程一经建立，其中各项系数 在整个求解过程中不断更新。 3 ）是否收敛判断：是指用迭代法求解代数方程是否收敛，即本次迭代计算所得之解与上一次迭代计算所得之解的偏差是否小于允许值。 ;（6） 解的分析 ;4 建立离散方程的常用方法：;(1) 泰勒级数展开法;将上两式相加可得;?根据导热问题的控制方程 ( 导热微分方程 );(2) 控制容积平衡法(热平衡法);对元体 (m,n). 根据能量守恒定律可知： ;说明： ① 上述分析与推导是在笛卡儿坐标系中进行的； ② 热平衡法概念清晰，过程简捷； ③ 热平衡法与建立微分方程的思路与过程一致，但不同的是前者是有限大小的元体，后者是微元体。 ;4-2 边界节点离散方程的建立 及代数方程的求解;如图所示 边界节点 (m,n) 只能代表半个元体，若边界上有向该元体传递的热流密度为 ，据能量守恒定律对该元体有： ;;(2) 外部角点;;(3) 内部角点;;讨论关于边界热流密度的三种情况： ;（a）平直边界;2 代数方程的求解方法 ;2 迭代法目前应用较多的是： ;设有一三元方程组： ;采用高斯——赛德尔迭代法的步骤： （1）将三元方程变形为迭式方程： ;;判断迭代是否收敛的准则：;说明： 1 ）对于一个代数方程组，若选用的迭代方式不合适，有可能导致发散，即称迭代过程发散； ;3 ）采用热平衡法导出差分方程时，若每一个方程都选用导出该方程中心节点的温度作为迭代变量，则上述条件必满足，迭代一定收敛。