热传导中积分近似法与数值法的比较分析.pdfVIP

热传导中积分近似法与数值法的比较分析 1 1, 2 1 2 1 刘少培 ，许彬 ，张敏 ，John C. Chai ，张钧波 1 南京理工大学动力工程学院，南京(210094) 2 南洋理工大学机械与航天学院，新加坡(639798) E-mail ：mz2455@163.com 摘 要：在介绍积分近似方法基本概念和步骤的基础上，在结构化网格中，对非稳态热传导 问题，同时进行数值解、精确解和积分近似解的计算。通过比较计算结果，三者之间得到令 人满足的一致性。这些充分证明积分近似方法在非稳态热传导问题求解中的实用性和可靠 性。 关键词：积分近似方法，数值方法，精确解析方法 积分近似法在求解偏微分方程中的应用可追溯到冯卡门(VonKarman) 与波尔豪森 (Pohlhausen)，他们用近似分析方法求解了流体力学中的动量边界层方程与能量边界层方程。 兰代尔(Landahl)把它用于生物物理领域求解分离浓缩物的扩散方程。默克(Merk)用这种方法 求解了二维稳态的熔化问题。从此以后，这种方法就应用于求解各类一维瞬态热传导问题， 熔化与凝固问题以及诸如海水中的冰在熔化过程中，以及聚合物在熔化及热压过程中的热量 传输和动量传输问题[1-3] 。 对于一定边界条件下的一维瞬态热传导边值问题，无论是线性的或非线性的，使用积分 近似法都是很简捷的，而且十分方便。尽管用积分法所得的结果是近似的，但当把由此得到 的很多近似结果与精确解的结果作比较时，即可发现，就工程应用观点来看，其准确度一般 已经达到令人满意的程度。本文先介绍一些基本概念，这些基本概念对于把积分法应用于求 解简单的半无限大物体瞬态热传导问题是有用的[4-5] 。 1. 积分近似法的基本概念及步骤 精确求解在特定边界条件及初始条件下，某一区域内的热传导微分方程，由此得到的解 对该区域内每一点都应满足，而用积分方法得到的解对该区域只是平均地满足。下面，我们 应用积分法求解具有特定边界条件、均匀初始条件和无热源的半无限大物体中一维的瞬态热 传导问题，概括地叙述一下用积分法分析问题的基本步骤。 （a ）将热传导微分方程，对称为热层的一表观厚度δ(t) 进行积分，可把偏微分方程中 有关空间变量的导数去掉。对热层厚度作如下定义：若从实际应用角度来看，超过某一厚度 就不再存在热流，则将此厚度定义为热层。因此，超过δ(t) ，初始温度分布就不再受影响。 由此得到的方程称为能量积分方程（也称热平衡积分）。 （b ）选某一合适的剖面作为热层内的温度分布。通常选某一多项式为剖面。经验已经 证明：所选的多项式高于四次后，解的精度不再有明显的改进。多项式中的系数可根据实 际的（有时还需补充几个推导得到）边界条件予以确定，并用热层厚度δ(t) 来表示。 （c ）把所得到的温度剖面代入能量积分方程，在进行简要的运算之后，即可得到关于 热层厚度δ(t) 以时间为自变量的常微分方程。这个微分方程的解，满足特定的初始条件， 并给出δ(t) ，它是时间的函数。 （d ）从第三步得知δ(t) 后，即可知道温度分布，它是时间与物体内位置的函数。 - 1 - 2. 算例分析 一半无限大物体的瞬态热传导问题，初始时该半无限大物体有均匀温度T ，在时间 1 [1] t 0 时，边界面维持为恒温T 。该问题的数学描述为 ： 0