A sub-cell spectral-element simulator for 2D rigid-body-fluid interaction problems

In this work, we propose an immersed-boundary-type simulation method for the two-way coupling problems between a rigid body and a fluid. In this simulation tool, the entire fluid-solid domain is treated as an incompressible fluid with non-uniform density and the no-slip boundary condition at the rigid body surface is enforced by the penalization method. The fluid solver is developed in use of the spectral element method for the spatial discretization and the mixed explicit/implicit scheme for the temporal discretization. An additional Lagrangian mesh is employed and attached to the rigid body in order to trace the rigid body and to perform the area integration over the solid domain. Besides, a so-called sub-cell scheme is developed to smooth the discontinuity at the fluid-solid interface. The validity and accuracy of the proposed simulation method were examined well by applying it to the sedimentation problems of circular, triangular, square, as well as elliptic cylinders in a channel. An accuracy of 2nd order was observed, probably due to the use of triangular Lagrangian elements and consequently a piecewise-linear approximation of the rigid body shape and also due to the 2nd-order interpolation involved in the sub-cell scheme.

[1]  Valentino Pediroda,et al.  Fictitious Domain approach with hp-finite element approximation for incompressible fluid flow , 2009, J. Comput. Phys..

[2]  Charles S. Peskin,et al.  Modeling prosthetic heart valves for numerical analysis of blood flow in the heart , 1980 .

[3]  Zhaosheng Yu,et al.  A direct-forcing fictitious domain method for particulate flows , 2007, J. Comput. Phys..

[4]  Lucy T. Zhang,et al.  Semi-implicit formulation of the immersed finite element method , 2011, Computational Mechanics.

[5]  B. Fornberg.  A numerical study of steady viscous flow past a circular cylinder , 1980, Journal of Fluid Mechanics.

[6]  N. Patankar,et al.  A fast computation technique for the direct numerical simulation of rigid particulate flows , 2005 .

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

[8]  Sourabh V. Apte,et al.  A variable-density fictitious domain method for particulate flows with broad range of particle-fluid density ratios , 2013, J. Comput. Phys..

[9]  Wing Kam Liu,et al.  Immersed finite element method for rigid body motions in the incompressible Navier–Stokes flow , 2008 .

[10]  A. Roshko.  On the development of turbulent wakes from vortex streets , 1953 .

[11]  Emmanuel Maitre,et al.  Convergence Analysis of a Penalization Method for the Three-Dimensional Motion of a Rigid Body in an Incompressible Viscous Fluid , 2010, SIAM J. Numer. Anal..

[12]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[13]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .

[14]  J. Happel,et al.  Low Reynolds number hydrodynamics: with special applications to particulate media , 1973 .

[15]  Petros Koumoutsakos,et al.  Simulations of single and multiple swimmers with non-divergence free deforming geometries , 2011, J. Comput. Phys..

[16]  Mathieu Martin,et al.  A numerical method for fully resolved simulation (FRS) of rigid particle-flow interactions in complex flows , 2009, J. Comput. Phys..

[17]  François Peters,et al.  A fictitious domain approach for the simulation of dense suspensions , 2014, J. Comput. Phys..

[18]  C. Peskin,et al.  A three-dimensional computational method for blood flow in the heart. 1. Immersed elastic fibers in a viscous incompressible fluid , 1989 .

[19]  S. Orszag,et al.  High-order splitting methods for the incompressible Navier-Stokes equations , 1991 .

[20]  Anthony Wachs,et al.  A DEM-DLM/FD method for direct numerical simulation of particulate flows: Sedimentation of polygonal isometric particles in a Newtonian fluid with collisions , 2009 .

[21]  David L. Finn,et al.  Falling Paper and Flying Business Cards , 2007 .

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

[23]  Timothy C. Warburton,et al.  Basis Functions for Triangular and Quadrilateral High-Order Elements , 1999, SIAM J. Sci. Comput..

[24]  Charles S. Peskin,et al.  Improved Volume Conservation in the Computation of Flows with Immersed Elastic Boundaries , 1993 .

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

[26]  Mathieu Coquerelle,et al.  ARTICLE IN PRESS Available online at www.sciencedirect.com Journal of Computational Physics xxx (2008) xxx–xxx , 2022 .

[27]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[28]  C. Williamson.  Defining a Universal and Continuous Strouhal-Reynolds Number Relationship for the Laminar Vortex She , 1988 .

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

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

[31]  Lucy T. Zhang.  Immersed finite element method for fluid-structure interactions , 2007 .

[32]  Lucy T. Zhang,et al.  Immersed finite element method , 2004 .