Ethan Kung

An integrated approach to patient-specific predictive modeling for single ventricle heart palliation

In patients with congenital heart disease and a single ventricle (SV), ventricular support of the circulation is inadequate, and staged palliative surgery (usually 3 stages) is needed for treatment. In the various palliative surgical stages individual differences in the circulation are important and patient-specific surgical planning is ideal. In this study, an integrated approach between clinicians and engineers has been developed, based on patient-specific multi-scale models, and is here applied to predict stage 2 surgical outcomes. This approach involves four distinct steps: (1) collection of pre-operative clinical data from a patient presenting for SV palliation, (2) construction of the pre-operative model, (3) creation of feasible virtual surgical options which couple a three-dimensional model of the surgical anatomy with a lumped parameter model (LPM) of the remainder of the circulation and (4) performance of post-operative simulations to aid clinical decision making. The pre-operative model is described, agreeing well with clinical flow tracings and mean pressures. Two surgical options (bi-directional Glenn and hemi-Fontan operations) are virtually performed and coupled to the pre-operative LPM, with the hemodynamics of both options reported. Results are validated against postoperative clinical data. Ultimately, this work represents the first patient-specific predictive modeling of stage 2 palliation using virtual surgery and closed-loop multi-scale modeling.

In Vitro Validation of Finite-Element Model of AAA Hemodynamics Incorporating Realistic Outlet Boundary Conditions

The purpose of this study is to validate numerical simulations of flow and pressure in an abdominal aortic aneurysm (AAA) using phase-contrast magnetic resonance imaging (PCMRI) and an in vitro phantom under physiological flow and pressure conditions. We constructed a two-outlet physical flow phantom based on patient imaging data of an AAA and developed a physical Windkessel model to use as outlet boundary conditions. We then acquired PCMRI data in the phantom while it operated under conditions mimicking a resting and a light exercise physiological state. Next, we performed in silico numerical simulations and compared experimentally measured velocities, flows, and pressures in the in vitro phantom to those computed in the in silico simulations. There was a high degree of agreement in all of the pressure and flow waveform shapes and magnitudes between the experimental measurements and simulated results. The average pressures and flow split difference between experiment and simulation were all within 2%. Velocity patterns showed good agreement between experimental measurements and simulated results, especially in the case of whole-cycle averaged comparisons. We demonstrated methods to perform in vitro phantom experiments with physiological flows and pressures, showing good agreement between numerically simulated and experimentally measured velocity fields and pressure waveforms in a complex patient-specific AAA geometry.

Moving Domain Computational Fluid Dynamics to Interface with an Embryonic Model of Cardiac Morphogenesis

Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis.

A Simulation Protocol for Exercise Physiology in Fontan Patients Using a Closed Loop Lumped-Parameter Model

BACKGROUND Reduced exercise capacity is nearly universal among Fontan patients, though its etiology is not yet fully understood. While previous computational studies have attempted to model Fontan exercise, they did not fully account for global physiologic mechanisms nor directly compare results against clinical and physiologic data. METHODS In this study, we developed a protocol to simulate Fontan lower-body exercise using a closed-loop lumped-parameter model describing the entire circulation. We analyzed clinical exercise data from a cohort of Fontan patients, incorporated previous clinical findings from literature, quantified a comprehensive list of physiological changes during exercise, translated them into a computational model of the Fontan circulation, and designed a general protocol to model Fontan exercise behavior. Using inputs of patient weight, height, and if available, patient-specific reference heart rate (HR) and oxygen consumption, this protocol enables the derivation of a full set of parameters necessary to model a typical Fontan patient of a given body-size over a range of physiologic exercise levels. RESULTS In light of previous literature data and clinical knowledge, the model successfully produced realistic trends in physiological parameters with exercise level. Applying this method retrospectively to a set of clinical Fontan exercise data, direct comparison between simulation results and clinical data demonstrated that the model successfully reproduced the average exercise response of a cohort of typical Fontan patients. CONCLUSION This work is intended to offer a foundation for future advances in modeling Fontan exercise, highlight the needs in clinical data collection, and provide clinicians with quantitative reference exercise physiologies for Fontan patients.

Hemodynamic effects of left pulmonary artery stenosis after superior cavopulmonary connection: a patient-specific multiscale modeling study.

OBJECTIVE Currently, no quantitative guidelines have been established for treatment of left pulmonary artery (LPA) stenosis. This study aims to quantify the effects of LPA stenosis on postoperative hemodynamics for single-ventricle patients undergoing stage II superior cavopulmonary connection (SCPC) surgery, using a multiscale computational approach. METHODS Image data from 6 patients were segmented to produce 3-dimensional models of the pulmonary arteries before stage II surgery. Pressure and flow measurements were used to tune a 0-dimensional model of the entire circulation. Postoperative geometries were generated through stage II virtual surgery; varying degrees of LPA stenosis were applied using mesh morphing and hemodynamics assessed through coupled 0-3-dimensional simulations. To relate metrics of stenosis to clinical classifications, pediatric cardiologists and surgeons ranked the degrees of stenosis in the models. The effects of LPA stenosis were assessed based on left-to-right pulmonary artery flow split ratios, mean pressure drop across the stenosis, cardiac pressure-volume loops, and other clinically relevant parameters. RESULTS Stenosis of >65% of the vessel diameter was required to produce a right pulmonary artery:LPA flow split <30%, and/or a mean pressure drop of >3.0 mm Hg, defined as clinically significant changes. CONCLUSIONS The effects of <65% stenosis on SCPC hemodynamics and physiology were minor and may not justify the increased complexity of adding LPA arterioplasty to the SCPC operation. However, in the longer term, pulmonary augmentation may affect outcomes of the Fontan completion surgery, as pulmonary artery distortion is a risk factor that may influence stage III physiology.

Computational Modeling of Pathophysiologic Responses to Exercise in Fontan Patients

Reduced exercise capacity is nearly universal among Fontan patients. Although many factors have emerged as possible contributors, the degree to which each impacts the overall hemodynamics is largely unknown. Computational modeling provides a means to test hypotheses of causes of exercise intolerance via precisely controlled virtual experiments and measurements. We quantified the physiological impacts of commonly encountered, clinically relevant dysfunctions introduced to the exercising Fontan system via a previously developed lumped-parameter model of Fontan exercise. Elevated pulmonary arterial pressure was observed in all cases of dysfunction, correlated with lowered cardiac output (CO), and often mediated by elevated atrial pressure. Pulmonary vascular resistance was not the most significant factor affecting exercise performance as measured by CO. In the absence of other dysfunctions, atrioventricular valve insufficiency alone had significant physiological impact, especially under exercise demands. The impact of isolated dysfunctions can be linearly summed to approximate the combined impact of several dysfunctions occurring in the same system. A single dominant cause of exercise intolerance was not identified, though several hypothesized dysfunctions each led to variable decreases in performance. Computational predictions of performance improvement associated with various interventions should be weighed against procedural risks and potential complications, contributing to improvements in routine patient management protocol.

Superior performance of continuous over pulsatile flow ventricular assist devices in the single ventricle circulation: A computational study

This study compares the physiological responses of systemic-to-pulmonary shunted single ventricle patients to pulsatile and continuous flow ventricular assist devices (VADs). Performance differences between pulsatile and continuous flow VADs have been clinically observed, but the underlying mechanism remains poorly understood. Six systemic-to-pulmonary shunted single ventricle patients (mean BSA=0.30m2) were computationally simulated using a lumped-parameter network tuned to match patient specific clinical data. A first set of simulations compared current clinical implementation of VADs in single ventricle patients. A second set modified pulsatile flow VAD settings with the goal to optimize cardiac output (CO). For all patients, the best-case continuous flow VAD CO was at least 0.99L/min greater than the optimized pulsatile flow VAD CO (p=0.001). The 25 and 50mL pulsatile flow VADs exhibited incomplete filling at higher heart rates that reduced CO as much as 9.7% and 37.3% below expectations respectively. Optimization of pulsatile flow VAD settings did not achieve statistically significant (p<0.05) improvement to CO. Results corroborate clinical experience that continuous flow VADs produce higher CO and superior ventricular unloading in single ventricle patients. Impaired filling leads to performance degradation of pulsatile flow VADs in the single ventricle circulation.

Does TCPC power loss really affect exercise capacity

To the Editor, We read with interest the article by Khiabani et al ,1 where the authors examine the relationship between power loss in the total cavopulmonary connection (TCPC) and clinical exercise testing. Using an indexed power loss, ‘iPL’, they report that higher iPL correlates with worse minute oxygen consumption and exercise work at anaerobic threshold. Based on this, the authors suggest that power loss in the TCPC could affect exercise performance in a patient with Fontan circulation. In the manuscript, the authors attempt to discover a correlation of the hydraulics of the Fontan circulation with exercise performance. In doing so, they use the unique parameter, iPL, instead of the unadjusted power loss. We believe that this approach is misleading and leads to the wrong conclusion. The term iPL was defined as:![Formula][1] where PL, ρ, Q … [1]: /embed/graphic-1.gif

Multiscale Modeling of Cardiovascular Flows

....................................................................................................................................................... 3

A hemi fontan operation performed by an engineer: Considerations on virtual surgery

Single ventricle circulation, characterized at birth by a rudimentary or absent left or right ventricle, presents a challenging and life-threatening physiological scenario. Surgical palliation aims to restore the balance between systemic and pulmonary blood flow and is staged, each of the three stages presenting the surgeon with different options:- Stage 1 (Norwood procedure) involves different types of shunting to source pulmonary blood flow, or recently a hybrid approach [1];- Stage 2 can involve a superior cavopulmonary connection (Glenn operation) or patching between the right atrium and the pulmonary arteries (Hemi Fontan operation [2]);- Stage 3 involves a total cavopulmonary connection with extracardiac conduit or lateral tunnel, or with novel alternatives such as the Y-graft [3].Copyright