Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries

The simulation of blood flow and pressure in arteries requires outflow boundary conditions that incorporate models of downstream domains. We previously described a coupled multidomain method to couple analytical models of the downstream domains with 3D numerical models of the upstream vasculature. This prior work either included pure resistance boundary conditions or impedance boundary conditions based on assumed periodicity of the solution. However, flow and pressure in arteries are not necessarily periodic in time due to heart rate variability, respiration, complex transitional flow or acute physiological changes. We present herein an approach for prescribing lumped parameter outflow boundary conditions that accommodate transient phenomena. We have applied this method to compute haemodynamic quantities in different physiologically relevant cardiovascular models, including patient-specific examples, to study non-periodic flow phenomena often observed in normal subjects and in patients with acquired or congenital cardiovascular disease. The relevance of using boundary conditions that accommodate transient phenomena compared with boundary conditions that assume periodicity of the solution is discussed

A finite strain nonlinear human mitral valve model with fluid structure interaction

A simulated human mitral valve under a physiological pressure loading is developed using a hybrid finite element immersed boundary method, which incorporates experimentally based constitutive laws in a three-dimensional fluid-structure interaction framework. A transversely isotropic material constitutive model is used for characterizing the mechanical behaviour of the mitral valve tissue based on recent mechanical tests of healthy human mitral leaflets. Our results show good agreement, in terms of the flow rate and the closing and opening configurations, with the measurements from the magnetic resonance images. The stresses in the anterior leaflet are found to be higher than those in the posterior leaflet, and concentrated around the annulus trigons and free edges of the valve leaflets. Those areas are located where the leaflet has the highest curvature. Effects of the chordae tendineae in the material model are studied and the results show that these chordae play an important role in providing a secondary orifice for the flow when valve opens. Although there are some discrepancies to be overcome in future works, our simulations show that the developed computational model is promising in mimicking the in vivo mitral valve dynamics and providing important information that are not obtainable by in vivo measurements. This article is protected by copyright. All rights reserved

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. With advancements in cardiac imaging, nonlinear mechanics and computational techniques, it is now possible to explore the mechanics of valve-heart interactions using anatomically and physiologically realistic models. 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 the cardiac function from diastole to systole, and investigate how myocardial active relaxation function affects the LV pump function. The results of the new model agree with in vivo measurements, and demonstrate that the diastolic filling pressure increases significantly with impaired myocardial active relaxation to maintain the normal cardiac output. The coupled model has the potential to advance fundamental knowledge of mechanisms underlying MV-LV interaction, and help in risk stratification and optimization of therapies for heart diseases.Comment: 25 pages, 6 figure

Estimating prognosis in patients with acute myocardial infarction using personalized computational heart models

Biomechanical computational models have potential prognostic utility in patients after an acute ST-segment–elevation myocardial infarction (STEMI). In a proof-of-concept study, we defined two groups (1) an acute STEMI group (n = 6, 83% male, age 54 ± 12 years) complicated by left ventricular (LV) systolic dysfunction; (2) an age- and sex- matched hyper-control group (n = 6, 83% male, age 46 ± 14 years), no prior history of cardiovascular disease and normal systolic blood pressure (SBP < 130 mmHg). Cardiac MRI was performed in the patients (2 days & 6 months post-STEMI) and the volunteers, and biomechanical heart models were synthesized for each subject. The candidate parameters included normalized active tension (ATnorm) and active tension at the resting sarcomere length (Treq, reflecting required contractility). Myocardial contractility was inversely determined from personalized heart models by matching CMR-imaged LV dynamics. Compared with controls, patients with recent STEMI exhibited increased LV wall active tension when normalized by SBP. We observed a linear relationship between Treq 2 days post-MI and global longitudinal strain 6 months later (r = 0.86; p = 0.03). Treq may be associated with changes in LV function in the longer term in STEMI patients complicated by LV dysfunction. Further studies seem warranted

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

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

Modelling the evolution of cerebral aneurysms: biomechanics, mechanobiology and multiscale modelling

Intracranial aneurysms (IAs) are abnormal dilatations of the cerebral vasculature. Computational modelling may shed light on the aetiology of the disease and lead to improved criteria to assist diagnostic decisions. We briefly review models of aneurysm evolution to date and present a novel fluid-solid-growth (FSG) framework for patient-specific modelling of IA evolution. We illustrate its application to 4 clinical cases depicting an IA. The section of arterial geometry containing the IA is removed and replaced with a cylindrical section: this represents an idealised section of healthy artery upon which IA evolution is simulated. The utilisation of patient-specific geometries enables G&R to be explicitly linked to physiologically realistic spatial distributions and magnitudes of haemodynamic stimuli. In this study, we investigate the hypothesis that elastin degradation is driven by locally low wall shear stress (WSS). In 3 out of 4 cases, the evolved model IA geometry is qualitatively similar to the corresponding in vivo IA geometry. This suggests some tentative support for the hypothesis that low WSS plays a role in the mechanobiology of IA evolution

A Rapid and Computationally Inexpensive Method to Virtually Implant Current and Next-Generation Stents into Subject-Specific Computational Fluid Dynamics Models

Computational modeling is often used to quantify hemodynamic alterations induced by stenting, but frequently uses simplified device or vascular representations. Based on a series of Boolean operations, we developed an efficient and robust method for assessing the influence of current and next-generation stents on local hemodynamics and vascular biomechanics quantified by computational fluid dynamics. Stent designs were parameterized to allow easy control over design features including the number, width and circumferential or longitudinal spacing of struts, as well as the implantation diameter and overall length. The approach allowed stents to be automatically regenerated for rapid analysis of the contribution of design features to resulting hemodynamic alterations. The applicability of the method was demonstrated with patient-specific models of a stented coronary artery bifurcation and basilar trunk aneurysm constructed from medical imaging data. In the coronary bifurcation, we analyzed the hemodynamic difference between closed-cell and open-cell stent geometries. We investigated the impact of decreased strut size in stents with a constant porosity for increasing flow stasis within the stented basilar aneurysm model. These examples demonstrate the current method can be used to investigate differences in stent performance in complex vascular beds for a variety of stenting procedures and clinical scenarios

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