Search

Mechanical forces are essential for proper growth and remodeling of the primitive pharyngeal arch arteries (PAAs) into the great vessels of the heart. Despite general acknowledgement of a hemodynamic-malformation link, the direct correlation between hemodynamics and PAA morphogenesis remains poorly understood. The elusiveness is largely due to difficulty in performing isolated hemodynamic perturbations and quantifying changes in-vivo. Previous in-vivo arch artery occlusion/ablation experiments either did not isolate the effects of hemodynamics, did not analyze the results in a 3D context or did not consider the effects of varying degrees of occlusion. Here, we overcome these limitations by combining minimally invasive occlusion experiments in the avian embryo with 3D anatomical models of development and in-silico testing of experimental phenomenon. We detail morphological and hemodynamic changes 24 hours post vessel occlusion. 3D anatomical models showed that occlusion geometries had more circular cross-sectional areas and more elongated arches than their control counterparts. Computational fluid dynamics revealed a marked change in wall shear stress-morphology trends. Instantaneous (in-silico) occlusion models provided mechanistic insights into the dynamic vessel adaptation process, predicting pressure-area trends for a number of experimental occlusion arches. We follow the propagation of small defects in a single embryo Hamburger Hamilton (HH) Stage 18 embryo to a more serious defect in an HH29 embryo. Results demonstrate that hemodynamic perturbation of the presumptive aortic arch, through varying degrees of vessel occlusion, overrides natural growth mechanisms and prevents it from becoming the dominant arch of the aorta., (© 2021. Biomedical Engineering Society.)

Image-based in silico modeling tools provide detailed velocity and particle deposition data. However, care must be taken when prescribing boundary conditions to model lung physiology in health or disease, such as in emphysema. In this study, the respiratory resistance and compliance were obtained by solving an inverse problem; a 0D global model based on healthy and emphysematous rat experimental data. Multi-scale CFD simulations were performed by solving the 3D Navier-Stokes equations in an MRI-derived rat geometry coupled to a 0D model. Particles with 0.95 μm diameter were tracked and their distribution in the lung was assessed. Seven 3D-0D simulations were performed: healthy, homogeneous, and five heterogeneous emphysema cases. Compliance (C) was significantly higher (p = 0.04) in the emphysematous rats (C = 0.37 ± 0.14 cm(3)/cmH2O) compared to the healthy rats (C = 0.25 ± 0.04 cm(3)/cmH2O), while the resistance remained unchanged (p = 0.83). There were increases in airflow, particle deposition in the 3D model, and particle delivery to the diseased regions for the heterogeneous cases compared to the homogeneous cases. The results highlight the importance of multi-scale numerical simulations to study airflow and particle distribution in healthy and diseased lungs. The effect of particle size and gravity were studied. Once available, these in silico predictions may be compared to experimental deposition data.

Coronary flow is different from the flow in other parts of the arterial system because it is influenced by the contraction and relaxation of the heart. To model coronary flow realistically, the compressive force of the heart acting on the coronary vessels needs to be included. In this study, we developed a method that predicts coronary flow and pressure of three-dimensional epicardial coronary arteries by considering models of the heart and arterial system and the interactions between the two models. For each coronary outlet, a lumped parameter coronary vascular bed model was assigned to represent the impedance of the downstream coronary vascular networks absent in the computational domain. The intramyocardial pressure was represented with either the left or right ventricular pressure depending on the location of the coronary arteries. The left and right ventricular pressure were solved from the lumped parameter heart models coupled to a closed loop system comprising a three-dimensional model of the aorta, three-element Windkessel models of the rest of the systemic circulation and the pulmonary circulation, and lumped parameter models for the left and right sides of the heart. The computed coronary flow and pressure and the aortic flow and pressure waveforms were realistic as compared to literature data.

Aortic flow and pressure result from the interactions between the heart and arterial system. In this work, we considered these interactions by utilizing a lumped parameter heart model as an inflow boundary condition for three-dimensional finite element simulations of aortic blood flow and vessel wall dynamics. The ventricular pressure-volume behavior of the lumped parameter heart model is approximated using a time varying elastance function scaled from a normalized elastance function. When the aortic valve is open, the coupled multidomain method is used to strongly couple the lumped parameter heart model and three-dimensional arterial models and compute ventricular volume, ventricular pressure, aortic flow, and aortic pressure. The shape of the velocity profiles of the inlet boundary and the outlet boundaries that experience retrograde flow are constrained to achieve a robust algorithm. When the aortic valve is closed, the inflow boundary condition is switched to a zero velocity Dirichlet condition. With this method, we obtain physiologically realistic aortic flow and pressure waveforms. We demonstrate this method in a patient-specific model of a normal human thoracic aorta under rest and exercise conditions and an aortic coarctation model under pre- and post-interventions.

Congenital heart defects with a single functional ventricle, such as hypoplastic left heart syndrome and tricuspid atresia, require a staged surgical approach to separate the systemic and pulmonary circulations. Ultimately, the venous or pulmonary side of the heart is bypassed by directly connecting the vena cava to the pulmonary arteries with a modified t-shaped junction. The Fontan procedure (total cavopulmonary connection, TCPC) completes this process of separation. To date, computational fluid dynamics (CFD) simulations in this low pressure, passive flow, intrathoracic system have neglected the presumed important effects of respiration on physiology and higher "stress" states such as with exercise have never been considered. We hypothesize that incorporating effects of respiration and exercise would provide more realistic estimates of TCPC performance. Time-dependent, 3D blood flow simulations are performed by a custom finite element solver for two patient-specific Fontan models with a novel respiration model, developed to generate physiologic time-varying flow conditions. Blood flow features, pressure, and energy efficiency are analyzed at rest and with increasing flow rates to simulate exercise conditions. The simulations produce realistic pressure and flow data, comparable to that measured by catheterization and echocardiography, and demonstrate substantial increases in energy dissipation (i.e. decreased performance) with exercise and respiration due to increasing intensity of small scale vortices in the flow. As would be expected, these changes are highly dependent on patient-specific anatomy and Fontan geometry. We propose that respiration and exercise should be incorporated into TCPC CFD simulations to provide increasingly realistic evaluations of TCPC performance.