二维瞬态传热偏微分方程的求解.docxVIP

二维瞬态传热偏微分方程的求解 1 典型数值模拟方法的建立 热传导广泛应用于各种材料的加工中。在硅铁粉中的硬化过程可以被认为是一个没有均匀厚度的小板的传热过程。此外，轴向同构金属铸造的定向硬化过程可以简化为二维截面值模型，以实现传热模拟。要掌握这些成型过程的最优工艺参数,关键是正确求解关于传热问题的偏微分方程,典型的数值方法有有限元、有限差分和有限容积法。MATLAB是一款具有精确数值计算功能和丰富图形处理函数的软件,传热偏微分方程的数值解可直观地以二维、三维图形方式显示在屏幕,之前有一些文章介绍用MATLAB实现有限差分方法求解传热问题,本文主要讨论一、二维热传导方程有限容积法的MATLAB实现。 2 热值的建模和参数 2.1 单元网格节点值解 有限容积法是一种通过代数方程组代替原有难以直接求解的传热偏微分方程的数值计算方法,与有限差分法、有限元方法非常相似,都是求解离散单元网格节点值从而获得传热问题的解。如图1所示,对于二维问题,“有限体积单元”是指以节点P为中心,由东(e)、西(w)、南(s)、北(n)方向界面所形成的封闭面积(图1中的虚线框),相邻的体积单元通过各自界面传递热流量(图1中flux),某个有限体积单元的流入热量和与之相邻有限体积单元的流出热量相等,所以它是一种能量守恒的计算方法,因此有限容积法被广泛应用于热量、动量传递问题的数值求解。 2.2 直方向的空间坐标