第4章 扩散方程的数值解法及其应用.pptVIP

传热与流体流动的数值计算*第1页，共30页，星期日，2025年，2月5日§4.1一维导热问题一、一维导热问题的通用控制方程坐标系空间变量x面积因子A(x)直角x1圆柱rr球rr2*第2页，共30页，星期日，2025年，2月5日一、一维导热问题的通用控制方程（续）*第3页，共30页，星期日，2025年，2月5日二、用控制容积积分法导出离散形式假设：（源项线性化）*第4页，共30页，星期日，2025年，2月5日二、用控制容积积分法导出离散形式（续）分段线性插值（中心差分）*第5页，共30页，星期日，2025年，2月5日三、界面上当量导热系数的确定方法根据界面上热流密度连续的原则：根据界面上当量导热系数的含义：调和平均法：*第6页，共30页，星期日，2025年，2月5日三、界面上当量导热系数的确定方法（续）算术平均法：调和平均法优于算术平均法。导热系数发生阶跃变化时的处理方法：(1)把阶跃面作为控制容积的分界面。(2)把阶跃面设置成一个节点的位置。（更准确）*第7页，共30页，星期日，2025年，2月5日四、一维非稳态导热方程及其离散化采用全隐格式离散：*第8页，共30页，星期日，2025年，2月5日四、一维非稳态导热方程及其离散化（续）令：得：其中：*第9页，共30页，星期日，2025年，2月5日§4.2多维非稳态导热问题的全隐格式一、三种正交坐标系中的全隐离散方程(1)直角坐标系非稳态项：*第10页，共30页，星期日，2025年，2月5日一、三种正交坐标系中的全隐离散方程（续）源项：扩散项：*第11页，共30页，星期日，2025年，2月5日一、三种正交坐标系中的全隐离散方程（续）整理得：其中：*第12页，共30页，星期日，2025年，2月5日一、三种正交坐标系中的全隐离散方程（续）(2)圆柱轴对称坐标系*第13页，共30页，星期日，2025年，2月5日一、三种正交坐标系中的全隐离散方程（续）(3)极坐标系*第14页，共30页，星期日，2025年，2月5日二、三种坐标系中全隐格式的通用离散形式三种坐标系的区别：x-yx-rq-r东西坐标xxq引入尺度系数SX南北坐标yrr引入名义半径R通用离散形式：*第15页，共30页，星期日，2025年，2月5日二、三种坐标系中全隐格式的通用离散形式（续）*第16页，共30页，星期日，2025年，2月5日二、三种坐标系中全隐格式的通用离散形式（续）*第17页，共30页，星期日，2025年，2月5日§4.3源项及边界条件的处理一、非常数源项的线性化处理说明：1、当源项为未知量的函数时，线性化处理比假定为常数更合理。2、线性化处理是建立代数方程所必须的。3、为了保证代数方程迭代求解收敛，要求SP0。4、SP绝对值的大小影响到迭代过程中温度的变化速度。*第18页，共30页，星期日，2025年，2月5日一、非常数源项的线性化处理（续）-举例*第19页，共30页，星期日，2025年，2月5日二、边界条件的处理1、补充边界节点代数方程的方法第2类BC：内节点法Dx=0，得第3类BC：得：由控制容积平衡法导出离散方程物理意义明确。对于不规则区域，可用区域扩充法使之成为规则区域。*第20页，共30页，星期日，2025年，2月5日二、边界条件的处理（续）*第21页，共30页，星期日，2025年，2月5日传热与流体流动的数值计算