An Immersed Boundary method with divergence-free velocity interpolation and force spreading

Pressure boundary conditions for immersed-boundary methods

Immersed boundary methods have seen an enormous increase in popularity over the past two decades, especially for problems involving complex moving/deforming boundaries. In most cases, the boundary conditions on the immersed body are enforced via forcing functions in the momentum equations, which in the case of fractional step methods may be problematic due to: i) creation of slip-errors resulting from the lack of explicitly enforcing boundary conditions on the (pseudo-)pressure on the immersed body; ii) coupling of the solution in the fluid and solid domains via the Poisson equation. Examples of fractional-step formulations that simultaneously enforce velocity and pressure boundary conditions have also been developed, but in most cases the standard Poisson equation is replaced by a more complex system which requires expensive iterative solvers. In this work we propose a new formulation to enforce appropriate boundary conditions on the pseudo-pressure as part of a fractional-step approach. The overall treatment is inspired by the ghost-fluid method typically utilized in two-phase flows. The main advantage of the algorithm is that a standard Poisson equation is solved, with all the modifications needed to enforce the boundary conditions being incorporated within the right-hand side. As a result, fast solvers based on trigonometric transformations can be utilized. We demonstrate the accuracy and robustness of the formulation for a series of problems with increasing complexity.

On the Lagrangian-Eulerian Coupling in the Immersed Finite Element/Difference Method

The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure system and a Lagrangian description of the deformations, stresses, and resultant forces of the immersed structure. Integral transforms with Dirac delta function kernels couple the Eulerian and Lagrangian variables, and in practice, discretizations of these integral transforms use regularized delta function kernels. Many different kernel functions have been proposed, but prior numerical work investigating the impact of the choice of kernel function on the accuracy of the methodology has often been limited to simplified test cases or Stokes flow conditions that may not reflect the method's performance in applications, particularly at intermediate-to-high Reynolds numbers, or under different loading conditions. This work systematically studies the effect of the choice of regularized delta function in several fluid-structure interaction benchmark tests using the immersed finite element/difference (IFED) method, which is an extension of the IB method that uses a finite element structural discretization combined with a Cartesian grid finite difference method for the incompressible Navier-Stokes equations. Whereas the conventional IB method spreads forces from the nodes of the structural mesh and interpolates velocities to those nodes, the IFED formulation evaluates the regularized delta function on a collection of interaction points that can be chosen to be denser than the nodes of the Lagrangian mesh. This opens the possibility of using structural discretizations with wide node spacings that would produce gaps in the Eulerian force in nodally coupled schemes (e.g., if the node spacing is comparable to or broader than the support of the regularized delta functions). Earlier work with this methodology suggested that such coarse structural meshes can yield improved accuracy for shear-dominated cases and, further, found that accuracy improves when the structural mesh spacing is increased. However, these results were limited to simple test cases that did not include substantial pressure loading on the structure. This study investigates the effect of varying the relative mesh widths of the Lagrangian and Eulerian discretizations in a broader range of tests. Our results indicate that kernels satisfying a commonly imposed even-odd condition require higher resolution to achieve similar accuracy as kernels that do not satisfy this condition. We also find that narrower kernels are more robust, in the sense that they yield results that are less sensitive to relative changes in the Eulerian and Lagrangian mesh spacings, and that structural meshes that are substantially coarser than the Cartesian grid can yield high accuracy for shear-dominated cases but not for cases with large normal forces. We verify our results in a large-scale FSI model of a bovine pericardial bioprosthetic heart valve in a pulse duplicator.

A fully resolved multiphysics model of gastric peristalsis and bolus emptying in the upper gastrointestinal tract

Over the past few decades, in silico modeling of organ systems has significantly furthered our understanding of their physiology and biomechanical function. In spite of the relative importance of the digestive system in normal functioning of the human body, there is a scarcity of high-fidelity models for the upper gastrointestinal tract including the esophagus and the stomach. In this work, we present a detailed numerical model of the upper gastrointestinal tract that not only accounts for the fiber architecture of the muscle walls, but also the multiphasic components they help transport during normal digestive function. Construction details for 3D models of representative stomach geometry are presented along with a simple strategy for assigning circular and longitudinal muscle fiber orientations for each layer. We developed a fully resolved model of the stomach to simulate gastric peristalsis by systematically activating muscle fibers embedded in the stomach. Following this, for the first time, we simulate gravity-driven bolus emptying into the stomach due to density differences between ingested contents and fluid contents of the stomach. Finally, we present a case of retrograde flow of fluid from the stomach into the esophagus, resembling the phenomenon of acid reflux. This detailed computational model of the upper gastrointestinal tract provides a foundation for future models to investigate the biomechanics of acid reflux and probe various strategies for gastric bypass surgeries to address the growing problem of obesity.

A novel interpolation-free sharp-interface immersed boundary method

This paper describes a novel 2^nd order direct forcing immersed boundary method designed for simulation of 2D and 3D incompressible flow problems with complex immersed boundaries. In this formulation, each cell cut by the immersed boundary (IB) is reshaped to conform to the shape of the IB. IBs are modeled as a series of 2D planes in 3D space that connect seamlessly at the edges of the cut cells, in a way that mimics conformal grid. IBs are represented in a continuous and consistent fashion from one cell to another, thus eliminating spatial pressure oscillations originating from inconsistent description of the IB as well as the traditional stair-step problem, leading to a more accurate resolution of the boundary layer. Boundary conditions are enforced at the exact location of the IB devoid of interpolation, which guarantees sound simulations even on grids with high aspect ratio, and enables simulations of flow packed with multiple IBs in close proximity. Boundary conditions for each phase across the IB are enforced independently, yielding a unique capability to solve flows with zero-thickness IBs. Simulations of a large number of 2D and 3D test cases confirm the prowess of the devised immersed boundary method in solving flows over multiple loosely/closely-packed IBs; stationary, moving and highly morphing IBs; as well as IBs with zero-thickness. Results show that predictions from the proposed scheme agree remarkably well with theoretical models and experimental data for both integral variables and local flow fields and they are often with less than 1% deviation from solutions obtained by conformal grid of similar resolution.

Projection-based model reduction for the immersed boundary method

Fluid-structure interactions are central to many biomolecular processes, and they impose a great challenge for computational and modeling methods. In this paper, we consider the immersed boundary method (IBM) for biofluid systems, and to alleviate the computational cost, we apply reduced-order techniques to eliminate the degrees of freedom associated with the large number of fluid variables. We show how reduced models can be derived using Petrov-Galerkin projection and subspaces that maintain the incompressibility condition. More importantly, the reduced-order model (ROM) is shown to preserve the Lyapunov stability. We also address the practical issue of computing coefficient matrices in the ROM using an interpolation technique. The efficiency and robustness of the proposed formulation are examined with test examples from various applications.

Stabilization approaches for the hyperelastic immersed boundary method for problems of large-deformation incompressible elasticity

The immersed boundary method is a mathematical framework for modeling fluid-structure interaction. This formulation describes the momentum, viscosity, and incompressibility of the fluid-structure system in Eulerian form, and it uses Lagrangian coordinates to describe the structural deformations, stresses, and resultant forces. Integral transforms with Dirac delta function kernels connect the Eulerian and Lagrangian frames. The fluid and the structure are both typically treated as incompressible materials. Upon discretization, however, the incompressibility of the structure is only maintained approximately. To obtain an immersed method for incompressible hyperelastic structures that is robust under large structural deformations, we introduce a volumetric energy in the solid region that stabilizes the formulation and improves the accuracy of the numerical scheme. This formulation augments the discrete Lagrange multiplier for the incompressibility constraint, thereby improving the original method's accuracy. This volumetric energy is incorporated by decomposing the strain energy into isochoric and dilatational components, as in standard solid mechanics formulations of nearly incompressible elasticity. We study the performance of the stabilized method using several quasi-static solid mechanics benchmarks, a dynamic fluid-structure interaction benchmark, and a detailed three-dimensional model of esophageal transport. The accuracy achieved by the stabilized immersed formulation is comparable to that of a stabilized finite element method for incompressible elasticity using similar numbers of structural degrees of freedom.

Modeling the mitral valve

This work is concerned with modeling and simulation of the mitral valve, one of the four valves in the human heart. The valve is composed of leaflets, the free edges of which are supported by a system of chordae, which themselves are anchored to the papillary muscles inside the left ventricle. First, we examine valve anatomy and present the results of original dissections. These display the gross anatomy and information on fiber structure of the mitral valve. Next, we build a model valve following a design-based methodology, meaning that we derive the model geometry and the forces that are needed to support a given load and construct the model accordingly. We incorporate information from the dissections to specify the fiber topology of this model. We assume the valve achieves mechanical equilibrium while supporting a static pressure load. The solution to the resulting differential equations determines the pressurized configuration of the valve model. To complete the model, we then specify a constitutive law based on a stress-strain relation consistent with experimental data that achieves the necessary forces computed in previous steps. Finally, using the immersed boundary method, we simulate the model valve in fluid in a computer test chamber. The model opens easily and closes without leak when driven by physiological pressures over multiple beats. Further, its closure is robust to driving pressures that lack atrial systole or are much lower or higher than normal.