Immersed boundary method

Last updated

In computational fluid dynamics, the immersed boundary method originally referred to an approach developed by Charles Peskin in 1972 to simulate fluid-structure (fiber) interactions. [1] Treating the coupling of the structure deformations and the fluid flow poses a number of challenging problems for numerical simulations (the elastic boundary changes the flow of the fluid and the fluid moves the elastic boundary simultaneously). In the immersed boundary method the fluid is represented in an Eulerian coordinate system and the structure is represented in Lagrangian coordinates. For Newtonian fluids governed by the Navier–Stokes equations, the fluid equations are

Contents

and if the flow is incompressible, we have the further condition that

The immersed structures are typically represented as a collection of one-dimensional fibers, denoted by . Each fiber can be viewed as a parametric curve where is the Lagrangian coordinate along the fiber and is time. The physics of the fiber is represented via a fiber force distribution function . Spring forces, bending resistance or any other type of behavior can be built into this term. The force exerted by the structure on the fluid is then interpolated as a source term in the momentum equation using

where is the Dirac δ function. The forcing can be extended to multiple dimensions to model elastic surfaces or three-dimensional solids. Assuming a massless structure, the elastic fiber moves with the local fluid velocity and can be interpolated via the delta function

where denotes the entire fluid domain. Discretization of these equations can be done by assuming an Eulerian grid on the fluid and a separate Lagrangian grid on the fiber. Approximations of the Delta distribution by smoother functions will allow us to interpolate between the two grids. Any existing fluid solver can be coupled to a solver for the fiber equations to solve the Immersed Boundary equations. Variants of this basic approach have been applied to simulate a wide variety of mechanical systems involving elastic structures which interact with fluid flows.

Since the original development of this method by Peskin, a variety of approaches have been developed to simulate flow over complicated immersed bodies on grids that do not conform to the surface of the body. These include methods such as the immersed interface method, the Cartesian grid method, the ghost fluid method and the cut-cell method. Mittal and Iaccarino [2] refer to all these (and other related) methods as Immersed Boundary Methods and provide various categorizations of these methods. From the point of view of implementation, they categorize immersed boundary methods into continuous forcing and discrete forcing methods. In the former, a force term is added to the continuous Navier-Stokes equations before discretization, whereas in the latter, the forcing is applied (explicitly or implicitly) to the discretized equations. Under this taxonomy, Peskin's original method is a continuous forcing method whereas Cartesian grid, cut-cell and the ghost-fluid methods are discrete forcing methods.

For simulations of viscoelastic fluids, curved fluid interfaces, microscopic biophysical systems (proteins in lipid bilayer membranes, swimmers), and engineered devices, further variants and extensions of the immersed boundary method have been developed, such as the Stochastic Immersed Boundary Methods of Atzberger, Kramer, and Peskin [3] [4] , Stochastic Eulerian Lagrangian Methods of Atzberger [5] [6] [7] , Massed Immersed Boundary Methods of Moria [8] , and Rotational Immersed Boundary Methods of Olson, Lim, Cortez [9] .

In general, for immersed boundary methods and related variants, there is an active research community that is still developing new techniques and related software implementations and incorporating related techniques into simulation packages and CAD engineering software.

See also

Software: Numerical codes

Notes

  1. Peskin, Charles S (1972-10-01). "Flow patterns around heart valves: A numerical method". Journal of Computational Physics. 10 (2): 252–271. Bibcode:1972JCoPh..10..252P. doi:10.1016/0021-9991(72)90065-4. ISSN   0021-9991.
  2. Mittal & Iaccarino 2005.
  3. Atzberger, Paul J. (2011). "Stochastic Eulerian Lagrangian Methods for Fluid Structure Interactions with Thermal Fluctuations". Journal of Computational Physics. 230 (8): 2821–2837. arXiv: 1009.5648 . Bibcode:2011JCoPh.230.2821A. doi:10.1016/j.jcp.2010.12.028. S2CID   6067032.
  4. Rower, David A.; Padidar, Misha; Atzberger, Paul J. (April 2022). "Surface fluctuating hydrodynamics methods for the drift-diffusion dynamics of particles and microstructures within curved fluid interfaces". Journal of Computational Physics. 455: 110994. arXiv: 1906.01146 . doi:10.1016/j.jcp.2022.110994.
  5. Atzberger, Paul J. (2011). "Stochastic Eulerian Lagrangian Methods for Fluid Structure Interactions with Thermal Fluctuations". Journal of Computational Physics. 230 (8): 2821–2837. arXiv: 1009.5648 . Bibcode:2011JCoPh.230.2821A. doi:10.1016/j.jcp.2010.12.028. S2CID   6067032.
  6. Atzberger, Paul (2013), "Incorporating Shear into Stochastic Eulerian Lagrangian Methods for Rheological Studies of Complex Fluids and Soft Materials", Physica D, 265: 57–70, arXiv: 2212.10651 , doi:10.1016/j.physd.2013.09.002
  7. Atzberger, Paul (2016). "Hydrodynamic Coupling of Particle Inclusions Embedded in Curved Lipid Bilayer Membranes". Soft Matter, The Royal Society of Chemistry. 12: 6685–6707. arXiv: 1601.06461 . doi:10.1039/C6SM00194G.
  8. Moria, Yoichiro; Peskin, Charles S. (2008). "Implicit Second-Order Immersed Boundary Methods with Boundary Mass". Computer Methods in Applied Mechanics and Engineering. 197 (25–28): 2049–2067. Bibcode:2008CMAME.197.2049M. doi:10.1016/j.cma.2007.05.028.
  9. Olson, S.; Lim, S.; Cortez, R. (2013). "Modeling the dynamics of an elastic rod with intrinsic curvature and twist using a regularized Stokes formulation". Journal of Computational Physics. 238. doi:10.1016/j.jcp.2012.12.026.

Related Research Articles

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

<span class="mw-page-title-main">Particle-in-cell</span> Mathematical technique used to solve a certain class of partial differential equations

In plasma physics, the particle-in-cell (PIC) method refers to a technique used to solve a certain class of partial differential equations. In this method, individual particles in a Lagrangian frame are tracked in continuous phase space, whereas moments of the distribution such as densities and currents are computed simultaneously on Eulerian (stationary) mesh points.

<span class="mw-page-title-main">Smoothed-particle hydrodynamics</span> Method of hydrodynamics simulation

Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan and Lucy in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method, and the resolution of the method can easily be adjusted with respect to variables such as density.

<span class="mw-page-title-main">Lagrangian and Eulerian specification of the flow field</span> Computational fluid dynamics tools

In classical field theories, the Lagrangian specification of the flow field is a way of looking at fluid motion where the observer follows an individual fluid parcel as it moves through space and time. Plotting the position of an individual parcel through time gives the pathline of the parcel. This can be visualized as sitting in a boat and drifting down a river.

<span class="mw-page-title-main">Fluid–structure interaction</span>

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow. Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat.

<span class="mw-page-title-main">Meshfree methods</span> Methods in numerical analysis not requiring knowledge of neighboring points

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

A stochastic simulation is a simulation of a system that has variables that can change stochastically (randomly) with individual probabilities.

<span class="mw-page-title-main">Bram van Leer</span> Dutch mathematician

Bram van Leer is Arthur B. Modine Emeritus Professor of aerospace engineering at the University of Michigan, in Ann Arbor. He specializes in Computational fluid dynamics (CFD), fluid dynamics, and numerical analysis. His most influential work lies in CFD, a field he helped modernize from 1970 onwards. An appraisal of his early work has been given by C. Hirsch (1979)

Stokesian dynamics is a solution technique for the Langevin equation, which is the relevant form of Newton's 2nd law for a Brownian particle. The method treats the suspended particles in a discrete sense while the continuum approximation remains valid for the surrounding fluid, i.e., the suspended particles are generally assumed to be significantly larger than the molecules of the solvent. The particles then interact through hydrodynamic forces transmitted via the continuum fluid, and when the particle Reynolds number is small, these forces are determined through the linear Stokes equations. In addition, the method can also resolve non-hydrodynamic forces, such as Brownian forces, arising from the fluctuating motion of the fluid, and interparticle or external forces. Stokesian Dynamics can thus be applied to a variety of problems, including sedimentation, diffusion and rheology, and it aims to provide the same level of understanding for multiphase particulate systems as molecular dynamics does for statistical properties of matter. For rigid particles of radius suspended in an incompressible Newtonian fluid of viscosity and density , the motion of the fluid is governed by the Navier–Stokes equations, while the motion of the particles is described by the coupled equation of motion:

<span class="mw-page-title-main">Volume of fluid method</span> Free-surface modelling technique

In computational fluid dynamics, the volume of fluid (VOF) method is a family of free-surface modelling techniques, i.e. numerical techniques for tracking and locating the free surface. They belong to the class of Eulerian methods which are characterized by a mesh that is either stationary or is moving in a certain prescribed manner to accommodate the evolving shape of the interface. As such, VOF methods are advection schemes capturing the shape and position of the interface, but are not standalone flow solving algorithms. The Navier–Stokes equations describing the motion of the flow have to be solved separately.

A phase-field model is a mathematical model for solving interfacial problems. It has mainly been applied to solidification dynamics, but it has also been applied to other situations such as viscous fingering, fracture mechanics, hydrogen embrittlement, and vesicle dynamics.

Vorticity confinement (VC), a physics-based computational fluid dynamics model analogous to shock capturing methods, was invented by Dr. John Steinhoff, professor at the University of Tennessee Space Institute, in the late 1980s to solve vortex dominated flows. It was first formulated to capture concentrated vortices shed from the wings, and later became popular in a wide range of research areas. During the 1990s and 2000s, it became widely used in the field of engineering.

Particle-laden flows refers to a class of two-phase fluid flow, in which one of the phases is continuously connected and the other phase is made up of small, immiscible, and typically dilute particles. Fine aerosol particles in air is an example of a particle-laden flow; the aerosols are the dispersed phase, and the air is the carrier phase.

In computational fluid dynamics, the Stochastic Eulerian Lagrangian Method (SELM) is an approach to capture essential features of fluid-structure interactions subject to thermal fluctuations while introducing approximations which facilitate analysis and the development of tractable numerical methods. SELM is a hybrid approach utilizing an Eulerian description for the continuum hydrodynamic fields and a Lagrangian description for elastic structures. Thermal fluctuations are introduced through stochastic driving fields. Approaches also are introduced for the stochastic fields of the SPDEs to obtain numerical methods taking into account the numerical discretization artifacts to maintain statistical principles, such as fluctuation-dissipation balance and other properties in statistical mechanics.

Model order reduction (MOR) is a technique for reducing the computational complexity of mathematical models in numerical simulations. As such it is closely related to the concept of metamodeling, with applications in all areas of mathematical modelling.

Quadrature-based moment methods (QBMM) are a class of computational fluid dynamics (CFD) methods for solving Kinetic theory and is optimal for simulating phases such as rarefied gases or dispersed phases of a multiphase flow. The smallest "particle" entities which are tracked may be molecules of a single phase or granular "particles" such as aerosols, droplets, bubbles, precipitates, powders, dust, soot, etc. Moments of the Boltzmann equation are solved to predict the phase behavior as a continuous (Eulerian) medium, and is applicable for arbitrary Knudsen number and arbitrary Stokes number . Source terms for collision models such as Bhatnagar-Gross-Krook (BGK) and models for evaporation, coalescence, breakage, and aggregation are also available. By retaining a quadrature approximation of a probability density function (PDF), a set of abscissas and weights retain the physical solution and allow for the construction of moments that generate a set of partial differential equations (PDE's). QBMM has shown promising preliminary results for modeling granular gases or dispersed phases within carrier fluids and offers an alternative to Lagrangian methods such as Discrete Particle Simulation (DPS). The Lattice Boltzmann Method (LBM) shares some strong similarities in concept, but it relies on fixed abscissas whereas quadrature-based methods are more adaptive. Additionally, the Navier–Stokes equations(N-S) can be derived from the moment method approach.

<span class="mw-page-title-main">Numerical modeling (geology)</span> Technique to solve geological problems by computational simulation

In geology, numerical modeling is a widely applied technique to tackle complex geological problems by computational simulation of geological scenarios.

Hybrid stochastic simulations are a sub-class of stochastic simulations. These simulations combine existing stochastic simulations with other stochastic simulations or algorithms. Generally they are used for physics and physics-related research. The goal of a hybrid stochastic simulation varies based on context, however they typically aim to either improve accuracy or reduce computational complexity. The first hybrid stochastic simulation was developed in 1985.

Lagrangian ocean analysis is a way of analysing ocean dynamics by computing the trajectories of virtual fluid particles, following the Lagrangian perspective of fluid flow, from a specified velocity field. Often, the Eulerian velocity field used as an input for Lagrangian ocean analysis has been computed using an ocean general circulation model (OGCM). Lagrangian techniques can be employed on a range of scales, from modelling the dispersal of biological matter within the Great Barrier Reef to global scales. Lagrangian ocean analysis has numerous applications, from modelling the diffusion of tracers, through the dispersal of aircraft debris and plastics, to determining the biological connectivity of ocean regions.

<span class="mw-page-title-main">Physics-informed neural networks</span> Technique to solve partial differential equations

Physics-informed neural networks (PINNs), also referred to as Theory-Trained Neural Networks (TTNs), are a type of universal function approximators that can embed the knowledge of any physical laws that govern a given data-set in the learning process, and can be described by partial differential equations (PDEs).They overcome the low data availability of some biological and engineering systems that makes most state-of-the-art machine learning techniques lack robustness, rendering them ineffective in these scenarios. The prior knowledge of general physical laws acts in the training of neural networks (NNs) as a regularization agent that limits the space of admissible solutions, increasing the correctness of the function approximation. This way, embedding this prior information into a neural network results in enhancing the information content of the available data, facilitating the learning algorithm to capture the right solution and to generalize well even with a low amount of training examples.

References