摘要
采用作者提出的基于区域分解、有限差分法与涡法杂交的数值方法，结合高阶隐式差分格式，和以修正的不完全ＬＵ分解为预处理器的共轭梯度法作求解器．系统地研究了雷诺数Ｒｅ＝１０００，旋转速度比ａ∈（０．５，３．２５）范围内，绕旋转圆柱从突然起动到充分发展，长时间内尾流旋涡结构和阻力、升力系数的变化规律．计算所得流动图案与实验流场显示符合很好．数值试验证实了临界状态的存在，并首次给出了临界状态时的旋涡结构特性。
A hybrid finite difference and vortex method, which is based on domain decomposition and proposed by the authors (1992), is used forcalculating the flow around a rotating circular cylinder at Reynoldsnumber Re=1 000 and the rotation velocity ratio a∈(0.5,3.25) respectively. Some improvments are made in solution process: ahigh order full implicit difference algorithm is adopted, and the Solvers is replaced by a very efficient modified incomplete LU decomposition conjugate gradient method (MILU-CG). The long-time fullydeveloped features about the variations of the vortex patterns in thewake , and drag, lift forces on the cylinder are given. The calculated streamline contour plots agreed well with the experimentalvisualized flow pictures. The existence of the critical state isnumerically shown for the first time, and the vortex patterns at corresponding states are given, too . The maximum lift to drag forceratio can be obtained near the critical state.
出处
《力学学报》
EI
CSCD
北大核心
1997年第1期8-16,共9页
Chinese Journal of Theoretical and Applied Mechanics
基金
江苏省自然科学基金
中国科学院力研究所LNM开放实验室资助
关键词
旋转圆柱
涡尾流
有限差分法
涡法
rotating circular cylinder, vortex wake flow, finite difference method,vortex method, preconditioned conjugate gradient method, incomplete LU decomposition
