Methods on applying stream-function restraints in differential quadrature modelling of two-dimensional flow
WANG Tong1,HE Tao1,3,CAO Shuyang2
1.College of Civil Engineering,Shanghai Normal University,Shanghai 201418,China;
2.State Key Laboratory for Disaster Reduction in Civil Engineering,Tongji University,Shanghai 200092,China;
3.School of Engineering,University of Birmingham,Birmingham B15 2TT,UK
Abstract:The 2D lid-driven cavity flow was simulated by applying the differential quadrature method to solve the stream function-vorticity equations.There were two boundary conditions,one Dirichlet and one Neumann,for the stream function equation at each solid boundary though the stream function equation was just second order.Analysis on this over- specified problem was carried out,based on which a new applying method was proposed: the Neumann condition was considered in calculating the vorticity at the boundary while only the Dirichlet condition was applied in the stream function equation.Validity of this method was verified by comparing its numerical results with benchmark data.Two other existing methods,the one-layer approach and the two-layer approach were shown as contrasts.Trial calculations indicate that the one-layer approach is sensitive to the parity of grid numbers and is not suitable for the present problem.Comparisons between the new method and the two-layer approach show that the former is not only more accurate but also more convenient to be used in practice for avoiding the over-specified problem.
[1] Bellman R E, Casti J. Differential quadrature and long-term integration [J]. Journal of Mathematical Analysis and Applications, 1971, 34(2): 235–238.
[2] Bert C W, Malik M. Differential quadrature method in computational mechanics: a review [J]. Applied Mechanics Reviews, 1996, 49(1): 1–28.
[3] Shu C. Differential quadrature and its application in engineering [M]. London: Springer, 2000.
[4] 程昌钧, 朱正佑. 微分求积方法及其在力学应用中的若干新进展 [J]. 上海大学学报(自然科学版), 2009, 15(6): 551–559.
CHENG Chang-jun, ZHU Zheng-you. Recent advances in differential quadrature method with applications to mechanics [J]. Journal of Shanghai University (Natural Science Edition), 2009, 15(6): 551–559.
[5] 李鸿晶, 王通. 结构地震反应分析的逐步微分积分方法 [J]. 力学学报, 2011, 43(2): 430–435.
LI Hong-jing, WANG Tong. A time-stepping method of seismic response analysis for structures using differential quadrature rule [J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(2): 430–435.
[6] 武吉梅, 陈媛, 王砚, 等. 基于微分求积法的印刷运动薄膜动力稳定性分析 [J]. 振动与冲击, 2015, 34(20): 57–60.
WU Ji-mei, CHEN Yuan, WANG Yan et al. Study of kinetic stability for membrane based on differential quadrature method [J]. Journal of Vibration and Shock, 2015, 34(20): 57–60.
[7] 董宇, 杨翊仁, 鲁丽. 基于微分求积法的轴向流作用下二维板复杂响应研究 [J]. 振动与冲击, 2015, 34(6): 46–51.
DONG Yu,YANG Yi-ren,LU Li. On the complicated responses of a two-dimensional plate in axial liquid flow by differential quadrature method [J]. Journal of Vibration and Shock, 2015, 34(6): 46–51.
[8] Ghia U, Ghia K N, Shin C T. High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method [J]. Journal of Computational Physics, 1982, 48(3): 387–411.
[9] Striz A G, Chen W. Application of the differential quadrature method to the driven cavity flow problem [J]. International Journal of Non-Linear Mechanics, 1994, 29(5): 665–670.
[10] Shu C, Wang L, Chew Y T. Comparative studies of three approaches for GDQ computation of incompressible Navier- Stokes equations in primitive variable form [J]. International Journal of Computational Fluid Dynamics, 2004, 18(5): 401– 412.
[11] Bruneau C H, Saad M. The 2D lid-driven cavity problem revisited [J]. Computers & Fluids, 2006, 35(3): 326–348.
[12] Marchi C H, Suero R, Araki L K. The lid-driven square cavity flow: numerical solution with a 1024 x 1024 grid [J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2009, 31(3): 186–198.
[13] Davis M B, Carey G F. Iterative solution of the stream function-vorticity equations using a multigrid solver with finite elements [J]. Communications in Numerical Methods in Engineering, 1993, 9(7): 587–594.
[14] Shu C, Xue H. Comparison of two approaches for implementing stream function boundary conditions in DQ simulation of natural convection in a square cavity [J]. International Journal of Heat and Fluid Flow, 1998, 19(1): 59–68.
[15] 刘剑, 王鑫伟. 基于微分求积法的逐步积分法与常用时间积分法的比较 [J]. 力学季刊, 2008, 29(2): 304–309.
LIU Jian, WANG Xin-wei. Comparisons of successive integration method based on differential quadrature with commonly used time integration schemes [J]. Chinese Quarterly of Mechanics, 2008, 29(2): 304–309.
[16] 汪芳宗, 廖小兵, 谢雄. 微分求积法的特性及其改进 [J]. 计算力学学报, 2015, 32(6): 765–771.
WANG Fang-zong, LIAO Xiao-bing, XIE Xiong. Characteristic of the differential quadrature method and its improvement [J]. Chinese Journal of Computational Mechanics, 2015, 32(6): 765–771.
[17] Ku H C, Taylor T D, Hirsh R S. Pseudospectral methods for solution of the incompressible Navier-Stokes equations [J]. Computers & Fluids, 1987, 15(2): 195–214.
[18] Wang X, Liu F, Wang X, Gan L. New approaches in application of differential quadrature method to fourth-order differential equations [J]. Communications in Numerical Methods in Engineering, 2005, 21(2): 61–71.