Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow-structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier-Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.

[1]  Peter McCorquodale,et al.  A Cartesian grid embedded boundary method for the heat equation on irregular domains , 2001 .

[2]  Tony W. H. Sheu,et al.  An effective explicit pressure gradient scheme implemented in the two-level non-staggered grids for incompressible Navier-Stokes equations , 2008, J. Comput. Phys..

[3]  M. Minion A Projection Method for Locally Refined Grids , 1996 .

[4]  Simon Tavener,et al.  Bifurcation for flow past a cylinder between parallel planes , 1995, Journal of Fluid Mechanics.

[5]  M. Provansal,et al.  Bénard-von Kármán instability: transient and forced regimes , 1987, Journal of Fluid Mechanics.

[6]  D. Young,et al.  Simulation of laminar vortex shedding flow past cylinders using a coupled BEM and FEM model , 2001 .

[7]  Luis Chacón,et al.  A fully implicit, nonlinear adaptive grid strategy , 2006, J. Comput. Phys..

[8]  S. Dennis,et al.  Numerical solutions for steady flow past a circular cylinder at Reynolds numbers up to 100 , 1970, Journal of Fluid Mechanics.

[9]  W. Shyy,et al.  Elafint: a Mixed Eulerian-Lagrangian Method for Fluid Flows with Complex and Moving Boundaries , 1996 .

[10]  Yih-Ferng Peng,et al.  Transition in a 2-D lid-driven cavity flow , 2003 .

[11]  A. Sau Hopf bifurcations in the wake of a square cylinder , 2009 .

[12]  T. Hsu,et al.  Three-dimensional evolution of vortical structures and associated flow bifurcations in the wake of two side-by-side square cylinders , 2007 .

[13]  A. Chorin Numerical Solution of the Navier-Stokes Equations* , 1989 .

[14]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[15]  Liang Cheng,et al.  A finite difference solution of the shear flow over a circular cylinder , 2000 .

[16]  L. Kaiktsis,et al.  Control of circular cylinder wakes using base mass transpiration , 2001 .

[17]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[18]  K. Sreenivasan,et al.  On the formation and suppression of vortex ‘shedding’ at low Reynolds numbers , 1990, Journal of Fluid Mechanics.

[19]  Gianluca Iaccarino,et al.  An approach to local refinement of structured grids , 2002 .

[20]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[21]  M. Berger,et al.  An Adaptive Version of the Immersed Boundary Method , 1999 .

[22]  Robert D. Russell,et al.  Anr-Adaptive Finite Element Method Based upon Moving Mesh PDEs , 1999 .

[23]  Chang Shu,et al.  A stencil adaptive algorithm for finite difference solution of incompressible viscous flows , 2006, J. Comput. Phys..

[24]  John B. Bell,et al.  Cartesian grid method for unsteady compressible flow in irregular regions , 1995 .

[25]  Pingwen Zhang,et al.  Moving mesh methods in multiple dimensions based on harmonic maps , 2001 .

[26]  Odd M. Faltinsen,et al.  A local directional ghost cell approach for incompressible viscous flow problems with irregular boundaries , 2008, J. Comput. Phys..

[27]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity , 2000 .

[28]  H. Sakamoto,et al.  Optimum Suppression of Fluid Forces Acting on a Circular Cylinder , 1994 .

[29]  Bin Teng,et al.  Numerical simulation of viscous flow past two circular cylinders of different diameters , 2005 .

[30]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[31]  C. Peskin The Fluid Dynamics of Heart Valves: Experimental, Theoretical, and Computational Methods , 1982 .

[32]  Hiroshi Sakamoto,et al.  An optimum suppression of fluid forces by controlling a shear layer separated from a square prism , 1991 .

[33]  W. Shyy,et al.  ELAFINT - A mixed Eulerian-Lagrangian method for fluid flows with complex and moving boundaries , 1994 .

[34]  A. Sau Growth of unsteady Hopf bifurcation modes and their swapping in a transitional coupled wake , 2008 .

[35]  Alexei M. Khokhlov,et al.  Fully Threaded Tree Algorithms for Adaptive Refinement Fluid Dynamics Simulations , 1997, astro-ph/9701194.

[36]  T. Sheu,et al.  Structural development of vortical flows around a square jet in cross-flow , 2004, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[37]  A. Nadim,et al.  An Embedding Method for Bluff Body Flows: Interactions of Two Side-by-Side Cylinder Wakes , 2003 .

[38]  O. C. Zienkiewicz,et al.  Adaptive techniques in the finite element method , 1988 .

[39]  A. S. Grove,et al.  An experimental investigation of the steady separated flow past a circular cylinder , 1964, Journal of Fluid Mechanics.

[40]  Charles Dalton,et al.  The Suppression of Lift on a Circular Cylinder due to Vortex Shedding at Moderate Reynolds Numbers , 2001 .

[41]  Phillip Colella,et al.  A cell-centered adaptive projection method for the incompressible Navier-Stokes equations in three dimensions , 2007, J. Comput. Phys..

[42]  Khodor Khadra,et al.  Fictitious domain approach for numerical modelling of Navier–Stokes equations , 2000 .

[43]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[44]  S. Balachandar,et al.  Direct Numerical Simulation of Flow Past Elliptic Cylinders , 1996 .

[45]  J. Ferziger,et al.  A ghost-cell immersed boundary method for flow in complex geometry , 2002 .

[46]  Margot Gerritsen,et al.  Designing an efficient solution strategy for fluid flows. 1. A stable high order finite difference scheme and sharp shock resolution for the Euler equations , 1996 .

[47]  Charles S. Peskin,et al.  Flow patterns around heart valves: a digital computer method for solving the equations of motion , 1973 .

[48]  M. D. Salas,et al.  Euler calculations for multielement airfoils using Cartesian grids , 1986 .

[49]  Daniel F. Martin,et al.  A Cell-Centered Adaptive Projection Method for the Incompressible Euler Equations , 2000 .

[50]  Boyce E. Griffith,et al.  An adaptive, formally second order accurate version of the immersed boundary method , 2007, J. Comput. Phys..

[51]  P. Colella,et al.  A Conservative Adaptive Projection Method for the Variable Density Incompressible Navier-Stokes Equations , 1998 .