浅析对称边界在数值计算中的应用.docxVIP

-1-

浅析对称边界在数值计算中的应用

许彬1，JohnC.Chai2，陈少松1，张敏1，王学德1

1南京理工大学动力工程学院，南京(210094)

2南洋理工大学机械与航天学院，新加坡(639798)

摘要：在结构化网格中，从一维拓展到三维，采用基元中心有限容积法和全隐时间格式求解非稳态热传导问题，并将数值解同带有误差函数的精确解进行比较，得到令人满意的一致结果，这充分显示了对称边界在数值计算中的重要作用。

关键词：对称边界，有限容积法，误差函数

对称边界的选取在数值计算中有着十分重要的意义。首先，选取适当的对称边界可以减少计算区域的网格数目和求解中的计算量；其次，可以使求解的问题简洁明了；最后，对称边界的选取可以使通过计算得到的图形更加直观生动，对问题有更深刻的认识[1-4]。

本文在给出对称边界的基本概念和关系式的前提下，通过一个经典的非稳态热传导算例，从一维拓展到三维，全面展示对称边界在数值计算中的精彩之处，以此为读者在数值计算中提供有益的参考。

1基本概念和关系式

对称边界的选取必须满足一个条件，就是该边界上的通量必须为零，也就是说某标量在该边界上为第二类边界条件，其法向变化率为零，数学表达式为，

(1)

上式中φ为任意标量。对于热传导方程有，

qn=?k(2)

对于N-S方程有，

(3)

下面我们给出一个热传导问题作为参考算例，同时将数值计算的结果同数学分析后的精确解进行比较，以此验证数值计算中对称边界选取的重要性。

2算例分析

2.1一维热传导算例

一个半无限大物体（0≤x∞)初始温度是T0，当时间t0时，在x=0的边界处的温度始终为零度。该问题的数学描述为[1]，

0x∞,t0(4)

T(x,t)=0x=0,t0(5)

-2-

T(x,t)=T0x≥0,t=0(6)

通过数学分析，可得该问题的精确解为，

在对半无限大空间内的非稳态热传导问题进行数值模拟时，初始温度取常数T0=20℃,其他参数都取1，单位为国际单位。为了模拟这个问题，将上下两个边界设为对称边界，左边界给定温度0℃,右边界为绝热边界。计算时采用的尺寸为1×1，网格数为10×10，时间步长取0.01s。

图1为t=0.03s时方形区域内温度分布情况，同时给出了数值解和精确解，虚线为精确解，背景云图为数值解。图2为三个不同位置温度随时间变化的曲线。从两图中可以看出数值解和精确解吻合的比较好。精确解略低于数值解，这是由于在模拟一维半无限大热传导问题时，一方面用矩形区域近似代替半无限大区域，另一方面精确解中带有误差函数，在进行数值计算（双精度的计算）时，会存在舍入误差。

18.8

18.817.516.315.013.812.511.310.08.87.56.35.03.82.51.3

图1t=0.03s（数值解/精确解）温度场分布图2不同位置温度随时间变化图

2.2二维热传导算例

将上述一维热传导算例推广到二维空间，即为二维热传导问题。问题描述如下，一个半无限大的角域，0≤x∞,0≤y∞,初始温度是T0，当时间t0时，在x=0与y=0的边界处温度始终为零度。

由数学分离变量法推导可得，这个问题的解可以表示为如下两个一维问题解的乘积：

（1）T1(x,t)为半无限大区域（0≤x∞)的解，初始温度F1(x)=T0，当时间t0时，x=0处的边界温度始终为零度；

（2）T2(y,t)为半无限大区域（0≤y∞)的解，初始温度F2