Daniel R. Einstein

Fluid–structure interaction models of the mitral valve: function in normal and pathological states

Successful mitral valve repair is dependent upon a full understanding of normal and abnormal mitral valve anatomy and function. Computational analysis is one such method that can be applied to simulate mitral valve function in order to analyse the roles of individual components and evaluate proposed surgical repair. We developed the first three-dimensional finite element computer model of the mitral valve including leaflets and chordae tendineae; however, one critical aspect that has been missing until the last few years was the evaluation of fluid flow, as coupled to the function of the mitral valve structure. We present here our latest results for normal function and specific pathological changes using a fluid–structure interaction model. Normal valve function was first assessed, followed by pathological material changes in collagen fibre volume fraction, fibre stiffness, fibre splay and isotropic stiffness. Leaflet and chordal stress and strain and papillary muscle force were determined. In addition, transmitral flow, time to leaflet closure and heart valve sound were assessed. Model predictions in the normal state agreed well with a wide range of available in vivo and in vitro data. Further, pathological material changes that preserved the anisotropy of the valve leaflets were found to preserve valve function. By contrast, material changes that altered the anisotropy of the valve were found to profoundly alter valve function. The addition of blood flow and an experimentally driven microstructural description of mitral tissue represent significant advances in computational studies of the mitral valve, which allow further insight to be gained. This work is another building block in the foundation of a computational framework to aid in the refinement and development of a truly non-invasive diagnostic evaluation of the mitral valve. Ultimately, it represents the basis for simulation of surgical repair of pathological valves in a clinical and educational setting.

A Computationally Efficient Formal Optimization of Regional Myocardial Contractility in a Sheep With Left Ventricular Aneurysm

A non-invasive method for estimating regional myocardial contractility in vivo would be of great value in the design and evaluation of new surgical and medical strategies to treat and/or prevent infarction-induced heart failure. As a first step towards developing such a method, an explicit finite element (FE) model-based formal optimization of regional myocardial contractility in a sheep with left ventricular (LV) aneurysm was performed using tagged magnetic resonance (MR) images and cardiac catheterization pressures. From the tagged MR images, 3-dimensional (3D) myocardial strains, LV volumes and geometry for the animal-specific 3D FE model of the LV were calculated, while the LV pressures provided physiological loading conditions. Active material parameters (T(max_B) and T(max_R)) in the non-infarcted myocardium adjacent to the aneurysm (borderzone) and in myocardium remote from the aneurysm were estimated by minimizing the errors between FE model-predicted and measured systolic strains and LV volumes using the successive response surface method for optimization. The significant depression in optimized T(max_B) relative to T(max_R) was confirmed by direct ex vivo force measurements from skinned fiber preparations. The optimized values of T(max_B) and T(max_R) were not overly sensitive to the passive material parameters specified. The computation time of less than 5 hours associated with our proposed method for estimating regional myocardial contractility in vivo makes it a potentially very useful clinical tool.

Comparative computational modeling of airflows and vapor dosimetry in the respiratory tracts of rat, monkey, and human

Computational fluid dynamics (CFD) models are useful for predicting site-specific dosimetry of airborne materials in the respiratory tract and elucidating the importance of species differences in anatomy, physiology, and breathing patterns. We improved the imaging and model development methods to the point where CFD models for the rat, monkey, and human now encompass airways from the nose or mouth to the lung. A total of 1272, 2172, and 135 pulmonary airways representing 17±7, 19±9, or 9±2 airway generations were included in the rat, monkey and human models, respectively. A CFD/physiologically based pharmacokinetic model previously developed for acrolein was adapted for these anatomically correct extended airway models. Model parameters were obtained from the literature or measured directly. Airflow and acrolein uptake patterns were determined under steady-state inhalation conditions to provide direct comparisons with prior data and nasal-only simulations. Results confirmed that regional uptake was sensitive to airway geometry, airflow rates, acrolein concentrations, air:tissue partition coefficients, tissue thickness, and the maximum rate of metabolism. Nasal extraction efficiencies were predicted to be greatest in the rat, followed by the monkey, and then the human. For both nasal and oral breathing modes in humans, higher uptake rates were predicted for lower tracheobronchial tissues than either the rat or monkey. These extended airway models provide a unique foundation for comparing material transport and site-specific tissue uptake across a significantly greater range of conducting airways in the rat, monkey, and human than prior CFD models.

Application of Magnetic Resonance (MR) Imaging for the Development and Validation of Computational Fluid Dynamic (CFD) Models of the Rat Respiratory System

Computational fluid dynamic (CFD) models of the respiratory system provide a quantitative basis for extrapolating the localized dose of inhaled materials and improving human health risk assessments based upon inhalation studies conducted in animals. Nevertheless, model development and validation have historically been tedious and time-consuming tasks. In recognition of this, we previously reported on the use of proton (1H) magnetic resonance (MR) imaging for visualizing nasal-sinus passages in the rat, and for speeding computational mesh generation. Here, the generation and refinement of meshes for rat nasal airways are described in more detail and simulated airflows are presented. To extend the CFD models to the complete respiratory tract, three-dimensional (3D) 1H MR imaging of rat pulmonary casts was also utilized to construct pulmonary airway meshes using procedures developed for the nasal airways. Furthermore, the feasibility of validating CFD predictions with MR was tested by imaging hyperpolarized 3He gas at physiological flow rates in a straight pipe with a diameter comparable to the rat trachea. Results from these diverse studies highlight the potential utility of MR imaging not only for speeding CFD development but also possibly for model validation.

Inverse parameter fitting of biological tissues: a response surface approach.

In this paper, we present the application of a semi-global inverse method for determining material parameters of biological tissues. The approach is based on the successive response surface method, and is illustrated by fitting constitutive parameters to two nonlinear anisotropic constitutive equations, one for aortic sinus and aortic wall, the other for aortic valve tissue. Material test data for the aortic sinus consisted of two independent orthogonal uniaxial tests. Material test data for the aortic valve was obtained from a dynamic inflation test. In each case, a numerical simulation of the experiment was performed and predictions were compared to the real data. For the uniaxial test simulation, the experimental targets were force at a measured displacement. For the inflation test, the experimental targets were the three-dimensional coordinates of material markers at a given pressure. For both sets of tissues, predictions with converged parameters showed excellent agreement with the data, and we found that the method was able to consistently identify model parameters. We believe the method will find wide application in biomedical material characterization and in diagnostic imaging.

Dor procedure for dyskinetic anteroapical myocardial infarction fails to improve contractility in the border zone

BACKGROUND Endoventricular patch plasty (Dor) is used to reduce left ventricular volume after myocardial infarction and subsequent left ventricular remodeling. METHODS AND RESULTS End-diastolic and end-systolic pressure-volume and Starling relationships were measured, and magnetic resonance images with noninvasive tags were used to calculate 3-dimensional myocardial strain in 6 sheep 2 weeks before and 2 and 6 weeks after the Dor procedure. These experimental results were previously reported. The imaging data from 1 sheep were incomplete. Animal specific finite element models were created from the remaining 5 animals using magnetic resonance images and left ventricular pressure obtained at early diastolic filling. Finite element models were optimized with 3-dimensional strain and used to determine systolic material properties, T(max,skinned-fiber), and diastolic and systolic stress in remote myocardium and border zone. Six weeks after the Dor procedure, end-diastolic and end-systolic stress in the border zone were substantially reduced. However, although there was a slight increase in T(max,skinned-fiber) in the border zone near the myocardial infarction at 6 weeks, the change was not significant. CONCLUSIONS The Dor procedure decreases end-diastolic and end-systolic stress but fails to improve contractility in the infarct border zone. Future work should focus on measures that will enhance border zone function alone or in combination with surgical remodeling.

High resolution lung airway cast segmentation with proper topology suitable for computational fluid dynamic simulations

Developing detailed lung airway models is an important step towards understanding the respiratory system. While modern imaging and airway casting approaches have dramatically improved the potential detail of such models, challenges have arisen in image processing as the demand for greater detail pushes the image processing approaches to their limits. Airway segmentations with proper topology have neither loops nor invalid voxel-to-voxel connections. Here we describe a new technique for segmenting airways with proper topology and apply the approach to an image volume generated by magnetic resonance imaging of a silicone cast created from an excised monkey lung.

Phase-Contrast MRI and CFD Modeling of Apparent 3He Gas Flow in Rat Pulmonary Airways

Phase-contrast (PC) magnetic resonance imaging (MRI) with hyperpolarized ³He is potentially useful for developing and testing patient-specific models of pulmonary airflow. One challenge, however, is that PC-MRI provides apparent values of local ³He velocity that not only depend on actual airflow but also on gas diffusion. This not only blurs laminar flow patterns in narrow airways but also introduces anomalous airflow structure that reflects gas-wall interactions. Here, both effects are predicted in a live rat using computational fluid dynamics (CFD), and for the first time, simulated patterns of apparent ³He gas velocity are compared with in vivo PC-MRI. Results show (1) that correlations (R²) between measured and simulated airflow patterns increase from 0.23 to 0.79 simply by accounting for apparent ³He transport, and (2) that remaining differences are mainly due to uncertain airway segmentation and partial volume effects stemming from relatively coarse MRI resolution. Higher-fidelity testing of pulmonary airflow predictions should therefore be possible with future imaging improvements.

Multiscale modeling in computational biomechanics

Biomechanics is broadly defined as the scientific discipline that investigates the effects of forces acting on and within biological structures. The realm of biomechanics includes the circulatory and respiratory systems, tissue mechanics and mechanotransduction, and the musculoskeletal system and motor control. As in many other biological phenomena, many spatial scales are crossed by biomechanics research: intracellular, multicellular, and extracellular matrices; and tissue, organ, and multiorgan systems. It is well established that the effect of forces at higher scales influence behavior at lower scales and that lower-scale properties influence higher-scale response. However, computational methods that incorporate these interactions in biomechanics are relatively rare. In general, computational models that include representation of multiple spatial or temporal scales are loosely defined as multiscale. The fact that multiscale modeling is not well defined lends the term to a variety of scenarios within the computational physiology community. In biomechanics, multiscale modeling may mean establishing a hierarchical link between the spatial and temporal scales, while the output of a larger-scale system is passed through a finely detailed representation at a lower scale (e.g., body-level movement simulations that provide net joint loading for tissue-level stress analysis). In reality, multiscale modeling may require more intricate representation of interactions among scales. A concurrent simulation strategy is inevitable to adequately represent nonlinear associations that have been known for decades [1].

Development of a rhesus monkey lung geometry model and application to particle deposition in comparison to humans

The exposure-dose-response characterization of an inhalation hazard established in an animal species needs to be translated to an equivalent characterization in humans relative to comparable doses or exposure scenarios. Here, the first geometry model of the conducting airways for rhesus monkeys is developed based upon CT images of the conducting airways of a 6-month-old male, rhesus monkey. An algorithm was developed for adding the alveolar region airways using published rhesus morphometric data. The resultant lung geometry model can be used in mechanistic particle or gaseous dosimetry models. Such dosimetry models require estimates of the upper respiratory tract volume of the animal and the functional residual capacity, as well as of the tidal volume and breathing frequency of the animal. The relationship of these variables to rhesus monkeys of differing body weights was established by synthesizing and modeling published data as well as modeling pulmonary function measurements on 121 rhesus control animals. Deposition patterns of particles up to 10 µm in size were examined for endotracheal and and up to 5 µm for spontaneous breathing in infant and young adult monkeys and compared to those for humans. Deposition fraction of respirable size particles was found to be higher in the conducting airways of infant and young adult rhesus monkeys compared to humans. Due to the filtering effect of the conducting airways, pulmonary deposition in rhesus monkeys was lower than that in humans. Future research areas are identified that would either allow replacing assumptions or improving the newly developed lung model.