偏积分方法与数值计算方法比较分析∫.pdfVIP

偏积分方法与数值计算方法比较分析 1 2 1 1 1 张敏 ，John C. Chai ，刘少培 ，张钧波 ，许彬 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天学院，新加坡 (639798) E-mail ：mz2455@163.com 摘 要：在介绍偏积分近似方法基本概念的基础上，在结构化网格中，对稳态热传导问题， 同时进行数值解、精确解和偏积分近似解的计算。通过比较计算结果，得到了令人满足的一 致性。这些充分证明偏积分近似方法的可靠性和实用性。 关键词：偏积分近似方法，数值方法，精确解 分离变量法是解析求解各种热传导问题最有效工具之一；而伽略金法又是近似分析或数 值分析求解这类问题几乎万能的手段之一。它们两者的结合，又可产生一种近似求解方法， 这就是偏积分法。 当一个多维热传导问题，在维与维之间可以分离时，偏积分方法在非齐次方向或非线性 方向，首先选择伽略金近似；然后进行积分。这样我们就可得到另一维方向变量的常微分方 程，进行求解最后得到偏积分方法的近似解[1-5] 。 在前面的文章中，已用伽略金近似法求解了二维稳态热传导问题，在那里所采用的试探 解T x, y 是由变量为 x 与 y 的函数形式业已确定的函数φ x , y 所构成。如果试探解 ( ) ( ) T x, y 的构造，对一个变量，比如说y，其剖面的函数形式是确定的，而对另一个变量， ( ) 其函数形式是待定的，比如说，变量x 的函数X x 是待定的，那么，仍可望得到一种较好 ( ) 的近似解。总之，对于变量为x , y 的二维问题，在变量y 方面应用伽略金法，其结果将得到 一个函数X x 的常微分方程。这个方法的一个优点是，在温度剖面的函数形式有一个方向 ( ) 在预先不能被选定的情况下，可根据问题的特点，从由此得到的常微分方程的解来确定。 1. 偏积分法的基本理论 [1] 在有限矩形区域内，稳态热传导问题具有如下一般的形式 ， 2 2 ∂ T ∂ T g + + 0 , 0 =x a , 0 y f x (1.1a) 2 2 ( ) ∂x ∂y k T 0 x 0, x a, y f x (1.1b) ( ) ∂T 0 y 0 (1.1c) y ∂ 将伽略金近似法应用于微分方程(1.1)，并对变量y 作偏积分可得， 2 2 f x