Search

Objective: We propose a novel method to obtainmap patient-specific blood velocity profiles (obtained from imaging data such as 2D flow MRI or 3D colour Doppler ultrasound) and map them to geometric vascular models suitable to perform CFD simulations of haemodynamics. We describe the implementation and utilisation of the method within an open-source computational hemodynamics simulation software (CRIMSON). Methods: tThe proposed method establishes point-wise correspondences between the contour of a fixed geometric model and time-varying contours containing the velocity image data, from which a continuous, smooth and cyclic deformation field is calculated. Our methodology is validated using synthetic data, and demonstrated using two different in-vivo aortic velocity datasets: a healthy subject with normal tricuspid valve and a patient with bicuspid aortic valve. Results: We compare the performance of our method with results obtained with the state-of-the-art Schwarz-Christoffel method, in terms of preservation of velocities and execution time. Our method is as accurate as the Schwarz-Christoffel method, while being over 8 times faster. The proposed method can preserve either the flow rate or the velocity field through the surface, and can cope with inconsistencies in motion and contour shape. Conclusions: Our results show that the method is as accurate as the Schwarz-Christoffel method in terms of maintaining the velocity distributions, while being more computationally efficient.Our mapping method can accurately preserve either the flow rate or the velocity field through the surface, and can cope with inconsistencies in motion and contour shape. Significance: The proposed method and its integration into the CRIMSON software enable a streamlined approach towards incorporating more patient-specific data in blood flow simulations.

A major challenge in constructing three dimensional patient specific hemodynamic models is the calibration of model parameters to match patient data on flow, pressure, wall motion, etc. acquired in the clinic. Current workflows are manual and time-consuming. This work presents a flexible computational framework for model parameter estimation in cardiovascular flows that relies on the following fundamental contributions. (i) A Reduced-Order Unscented Kalman Filter (ROUKF) model for data assimilation for wall material and simple lumped parameter network (LPN) boundary condition model parameters. (ii) A constrained least squares augmentation (ROUKF-CLS) for more complex LPNs. (iii) A “Netlist” implementation, supporting easy filtering of parameters in such complex LPNs. The ROUKF algorithm is demonstrated using non-invasive patient-specific data on anatomy, flow and pressure from a healthy volunteer. The ROUKF-CLS algorithm is demonstrated using synthetic data on a coronary LPN. The methods described in this paper have been implemented as part of the CRIMSON hemodynamics software package.

In this work, we describe the CRIMSON (CardiovasculaR Integrated Modelling and SimulatiON) software environment. CRIMSON provides a powerful, customizable and user-friendly system for performing three-dimensional and reduced-order computational haemodynamics studies via a pipeline which involves: 1) segmenting vascular structures from medical images; 2) constructing analytic arterial and venous geometric models; 3) performing finite element mesh generation; 4) designing, and 5) applying boundary conditions; 6) running incompressible Navier-Stokes simulations of blood flow with fluid-structure interaction capabilities; and 7) post-processing and visualizing the results, including velocity, pressure and wall shear stress fields. A key aim of CRIMSON is to create a software environment that makes powerful computational haemodynamics tools accessible to a wide audience, including clinicians and students, both within our research laboratories and throughout the community. The overall philosophy is to leverage best-in-class open source standards for medical image processing, parallel flow computation, geometric solid modelling, data assimilation, and mesh generation. It is actively used by researchers in Europe, North and South America, Asia, and Australia. It has been applied to numerous clinical problems; we illustrate applications of CRIMSON to real-world problems using examples ranging from pre-operative surgical planning to medical device design optimization., Author summary We provide the first full presentation in the literature of CRIMSON, the Cardiovascular Integrated Modelling and Simulation Package. CRIMSON consists of a graphical user interface desktop computer program for creating geometric models of blood vessels from medical imaging scans, specifying parameters such as the stiffness of the artery walls, the resistance of connected vessels which are not visible on the scans, and determining the appropriate parameters for all aspects of the model. CRIMSON additionally consists of the Flowsolver, a high-performance computing package which simulates the flow of blood through the models created in the graphical user interface. Whilst several packages which can simulate blood flow exist, most target general fluid simulations, and this lack of specialisation means that blood flow simulation is harder to perform, and can require ad hoc (and potentially scientifically-limiting) workflow decisions. CRIMSON’s specialisation deals with these problems, as well as presenting a number of unique features which are unavailable elsewhere.

Objectives Aortic arch repair remains a major surgical challenge. Multiple manufacturers are developing branched endografts for Zone 0 endovascular repair, extending the armamentarium for minimally invasive treatment of aortic arch pathologies. We hypothesize that the design of the Zone 0 endograft has a significant impact on the postoperative haemodynamic performance, particularly in the cervical arteries. The goal of our study was to compare the postoperative haemodynamic performance of different Zone 0 endograft designs. Methods Patient-specific, clinically validated, computational fluid dynamics simulations were performed in a 71-year-old woman with a 6.5-cm saccular aortic arch aneurysm. Additionally, 4 endovascular repair scenarios using different endograft designs were created. Haemodynamic performance was evaluated by calculation of postoperative changes in blood flow and platelet activation potential (PLAP) in the cervical arteries. Results Preoperative cervical blood flow and mean PLAP were 1080 ml/min and 151.75, respectively. Cervical blood flow decreased and PLAP increased following endovascular repair in all scenarios. Endografts with 2 antegrade inner branches performed better compared to single-branch endografts. Scenario 3 performed the worst with a decrease in the total cervical blood flow of 4.8%, a decrease in the left hemisphere flow of 6.7% and an increase in the mean PLAP of 74.3%. Conclusions Endograft design has a significant impact on haemodynamic performance following Zone 0 endovascular repair, potentially affecting cerebral blood flow during follow-up. Our results demonstrate the use of computational modelling for virtual testing of therapeutic interventions and underline the need to monitor the long-term outcomes in this cohort of patients.

International audience; In this paper, we perform a verification study of the Coupled-Momentum Method (CMM), a 3D fluid-structure interaction (FSI) model which uses a thin linear elastic membrane and linear kinematics to describe the mechanical behavior of the vessel wall. The verification of this model is done using Womersley's deformable wall analytical solution for pulsatile flow in a semi-infinite cylindrical vessel. This solution is, under certain premises, the analytical solution of the CMM and can thus be used for model verification. For the numerical solution, we employ an impedance boundary condition to define a reflection-free outflow boundary condition and thus mimic the physics of the analytical solution, which is defined on a semi-infinite domain. We first provide a rigorous derivation of Womersley's deformable wall theory via scale analysis. We then illustrate different characteristics of the analytical solution such as space-time wave periodicity and attenuation. Finally, we present the verification tests comparing the CMM with Womersley's theory.

This paper presents a new mathematical model of the dynamic control of coronary resistance. It shows a remarkable ability to predict coronary flow in an exercising patient, which would otherwise be impossible, and provides new insight into the purpose and action of the coronary flow control systems. It is applicable as part of a controlled boundary condition at the coronary outlets of a three-dimensional Navier-Stokes simulation of hemodynamics., This work presents a mathematical model of the metabolic feedback and adrenergic feedforward control of coronary blood flow that occur during variations in the cardiac workload. It is based on the physiological observations that coronary blood flow closely follows myocardial oxygen demand, that myocardial oxygen debts are repaid, and that control oscillations occur when the system is perturbed and so are phenomenological in nature. Using clinical data, we demonstrate that the model can provide patient-specific estimates of coronary blood flow changes between rest and exercise, requiring only the patient's heart rate and peak aortic pressure as input. The model can be used in zero-dimensional lumped parameter network studies or as a boundary condition for three-dimensional multidomain Navier-Stokes blood flow simulations. For the first time, this model provides feedback control of the coronary vascular resistance, which can be used to enhance the physiological accuracy of any hemodynamic simulation, which includes both a heart model and coronary arteries. This has particular relevance to patient-specific simulation for which heart rate and aortic pressure recordings are available. In addition to providing a simulation tool, under our assumptions, the derivation of our model shows that β-feedforward control of the coronary microvascular resistance is a mathematical necessity and that the metabolic feedback control must be dependent on two error signals: the historical myocardial oxygen debt, and the instantaneous myocardial oxygen deficit.

Numerical simulations of cardiovascular mass transport pose significant challenges due to the wide range of P\'eclet numbers and backflow at Neumann boundaries. In this paper we present and discuss several numerical tools to address these challenges in the context of a stabilized finite element computational framework. To overcome numerical instabilities when backflow occurs at Neumann boundaries, we propose an approach based on the prescription of the total flux. In addition, we introduce a "consistent flux" outflow boundary condition and demonstrate its superior performance over the traditional zero diffusive flux boundary condition. Lastly, we discuss discontinuity capturing (DC) stabilization techniques to address the well-known oscillatory behavior of the solution near the concentration front in advection-dominated flows.We present numerical examples in both idealized and patient-specific geometries to demonstrate the efficacy of the proposed procedures. The three contributions dis-cussed in this paper enable to successfully address commonly found challenges when simulating mass transport processes in cardiovascular flows.

For babies born with hypoplastic left heart syndrome, several open-heart surgeries are required. During Stage I, a Norwood procedure is performed to construct an appropriate circulation to both the systemic and the pulmonary arteries. The pulmonary arteries receive flow from the systemic circulation, often by using a Blalock-Taussig (BT) shunt between the innominate artery and the right pulmonary artery. This procedure causes significantly disturbed flow in the pulmonary arteries. In this study, we use computational hemodynamic simulations to demonstrate its capacity for examining the properties of the flow through and near the BT shunt. Initially, we construct a computational model which produces blood flow and pressure measurements matching the clinical magnetic resonance imaging (MRI) and catheterization data. Achieving this required us to determine the level of BT shunt occlusion; because the occlusion is below the MRI resolution, this information is difficult to recover without the aid of computational simulations. We determined that the shunt had undergone an effective diameter reduction of 22% since the time of surgery. Using the resulting geometric model, we show that we can computationally reproduce the clinical data. We then replace the BT shunt by with a hypothetical alternative shunt design with a flare at the distal end. Investigation of the impact of the shunt design reveals that the flare can increase pulmonary pressure by as much as 7%, and flow by as much as 9% in the main pulmonary branches, which may be beneficial to the pulmonary circulation.

Computational fluid dynamics (CFD) provides a noninvasive method to functionally assess aortic hemodynamics. The thoracic aorta has an anatomically complex inlet comprising of the aortic valve and root, which is highly prone to different morphologies and pathologies. We investigated the effect of using patient-specific (PS) inflow velocity profiles compared to idealized profiles based on the patient's flow waveform. A healthy 31 yo with a normally functioning tricuspid aortic valve (subject A), and a 52 yo with a bicuspid aortic valve (BAV), aortic valvular stenosis, and dilated ascending aorta (subject B) were studied. Subjects underwent MR angiography to image and reconstruct three-dimensional (3D) geometric models of the thoracic aorta. Flow-magnetic resonance imaging (MRI) was acquired above the aortic valve and used to extract the patient-specific velocity profiles. Subject B's eccentric asymmetrical inflow profile led to highly complex velocity patterns, which were not replicated by the idealized velocity profiles. Despite having identical flow rates, the idealized inflow profiles displayed significantly different peak and radial velocities. Subject A's results showed some similarity between PS and parabolic inflow profiles; however, other parameters such as Flowasymmetry were significantly different. Idealized inflow velocity profiles significantly alter velocity patterns and produce inaccurate hemodynamic assessments in the thoracic aorta. The complex structure of the aortic valve and its predisposition to pathological change means the inflow into the thoracic aorta can be highly variable. CFD analysis of the thoracic aorta needs to utilize fully PS inflow boundary conditions in order to produce truly meaningful results.

We present the possible utility of computational fluid dynamics in the assessment of thrombus formation and virtual surgical planning illustrated in a patient with aortic thrombus in a kinked ascending aortic graft following thoracic endovascular aortic repair.A patient-specific three-dimensional model was built from computed tomography. Additionally, we modeled 3 virtual aortic interventions to assess their effect on thrombosis potential: (1) open surgical repair, (2) conformable endografting, and (3) single-branched endografting. Flow waveforms were extracted from echocardiography and used for the simulations. We used the computational index termed platelet activation potential (PLAP) representing accumulated shear rates of fluid particles within a fluid domain to assess thrombosis potential.The baseline model revealed high PLAP in the entire arch (119.8 ± 42.5), with significantly larger PLAP at the thrombus location (125.4 ± 41.2, p0.001). Surgical repair showed a 37% PLAP reduction at the thrombus location (78.6 ± 25.3, p0.001) and a 24% reduction in the arch (91.6 ± 28.9, p0.001). Single-branched endografting reduced PLAP in the thrombus region by 20% (99.7 ± 24.6, p0.001) and by 14% in the arch (103.8 ± 26.1, p0.001), whereas a more conformable endograft did not have a profound effect, resulting in a modest 4% PLAP increase (130.6 ± 43.7, p0.001) in the thrombus region relative to the baseline case.Regions of high PLAP were associated with aortic thrombus. Aortic repair resolved pathologic flow patterns, reducing PLAP. Branched endografting also relieved complex flow patterns reducing PLAP. Computational fluid dynamics may assist in the prediction of aortic thrombus formation in hemodynamically complex cases and help guide repair strategies.