MILU-CG方法和绕旋转圆柱流动的数值研究  

MILU-CG Method and the Numerical Study on the Flow around a Rotating Circular Cylinder

作　　者：凌国平[1,2] 凌国灿[1,2] 

机构地区：[1]苏州大学数学系 [2]中国科学院力学研究所

出　　处：《应用数学和力学》1998年第8期731-739,共9页Applied Mathematics and Mechanics

基　　金：江苏省自然科学基金和中国科学院力学研究所非线性力学开放实验室资助

摘　　要：本文采用以修正的不完全ＬＵ分解作预处理器的共轭梯度法（ＭＩＬＵ＿ＣＧ），结合高阶隐式差分格式，改进了作者（１９９２）提出的基于区域分解、有限差分法与涡法杂交的数值方法（ＨＤＶ）·系统地研究了雷诺数Ｒｅ＝１０００，２００，旋转速度比α∈（０５，３２５）范围内，绕旋转圆柱从突然起动到充分发展，长时间内尾流旋涡结构和阻力、升力系数的变化规律·计算所得流线与实验流场显示相比，完全吻合·首次揭示了临界状态时的旋涡结构特性。 hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number Re=1000,200 and the angular to recti-linear speed ratio α∈(05,325) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

关 键 词：旋转圆柱 涡结构 MILU-CG法 非定常流动 

分 类 号：O357.1[理学—流体力学]