In vitro validation of finite-element model of AAA hemodynamics incorporating realistic outlet boundary conditions

Hemodynamic Analysis in Aortic Dilatation after Arterial Switch Operation for Patients with Transposition of Great Arteries Using Computational Fluid Dynamics

After an arterial switch operation for complete transposition of the great arteries, neo-aortic root dilatation occurs, with unclear hemodynamic effects. This study analyzes three groups (severe dilation, mild dilation, and normal) using computational fluid dynamics (CFD) on cardiac CT scans. Aortic arch angles in severe (median 72.3, range: 68.5-77.2) and mild dilation (76.6, 71.1-85.2) groups are significantly smaller than the normal group (97.3, 87.4-99.0). In the normal and mild dilatation groups, Wall Shear Stress (WSS) exhibits a consistent pattern: it is lowest at the aortic root, gradually increases until just before the bend in the aortic arch, peaks, and then subsequently decreases. However, severe dilation shows disrupted WSS patterns, notably lower in the distal ascending aorta, attributed to local recirculation. This unique WSS pattern observed in severely dilated patients, especially in the transverse aorta. CFD plays an essential role in comprehensively studying the pathophysiology underlying aortic dilation in this population.

Integrated particle image velocimetry and fluid-structure interaction analysis for patient-specific abdominal aortic aneurysm studies

BACKGROUND

Understanding the hemodynamics of an abdominal aortic aneurysm (AAA) is crucial for risk assessment and treatment planning. This study introduces a low-cost, patient-specific in vitro AAA model to investigate hemodynamics using particle image velocimetry (PIV) and flow-simulating circuit, validated through fluid-structure interaction (FSI) simulations.

METHODS

In this study, 3D printing was employed to manufacture a flexible patient-specific AAA phantom using a lost-core casting technique. A pulsatile flow circuit was constructed using off-the-shelf components. A particle image velocimetry (PIV) setup was built using an affordable laser source and global shutter camera, and finally, the flow field inside the AAA was analyzed using open-source software. Fluid-structure interaction (FSI) simulations were performed to enhance our understanding of the flow field, and the results were validated by PIV analysis. Both steady-state and transient flow conditions were investigated.

RESULTS

Our experimental setup replicated physiological conditions, analyzing arterial wall deformations and flow characteristics within the aneurysm. Under constant flow, peak wall deformations and flow velocities showed deviations within - 12% to + 27% and - 7% to + 5%, respectively, compared to FSI simulations. Pulsatile flow conditions further demonstrated a strong correlation (Pearson coefficient 0.85) in flow velocities and vectors throughout the cardiac cycle. Transient phenomena, particularly the formation and progression of vortex structures during systole, were consistently depicted between experimental and numerical models.

CONCLUSIONS

By bridging high-fidelity experimental observations with comprehensive computational analyses, this study underscores the potential of integrated methodologies in enhancing our understanding of AAA pathophysiology. The convergence of realistic AAA phantoms, precise PIV measurements at affordable cost point, and validated FSI models heralds a new paradigm in vascular research, with significant implications for personalized medicine and bioengineering innovations.

Deep-learning-based image segmentation for image-based computational hemodynamic analysis of abdominal aortic aneurysms: a comparison study

Computational hemodynamics is increasingly being used to quantify hemodynamic characteristics in and around abdominal aortic aneurysms (AAA) in a patient-specific fashion. However, the time-consuming manual annotation hinders the clinical translation of computational hemodynamic analysis. Thus, we investigate the feasibility of using deep-learning-based image segmentation methods to reduce the time required for manual segmentation. Two of the latest deep-learning-based image segmentation methods, ARU-Net and CACU-Net, were used to test the feasibility of automated computer model creation for computational hemodynamic analysis. Morphological features and hemodynamic metrics of 30 computed tomography angiography (CTA) scans were compared between pre-dictions and manual models. The DICE score for both networks was 0.916, and the correlation value was above 0.95, indicating their ability to generate models comparable to human segmentation. The Bland-Altman analysis shows a good agreement between deep learning and manual segmentation results. Compared with manual (computational hemodynamics) model recreation, the time for automated computer model generation was significantly reduced (from ∼2 h to ∼10 min). Automated image segmentation can significantly reduce time expenses on the recreation of patient-specific AAA models. Moreover, our study showed that both CACU-Net and ARU-Net could accomplish AAA segmentation, and CACU-Net outperformed ARU-Net in terms of accuracy and time-saving.

A Protocol for Coupling Volumetrically Dynamic In-Vitro Experiments to Numerical Physiology Simulation for a Hybrid Cardiovascular Model

OBJECTIVE

The Physiology Simulation Coupled Experiment (PSCOPE) is a hybrid modeling framework that enables a physical fluid experiment to operate in the context of a closed-loop computational simulation of cardiovascular physiology. Previous PSCOPE methods coupled rigid experiments to a lumped parameter network (LPN) of physiology but are incompatible with volumetrically dynamic experiments where fluid volume varies periodically. We address this limitation by introducing a method capable of coupling multi-branch and volumetrically dynamic in-vitro experiments to an LPN.

METHODS

Our proposed method utilizes an iterative weighted-averaging algorithm to identify the unique solution waveforms for a given PSCOPE model. We confirm the accuracy of these PSCOPE solutions by integrating mathematical surrogates of in-vitro experiments directly into the LPN to derive reference solutions, which serve as the gold standard to validate the solutions obtained from using our proposed method to couple the same mathematical surrogates to the LPN. Finally, we illustrate a practical application of our PSCOPE method by coupling an in-vitro renal circulation experiment to the LPN.

RESULTS

Compared to the reference solution, the normalized root mean square error of the flow and pressure waveforms were 0.001%∼0.55%, demonstrating the accuracy of the coupling method.

CONCLUSION

We successfully coupled the in-vitro experiment to the LPN, demonstrating the real-world performance within the constraints of sensor and actuation limitations in the physical experiment.

SIGNIFICANCE

This study introduces a PSCOPE method that can be used to investigate medical devices and anatomies that exhibit periodic volume changes, expanding the utility of the hybrid framework.

Validation of the Reduced Unified Continuum Formulation Against In Vitro 4D-Flow MRI

We previously introduced and verified the reduced unified continuum formulation for vascular fluid-structure interaction (FSI) against Womersley's deformable wall theory. Our present work seeks to investigate its performance in a patient-specific aortic setting in which assumptions of idealized geometries and velocity profiles are invalid. Specifically, we leveraged 2D magnetic resonance imaging (MRI) and 4D-flow MRI to extract high-resolution anatomical and hemodynamic information from an in vitro flow circuit embedding a compliant 3D-printed aortic phantom. To accurately reflect experimental conditions, we numerically implemented viscoelastic external tissue support, vascular tissue prestressing, and skew boundary conditions enabling in-plane vascular motion at each inlet and outlet. Validation of our formulation is achieved through close quantitative agreement in pressures, lumen area changes, pulse wave velocity, and early systolic velocities, as well as qualitative agreement in late systolic flow structures. Our validated suite of FSI techniques offers a computationally efficient approach for numerical simulation of vascular hemodynamics. This study is among the first to validate a cardiovascular FSI formulation against an in vitro flow circuit involving a compliant vascular phantom of complex patient-specific anatomy.

Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases

Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.

A modified method of computed fluid dynamics simulation in abdominal aorta and visceral arteries

PURPOSE

The flow velocity of visceral arteries was measured by 2D PCMRI to produce the patient-specific flow BC imposed on the outlets of visceral arteries in CFD simulation. This modified method aimed to improve the CFD accuracy in the abdominal aorta and visceral arteries.

METHODS

A volunteer underwent non-contrast-enhanced MRA to scan the abdominal aorta and visceral arteries, and 2D PCMRI to obtain the flow velocity of the aforementioned vessels. The three-dimensional geometric model was reconstructed using the MRI scan data of the abdominal aorta and visceral arteries. The flow waveforms measured by 2D PCMRI were processed and then imposed on the aortic inlet and the outlets of all visceral arteries as the flow BC. The RCR parameters of the three elements Windkessel model were modulated and imposed on the aortic outlet. CFD simulation was run in the open-source software: svSolver. The same volunteer underwent 4D flow MRI to compare the flow field with those extracted from CFD results.

RESULTS

Four specific time points in a cardiac cycle and three cross-sectional planes of aorta were selected to analyze the flow field, pressure and wall shear stress (WSS) from CFD. The flow waveforms and streamlines of CFD agreed with those of 4D flow MRI. The pressure waveforms, pressure distribution and WSS distribution from CFD conformed with the physiological condition of human body.

CONCLUSION

These results suggest this modified CFD method may yield reasonable flow field, pressure and WSS in the abdominal aorta and visceral arteries.

Design, Development, and Temporal Evaluation of a Magnetic Resonance Imaging-Compatible In Vitro Circulation Model Using a Compliant Abdominal Aortic Aneurysm Phantom

Biomechanical characterization of abdominal aortic aneurysms (AAAs) has become commonplace in rupture risk assessment studies. However, its translation to the clinic has been greatly limited due to the complexity associated with its tools and their implementation. The unattainability of patient-specific tissue properties leads to the use of generalized population-averaged material models in finite element analyses, which adds a degree of uncertainty to the wall mechanics quantification. In addition, computational fluid dynamics modeling of AAA typically lacks the patient-specific inflow and outflow boundary conditions that should be obtained by nonstandard of care clinical imaging. An alternative approach for analyzing AAA flow and sac volume changes is to conduct in vitro experiments in a controlled laboratory environment. In this study, we designed, built, and characterized quantitatively a benchtop flow loop using a deformable AAA silicone phantom representative of a patient-specific geometry. The impedance modules, which are essential components of the flow loop, were fine-tuned to ensure typical intraluminal pressure conditions within the AAA sac. The phantom was imaged with a magnetic resonance imaging (MRI) scanner to acquire time-resolved images of the moving wall and the velocity field inside the sac. Temporal AAA sac volume changes lead to a corresponding variation in compliance throughout the cardiac cycle. The primary outcome of this work was the design optimization of the impedance elements, the quantitative characterization of the resistive and capacitive attributes of a compliant AAA phantom, and the exemplary use of MRI for flow visualization and quantification of the deformed AAA geometry.

Multiscale Computational Fluid Dynamics Modeling for Personalized Liver Cancer Radioembolization Dosimetry

Yttrium-90 (90Y) radioembolization is a minimally invasive procedure increasingly used for advanced liver cancer treatment. In this method, radioactive microspheres are injected into the hepatic arterial bloodstream to target, irradiate, and kill cancer cells. Accurate and precise treatment planning can lead to more efficient and safer treatment by delivering a higher radiation dose to the tumor while minimizing the exposure of the surrounding liver parenchyma. Treatment planning primarily relies on the estimated radiation dose delivered to tissue. However, current methods used to estimate the dose are based on simplified assumptions that make the dosimetry results unreliable. In this work, we present a computational model to predict the radiation dose from the 90Y activity in different liver segments to provide a more realistic and personalized dosimetry. Computational fluid dynamics (CFD) simulations were performed in a 3D hepatic arterial tree model segmented from cone-beam CT angiographic data obtained from a patient with hepatocellular carcinoma (HCC). The microsphere trajectories were predicted from the velocity field. 90Y dose distribution was then calculated from the volumetric distribution of the microspheres. Two injection locations were considered for the microsphere administration, a lobar and a selective injection. Results showed that 22% and 82% of the microspheres were delivered to the tumor, after each injection, respectively, and the combination of both injections ultimately delivered 49% of the total administered 90Y microspheres to the tumor. Results also illustrated the nonhomogeneous distribution of microspheres between liver segments, indicating the importance of developing patient-specific dosimetry methods for effective radioembolization treatment.

Roudsari B, Vu CT, Roncali E. Computational Modeling of the Liver Arterial Blood Flow for Microsphere Therapy: Effect of Boundary Conditions

Transarterial embolization is a minimally invasive treatment for advanced liver cancer using microspheres loaded with a chemotherapeutic drug or radioactive yttrium-90 (90Y) that are injected into the hepatic arterial tree through a catheter. For personalized treatment, the microsphere distribution in the liver should be optimized through the injection volume and location. Computational fluid dynamics (CFD) simulations of the blood flow in the hepatic artery can help estimate this distribution if carefully parameterized. An important aspect is the choice of the boundary conditions imposed at the inlet and outlets of the computational domain. In this study, the effect of boundary conditions on the hepatic arterial tree hemodynamics was investigated. The outlet boundary conditions were modeled with three-element Windkessel circuits, representative of the downstream vasculature resistance. Results demonstrated that the downstream vasculature resistance affected the hepatic artery hemodynamics such as the velocity field, the pressure field and the blood flow streamline trajectories. Moreover, the number of microspheres received by the tumor significantly changed (more than 10% of the total injected microspheres) with downstream resistance variations. These findings suggest that patient-specific boundary conditions should be used in order to achieve a more accurate drug distribution estimation with CFD in transarterial embolization treatment planning.

Computational Hemodynamic Modeling of Arterial Aneurysms: A Mini-Review

Arterial aneurysms are pathological dilations of blood vessels, which can be of clinical concern due to thrombosis, dissection, or rupture. Aneurysms can form throughout the arterial system, including intracranial, thoracic, abdominal, visceral, peripheral, or coronary arteries. Currently, aneurysm diameter and expansion rates are the most commonly used metrics to assess rupture risk. Surgical or endovascular interventions are clinical treatment options, but are invasive and associated with risk for the patient. For aneurysms in locations where thrombosis is the primary concern, diameter is also used to determine the level of therapeutic anticoagulation, a treatment that increases the possibility of internal bleeding. Since simple diameter is often insufficient to reliably determine rupture and thrombosis risk, computational hemodynamic simulations are being developed to help assess when an intervention is warranted. Created from subject-specific data, computational models have the potential to be used to predict growth, dissection, rupture, and thrombus-formation risk based on hemodynamic parameters, including wall shear stress, oscillatory shear index, residence time, and anomalous blood flow patterns. Generally, endothelial damage and flow stagnation within aneurysms can lead to coagulation, inflammation, and the release of proteases, which alter extracellular matrix composition, increasing risk of rupture. In this review, we highlight recent work that investigates aneurysm geometry, model parameter assumptions, and other specific considerations that influence computational aneurysm simulations. By highlighting modeling validation and verification approaches, we hope to inspire future computational efforts aimed at improving our understanding of aneurysm pathology and treatment risk stratification.

An algorithm for coupling multibranch in vitro experiment to numerical physiology simulation for a hybrid cardiovascular model

The hybrid cardiovascular modeling approach integrates an in vitro experiment with a computational lumped-parameter simulation, enabling direct physical testing of medical devices in the context of closed-loop physiology. The interface between the in vitro and computational domains is essential for properly capturing the dynamic interactions of the two. To this end, we developed an iterative algorithm capable of coupling an in vitro experiment containing multiple branches to a lumped-parameter physiology simulation. This algorithm identifies the unique flow waveform solution for each branch of the experiment using an iterative Broyden's approach. For the purpose of algorithm testing, we first used mathematical surrogates to represent the in vitro experiments and demonstrated five scenarios where the in vitro surrogates are coupled to the computational physiology of a Fontan patient. This testing approach allows validation of the coupling result accuracy as the mathematical surrogates can be directly integrated into the computational simulation to obtain the "true solution" of the coupled system. Our algorithm successfully identified the solution flow waveforms in all test scenarios with results matching the true solutions with high accuracy. In all test cases, the number of iterations to achieve the desired convergence criteria was less than 130. To emulate realistic in vitro experiments in which noise contaminates the measurements, we perturbed the surrogate models by adding random noise. The convergence tolerance achievable with the coupling algorithm remained below the magnitudes of the added noise in all cases. Finally, we used this algorithm to couple a physical experiment to the computational physiology model to demonstrate its real-world applicability.

Anthropomorphic and biomechanical mockup for abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) is an asymptomatic condition due to the dilation of abdominal aorta along with progressive wall degeneration, where rupture of AAA is life-threatening. Failures of AAA endovascular repair (EVAR) reflect our inadequate knowledge about the complex interaction between the aortic wall and medical devices. In this regard, we are presenting a hydrogel-based anthropomorphic mockup (AMM) to better understand the biomechanical constraints during EVAR. By adjusting the cryogenic treatments, we tailored the hydrogel to mimic the mechanical behavior of human AAA wall, thrombus and abdominal fat. A specific molding sequence and a pressurizing system were designed to reproduce the geometrical and diseased characteristics of AAA. A mechanically, anatomically and pathologically realistic AMM for AAA was developed for the first time, EVAR experiments were then performed with and without the surrounding fat. Substantial displacements of the aortic centerlines and vessel expansion were observed in the case without surrounding fat, revealing an essential framework created by the surrounding fat to account for the interactions with medical devices. In conclusion, the importance to consider surrounding tissue for the global deformation of AAA during EVAR was highlighted. Furthermore, potential use of this AMM for medical training was also suggested.

The Implementation of an Adjustable Afterload Module for Ex Situ Heart Perfusion

PURPOSE

Windkessel impedance analysis has proven to be an effective technique for instituting artificial afterload on ex situ hearts. Traditional fixed parameter afterload modules, however, are unable to handle the changing contractile conditions associated with prolonged ex situ heart perfusion. In this paper, an adjustable afterload module is described comprising of three fully adjustable sub-components: a systemic resistor, a proximal resistor and a compliance chamber.

METHODS

Using a centrifugal pump, the systemic resistor and compliance chamber were subjected to testing across their operating ranges, whereby the predictability of resistance and compliance values was evaluated. The components were then assembled, and the full module tested on three separate porcine hearts perfused for 6 h with success defined by the ability to maintain physiological systolic and diastolic aortic pressures across flow rate variability.

RESULTS

For both the systemic resistor and compliance chamber, experimental measurements agreed with their theoretical equivalents, with coefficients of determination of 0.99 and 0.97 for the systemic resistor and compliance chamber, respectively. During ex situ perfusion, overall 95% confidence intervals demonstrate that physiological systolic (95-96.21 mmHg) and diastolic (26.8-28.8 mmHg) pressures were successfully maintained, despite large variability in aortic flow. Left ventricular contractile parameters, were found to be in line with those in previous studies, suggesting the afterload module has no detrimental impact on functional preservation.

CONCLUSIONS

We conclude that due to the demonstrable control of our afterload module, we can maintain physiological aortic pressures in a passive afterload working mode across prolonged perfusion periods, enabling effective perfusion regardless of contractile performance.

Biomechanical Investigation of Disturbed Hemodynamics-Induced Tissue Degeneration in Abdominal Aortic Aneurysms Using Computational and Experimental Techniques

Abdominal aortic aneurysm (AAA) is the dilatation of the aorta beyond 50% of the normal vessel diameter. It is reported that 4-8% of men and 0.5-1% of women above 50 years of age bear an AAA and it accounts for ~15,000 deaths per year in the United States alone. If left untreated, AAA might gradually expand until rupture; the most catastrophic complication of the aneurysmal disease that is accompanied by a striking overall mortality of 80%. The precise mechanisms leading to AAA rupture remains unclear. Therefore, characterization of disturbed hemodynamics within AAAs will help to understand the mechanobiological development of the condition which will contribute to novel therapies for the condition. Due to geometrical complexities, it is challenging to directly quantify disturbed flows for AAAs clinically. Two other approaches for this investigation are computational modeling and experimental flow measurement. In computational modeling, the problem is first defined mathematically, and the solution is approximated with numerical techniques to get characteristics of flow. In experimental flow measurement, once the setup providing physiological flow pattern in a phantom geometry is constructed, velocity measurement system such as particle image velocimetry (PIV) enables characterization of the flow. We witness increasing number of applications of these complimentary approaches for AAA investigations in recent years. In this paper, we outline the details of computational modeling procedures and experimental settings and summarize important findings from recent studies, which will help researchers for AAA investigations and rupture mechanics.

4D Flow MRI Estimation of Boundary Conditions for Patient Specific Cardiovascular Simulation

Accurate image based cardiovascular simulations require patient specific boundary conditions (BCs) for inlets, outlets and vessel wall mechanical properties. While inlet BCs are typically determined non-invasively, invasive pressure catheterization is often used to determine patient specific outlet BCs and vessel wall mechanical properties. A method using 4D Flow MRI to non-invasively determine both patient specific outlet BCs and vessel wall mechanical properties is presented and results for both in vitro validation with a latex tube and an in vivo pulmonary artery stenosis (PAS) stent intervention are presented. For in vitro validation, acceptable agreement is found between simulation using BCs from 4D Flow MRI and benchtop measurements. For the PAS virtual intervention, simulation correctly predicts flow distribution with 9% error compared to MRI. Using 4D Flow MRI to noninvasively determine patient specific BCs increases the ability to use image based simulations as pressure catheterization is not always performed.

A Hybrid Experimental-Computational Modeling Framework for Cardiovascular Device Testing

Significant advances in biomedical science often leverage powerful computational and experimental modeling platforms. We present a framework named physiology simulation coupled experiment ("PSCOPE") that can capitalize on the strengths of both types of platforms in a single hybrid model. PSCOPE uses an iterative method to couple an in vitro mock circuit to a lumped-parameter numerical simulation of physiology, obtaining closed-loop feedback between the two. We first compared the results of Fontan graft obstruction scenarios modeled using both PSCOPE and an established multiscale computational fluid dynamics method; the normalized root-mean-square error values of important physiologic parameters were between 0.1% and 2.1%, confirming the fidelity of the PSCOPE framework. Next, we demonstrate an example application of PSCOPE to model a scenario beyond the current capabilities of multiscale computational methods-the implantation of a Jarvik 2000 blood pump for cavopulmonary support in the single-ventricle circulation; we found that the commercial Jarvik 2000 controller can be modified to produce a suitable rotor speed for augmenting cardiac output by approximately 20% while maintaining blood pressures within safe ranges. The unified modeling framework enables a testing environment which simultaneously operates a medical device and performs computational simulations of the resulting physiology, providing a tool for physically testing medical devices with simulated physiologic feedback.

[A correction method for calculating resistance of Westerhof 's resistor]

The objective of the mock circulatory system (MCS) is to construct the characteristics of cardiovascular hemodynamics. Westerhof 's resistor that often regarded as the laminar flow resistance in the MCS, is commonly used to simulate the peripheral resistance of the cardiovascular system. However, the theoretical calculation value of fluid resistance of the Westerhof 's resistor shows distinguished difference with the actual needed value. If the theoretical resistance is regarded as the actual needed one and be used directly in the experiment, the experimental accuracy would not be acceptable. In order to improve the accuracy, an effective correction method for calculating the resistance of Westerhof 's resistor was proposed in this paper. Simulation software was also developed to compute accurately the capillary number, total length and resistance. The results demonstrate the proposed method is able to reduce the difficulty and complexity of the design of the resistor, which would obviously increase the manufactured precision of the Westerhof 's resistor. Simulation software would provide great support to the construction of various MCSs.

Use of the FDA nozzle model to illustrate validation techniques in computational fluid dynamics (CFD) simulations

A "credible" computational fluid dynamics (CFD) model has the potential to provide a meaningful evaluation of safety in medical devices. One major challenge in establishing "model credibility" is to determine the required degree of similarity between the model and experimental results for the model to be considered sufficiently validated. This study proposes a "threshold-based" validation approach that provides a well-defined acceptance criteria, which is a function of how close the simulation and experimental results are to the safety threshold, for establishing the model validity. The validation criteria developed following the threshold approach is not only a function of Comparison Error, E (which is the difference between experiments and simulations) but also takes in to account the risk to patient safety because of E. The method is applicable for scenarios in which a safety threshold can be clearly defined (e.g., the viscous shear-stress threshold for hemolysis in blood contacting devices). The applicability of the new validation approach was tested on the FDA nozzle geometry. The context of use (COU) was to evaluate if the instantaneous viscous shear stress in the nozzle geometry at Reynolds numbers (Re) of 3500 and 6500 was below the commonly accepted threshold for hemolysis. The CFD results ("S") of velocity and viscous shear stress were compared with inter-laboratory experimental measurements ("D"). The uncertainties in the CFD and experimental results due to input parameter uncertainties were quantified following the ASME V&V 20 standard. The CFD models for both Re = 3500 and 6500 could not be sufficiently validated by performing a direct comparison between CFD and experimental results using the Student's t-test. However, following the threshold-based approach, a Student's t-test comparing |S-D| and |Threshold-S| showed that relative to the threshold, the CFD and experimental datasets for Re = 3500 were statistically similar and the model could be considered sufficiently validated for the COU. However, for Re = 6500, at certain locations where the shear stress is close the hemolysis threshold, the CFD model could not be considered sufficiently validated for the COU. Our analysis showed that the model could be sufficiently validated either by reducing the uncertainties in experiments, simulations, and the threshold or by increasing the sample size for the experiments and simulations. The threshold approach can be applied to all types of computational models and provides an objective way of determining model credibility and for evaluating medical devices.

A Real-Time Programmable Pulsatile Flow Pump for In Vitro Cardiovascular Experimentation

Benchtop in vitro experiments are valuable tools for investigating the cardiovascular system and testing medical devices. Accurate reproduction of the physiologic flow waveforms at various anatomic locations is an important component of these experimental methods. This study discusses the design, construction, and testing of a low-cost and fully programmable pulsatile flow pump capable of continuously producing unlimited cycles of physiologic waveforms. It consists of a gear pump actuated by an AC servomotor and a feedback algorithm to achieve highly accurate reproduction of flow waveforms for flow rates up to 300 ml/s across a range of loading conditions. The iterative feedback algorithm uses the flow error values in one iteration to modify the motor control waveform for the next iteration to better match the desired flow. Within four to seven iterations of feedback, the pump replicated desired physiologic flow waveforms to within 2% normalized RMS error (for flow rates above 20 mL/s) under varying downstream impedances. This pump device is significantly more affordable (∼10% of the cost) than current commercial options. More importantly, the pump can be controlled via common scientific software and thus easily implemented into large automation frameworks.

A High Performance Pulsatile Pump for Aortic Flow Experiments in 3-Dimensional Models

Aortic pathologies such as coarctation, dissection, and aneurysm represent a particularly emergent class of cardiovascular diseases. Computational simulations of aortic flows are growing increasingly important as tools for gaining understanding of these pathologies, as well as for planning their surgical repair. In vitro experiments are required to validate the simulations against real world data, and the experiments require a pulsatile flow pump system that can provide physiologic flow conditions characteristic of the aorta. We designed a newly capable piston-based pulsatile flow pump system that can generate high volume flow rates (850 mL/s), replicate physiologic waveforms, and pump high viscosity fluids against large impedances. The system is also compatible with a broad range of fluid types, and is operable in magnetic resonance imaging environments. Performance of the system was validated using image processing-based analysis of piston motion as well as particle image velocimetry. The new system represents a more capable pumping solution for aortic flow experiments than other available designs, and can be manufactured at a relatively low cost.

Computational modeling and engineering in pediatric and congenital heart disease

PURPOSE OF REVIEW

Recent methodological advances in computational simulations are enabling increasingly realistic simulations of hemodynamics and physiology, driving increased clinical utility. We review recent developments in the use of computational simulations in pediatric and congenital heart disease, describe the clinical impact in modeling in single-ventricle patients, and provide an overview of emerging areas.

RECENT FINDINGS

Multiscale modeling combining patient-specific hemodynamics with reduced order (i.e., mathematically and computationally simplified) circulatory models has become the de-facto standard for modeling local hemodynamics and 'global' circulatory physiology. We review recent advances that have enabled faster solutions, discuss new methods (e.g., fluid structure interaction and uncertainty quantification), which lend realism both computationally and clinically to results, highlight novel computationally derived surgical methods for single-ventricle patients, and discuss areas in which modeling has begun to exert its influence including Kawasaki disease, fetal circulation, tetralogy of Fallot (and pulmonary tree), and circulatory support.

SUMMARY

Computational modeling is emerging as a crucial tool for clinical decision-making and evaluation of novel surgical methods and interventions in pediatric cardiology and beyond. Continued development of modeling methods, with an eye towards clinical needs, will enable clinical adoption in a wide range of pediatric and congenital heart diseases.

Examination of near-wall hemodynamic parameters in the renal bridging stent of various stent graft configurations for repairing visceral branched aortic aneurysms

OBJECTIVE

This study examined the flow behavior of four stent graft configurations for endovascular repair of complex aneurysms of the descending aorta.

METHODS

Computational fluid dynamics models with transient boundary conditions and rigid wall simplifying assumptions were developed and used with four distinct geometries to compare various near-wall hemodynamic parameters.

RESULTS

Graphic plots for time-averaged wall shear stress, oscillating shear index, and relative residence time were presented and compared among the four stent graft configurations of interest.

CONCLUSIONS

Abrupt 90° and 180° changes in stent geometry (particularly in the side branches) cause a high momentum change and thus increased flow separation and mixing, which has significant implications in blood flow characteristics near the wall. By comparison, longer bridging stents provide more gradual changes in momentum, thus allowing blood flow to develop before reaching the target vessel.

Effects of intraluminal thrombus on patient-specific abdominal aortic aneurysm hemodynamics via stereoscopic particle image velocity and computational fluid dynamics modeling

The pathology of the human abdominal aortic aneurysm (AAA) and its relationship to the later complication of intraluminal thrombus (ILT) formation remains unclear. The hemodynamics in the diseased abdominal aorta are hypothesized to be a key contributor to the formation and growth of ILT. The objective of this investigation is to establish a reliable 3D flow visualization method with corresponding validation tests with high confidence in order to provide insight into the basic hemodynamic features for a better understanding of hemodynamics in AAA pathology and seek potential treatment for AAA diseases. A stereoscopic particle image velocity (PIV) experiment was conducted using transparent patient-specific experimental AAA models (with and without ILT) at three axial planes. Results show that before ILT formation, a 3D vortex was generated in the AAA phantom. This geometry-related vortex was not observed after the formation of ILT, indicating its possible role in the subsequent appearance of ILT in this patient. It may indicate that a longer residence time of recirculated blood flow in the aortic lumen due to this vortex caused sufficient shear-induced platelet activation to develop ILT and maintain uniform flow conditions. Additionally, two computational fluid dynamics (CFD) modeling codes (Fluent and an in-house cardiovascular CFD code) were compared with the two-dimensional, three-component velocity stereoscopic PIV data. Results showed that correlation coefficients of the out-of-plane velocity data between PIV and both CFD methods are greater than 0.85, demonstrating good quantitative agreement. The stereoscopic PIV study can be utilized as test case templates for ongoing efforts in cardiovascular CFD solver development. Likewise, it is envisaged that the patient-specific data may provide a benchmark for further studying hemodynamics of actual AAA, ILT, and their convolution effects under physiological conditions for clinical applications.

Computational Models of Aortic Coarctation in Hypoplastic Left Heart Syndrome: Considerations on Validation of a Detailed 3D model

Background

Reliability of computational models for cardiovascular investigations strongly depends on their validation against physical data. This study aims to experimentally validate a computational model of complex congenital heart disease (i.e., surgically palliated hypoplastic left heart syndrome with aortic coarctation) thus demonstrating that hemodynamic information can be reliably extrapolated from the model for clinically meaningful investigations.

Materials and methods

A patient-specific aortic arch model was tested in a mock circulatory system and the same flow conditions were re-created in silico, by setting an appropriate lumped parameter network (LPN) attached to the same three-dimensional (3D) aortic model (i.e., multi-scale approach). The model included a modified Blalock-Taussig shunt and coarctation of the aorta. Different flow regimes were tested as well as the impact of uncertainty in viscosity.

Results

Computational flow and pressure results were in good agreement with the experimental signals, both qualitatively, in terms of the shape of the waveforms, and quantitatively (mean aortic pressure 62.3 vs. 65.1 mmHg, 4.8% difference; mean aortic flow 28.0 vs. 28.4% inlet flow, 1.4% difference; coarctation pressure drop 30.0 vs. 33.5 mmHg, 10.4% difference), proving the reliability of the numerical approach. It was observed that substantial changes in fluid viscosity or using a turbulent model in the numerical simulations did not significantly affect flows and pressures of the investigated physiology. Results highlighted how the non-linear fluid dynamic phenomena occurring in vitro must be properly described to ensure satisfactory agreement.

Conclusions

This study presents methodological considerations for using experimental data to preliminarily set up a computational model, and then simulate a complex congenital physiology using a multi-scale approach.

In Vitro Validation of Patient-Specific Hemodynamic Simulations in Coronary Aneurysms Caused by Kawasaki Disease

To perform experimental validation of computational fluid dynamics (CFD) applied to patient specific coronary aneurysm anatomy of Kawasaki disease. We quantified hemodynamics in a patient-specific coronary artery aneurysm physical phantom under physiologic rest and exercise flow conditions. Using phase contrast MRI (PCMRI), we acquired 3-component flow velocity at two slice locations in the aneurysms. We then performed numerical simulations with the same geometry and inflow conditions, and performed qualitative and quantitative comparisons of velocities between experimental measurements and simulation results. We observed excellent qualitative agreement in flow pattern features. The quantitative spatially and temporally varying differences in velocity between PCMRI and CFD were proportional to the flow velocity. As a result, the percent discrepancy between simulation and experiment was relatively constant regardless of flow velocity variations. Through 1D and 2D quantitative comparisons, we found a 5-17% difference between measured and simulated velocities. Additional analysis assessed wall shear stress differences between deformable and rigid wall simulations. This study demonstrated that CFD produced good qualitative and quantitative predictions of velocities in a realistic coronary aneurysm anatomy under physiological flow conditions. The results provide insights on factors that may influence the level of agreement, and a set of in vitro experimental data that can be used by others to compare against CFD simulation results. The findings of this study increase confidence in the use of CFD for investigating hemodynamics in the specialized anatomy of coronary aneurysms. This provides a basis for future hemodynamics studies in patient-specific models of Kawasaki disease.

Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains

Modeling blood flow in larger vessels using lattice-Boltzmann methods comes with a challenging set of constraints: a complex geometry with walls and inlets and outlets at arbitrary orientations with respect to the lattice, intermediate Reynolds (Re) number, and unsteady flow. Simple bounce-back is one of the most commonly used, simplest, and most computationally efficient boundary conditions, but many others have been proposed. We implement three other methods applicable to complex geometries [Guo, Zheng, and Shi, Phys. Fluids 14, 2007 (2002); Bouzidi, Firdaouss, and Lallemand, Phys. Fluids 13, 3452 (2001); Junk and Yang, Phys. Rev. E 72, 066701 (2005)] in our open-source application hemelb. We use these to simulate Poiseuille and Womersley flows in a cylindrical pipe with an arbitrary orientation at physiologically relevant Re number (1-300) and Womersley (4-12) numbers and steady flow in a curved pipe at relevant Dean number (100-200) and compare the accuracy to analytical solutions. We find that both the Bouzidi-Firdaouss-Lallemand (BFL) and Guo-Zheng-Shi (GZS) methods give second-order convergence in space while simple bounce-back degrades to first order. The BFL method appears to perform better than GZS in unsteady flows and is significantly less computationally expensive. The Junk-Yang method shows poor stability at larger Re number and so cannot be recommended here. The choice of collision operator (lattice Bhatnagar-Gross-Krook vs multiple relaxation time) and velocity set (D3Q15 vs D3Q19 vs D3Q27) does not significantly affect the accuracy in the problems studied.

Modeling single ventricle physiology: review of engineering tools to study first stage palliation of hypoplastic left heart syndrome

First stage palliation of hypoplastic left heart syndrome, i.e., the Norwood operation, results in a complex physiological arrangement, involving different shunting options (modified Blalock-Taussig, RV-PA conduit, central shunt from the ascending aorta) and enlargement of the hypoplastic ascending aorta. Engineering techniques, both computational and experimental, can aid in the understanding of the Norwood physiology and their correct implementation can potentially lead to refinement of the decision-making process, by means of patient-specific simulations. This paper presents some of the available tools that can corroborate clinical evidence by providing detailed insight into the fluid dynamics of the Norwood circulation as well as alternative surgical scenarios (i.e., virtual surgery). Patient-specific anatomies can be manufactured by means of rapid prototyping and such models can be inserted in experimental set-ups (mock circulatory loops) that can provide a valuable source of validation data as well as hydrodynamic information. Such models can be tuned to respond to differing the patient physiologies. Experimental set-ups can also be compatible with visualization techniques, like particle image velocimetry and cardiovascular magnetic resonance, further adding to the knowledge of the local fluid dynamics. Multi-scale computational models include detailed three-dimensional (3D) anatomical information coupled to a lumped parameter network representing the remainder of the circulation. These models output both overall hemodynamic parameters while also enabling to investigate the local fluid dynamics of the aortic arch or the shunt. As an alternative, pure lumped parameter models can also be employed to model Stage 1 palliation, taking advantage of a much lower computational cost, albeit missing the 3D anatomical component. Finally, analytical techniques, such as wave intensity analysis, can be employed to study the Norwood physiology, providing a mechanistic perspective on the ventriculo-arterial coupling for this specific surgical scenario.

Simulation based planning of surgical interventions in pediatric cardiology

Hemodynamics plays an essential role in the progression and treatment of cardiovascular disease. However, while medical imaging provides increasingly detailed anatomical information, clinicians often have limited access to hemodynamic data that may be crucial to patient risk assessment and treatment planning. Computational simulations can now provide detailed hemodynamic data to augment clinical knowledge in both adult and pediatric applications. There is a particular need for simulation tools in pediatric cardiology, due to the wide variation in anatomy and physiology in congenital heart disease patients, necessitating individualized treatment plans. Despite great strides in medical imaging, enabling extraction of flow information from magnetic resonance and ultrasound imaging, simulations offer predictive capabilities that imaging alone cannot provide. Patient specific simulations can be used for in silico testing of new surgical designs, treatment planning, device testing, and patient risk stratification. Furthermore, simulations can be performed at no direct risk to the patient. In this paper, we outline the current state of the art in methods for cardiovascular blood flow simulation and virtual surgery. We then step through pressing challenges in the field, including multiscale modeling, boundary condition selection, optimization, and uncertainty quantification. Finally, we summarize simulation results of two representative examples from pediatric cardiology: single ventricle physiology, and coronary aneurysms caused by Kawasaki disease. These examples illustrate the potential impact of computational modeling tools in the clinical setting.

Size-dependent predilections of cardiogenic embolic transport

While it is intuitively clear that aortic anatomy and embolus size could be important determinants for cardiogenic embolic stroke risk and stroke location, few data exist confirming or characterizing this hypothesis. The objective of this study is to use medical imaging and computational modeling to better understand if aortic anatomy and embolus size influence predilections for cardiogenic embolic transport and right vs. left hemisphere propensity. Anatomically accurate models of the human aorta and branch arteries to the head were reconstructed from computed tomography (CT) angiography of 10 patients. Blood flow was modeled by the Navier-Stokes equations using a well-validated flow solver with physiologic inflow and boundary conditions. Embolic particulate was released from the aortic root and tracked through the common carotid and vertebral arteries for a range of particle sizes. Cardiogenic emboli reaching the carotid and vertebral arteries appeared to have a strong size-destination relationship that varied markedly from expectations based on blood distribution. Observed trends were robust to modeling parameters. A patient's aortic anatomy appeared to significantly influence the probability a cardiogenic particle becomes embolic to the head. Right hemisphere propensity appeared dominant for cardiogenic emboli, which has been confirmed clinically. The predilections discovered through this modeling could represent an important mechanism underlying cardiogenic embolic stroke etiology.

New insights into the understanding of flow dynamics in an in vitro model for abdominal aortic aneurysms

An in vitro dynamics set-up of the flow in a compliant abdominal aortic aneurysm (AAA) model with an anterior posterior asymmetry, aorto-iliac bifurcation, and physiological inlet flow rate and outlet pressure waveforms was developed. The aims were first to show that the structural mechanical behavior of the used material to mimic the AAA wall was similar to this of patients with AAA and then to study the influence of the aorto-iliac bifurcation presence and to study the influence of the imbalanced flow rate in the iliac branches on the AAA flow field. 3D visualizations, never performed in the literature, have clearly put into evidence the development of a vortex ring generated at the AAA proximal neck during the decelerating phase of flow rate, which detaches and progresses downstream during the cardiac cycle, impinges on the anterior wall in the distal AAA region, breaks up, and separates into two vortices of which one rolls on upstream along the anterior wall. 2D particle image velocimetry measurements, swirling strength and enstrophy calculations allowed quantification of the vorticity, vortex trajectory and energy for the different geometrical and hydrodynamical conditions. The main results show that the instant and the intensity of the vortex ring impingement depend on the presence of the aorto-iliac bifurcation, with higher intensity, by about 90%, for an AAA without bifurcation. The imbalance of the flow rates into the iliac branches induces different propagation velocities of the vortex ring and lowers the intensity of the vortex impact by about 60%. The potential influence of the AAA dynamics is discussed in terms of AAA remodeling and rupture.

Virtual and real bench testing of a new percutaneous valve device: a case study

AIMS

To validate patient-specific computational testing of a second-generation device for percutaneous pulmonary valve implantation (PPVI), against realistic in vitro data.

METHODS AND RESULTS

Tests were initially carried out in a simple loading mode, performing a compliance test on a rapid prototyped cylinder. This model was reproduced computationally and validated against the experimental data. A second-generation PPVI stent-graft, with no valve mounted, was then deployed in a simplified cylindrical geometry, measuring its displacement when subjected to a pressure pulse. Experimental and computational measurements were in good agreement. Finally, having selected a patient regarded as unsuitable for first-generation PPVI, but potentially suitable for a second-generation device, the stent-graft was studied in the rapidly prototyped patient-specific right ventricular outflow tract (RVOT). Stent positioning and radial displacements with pulsatile flow were observed in a mock circuit using fluoroscopy imaging. Stent deformation and anchoring were measured both in vitro and computationally. Both tests indicated that the stent was well anchored in the RVOT, especially in the distal position, and its central region was rounded, ensuring, were a valve present, optimal valve function.

CONCLUSION

We suggest that an experimentally validated computational model can be used for preclinical device characterisation and patient selection.

Characterization of the transport topology in patient-specific abdominal aortic aneurysm models

Abdominal aortic aneurysm (AAA) is characterized by disturbed blood flow patterns that are hypothesized to contribute to disease progression. The transport topology in six patient-specific abdominal aortic aneurysms was studied. Velocity data were obtained by image-based computational fluid dynamics modeling, with magnetic resonance imaging providing the necessary simulation parameters. Finite-time Lyapunov exponent (FTLE) fields were computed from the velocity data, and used to identify Lagrangian coherent structures (LCS). The combination of FTLE fields and LCS was used to characterize topological flow features such as separation zones, vortex transport, mixing regions, and flow impingement. These measures offer a novel perspective into AAA flow. It was observed that all aneurysms exhibited coherent vortex formation at the proximal segment of the aneurysm. The evolution of the systolic vortex strongly influences the flow topology in the aneurysm. It was difficult to predict the vortex dynamics from the aneurysm morphology, motivating the application of image-based flow modeling.

Image-based modeling of hemodynamics in coronary artery aneurysms caused by Kawasaki disease

Kawasaki Disease (KD) is the leading cause of acquired pediatric heart disease. A subset of KD patients develops aneurysms in the coronary arteries, leading to increased risk of thrombosis and myocardial infarction. Currently, there are limited clinical data to guide the management of these patients, and the hemodynamic effects of these aneurysms are unknown. We applied patient-specific modeling to systematically quantify hemodynamics and wall shear stress in coronary arteries with aneurysms caused by KD. We modeled the hemodynamics in the aneurysms using anatomic data obtained by multi-detector computed tomography (CT) in a 10-year-old male subject who suffered KD at age 3 years. The altered hemodynamics were compared to that of a reconstructed normal coronary anatomy using our subject as the model. Computer simulations using a robust finite element framework were used to quantify time-varying shear stresses and particle trajectories in the coronary arteries. We accounted for the cardiac contractility and the microcirculation using physiologic downstream boundary conditions. The presence of aneurysms in the proximal coronary artery leads to flow recirculation, reduced wall shear stress within the aneurysm, and high wall shear stress gradients at the neck of the aneurysm. The wall shear stress in the KD subject (2.95-3.81 dynes/sq cm) was an order of magnitude lower than the normal control model (17.10-27.15 dynes/sq cm). Particle residence times were significantly higher, taking 5 cardiac cycles to fully clear from the aneurysmal regions in the KD subject compared to only 1.3 cardiac cycles from the corresponding regions of the normal model. In this novel quantitative study of hemodynamics in coronary aneurysms caused by KD, we documented markedly abnormal flow patterns that are associated with increased risk of thrombosis. This methodology has the potential to provide further insights into the effects of aneurysms in KD and to help risk stratify patients for appropriate medical and surgical interventions.

Simbios: an NIH national center for physics-based simulation of biological structures

Physics-based simulation provides a powerful framework for understanding biological form and function. Simulations can be used by biologists to study macromolecular assemblies and by clinicians to design treatments for diseases. Simulations help biomedical researchers understand the physical constraints on biological systems as they engineer novel drugs, synthetic tissues, medical devices, and surgical interventions. Although individual biomedical investigators make outstanding contributions to physics-based simulation, the field has been fragmented. Applications are typically limited to a single physical scale, and individual investigators usually must create their own software. These conditions created a major barrier to advancing simulation capabilities. In 2004, we established a National Center for Physics-Based Simulation of Biological Structures (Simbios) to help integrate the field and accelerate biomedical research. In 6 years, Simbios has become a vibrant national center, with collaborators in 16 states and eight countries. Simbios focuses on problems at both the molecular scale and the organismal level, with a long-term goal of uniting these in accurate multiscale simulations.

In vivo validation of numerical prediction for turbulence intensity in an aortic coarctation

This paper compares numerical predictions of turbulence intensity with in vivo measurement. Magnetic resonance imaging (MRI) was carried out on a 60-year-old female with a restenosed aortic coarctation. Time-resolved three-directional phase-contrast (PC) MRI data was acquired to enable turbulence intensity estimation. A contrast-enhanced MR angiography (MRA) and a time-resolved 2D PCMRI measurement were also performed to acquire data needed to perform subsequent image-based computational fluid dynamics (CFD) modeling. A 3D model of the aortic coarctation and surrounding vasculature was constructed from the MRA data, and physiologic boundary conditions were modeled to match 2D PCMRI and pressure pulse measurements. Blood flow velocity data was subsequently obtained by numerical simulation. Turbulent kinetic energy (TKE) was computed from the resulting CFD data. Results indicate relative agreement (error ≈10%) between the in vivo measurements and the CFD predictions of TKE. The discrepancies in modeled vs. measured TKE values were within expectations due to modeling and measurement errors.

In vitro validation of finite element analysis of blood flow in deformable models

The purpose of this article is to validate numerical simulations of flow and pressure incorporating deformable walls using in vitro flow phantoms under physiological flow and pressure conditions. We constructed two deformable flow phantoms mimicking a normal and a restricted thoracic aorta, and used a Windkessel model at the outlet boundary. We acquired flow and pressure data in the phantom while it operated under physiological conditions. Next, in silico numerical simulations were performed, and velocities, flows, and pressures in the in silico simulations were compared to those measured in the in vitro phantoms. The experimental measurements and simulated results of pressure and flow waveform shapes and magnitudes compared favorably at all of the different measurement locations in the two deformable phantoms. The average difference between measured and simulated flow and pressure was approximately 3.5 cc/s (13% of mean) and 1.5 mmHg (1.8% of mean), respectively. Velocity patterns also showed good qualitative agreement between experiment and simulation especially in regions with less complex flow patterns. We demonstrated the capabilities of numerical simulations incorporating deformable walls to capture both the vessel wall motion and wave propagation by accurately predicting the changes in the flow and pressure waveforms at various locations down the length of the deformable flow phantoms.