Computational fluid dynamics (CFD) tools have been extensively applied to study the hemodynamics in the total cavopulmonary connection (TCPC) in patients with only a single functioning ventricle. Without the contraction of a sub-pulmonary ventricle, pulsatility of flow through this connection is low and variable across patients, which is usually neglected in most numerical modeling studies. Recent studies suggest that such pulsatility can be non-negligible and can be important in hemodynamic predictions. The goal of this work is to compare the results of an in-house numerical methodology for simulating pulsatile TCPC flow with experimental results. Digital particle image velocimetry (DPIV) was acquired on TCPC in vitro models to evaluate the capability of the CFD tool in predicting pulsatile TCPC flow fields. In vitro hemodynamic measurements were used to compare the numerical prediction of power loss across the connection. The results demonstrated the complexity of the pulsatile TCPC flow fields and the validity of the numerical approach in simulating pulsatile TCPC flow dynamics in both idealized and complex patient specific models.
Total cavopulmonary connection is the result of a series of palliative surgical repairs performed on patients with single ventricle heart defects. The resulting anatomy has complex and unsteady hemodynamics characterized by flow mixing and flow separation. Although varying degrees of flow pulsatility have been observed in vivo, non-pulsatile (time-averaged) boundary conditions have traditionally been assumed in hemodynamic modeling, and only recently have pulsatile conditions been incorporated without completely characterizing their effect or importance. In this study, 3D numerical simulations with both pulsatile and non-pulsatile boundary conditions were performed for 24 patients with different anatomies and flow boundary conditions from Georgia Tech database. Flow structures, energy dissipation rates and pressure drops were compared under rest and simulated exercise conditions. It was found that flow pulsatility is the primary factor in determining the appropriate choice of boundary conditions, whereas the anatomic configuration and cardiac output had secondary effects. Results show that the hemodynamics can be strongly influenced by the presence of pulsatile flow. However, there was a minimum pulsatility threshold, identified by defining a weighted pulsatility index (wPI), above which the influence was significant. It was shown that when wPI < 30%, the relative error in hemodynamic predictions using time-averaged boundary conditions was less than 10% compared to pulsatile simulations. In addition, when wPI < 50, the relative error was less than 20%. A correlation was introduced to relate wPI to the relative error in predicting the flow metrics with non-pulsatile flow conditions.
Neuromusculoskeletal models solve the basic problem of determining how the body moves under the influence of external and internal forces. Existing biomechanical modeling programs often emphasize dynamics with the goal of finding a feed-forward neural program to replicate experimental data or of estimating force contributions or individual muscles. The computation of rigid-body dynamics, muscle forces, and activation of the muscles are often performed separately. We have developed an intrinsically forward computational platform (Neuromechanic, www.neuromechanic.com) that explicitly represents the interdependencies among rigid body dynamics, frictional contact, muscle mechanics, and neural control modules. This formulation has significant advantages for optimization and forward simulation, particularly with application to neural controllers with feedback or regulatory features. Explicit inclusion of all state dependencies allows calculation of system derivatives with respect to kinematic states and muscle and neural control states, thus affording a wealth of analytical tools, including linearization, stability analyses and calculation of initial conditions for forward simulations. In this review, we describe our algorithm for generating state equations and explain how they may be used in integration, linearization, and stability analysis tools to provide structural insights into the neural control of movement.
This study aims to investigate the capability of smoothed particle hydrodynamics (SPH), a fully Lagrangian mesh-free method, to simulate the bulk blood flow dynamics in two realistic left ventricular (LV) models. Three dimensional geometries and motion of the LV, proximal left atrium and aortic root are extracted from cardiac magnetic resonance imaging and multi-slice computed tomography imaging data. SPH simulation results are analyzed and compared with those obtained using a traditional finite volume-based numerical method, and to in vivo phase contrast magnetic resonance imaging and echocardiography data, in terms of the large-scale blood flow phenomena usually clinically measured. A quantitative comparison of the velocity fields and global flow parameters between the in silico models and the in vivo data shows a reasonable agreement, given the inherent uncertainties and limitations in the modeling and imaging techniques. The results indicate the capability of SPH as a promising tool for predicting clinically relevant large-scale LV flow information.
In this study, we present a fully-coupled fluid-structure interaction (FSI) framework that combines smoothed particle hydrodynamics (SPH) and nonlinear finite element (FE) method to investigate the coupled aortic and mitral valves structural response and the bulk intraventricular hemodynamics in a realistic left ventricle (LV) model during the entire cardiac cycle. The FSI model incorporates valve structures that consider native asymmetric leaflet geometries, anisotropic hyperelastic material models and human material properties.
Comparison of FSI results with subject-specific echocardiography data demonstrates that the SPH-FE approach is able to quantitatively predict the opening and closing times of the valves, the mitral leaflet opening and closing angles, and the large-scale intraventricular flow phenomena with a reasonable agreement. Moreover, comparison of FSI results with a LV model without valves reveals substantial differences in the flow field. Peak systolic velocities obtained from the FSI model and the LV model without valves are 2.56 m/s and 1.16 m/s, respectively, compared to the Doppler echo data of 2.17 m/s. The proposed SPH-FE FSI framework represents a further step towards modeling patient-specific coupled LV-valve dynamics, and has the potential to improve our understanding of cardiovascular physiology and to support professionals in clinical decision-making.
The ultimate goal of Fontan surgical planning is to provide additional insights into the clinical decision-making process. In its current state, surgical planning offers an accurate hemodynamic assessment of the pre-operative condition, provides anatomical constraints for potential surgical options, and produces decent post-operative predictions if boundary conditions are similar enough between the pre-operative and post-operative states. Moving forward, validation with post-operative data is a necessary step in order to assess the accuracy of surgical planning and determine which methodological improvements are needed. Future efforts to automate the surgical planning process will reduce the individual expertise needed and encourage use in the clinic by clinicians. As post-operative physiologic predictions improve, Fontan surgical planning will become an more effective tool to accurately model patient-specific hemodynamics.
The complexity of morphogenesis poses a fundamental challenge to understanding the mechanisms governing the formation of biological patterns and structures. Over the past century, numerous processes have been identified as critically contributing to morphogenetic events, but the interplay between the various components and aspects of pattern formation have been much harder to grasp. The combination of traditional biology with mathematical and computational methods has had a profound effect on our current understanding of morphogenesis and led to significant insights and advancements in the field. In particular, the theoretical concepts of reaction-diffusion systems and positional information, proposed by Alan Turing and Lewis Wolpert, respectively, dramatically influenced our general view of morphogenesis, although typically in isolation from one another. In recent years, agent-based modeling has been emerging as a consolidation and implementation of the two theories within a single framework. Agent-based models (ABMs) are unique in their ability to integrate combinations of heterogeneous processes and investigate their respective dynamics, especially in the context of spatial phenomena. In this review, we highlight the benefits and technical challenges associated with ABMs as tools for examining morphogenetic events. These models display unparalleled flexibility for studying various morphogenetic phenomena at multiple levels and have the important advantage of informing future experimental work, including the targeted engineering of tissues and organs.
Central nervous system (CNS) tissue motion of the brain occurs over 30 million cardiac cycles per year due to intracranial pressure differences caused by the pulsatile blood flow and cerebrospinal fluid (CSF) motion within the intracranial space. This motion has been found to be elevated in type 1 Chiari malformation. The impact of CNS tissue motion on CSF dynamics was assessed using a moving-boundary computational fluid dynamics (CFD) model of the cervical-medullary junction (CMJ). The cerebellar tonsils and spinal cord were modeled as rigid surfaces moving in the caudocranial direction over the cardiac cycle. The CFD boundary conditions were based on in vivo MR imaging of a 35-year old female Chiari malformation patient with ~150–300 µm motion of the cerebellar tonsils and spinal cord, respectively. Results showed that tissue motion increased CSF pressure dissociation across the CMJ and peak velocities up to 120 and 60%, respectively. Alterations in CSF dynamics were most pronounced near the CMJ and during peak tonsillar velocity. These results show a small CNS tissue motion at the CMJ can alter CSF dynamics for a portion of the cardiac cycle and demonstrate the utility of CFD modeling coupled with MR imaging to help understand CSF dynamics.
Cardiovascular simulations have great potential as a clinical tool for planning and evaluating patient-specific treatment strategies for those suffering from congenital heart diseases, specifically Fontan patients. However, several bottlenecks have delayed wider deployment of the simulations for clinical use; the main obstacle is simulation cost. Currently, time-averaged clinical flow measurements are utilized as numerical boundary conditions (BCs) in order to reduce the computational power and time needed to offer surgical planning within a clinical time frame. Nevertheless, pulsatile blood flow is observed in vivo, and its significant impact on numerical simulations has been demonstrated. Therefore, it is imperative to carry out a comprehensive study analyzing the sensitivity of using time-averaged BCs. In this study, sensitivity is evaluated based on the discrepancies between hemodynamic metrics calculated using time-averaged and pulsatile BCs; smaller discrepancies indicate less sensitivity. The current study incorporates a comparison between 3D patient-specific CFD simulations using both the time-averaged and pulsatile BCs for 101 Fontan patients. The sensitivity analysis involves two clinically important hemodynamic metrics: hepatic flow distribution (HFD) and indexed power loss (iPL). Paired demographic group comparisons revealed that HFD sensitivity is significantly different between single and bilateral superior vena cava cohorts but no other demographic discrepancies were observed for HFD or iPL. Multivariate regression analyses show that the best predictors for sensitivity involve flow pulsatilities, time-averaged flow rates, and geometric characteristics of the Fontan connection. These predictors provide patient-specific guidelines to determine the effectiveness of analyzing patient-specific surgical options with time-averaged BCs within a clinical time frame.