西 安 交 通  大  学  学  报

Journal of Xi'an Jiaotong University

采用算子分裂法求解不可压黏性流动
刘海湖¹，李开泰²，苏剑²，王尚锦¹
(1.西安交通大学能源与动力工程学院，710049，西安；2.西安交通大学理学院，710049，西安)

摘要：采用算子分裂算法求解Dirichlet边界条件的不可压黏性流动，通过时间离散，在每个时间步把Navier-Stokes方程分解成两个广义Stokes问题和一个非线性问题，分别采用共轭梯度算法求解这两类问题，从而逐个解决了不可压缩性和非线性两大数值困难，同时计算了一个存在解析解的广义Stokes问题和顶盖驱动流问题.计算结果表明，该算法求解不可压流动是可行的，并且具有精度高、稳定性好、收敛速度快的特点．
关键词：Stokes问题；算子分裂；共轭梯度法；顶盖驱动流；Dirichlet边界；有限元
中图分类号：TP391.7文献标识码：A文章编号：0253-987X(2006)12-1470-03