4,531 research outputs found

    Immersed boundary-finite element model of fluid-structure interaction in the aortic root

    Get PDF
    It has long been recognized that aortic root elasticity helps to ensure efficient aortic valve closure, but our understanding of the functional importance of the elasticity and geometry of the aortic root continues to evolve as increasingly detailed in vivo imaging data become available. Herein, we describe fluid-structure interaction models of the aortic root, including the aortic valve leaflets, the sinuses of Valsalva, the aortic annulus, and the sinotubular junction, that employ a version of Peskin's immersed boundary (IB) method with a finite element (FE) description of the structural elasticity. We develop both an idealized model of the root with three-fold symmetry of the aortic sinuses and valve leaflets, and a more realistic model that accounts for the differences in the sizes of the left, right, and noncoronary sinuses and corresponding valve cusps. As in earlier work, we use fiber-based models of the valve leaflets, but this study extends earlier IB models of the aortic root by employing incompressible hyperelastic models of the mechanics of the sinuses and ascending aorta using a constitutive law fit to experimental data from human aortic root tissue. In vivo pressure loading is accounted for by a backwards displacement method that determines the unloaded configurations of the root models. Our models yield realistic cardiac output at physiological pressures, with low transvalvular pressure differences during forward flow, minimal regurgitation during valve closure, and realistic pressure loads when the valve is closed during diastole. Further, results from high-resolution computations demonstrate that IB models of the aortic valve are able to produce essentially grid-converged dynamics at practical grid spacings for the high-Reynolds number flows of the aortic root

    Hybrid finite difference/finite element immersed boundary method

    Get PDF
    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

    Quasi-static imaged-based immersed boundary-finite element model of human left ventricle in diastole

    Get PDF
    SUMMARY: Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid–structure interaction simulations of the heart

    Three embedded techniques for finite element heat flow problem with embedded discontinuities

    Get PDF
    The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-017-1382-7The present paper explores the solution of a heat conduction problem considering discontinuities embedded within the mesh and aligned at arbitrary angles with respect to the mesh edges. Three alternative approaches are proposed as solutions to the problem. The difference between these approaches compared to alternatives, such as the eXtended Finite Element Method (X-FEM), is that the current proposal attempts to preserve the global matrix graph in order to improve performance. The first two alternatives comprise an enrichment of the Finite Element (FE) space obtained through the addition of some new local degrees of freedom to allow capturing discontinuities within the element. The new degrees of freedom are statically condensed prior to assembly, so that the graph of the final system is not changed. The third approach is based on the use of modified FE-shape functions that substitute the standard ones on the cut elements. The imposition of both Neumann and Dirichlet boundary conditions is considered at the embedded interface. The results of all the proposed methods are then compared with a reference solution obtained using the standard FE on a mesh containing the actual discontinuity.Peer ReviewedPostprint (author's final draft

    Spatially Adaptive Stochastic Methods for Fluid-Structure Interactions Subject to Thermal Fluctuations in Domains with Complex Geometries

    Full text link
    We develop stochastic mixed finite element methods for spatially adaptive simulations of fluid-structure interactions when subject to thermal fluctuations. To account for thermal fluctuations, we introduce a discrete fluctuation-dissipation balance condition to develop compatible stochastic driving fields for our discretization. We perform analysis that shows our condition is sufficient to ensure results consistent with statistical mechanics. We show the Gibbs-Boltzmann distribution is invariant under the stochastic dynamics of the semi-discretization. To generate efficiently the required stochastic driving fields, we develop a Gibbs sampler based on iterative methods and multigrid to generate fields with O(N)O(N) computational complexity. Our stochastic methods provide an alternative to uniform discretizations on periodic domains that rely on Fast Fourier Transforms. To demonstrate in practice our stochastic computational methods, we investigate within channel geometries having internal obstacles and no-slip walls how the mobility/diffusivity of particles depends on location. Our methods extend the applicability of fluctuating hydrodynamic approaches by allowing for spatially adaptive resolution of the mechanics and for domains that have complex geometries relevant in many applications

    A fully resolved active musculo-mechanical model for esophageal transport

    Full text link
    Esophageal transport is a physiological process that mechanically transports an ingested food bolus from the pharynx to the stomach via the esophagus, a multi-layered muscular tube. This process involves interactions between the bolus, the esophagus, and the neurally coordinated activation of the esophageal muscles. In this work, we use an immersed boundary (IB) approach to simulate peristaltic transport in the esophagus. The bolus is treated as a viscous fluid that is actively transported by the muscular esophagus, which is modeled as an actively contracting, fiber-reinforced tube. A simplified version of our model is verified by comparison to an analytic solution to the tube dilation problem. Three different complex models of the multi-layered esophagus, which differ in their activation patterns and the layouts of the mucosal layers, are then extensively tested. To our knowledge, these simulations are the first of their kind to incorporate the bolus, the multi-layered esophagus tube, and muscle activation into an integrated model. Consistent with experimental observations, our simulations capture the pressure peak generated by the muscle activation pulse that travels along the bolus tail. These fully resolved simulations provide new insights into roles of the mucosal layers during bolus transport. In addition, the information on pressure and the kinematics of the esophageal wall due to the coordination of muscle activation is provided, which may help relate clinical data from manometry and ultrasound images to the underlying esophageal motor function
    • …
    corecore