Hybrid finite difference/finite element immersed boundary method

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach employs a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach

The Smooth Forcing Extension Method: A High-Order Technique for Solving Elliptic Equations on Complex Domains

High-order numerical methods for solving elliptic equations over arbitrary domains typically require specialized machinery, such as high-quality conforming grids for finite elements method, and quadrature rules for boundary integral methods. These tools make it difficult to apply these techniques to higher dimensions. In contrast, fixed Cartesian grid methods, such as the immersed boundary (IB) method, are easy to apply and generalize, but typically are low-order accurate. In this study, we introduce the Smooth Forcing Extension (SFE) method, a fixed Cartesian grid technique that builds on the insights of the IB method, and allows one to obtain arbitrary orders of accuracy. Our approach relies on a novel Fourier continuation method to compute extensions of the inhomogeneous terms to any desired regularity. This is combined with the highly accurate Non-Uniform Fast Fourier Transform for interpolation operations to yield a fast and robust method. Numerical tests confirm that the technique performs precisely as expected on one-dimensional test problems. In higher dimensions, the performance is even better, in some cases yielding sub-geometric convergence. We also demonstrate how this technique can be applied to solving parabolic problems and for computing the eigenvalues of elliptic operators on general domains, in the process illustrating its stability and amenability to generalization

A sharp interface method for an immersed viscoelastic solid

The immersed boundary-finite element method (IBFE) is an approach to describing the dynamics of an elastic structure immersed in an incompressible viscous fluid. In this formulation, there are discontinuities in the pressure and viscous stress at fluid-structure interfaces. The standard immersed boundary approach, which connects the Lagrangian and Eulerian variables via integral transforms with regularized Dirac delta function kernels, smooths out these discontinuities, which generally leads to low order accuracy. This paper describes an approach to accurately resolve pressure discontinuities for these types of formulations, in which the solid may undergo large deformations. Our strategy is to decompose the physical pressure field into a sum of two pressure-like fields, one defined on the entire computational domain, which includes both the fluid and solid subregions, and one defined only on the solid subregion. Each of these fields is continuous on its domain of definition, which enables high accuracy via standard discretization methods without sacrificing sharp resolution of the pressure discontinuity. Numerical tests demonstrate that this method improves rates of convergence for displacements, velocities, stresses, and pressures, as compared to the conventional IBFE method. Further, it produces much smaller errors at reasonable numbers of degrees of freedom. The performance of this method is tested on several cases with analytic solutions, a nontrivial benchmark problem of incompressible solid mechanics, and an example involving a thick, actively contracting torus

Dynamic finite-strain modelling of the human left ventricle in health and disease using an immersed boundary-finite element method

Detailed models of the biomechanics of the heart are important both for developing improved interventions for patients with heart disease and also for patient risk stratification and treatment planning. For instance, stress distributions in the heart affect cardiac remodelling, but such distributions are not presently accessible in patients. Biomechanical models of the heart offer detailed three-dimensional deformation, stress and strain fields that can supplement conventional clinical data. In this work, we introduce dynamic computational models of the human left ventricle (LV) that are derived from clinical imaging data obtained from a healthy subject and from a patient with a myocardial infarction (MI). Both models incorporate a detailed invariant-based orthotropic description of the passive elasticity of the ventricular myocardium along with a detailed biophysical model of active tension generation in the ventricular muscle. These constitutive models are employed within a dynamic simulation framework that accounts for the inertia of the ventricular muscle and the blood that is based on an immersed boundary (IB) method with a finite element description of the structural mechanics. The geometry of the models is based on data obtained non-invasively by cardiac magnetic resonance (CMR). CMR imaging data are also used to estimate the parameters of the passive and active constitutive models, which are determined so that the simulated end-diastolic and end-systolic volumes agree with the corresponding volumes determined from the CMR imaging studies. Using these models, we simulate LV dynamics from end-diastole to end-systole. The results of our simulations are shown to be in good agreement with subject-specific CMR-derived strain measurements and also with earlier clinical studies on human LV strain distributions

Immersed Methods for Fluid-Structure Interaction

Fluid-structure interaction is ubiquitous in nature and occurs at all biological scales. Immersed methods provide mathematical and computational frameworks for modeling fluid-structure systems. These methods, which typically use an Eulerian description of the fluid and a Lagrangian description of the structure, can treat thin immersed boundaries and volumetric bodies, and they can model structures that are flexible or rigid or that move with prescribed deformational kinematics. Immersed formulations do not require body-fitted discretizations and thereby avoid the frequent grid regeneration that can otherwise be required for models involving large deformations and displacements. This article reviews immersed methods for both elastic structures and structures with prescribed kinematics. It considers formulations using integral operators to connect the Eulerian and Lagrangian frames and methods that directly apply jump conditions along fluid-structure interfaces. Benchmark problems demonstrate the effectiveness of these methods, and selected applications at Reynolds numbers up to approximately 20,000 highlight their impact in biological and biomedical modeling and simulation

A coupled mitral valve - left ventricle model with fluid-structure interaction

Understanding the interaction between the valves and walls of the heart is important in assessing and subsequently treating heart dysfunction. This study presents an integrated model of the mitral valve (MV) coupled to the left ventricle (LV), with the geometry derived from in vivo clinical magnetic resonance images. Numerical simulations using this coupled MV–LV model are developed using an immersed boundary/finite element method. The model incorporates detailed valvular features, left ventricular contraction, nonlinear soft tissue mechanics, and fluid-mediated interactions between the MV and LV wall. We use the model to simulate cardiac function from diastole to systole. Numerically predicted LV pump function agrees well with in vivo data of the imaged healthy volunteer, including the peak aortic flow rate, the systolic ejection duration, and the LV ejection fraction. In vivo MV dynamics are qualitatively captured. We further demonstrate that the diastolic filling pressure increases significantly with impaired myocardial active relaxation to maintain a normal cardiac output. This is consistent with clinical observations. The coupled model has the potential to advance our fundamental knowledge of mechanisms underlying MV–LV interaction, and help in risk stratification and optimisation of therapies for heart diseases

An immersed interface-lattice Boltzmann method for fluid-structure interaction

An immersed interface-lattice Boltzmann method (II-LBM) is developed for modelling fluid-structure systems. The key element of this approach is the determination of the jump conditions that are satisfied by the distribution functions within the framework of the lattice Boltzmann method when forces are imposed along a surface immersed in an incompressible fluid. In this initial II-LBM, the discontinuity related to the normal portion of the interfacial force is sharply resolved by imposing the relevant jump conditions using an approach that is analogous to imposing the corresponding pressure jump condition in the incompressible Navier-Stokes equations. We show that the jump conditions for the distribution functions are the same in both single-relaxation-time and multi-relaxation-time LBM formulations. Tangential forces are treated using the immersed boundary-lattice Boltzmann method (IB-LBM). The performance of the II-LBM method is compared to both the direct forcing IB-LBM for rigid-body fluid-structure interaction, and the classical IB-LBM for elastic interfaces. Higher order accuracy is observed with the II-LBM as compared to the IB-LBM for selected benchmark problems. Because the jump conditions of the distribution function also satisfy the continuity of the velocity field across the interface, the error in the velocity field is much smaller for the II-LBM than the IB-LBM. The II-LBM is also demonstrated to provide superior volume conservation when simulating flexible boundaries