从数值解浅析热传导方程精确解级数项选取.pdfVIP

从数值解浅析热传导方程精确解级数项的选取 1 2 2 黄庆宏 ，张敏 ，许彬 1 南京师范大学动力工程学院，江苏南京（210042 ） 2 南京理工大学动力工程学院，江苏南京（210094 ） 摘 要：用热传导方程的精确解检验数值求解过程的正确性是数值计算方法中一种重要的手 段。反之，用数值解来分析精确解的完整性也是一种有意义的方法。本文用基元有限容积法 和二阶迎风格式，在结构化网格中数值求解热传导方程，同时给出三个例子进行数值解和精 确解的比较，以此来说明某些教科书中古典导热方程精确求解过程的不足，并且分析阐述两 者之间的内在辨证关系和它们的互补性。 关键词：数值解，精确解，热传导方程 1. 引言 固体中的热传导问题是传热学领域中最古老、最经典和最成熟的分支之一。经过无数前 人的努力，对于导热方程在各种边界条件下的精确解形成了一套系统的求解方法，这使得它 [1-3] 们成为检验今天数值求解精确度的一块基石 。 随着数字化时代的到来，数值计算传热过程日臻完善。当今天用数值解来回顾古老的精 确解时，我们会发现一些有趣的现象，并使我们对传热研究中这块完美的领域有一个更深刻 [4-5] 的认识。本文基于这一点给出几个例子，以示之真 。 2. 基本导热方程 对一个标量物理变量φ，一般的稳态扩散方程可写成 1 ∂ ⎛⎜Γ ⋅F (x ) ∂φ⎞⎟+S 0 (1) f (x ) ∂x ⎜ φ i ∂x ⎟ φ i i ⎝ i ⎠ 其中，f (x ) 和F (x ) 是几何形状因子，S 是单位体积中的净源项，Γ 是对应于变量 i i φ φ φ 的扩散系数。根据方程(1)，在笛卡尔直角坐标系和圆柱坐标系中，导热方程分别为： ∂ ∂⎛ T ⎞ ∂ ∂⎛ T ⎞ ∂ ∂⎛ T ⎞ ⎜k ⎟+ ⎜k ⎟+ ⎜k ⎟+ST 0 (2) x x y y z z ∂ ∂⎝ ⎠ ∂ ∂⎝ ⎠ ∂ ∂⎝ ⎠ 1 ∂ ∂T 1 ∂ ∂T ∂ ∂T ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ⎜kr ⎟+ 2 ⎜k ⎟+ ⎜k ⎟+ST 0 (3) r r r r θ θ z z ∂ ⎝ ∂ ⎠ ∂ ⎝ ∂ ⎠ ∂ ∂⎝ ⎠ 其中T 为温度，k 为导热系数。在结构化网格中，控制方程(2-3)的离散方程之间存在着 不同的差异。 3. 算例 3.1 二维笛卡尔直角坐标系中复合边界条件的导热 问题 图 1 给出了这个问题的几何描述。顶边和左边界为绝热 边界条件，底边温度 T1 已知，右边界为第三类对流边界条件， z 方向视为无穷大。常物性、无内热源的控制方程及边界条 图1 几何形状和网格 件如下： -1-