4-1 导热问题数值解法(课件1).pptVIP

第四章 导热问题的数值解法;; 分析法：能获得所研究问题的精确解，可以为实验和数值计算提供比较依据；局限性很大，对复杂的问题无法求解；分析解具有普遍性，各种情况的影响清晰可见。;分析解法与数值解法的异同点： 相同点：根本目的相同，即确定： ① t=f( x, y, z,τ) ； ②热流量。 不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。 ; 数值求解的步骤 ;二维矩形域内稳态无内热源，常物性的导热问题;基本概念：网格线、步长、节点(内、外)、界面线、控制容积;三、建立节点物理量的代数方程（离散方程） ;四、设立迭代初场;五、求解代数方程组 ;2 ）非线性代数方程组：代数方程一经建立，其中各项系数在整个求解过程中不断更新。 3 ）是否收敛判断：是指用迭代法求解代数方程是否收敛，即本次迭代计算所得之解与上一次迭代计算所得之解的偏差是否小于允许值。 ;六、解的分析 ;一、建立控制方程及定解条件 ;二、区域离散化（确立节点）;三、建立节点物理量的代数方程（离散方程） ;三、建立节点物理量的代数方程（离散方程） ;中心差分;三、建立节点物理量的代数方程（离散方程） ;三、建立节点物理量的代数方程（离散方程） ;三、建立节点物理量的代数方程（离散方程） ;三、建立节点物理量的代数方程（离散方程） ;三、建立节点物理量的代数方程（离散方程） ;（3）边界节点的有限差分方程;从所有方向流入控制体的总热量＋ 控制体内热源生成热＝ 0;从所有方向流入控制体的总热量＋ 控制体内热源生成热＝ 0;从所有方向流入控制体的总热量＋ 控制体内热源生成热＝ 0;边界节点离散方程中的两个问题：;建立节点离散方程的泰勒级数法与热平衡法的比较： 泰勒级数法属于纯数学方法，而热平衡法基于能量守恒原理，物理概念明确，且推导过程简捷； 泰勒级数法对于建立边界节点的离散方程较困难； 当导热物体物性或内热源不均匀时，泰勒级数法不适用，而热平衡法能够方便处理。;例1： 对如图所示的圆截面直肋的一维稳态、无内热源、常物性导热问题，试分别列出内节点m和端部节点M的离散方程式。已知圆截面直径为d。;四、设立迭代初场;§4-2 稳态导热问题的数值解法;§4-2 稳态导热问题的数值解法;判断迭代是否收敛的准则：;§4-2 稳态导热问题的数值解法;例4-1：利用高斯-赛德尔迭代法求解下列方程：;假设t1、t2、t3的初始值均取零，迭代值为：;如此经过七次迭代后，在四位有效数字内得到了与精确解一致的结果。;假设t1、t2、t3的初始值均取零，迭代值为（经三次迭代）：;憋泰己咖茎蹿蹬肄组撬纲汁怂秸宛氖凿鉴殃顾淘文肤澜暗剩砌扇艺礼疹揭4-1 导热问题的数值解法(课件1)4-1 导热问题的数值解法(课件1);例2 ：某方形物体，导热系数为常数，已知各边界温度如图所示，试用高斯-塞德尔迭代法求其内部节点1、2、3、4点的温度。;解：1）列节点方程（内节点） 1： 2： 3： 4：;迭代次数;§4-2 稳态导热问题的数值解法;例3：如图所示，一等截面之类，高H=45mm，厚δ=10mm，肋根温度t0=100℃，流体温度tf=20℃，表面传热系数h=25W/(m2·K),肋片导热系数λ=50W/(m·K),设肋端绝热。将该肋片等分成4个节点。试列出节点2,3,4的离散方程式，并计算其温度。;解：这是一个一维稳态无内热源、常物性的导热问题。利用热平衡法列节点的离散方程。;利用迭代法解得：;怖规切挞魄竹掂佑惰揭雀挚示嗜掐析燎匪柏域院扔良特袱灭淀熟绝澈鸭轮4-1 导热问题的数值解法(课件1)4-1 导热问题的数值解法(课件1);瞩头泄绍翁囚酿找柳糖春藤楷谁拽闺离胎绦开崭厉适绍抬尤潭裔懂意葡梳4-1 导热问题的数值解法(课件1)4-1 导热问题的数值解法(课件1); 4-3 非稳态导热问题的数值解法;tf h;二、区域离散化（确立节点）;向前差分格式;控制方程离散化：;其中： ;从所有方向流入控制体的总热量 ＝ 控制体内能的增量;内节点 n;(2) 热平衡法; 4-3 非稳态导热问题的数值解法;显示格式存在稳定性问题：如果节点 tn(i) 前面的系数小于零，则数值解出现不稳定的震荡结果。;隐式格式;2．讨论一维导热问题显式差分格式稳定性限制的物理意义 由节点显式差分方程可以看出， 前的系数可能是负值。当系数为负时， 时刻节点温度越高 ， 时刻的温度却越低，违背了热力学第二定律。;表面温度在不同 Δτ 时的计算值;骗襟哺琴允乞悉尿埃骗盯服陀宦酞挎阂灰步秉灯摆与彝船酷甩驮至笛侄翼4-1 导热问题的数值解法(课件1)4-1 导热问题的数值解法(课件1); 4-3 非稳态导热问题的数值解法