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

浅析对称边界在数值计算中的应用 1 2 1 1 1 许彬 ，John C. Chai ，陈少松 ，张敏 ，王学德 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天学院，新加坡 (639798) E-mail ：norkistar@126.com 摘 要：在结构化网格中，从一维拓展到三维，采用基元中心有限容积法和全隐时间格式求 解非稳态热传导问题，并将数值解同带有误差函数的精确解进行比较，得到令人满意的一致 结果，这充分显示了对称边界在数值计算中的重要作用。 关键词：对称边界，有限容积法，误差函数 对称边界的选取在数值计算中有着十分重要的意义。首先，选取适当的对称边界可以减 少计算区域的网格数目和求解中的计算量；其次，可以使求解的问题简洁明了；最后，对称 [1-4] 边界的选取可以使通过计算得到的图形更加直观生动，对问题有更深刻的认识 。 本文在给出对称边界的基本概念和关系式的前提下，通过一个经典的非稳态热传导算 例，从一维拓展到三维，全面展示对称边界在数值计算中的精彩之处，以此为读者在数值计 算中提供有益的参考。 1 基本概念和关系式 对称边界的选取必须满足一个条件，就是该边界上的通量必须为零，也就是说某标量在 该边界上为第二类边界条件，其法向变化率为零，数学表达式为， ∂φ 0 (1) ∂n 上式中φ 为任意标量。对于热传导方程有， ∂T qn −k 0 (2) ∂n 对于N-S 方程有， ∂u v w ∂ ∂ 0 (3) ∂x ∂y ∂z 下面我们给出一个热传导问题作为参考算例，同时将数值计算的结果同数学分析后的 精确解进行比较，以此验证数值计算中对称边界选取的重要性。 2 算例分析 2.1 一维热传导算例 一个半无限大物体（0 ≤x ∞）初始温度是T ，当时间t 0 时，在x 0 的边界处 0 的温度始终为零度。该问题的数学描述为[1]， 2 ∂ T ( x, t) 1 ∂T ( x, t) 0 = ∞, 0 (4) x t ∂x 2 α ∂t T (x ,t) 0 x 0,t 0 (5) - 1 -