A fast immersed interface method for solving Stokes flows on irregular domains

We present a fast immersed interface method for solving the steady Stokes flows involving the rigid boundaries. The immersed rigid boundary is represented by a set of Lagrangian control points. In order to enforce the prescribed velocity at the rigid boundary, singular forces at the rigid boundary are applied on the fluid. The forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are approximated using the cubic splines. The strength of singular forces is determined by solving a small system of equations via the GMRES method. The Stokes equations are discretized using finite difference method with the incorporation of jump conditions on a staggered Cartesian grid and solved by the conjugate gradient Uzawa-type method. Numerical results demonstrate the accuracy and ability of the proposed method to simulate Stokes flows on irregular domains.

[1]  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 .

[2]  Cheng Wang,et al.  A Fast Finite Differenc Method For Solving Navier-Stokes Equations on Irregular Domains , 2003 .

[3]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[4]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[5]  D. L. Young,et al.  Method of fundamental solutions for multidimensional Stokes equations by the dual-potential formulation , 2006 .

[6]  Christoph Börgers,et al.  Domain imbedding methods for the Stokes equations , 1990 .

[7]  Vita Rutka,et al.  A staggered grid‐based explicit jump immersed interface method for two‐dimensional Stokes flows , 2008 .

[8]  Ronald Fedkiw,et al.  The immersed interface method. Numerical solutions of PDEs involving interfaces and irregular domains , 2007, Math. Comput..

[9]  M. Uhlmann An immersed boundary method with direct forcing for the simulation of particulate flows , 2005, 1809.08170.

[10]  E. Kim,et al.  A mixed galerkin method for computing the flow between eccentric rotating cylinders , 1998 .

[11]  A multigrid pseudospectral method for steady flow computation , 2003 .

[12]  Z. J. Wang,et al.  A cartesian grid method for modeling multiple moving objects in 2D incompressible viscous flow , 2003 .

[13]  Andreas Wiegmann,et al.  The Explicit-Jump Immersed Interface Method: Finite Difference Methods for PDEs with Piecewise Smooth Solutions , 2000, SIAM J. Numer. Anal..

[14]  Jeng-Tzong Chen,et al.  A new method for Stokes problems with circular boundaries using degenerate kernel and Fourier series , 2008 .

[15]  C. Alves,et al.  Density results using Stokeslets and a method of fundamental solutions for the Stokes equations , 2004 .

[16]  D. L. Young,et al.  Short Note: The method of fundamental solutions for 2D and 3D Stokes problems , 2006 .

[17]  C. Peskin The immersed boundary method , 2002, Acta Numerica.

[18]  J. Damasceno,et al.  Numerical simulation of two-dimensional flows over a circular cylinder using the immersed boundary method , 2003 .

[19]  Zhilin Li,et al.  A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow , 2008, Adv. Comput. Math..

[20]  D. Zorin,et al.  A fast solver for the Stokes equations with distributed forces in complex geometries , 2004 .

[21]  Randall J. LeVeque,et al.  Immersed Interface Methods for Stokes Flow with Elastic Boundaries or Surface Tension , 1997, SIAM J. Sci. Comput..

[22]  Zhilin Li,et al.  The immersed interface method for the Navier-Stokes equations with singular forces , 2001 .

[23]  Gregory H. Wannier,et al.  A contribution to the hydrodynamics of lubrication , 1950 .

[24]  Aaron L. Fogelson,et al.  Computational Methods for Continuum Models of Platelet Aggregation , 1999 .

[25]  Cornelis W. Oosterlee,et al.  Multigrid Methods for the Stokes System , 2006, Computing in Science & Engineering.

[26]  A multigrid finite difference approach to steady flow between eccentric rotating cylinders , 2000 .

[27]  H. Elman Multigrid and Krylov subspace methods for the discrete Stokes equations , 1994 .

[28]  D. L. Young,et al.  The method of fundamental solutions for Stokes flow in a rectangular cavity with cylinders , 2005 .

[29]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[30]  Aaron L. Fogelson,et al.  Continuum models of platelet aggregation: formulation and mechanical properties , 1992 .

[31]  Eric Yu Tau Numerical solution of the steady Stokes equations , 1992 .

[32]  Dongho Shin,et al.  Fast solvers for finite difference approximations for the stokes and navier-stokes equations , 1996, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[33]  Christoph BiJrgers Domain Imbedding Methods for the Stokes Equations , 1990 .

[34]  L. Gray,et al.  The boundary contour method for two-dimensional Stokes flow and incompressible elastic materials , 2002 .

[35]  H. Elman Preconditioners for saddle point problems arising in computational fluid dynamics , 2002 .

[36]  D. Young,et al.  Method of Fundamental Solutions for Three-Dimensional Stokes Flow in Exterior Field , 2006 .

[37]  Randall J. LeVeque,et al.  An Immersed Interface Method for Incompressible Navier-Stokes Equations , 2003, SIAM J. Sci. Comput..

[38]  Z. Jane Wang,et al.  An immersed interface method for simulating the interaction of a fluid with moving boundaries , 2006, J. Comput. Phys..

[39]  J. Spurk Boundary Layer Theory , 2019, Fluid Mechanics.

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

[41]  Boo Cheong Khoo,et al.  An immersed interface method for viscous incompressible flows involving rigid and flexible boundaries , 2006, J. Comput. Phys..

[42]  Jörg Peters,et al.  Fast Iterative Solvers for Discrete Stokes Equations , 2005, SIAM J. Sci. Comput..

[43]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[44]  Ahmed H. Sameh,et al.  An Efficient Iterative Method for the Generalized Stokes Problem , 1998, SIAM J. Sci. Comput..

[45]  Yijun Liu,et al.  A new fast multipole boundary element method for solving 2-D Stokes flow problems based on a dual BIE formulation , 2008 .

[46]  D. L. Young,et al.  The method of fundamental solutions for inverse 2D Stokes problems , 2005 .

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

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

[49]  R. Rivlin,et al.  Flow of a Newtonian fluid between eccentric rotating cylinders: Inertial effects , 1976 .

[50]  P. Swarztrauber,et al.  Efficient Fortran subprograms for the solution of separable elliptic partial differential equations , 1979 .