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

凌国平; 凌国灿

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

凌国平¹,
凌国灿²

1.
苏州大学数学系, 苏州 215006;
2.
中国科学院力学研究所, 北京 100080

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

详细信息

中图分类号: O241.3;O357
计量
- 文章访问数: 2730
- HTML全文浏览量: 147
- PDF下载量: 526
- 被引次数: 0
出版历程
- 收稿日期: 1996-01-31
- 修回日期: 1997-12-29
- 刊出日期: 1998-08-15

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

Ling Guoping¹,
Ling Guocan²

1.
Department of Mathematics, Suzhou University, Suzhou 215006, P. R. China;
2.
Institute of Mechanics, Academia Sinica, Beijing 100080, P. R. China

摘要

摘要: 本文采用以修正的不完全LU分解作预处理器的共轭梯度法(MILU-CG),结合高阶隐式差分格式,改进了作者(1992)提出的基于区域分解、有限差分法与涡法杂交的数值方法(HDV).系统地研究了雷诺数Re=1000,200,旋转速度比α∈(0.5,3.25)范围内,绕旋转圆柱从突然起动到充分发展,长时间内尾流旋涡结构和阻力、升力系数的变化规律.计算所得流线与实验流场显示相比,完全吻合.首次揭示了临界状态时的旋涡结构特性,并指出最佳升阻比就在该状态附近得到.
- 旋转圆柱 /
- 涡结构 /
- 有限差分法 /
- 涡法 /
- 预处理共轭梯度法 /
- 不完全LU分解
Abstract: A 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(MIL U-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 rectilinear speed ratio A I(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 pat terns at the states are given for the first time.The maximum lift to drag force ratio can be obtained nearby the critical states.
- rotating circular cylinder /
- vortex pattern /
- finite difference method /
- vortex method /
- preconditioned conjugate gradient method /
- incomplete LU decomposition

HTML全文

参考文献(12)

[1]	M.Coutanceau and C.Menard,Influence of rotation on the near-wake development behind an impulsively started circular cylinder,J.Fluid Mech.,158(1985),399-446.
[2]	H.M.Badr and S.C.R.Dennis,Time-dependent viscous flow past an impulsively started rotating and translating circular cylinder,J.Fluid Mech.,158(1985),447-488.
[3]	H.M.Badr,M.Coutanceau and S.C.R.Dennis,et al,Unsteady flow past a rotating circular cylinder at Reynolds numbers 10³ and 10⁴,J.Fluid Mech.,220(1990),459-484.
[4]	C.C.Chang and R.L.Chern,Vortex shedding from an impulsively started rotating and translating circular cylinder,J.Fluid Mech.,233(1991),265-298.
[5]	Y.M.Chen,Y.R.Ou and A.J.Pearlstein,Development to the wake behind a circular cylinder impulsively started into rotating and rectilinear motion,J.Fluid Mech.,253(1993),449-484.
[6]	Y.R.Ou,Mathematical modeling and numerical simulation in external flow control,Flow Control,Springer-Verlag(1995),219)255.
[7]	Y.T.Chew,M.Cheng and S.C.Luo,A numerical study of flow past a rotating circular cylinder using a hybrid vortex scheme,J.Fluid Mech.,299(1995),35-71.
[8]	凌国灿、凌国平、王运平,钝体分离旋涡流动的区域分解,杂交数值模拟——Ⅰ.理论方法及其应用,中国科学,A辑,(1)(1992),70-79.
[9]	凌国灿、凌国平、顾启扬,高雷诺数下绕钝体分离流动的一种新数值算法,华中理工大学学报,20(4)(1992),87-93.
[10]	Ling Guoping and Ling Guocan,A numerical study of the vortex motion in oscillating flow around a circular cylinder at low and middle Kcnumbers,Acta Mechanica Sinica,10(3)(1994),212-219.
[11]	吕涛、石济民、林振宝,《区域分解算法》,北京,科学出版社(1992),120-143.
[12]	石济民、林振宝、岳兴业,MILU-CG方法在三维三相黑油模型中的应用,石油勘探与开发,21(6)(1994),52-58.