线性椭圆型方程的数值网格生成方法.pdfVIP

　 辑第 15 卷第 2 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 15, . 2 A Ser A V o l N o 　　2000 年 6 月　　　　　　　JOU RNAL O F H YDROD YNAM ICS　　　　　June, 2000 文章编号:(2000) 线性椭圆型方程的数值网格 生成方法 李传亮, 孔祥言 ( 中国科学技术大学力学系, 安徽合肥 230027) 　　摘　要: 　用椭圆型方程生成的数值网格是一种贴体网格; 传统的网格生成方法需求解一组 非线性的椭圆型偏微分方程; 非线性椭圆型方程的求解非常困难且精度受到影响, 本文对此进行 了深入的理论研究, 提出了“用求解线性的椭圆型方程直接生成贴体网格”的新方法; 新方法简单 且易于实现, 并从理论上克服了传统方法存在的不足, 因此, 它对计算流体力学的发展具有一定的 理论意义和实用价值。 关　键　词: 　椭圆型方程; 数值计算; 贴体网格 中图分类号: 　 351. 2　　　文献标识码: O A 1　前言 ( ) 自 1974 年 . . 等人利用椭圆型方程 主要是 方程和 方程 J F T hom p son L ap lace Po isson 作为建立贴体坐标系的网格生成方法以来, 这种方法得到了迅速的发展。它不仅能处理二维或 三维问题, 而且还能处理定常和非定常问题。用椭圆型方程生成的贴体网格质量很高, 而且计 ( 算时间增加不多 因为求解L ap lace 方程或 Po isson 方程的差分迭代过程可以与求解物理方 ) 程的差分迭代过程一起进行 。目前, 椭圆型方程方法已成为一种重要的数值网格生成方法, 它 在计算流体力学中有着广泛的应用。但是, 它仍存在一个缺陷: 它将线性的椭圆型方程 ( ) L ap lace 方程方法或 Po isson 方程 首先变换成非线性的椭圆型方程, 然后求解生成贴体网 格[ 1, 2 ]。整个求解过程十分复杂。本文对此进行了研究, 从理论上提出了用线性的方程直接生成 网格系统的新方法, 并用计算实例说明了这种方法的优越性。这无疑会给数值计算工作带来许 多方便。 2　现用方法介绍 对物理平面上的区域D ( 图 1 (a) ) , 用贴体网格进行离散。所谓贴体网格, 就是一种曲线网 格, 它使物理平面上所有的边界均成为网格线, 内部网格线也与边界形状和汇源条件相一致 ( 图 1 ( ) ) 。生成这样的贴体网格系统是通过将D ( a 变换到计算平面上的区域D 来实现的 图 收稿日期: ( ) ( ) 作者简介: 　李传亮 1962～ , 男, 副教授, 博士; 孔祥言 1932～ , 男, 教授, 博导。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 李传亮: 椭圆型方程的数值网