A fast moving least squares approximation with adaptive Lagrangian mesh refinement for large scale immersed boundary simulations

Abstract In this paper we propose and test the validity of simple and easy-to-implement algorithms within the immersed boundary framework geared towards large scale simulations involving thousands of deformable bodies in highly turbulent flows. First, we introduce a fast moving least squares (fast-MLS) approximation technique with which we speed up the process of building transfer functions during the simulations which leads to considerable reductions in computational time. We compare the accuracy of the fast-MLS against the exact moving least squares (MLS) for the standard problem of uniform flow over a sphere. In order to overcome the restrictions set by the resolution coupling of the Lagrangian and Eulerian meshes in this particular immersed boundary method, we present an adaptive Lagrangian mesh refinement procedure that is capable of drastically reducing the number of required nodes of the basic Lagrangian mesh when the immersed boundaries can move and deform. Finally, a coarse-grained collision detection algorithm is presented which can detect collision events between several Lagrangian markers residing on separate complex geometries with minimal computational overhead.

[1]  Gretar Tryggvason,et al.  Multiscale computations of thin films in multiphase flows , 2010 .

[2]  Gui-Rong Liu,et al.  An Introduction to Meshfree Methods and Their Programming , 2005 .

[3]  D Dries Hegen,et al.  Element-free Galerkin methods in combination with finite element approaches , 1996 .

[4]  D. Cohen-Or,et al.  Robust moving least-squares fitting with sharp features , 2005, ACM Trans. Graph..

[5]  Lu De-tang,et al.  Curve and Surface Fitting Based on Moving Least-Squares Methods , 2004 .

[6]  Scott Schaefer,et al.  Image deformation using moving least squares , 2006, ACM Trans. Graph..

[7]  Hyun Gyu Kim,et al.  Analysis of thin beams, using the meshless local Petrov–Galerkin method, with generalized moving least squares interpolations , 1999 .

[8]  W. Breugem,et al.  Turbulent channel flow of dense suspensions of neutrally buoyant spheres , 2014, Journal of Fluid Mechanics.

[9]  G. Iaccarino,et al.  Immersed boundary technique for turbulent flow simulations , 2003 .

[10]  Giuseppe Pascazio,et al.  A moving-least-squares immersed boundary method for simulating the fluid-structure interaction of elastic bodies with arbitrary thickness , 2016, J. Comput. Phys..

[11]  Eitan Grinspun,et al.  Multimaterial mesh-based surface tracking , 2014, ACM Trans. Graph..

[12]  Guirong Liu,et al.  Smoothed Particle Hydrodynamics: A Meshfree Particle Method , 2003 .

[13]  Ravi Krishna Kolluri,et al.  Provably good moving least squares , 2005, SIGGRAPH Courses.

[14]  Ted Belytschko,et al.  Enforcement of essential boundary conditions in meshless approximations using finite elements , 1996 .

[15]  Eric Loth,et al.  Quasi-steady shape and drag of deformable bubbles and drops , 2008 .

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

[17]  Leif Kobbelt,et al.  A survey of point-based techniques in computer graphics , 2004, Comput. Graph..

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

[19]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[20]  Markus Uhlmann,et al.  Sedimentation of a dilute suspension of rigid spheres at intermediate Galileo numbers: the effect of clustering upon the particle motion , 2014, Journal of Fluid Mechanics.

[21]  Vamsi Spandan,et al.  A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes , 2016, J. Comput. Phys..

[22]  A. Prosperetti Life and death by boundary conditions , 2015, Journal of Fluid Mechanics.

[23]  T. Belytschko,et al.  Fracture and crack growth by element free Galerkin methods , 1994 .

[24]  Farid F. Abraham,et al.  Functional Dependence of Drag Coefficient of a Sphere on Reynolds Number , 1970 .

[25]  Earl H. Dowell,et al.  Modeling of Fluid-Structure Interaction , 2001 .

[26]  D. Lohse,et al.  Physical mechanisms governing drag reduction in turbulent Taylor–Couette flow with finite-size deformable bubbles , 2018, Journal of Fluid Mechanics.

[27]  Jochen Fröhlich,et al.  Imposing the free-slip condition with a continuous forcing immersed boundary method , 2015, J. Comput. Phys..

[28]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[29]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .

[30]  Jochen Fröhlich,et al.  An immersed boundary method for the simulation of bubbles with varying shape , 2016, J. Comput. Phys..

[31]  Elias Balaras,et al.  A moving-least-squares reconstruction for embedded-boundary formulations , 2009, J. Comput. Phys..

[32]  L. Brandt,et al.  Sedimentation of finite-size spheres in quiescent and turbulent environments , 2015, Journal of Fluid Mechanics.

[33]  S. Mortazavi,et al.  Drop breakup and drop pair coalescence using front-tracking method in three dimensions , 2018 .

[34]  A. Stuart,et al.  Algorithms for particle-field simulations with collisions , 2001 .