95 research outputs found
Fluid-structure interaction simulation of prosthetic aortic valves : comparison between immersed boundary and arbitrary Lagrangian-Eulerian techniques for the mesh representation
In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations' outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results
Immersed boundary-finite element model of fluid-structure interaction in the aortic root
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
Numerical Investigation of the Performance of Three Hinge Designs of Bileaflet Mechanical Heart Valves
Thromboembolic complications (TECs) of bileaflet mechanical heart valves (BMHVs) are believed to be due to the nonphysiologic mechanical stresses imposed on blood elements by the hinge flows. Relating hinge flow features to design features is, therefore, essential to ultimately design BMHVs with lower TEC rates. This study aims at simulating the pulsatile three-dimensional hinge flows of three BMHVs and estimating the TEC potential associated with each hinge design. Hinge geometries are constructed from micro-computed tomography scans of BMHVs. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Leaflet motion and flow boundary conditions are extracted from fluid–structure-interaction simulations of BMHV bulk flow. The numerical results are analyzed using a particle-tracking approach coupled with existing blood damage models. The gap width and, more importantly, the shape of the recess and leaflet are found to impact the flow distribution and TEC potential. Smooth, streamlined surfaces appear to be more favorable than sharp corners or sudden shape transitions. The developed framework will enable pragmatic and cost-efficient preclinical evaluation of BMHV prototypes prior to valve manufacturing. Application to a wide range of hinges with varying design parameters will eventually help in determining the optimal hinge design
Simulation of the Three-Dimensional Hinge Flow Fields of a Bileaflet Mechanical Heart Valve Under Aortic Conditions
Thromboembolic complications of bileaflet mechanical heart valves (BMHV) are believed to be due to detrimental stresses imposed on blood elements by the hinge flows. Characterization of these flows is thus crucial to identify the underlying causes for complications. In this study, we conduct three-dimensional pulsatile flow simulations through the hinge of a BMHV under aortic conditions. Hinge and leaflet geometries are reconstructed from the Micro-Computed Tomography scans of a BMHV. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Physiologic flow boundary conditions and leaflet motion are extracted from the Fluid–Structure Interaction simulations of the bulk of the flow through a BMHV. Calculations reveal the presence, throughout the cardiac cycle, of flow patterns known to be detrimental to blood elements. Flow fields are characterized by: (1) complex systolic flows, with rotating structures and slow reverse flow pattern, and (2) two strong diastolic leakage jets accompanied by fast reverse flow at the hinge bottom. Elevated shear stresses, up to 1920 dyn/cm2 during systole and 6115 dyn/cm2 during diastole, are reported. This study underscores the need to conduct three-dimensional simulations throughout the cardiac cycle to fully characterize the complexity and thromboembolic potential of the hinge flows
The fish tail motion forms an attached leading edge vortex
We propose an approach to falsification of oscillation properties of
parametric biological models, based on the recently developed techniques for
testing continuous and hybrid systems. In this approach, an oscillation
property can be specified using a hybrid automaton, which is then used to guide
the exploration in the state and input spaces to search for the behaviors that
do not satisfy the property. We illustrate the approach with the Laub-Loomis
model for spontaneous oscillations during the aggregation stage of
Dictyostelium.Comment: In Proceedings HSB 2013, arXiv:1308.572
Breaking the symmetry of a wavy channel alters the route to chaotic flow
We numerically explore the two-dimensional, incompressible, isothermal flow through a wavy channel, with a focus on how the channel geometry affects the routes to chaos at Reynolds numbers between 150 and 1000. We find that (i) the period-doubling route arises in a symmetric channel, (ii) the Ruelle-Takens-Newhouse route arises in an asymmetric channel, and (iii) the type-II intermittency route arises in both asymmetric and semiwavy channels. We also find that the flow through the semiwavy channel evolves from a quasiperiodic torus to an unstable invariant set (chaotic saddle), before eventually settling on a period-1 limit-cycle attractor. This study reveals that laminar channel flow at elevated Reynolds numbers can exhibit a variety of nonlinear dynamics. Specifically, it highlights how breaking the symmetry of a wavy channel can not only influence the critical Reynolds number at which chaos emerges, but also diversify the types of bifurcation encountered en route to chaos itself
- …