An efficient fluid-solid coupling algorithm for single-phase flows

We present a simple and efficient fluid-solid coupling method in two and three spatial dimensions. In particular, we consider the numerical approximation of the Navier-Stokes equations on irregular domains and propose a novel approach for solving the Hodge projection step on arbitrary shaped domains. This method is straightforward to implement and leads to a symmetric positive definite linear system for both the projection step and for the implicit treatment of the viscosity. We demonstrate the accuracy of our method in the L^1 and L^~ norms and present its removing the errors associated with the conventional rasterization-type discretizations. We apply this method to the simulation of a flow past a cylinder in two spatial dimensions and show that our method can reproduce the known stable and unstable regimes as well as correct lift and drag forces. We also apply this method to the simulation of a flow past a sphere in three spatial dimensions at low and moderate Reynolds number to reproduce the known steady axisymmetric and non-axisymmetric flow regimes. We further apply this algorithm to the coupling of flows with moving rigid bodies.

[1]  William W. Willmarth,et al.  Some experimental results on sphere and disk drag , 1971 .

[2]  C. Peskin Flow patterns around heart valves: A numerical method , 1972 .

[3]  A. Chorin A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[4]  R. Fedkiw,et al.  A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem , 2005 .

[5]  R. Bouard,et al.  Experimental determination of the main features of the viscous flow in the wake of a circular cylinder in uniform translation. Part 1. Steady flow , 1977, Journal of Fluid Mechanics.

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

[7]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[8]  C. Hirsch,et al.  Introduction to “Towards the Ultimate Conservative Difference Scheme. V. A Second-Order Sequel to Godunov's Method” , 1997 .

[9]  Kazufumi Ito,et al.  A well-conditioned augmented system for solving Navier-Stokes equations in irregular domains , 2009, J. Comput. Phys..

[10]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[11]  Chaoqun Liu,et al.  Preconditioned Multigrid Methods for Unsteady Incompressible Flows , 1997 .

[12]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[13]  G. Shortley,et al.  The Numerical Solution of Laplace's Equation , 1938 .

[14]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[15]  Frédéric Gibou,et al.  A second order accurate projection method for the incompressible Navier-Stokes equations on non-graded adaptive grids , 2006, J. Comput. Phys..

[16]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[17]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[18]  Frédéric Gibou,et al.  Geometric integration over irregular domains with application to level-set methods , 2007, J. Comput. Phys..

[19]  Frédéric Gibou,et al.  A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids , 2006, J. Comput. Phys..

[20]  M. Braza,et al.  Numerical study and physical analysis of the pressure and velocity fields in the near wake of a circular cylinder , 1986, Journal of Fluid Mechanics.

[21]  V. C. Patel,et al.  Flow past a sphere up to a Reynolds number of 300 , 1999, Journal of Fluid Mechanics.

[22]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[23]  George Em Karniadakis,et al.  A semi-Lagrangian high-order method for Navier-Stokes equations , 2001 .

[24]  D. Calhoun A Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity Equations in Irregular Regions , 2002 .

[25]  T. Aslam A partial differential equation approach to multidimensional extrapolation , 2004 .

[26]  Antony Jameson,et al.  A new implicit algorithm with multigrid for unsteady incompressible flow calculations , 1995 .

[27]  Robert Bridson,et al.  A fast variational framework for accurate solid-fluid coupling , 2007, SIGGRAPH 2007.