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

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

本文采用以修正的不完全ＬＵ分解作预处理器的共轭梯度法（ＭＩＬＵ＿ＣＧ），结合高阶隐式差分格式，改进了作者（１９９２）提出的基于区域分解、有限差分法与涡法杂交的数值方法（ＨＤＶ）·系统地研究了雷诺数Ｒｅ＝１０００，２００，旋转速度比α∈（０５，３２５）范围内，绕旋转圆柱从突然起动到充分发展，长时间内尾流旋涡结构和阻力、升力系数的变化规律·计算所得流线与实验流场显示相比，完全吻合·首次揭示了临界状态时的旋涡结构特性。

 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.