Patient-specific CFD models for intraventricular flow analysis from 3D ultrasound imaging: Comparison of three clinical cases.

BACKGROUND: As the intracardiac flow field is affected by changes in shape and motility of the heart, intraventricular flow features can provide diagnostic indications. Ventricular flow patterns differ depending on the cardiac condition and the exploration of different clinical cases can provide insights into how flow fields alter in different pathologies. METHODS: In this study, we applied a patient-specific computational fluid dynamics model of the left ventricle and mitral valve, with prescribed moving boundaries based on transesophageal ultrasound images for three cardiac pathologies, to verify the abnormal flow patterns in impaired hearts. One case (P1) had normal ejection fraction but low stroke volume and cardiac output, P2 showed low stroke volume and reduced ejection fraction, P3 had a dilated ventricle and reduced ejection fraction. RESULTS: The shape of the ventricle and mitral valve, together with the pathology influence the flow field in the left ventricle, leading to distinct flow features. Of particular interest is the pattern of the vortex formation and evolution, influenced by the valvular orifice and the ventricular shape. The base-to-apex pressure difference of maximum 2mmHg is consistent with reported data. CONCLUSION: We used a CFD model with prescribed boundary motion to describe the intraventricular flow field in three patients with impaired diastolic function. The calculated intraventricular flow dynamics are consistent with the diagnostic patient records and highlight the differences between the different cases. The integration of clinical images and computational techniques, therefore, allows for a deeper investigation intraventricular hemodynamics in patho-physiology.

Patient-specific CFD simulation of intraventricular haemodynamics based on 3D ultrasound imaging.

BACKGROUND: The goal of this paper is to present a computational fluid dynamic (CFD) model with moving boundaries to study the intraventricular flows in a patient-specific framework. Starting from the segmentation of real-time transesophageal echocardiographic images, a CFD model including the complete left ventricle and the moving 3D mitral valve was realized. Their motion, known as a function of time from the segmented ultrasound images, was imposed as a boundary condition in an Arbitrary Lagrangian-Eulerian framework. RESULTS: The model allowed for a realistic description of the displacement of the structures of interest and for an effective analysis of the intraventricular flows throughout the cardiac cycle. The model provides detailed intraventricular flow features, and highlights the importance of the 3D valve apparatus for the vortex dynamics and apical flow. CONCLUSIONS: The proposed method could describe the haemodynamics of the left ventricle during the cardiac cycle. The methodology might therefore be of particular importance in patient treatment planning to assess the impact of mitral valve treatment on intraventricular flow dynamics.

Unstructured hexahedral mesh generation of complex vascular trees using a multi-block grid-based approach.

The trend towards realistic numerical models of (pathologic) patient-specific vascular structures brings along larger computational domains and more complex geometries, increasing both the computation time and the operator time. Hexahedral grids effectively lower the computational run time and the required computational infrastructure, but at high cost in terms of operator time and minimal cell quality, especially when the computational analyses are targeting complex geometries such as aneurysm necks, severe stenoses and bifurcations. Moreover, such grids generally do not allow local refinements. As an attempt to overcome these limitations, a novel approach to hexahedral meshing is proposed in this paper, which combines the automated generation of multi-block structures with a grid-based method. The robustness of the novel approach is tested on common complex geometries, such as tree-like structures (including trifurcations), stenoses, and aneurysms. Additionally, the performance of the generated grid is assessed using two numerical examples. In the first example, a grid sensitivity analysis is performed for blood flow simulated in an abdominal mouse aorta and compared to tetrahedral grids with a prismatic boundary layer. In the second example, the fluid-structure interaction in a model of an aorta with aortic coarctation is simulated and the effect of local grid refinement is analyzed.

Plasma gasification of refuse derived fuel in a single-stage system using different gasifying agents.

The renewable evolution in the energy industry and the depletion of natural resources are putting pressure on the waste industry to shift towards flexible treatment technologies with efficient materials and/or energy recovery. In this context, a thermochemical conversion method of recent interest is plasma gasification, which is capable of producing syngas from a wide variety of waste streams. The produced syngas can be valorized for both energetic (heat and/or electricity) and chemical (ammonia, hydrogen or liquid hydrocarbons) end-purposes. This paper evaluates the performance of experiments on a single-stage plasma gasification system for the treatment of refuse-derived fuel (RDF) from excavated waste. A comparative analysis of the syngas characteristics and process yields was done for seven cases with different types of gasifying agents (CO2+O2, H2O, CO2+H2O and O2+H2O). The syngas compositions were compared to the thermodynamic equilibrium compositions and the performance of the single-stage plasma gasification of RDF was compared to that of similar experiments with biomass and to the performance of a two-stage plasma gasification process with RDF. The temperature range of the experiment was from 1400 to 1600 K and for all cases, a medium calorific value syngas was produced with lower heating values up to 10.9 MJ/Nm(3), low levels of tar, high levels of CO and H2 and which composition was in good agreement to the equilibrium composition. The carbon conversion efficiency ranged from 80% to 100% and maximum cold gas efficiency and mechanical gasification efficiency of respectively 56% and 95%, were registered. Overall, the treatment of RDF proved to be less performant than that of biomass in the same system. Compared to a two-stage plasma gasification system, the produced syngas from the single-stage reactor showed more favourable characteristics, while the recovery of the solid residue as a vitrified slag is an advantage of the two-stage set-up.

Influence of valve size, orientation and downstream geometry of an aortic BMHV on leaflet motion and clinically used valve performance parameters.

The aim of this study was to reconcile some of our own previous work and the work of others to generate a physiologically realistic numerical simulation environment that allows to virtually assess the performance of BMHVs. The model incorporates: (i) a left ventricular deformable model to generate a physiological inflow to the aortic valve; (ii) a patient-specific aortic geometry (root, arch and descending aorta); (iii) physiological pressure and flow boundary conditions. We particularly studied the influence of downstream geometry, valve size and orientation on leaflet kinematics and functional indices used in clinical routine. Compared to the straight tube geometry, the patient-specific aorta leads to a significant asynchronous movement of the valve, especially during the closing of the valve. The anterior leaflet starts to close first, impacts the casing at the closed position and remains in this position. At the same time, the posterior leaflet impacts the pivoting mechanisms at the fully open position. At the end of systole, this leaflet subsequently accelerates to the closed position, impacting the casing with an angular velocity of approximately -477 rad/s. The valve size greatly influences the transvalvular pressure gradient (TPG), but does not change the overall leaflet kinematics. This is in contrast to changes in valve orientation, where changing valve orientation induces large differences in leaflet kinematics, but the TPG remains approximately the same.

Validation of a numerical FSI simulation of an aortic BMHV by in vitro PIV experiments.

In this paper, a validation of a recently developed fluid-structure interaction (FSI) coupling algorithm to simulate numerically the dynamics of an aortic bileaflet mechanical heart valve (BMHV) is performed. This validation is done by comparing the numerical simulation results with in vitro experiments. For the in vitro experiments, the leaflet kinematics and flow fields are obtained via the particle image velocimetry (PIV) technique. Subsequently, the same case is numerically simulated by the coupling algorithm and the resulting leaflet kinematics and flow fields are obtained. Finally, the results are compared, revealing great similarity in leaflet motion and flow fields between the numerical simulation and the experimental test. Therefore, it is concluded that the developed algorithm is able to capture very accurately all the major leaflet kinematics and dynamics and can be used to study and optimize the design of BMHVs.

Haemodynamic impact of stent-vessel (mal)apposition following carotid artery stenting: mind the gaps!

Carotid artery stenting (CAS) has emerged as a minimally invasive alternative to endarterectomy but its use in clinical treatment is limited due to the post-stenting complications. Haemodynamic actors, related to blood flow in the stented vessel, have been suggested to play a role in the endothelium response to stenting, including adverse reactions such as in-stent restenosis and late thrombosis. Accessing the flow-related shear forces acting on the endothelium in vivo requires space and time resolutions which are currently not achievable with non-invasive clinical imaging techniques but can be obtained from image-based computational analysis. In this study, we present a framework for accurate determination of the wall shear stress (WSS) in a mildly stenosed carotid artery after the implantation of a stent, resembling the commercially available Acculink (Abbott Laboratories, Abbott Park, Illinois, USA). Starting from angiographic CT images of the vessel lumen and a micro-CT scan of the stent, a finite element analysis is carried out in order to deploy the stent in the vessel, reproducing CAS in silico. Then, based on the post-stenting anatomy, the vessel is perfused using a set of boundary conditions: total pressure is applied at the inlet, and impedances that are assumed to be insensitive to the presence of the stent are imposed at the outlets. Evaluation of the CAS outcome from a geometrical and haemodynamic perspective shows the presence of atheroprone regions (low time-average WSS, high relative residence time) colocalised with stent malapposition and stent strut interconnections. Stent struts remain unapposed in the ostium of the external carotid artery disturbing the flow and generating abnormal shear forces, which could trigger thromboembolic events.

FSI simulation of asymmetric mitral valve dynamics during diastolic filling.

In this article, we present a fluid-structure interaction algorithm accounting for the mutual interaction between two rigid bodies. The algorithm was used to perform a numerical simulation of mitral valve (MV) dynamics during diastolic filling. In numerical simulations of intraventricular flow and MV motion, the asymmetry of the leaflets is often neglected. In this study the MV was rendered as two rigid, asymmetric leaflets. The 2D simulations incorporated the dynamic interaction of blood flow and leaflet motion and an imposed subject-specific, transient left ventricular wall movement obtained from ultrasound recordings. By including the full Jacobian matrix in the algorithm, the speed of the simulation was enhanced by more than 20% compared to using a diagonal Jacobian matrix. Furthermore, our results indicate that important features of the flow field may not be predicted by the use of symmetric leaflets or in the absence of an adequate model for the left atrium.

An assessment of ductus venosus tapering and wave transmission from the fetal heart.

Pressure and flow pulsations in the fetal heart propagate through the precordial vein and the ductus venosus (DV) but are normally not transmitted into the umbilical vein (UV). Pulsations in the umbilical vein do occur, however, in early pregnancy and in pathological conditions. Such transmission into the umbilical vein is not well understood. In particular, the effect of the impedance changes in the DV due to its tapered geometry is not known. This paper presents a mathematical model that we developed to study the transmission of pulsations, originating in the fetal heart, through the DV to the umbilical vein. In our model, the tapered geometry of the DV was found to be of minor importance and the only effective reflection site in the DV appears to be at the DV inlet. Differences between the DV inlet and outlet flow were also found to be minor for medium to large umbilical vein-DV diameter ratios. Finally, the results of a previously proposed lumped model were found to agree well with the present model of the DV-umbilical vein bifurcation.

Optimisation of solute transport in dialysers using a three-dimensional finite volume model.

Dialyser manufacturers only provide limited information about mass removal under well-defined flow and solute conditions in commercially available dialysers for hemodialysis. This computational study aimed at assessing the solute transport efficiency in a dialyser for different geometries (fiber lengths and diameters). A three-dimensional finite volume model of a single fiber in a high flux polysulphone dialyser (Fresenius F60) was developed. Different equations describe blood and dialysate flow (Navier-Stokes), radial filtration flow (Darcy) and solute transport (convection-diffusion). Fluid and membrane properties were derived from in vitro and in vivo tests as well as from literature data. Urea (MW60) was used as marker to simulate small molecule removal, while middle molecule transport was modelled using vitamin B12 (MW1355) and inulin (MW5200). Keeping the fluid velocity in a single fiber constant, fiber diameter and length were changed in a wide range for evaluation of solute removal efficiency. Clearances were found enhanced by 13% (urea), 50% (vitamin B12) and 89% (inulin) for a fiber twice as long as a standard one and by 5.5% (vitamin B12) and 21% (inulin) for a fiber diameter of 150 mum instead of 200 mum. The impact of fiber dimensions was more pronounced for the middle molecules compared to urea.

Wave intensity analysis of left ventricular filling.

Wave intensity analysis (WIA) is a powerful technique to study pressure and flow velocity waves in the time domain in vascular networks. The method is based on the analysis of energy transported by the wave through computation of the wave intensity dI = dPdU, where dP and dU denote pressure and flow velocity changes per time interval, respectively. In this study we propose an analytical modification to the WIA so that it can be used to study waves in conditions of time varying elastic properties, such as the left ventricle (LV) during diastole. The approach is first analytically elaborated for a one-dimensional elastic tube-model of the left ventricle with a time-dependent pressure-area relationship. Data obtained with a validated quasi-three dimensional axi-symmetrical model of the left ventricle are employed to demonstrate this new approach. Along the base-apex axis close to the base wave intensity curves are obtained, both using the standard method and the newly proposed modified method. The main difference between the standard and modified wave intensity pattern occurs immediately after the opening of the mitral valve. Where the standard WIA shows a backward expansion wave, the modified analysis shows a forward compression wave. The proposed modification needs to be taken into account when studying left ventricular relaxation, as it affects the wave type.

Validation of a fluid-structure interaction model of a heart valve using the dynamic mesh method in fluent.

Simulations of coupled problems such as fluid-structure interaction (FSI) are becoming more and more important for engineering purposes. This is particularly true when modeling the aortic valve, where the FSI between the blood and the valve determines the valve movement and the valvular hemodynamics. Nevertheless only a few studies are focusing on the opening and closing behavior during the ejection phase (systole). In this paper, we present the validation of a FSI model using the dynamic mesh method of Fluent for the two-dimensional (2D) simulation of mechanical heart valves during the ejection phase of the cardiac cycle. The FSI model is successfully validated by comparing simulation results to experimental data obtained from in vitro studies using a CCD camera.

Validation of the coupling of magnetic resonance imaging velocity measurements with computational fluid dynamics in a U bend.

Magnetic resonance imaging (MRI) can be used in vivo in combination with computational fluid dynamics (CFD) to derive velocity profiles in space and time and accordingly, pressure drop and wall shear stress distribution in natural or artificial vessel segments. These hemodynamic data are difficult or impossible to acquire directly in vivo. Therefore, research has been performed combining MRI and CFD for flow simulations in flow phantoms, such as bends or anastomoses, and even in human vessels such as the aorta, the carotid, and the abdominal bifurcation. There is, however, no unanimity concerning the use of MRI velocity measurements as input for the inflow boundary condition of a CFD simulation. In this study, different input possibilities for the inflow boundary conditions are compared. MRI measurements of steady and pulsatile flow were performed on a U bend phantom, representing the aorta geometry. PAMFLOW (ESI Software, Krimpen aan den Ussel, The Netherlands), an industrial CFD software package, was used to solve the Navier-Stokes equations for incompressible flow. Three main parameters were found to influence the choice of an inflow boundary condition type. First, the flow rate through a vessel should be exact, since it proves to be a determining factor for the accuracy of the velocity profile. The other decisive parameters are the physiology of the flow profile and the required computer processing unit time. Our comparative study indicates that the best way to handle an inflow boundary condition is to use the velocities measured by MRI at the inflow plane as being fixed velocities. However, before using these MRI velocity data, they first should be corrected for the partial volume effect by filtering and second scaled in order to obtain the correct flow rate. This implies that a reliable flow rate measurement absolutely is needed for CFD calculations based on MRI velocity measurements.

Diastolic filling and pressure imaging: taking advantage of the information in a colour M-mode Doppler image.

The ability to derive non-invasively information on left ventricular diastolic function on one hand and pressure gradients on the other hand, makes Doppler ultrasound a very attractive tool in clinical practice. However, the limitations of the standard Doppler approaches in differentiating between normal and pseudonormal filling patterns, together with the limitations of the simplified Bernoulli equation for assessing pressure gradients, are well described. In this manuscript the role of colour M-mode Doppler echocardiography as a tool that can overcome these limitations is discussed. Relevant key concepts of the haemodynamics of left ventricular filling and its relationship with colour M-mode Doppler echocardiography are introduced.

Computer simulation of intraventricular flow and pressure gradients during diastole.

A two-dimensional axisymmetric computer model is developed for the simulation of the filling flow in the left ventricle (LV). The computed results show that vortices are formed during the acceleration phases of the filling waves. During the deceleration phases these are amplified and convected into the ventricle. The ratio of the maximal blood velocity at the mitral valve (peak E velocity) to the flow wave propagation velocity (WPV) of the filling wave is larger than 1. This hemodynamic behavior is also observed in experiments in vitro (Steen and Steen, 1994, Cardiovasc. Res., 28, pp. 1821-1827) and in measurements in vivo with color M-mode Doppler echocardiography (Stugaard et al., 1994, J. Am. Coll. Cardiol., 24, 663-670). Computed intraventricular pressure profiles are similar to observed profiles in a dog heart (Courtois et al., 1988, Circulation, 78, pp. 661-671). The long-term goal of the computer model is to study the predictive value of noninvasive parameters (e.g., velocities measured with Doppler echocardiography) on invasive parameters (e.g., pressures, stiffness of cardiac wall, time constant of relaxation). Here, we show that higher LV stiffness results in a smaller WPV for a given peak E velocity. This result may indicate an inverse relationship between WPV and LV stiffness, suggesting that WPV may be an important noninvasive index to assess LV diastolic stiffness, LV diastolic pressure and thus atrial pressure (preload).

A three-dimensional analysis of flow in the pivot regions of an ATS bileaflet valve.

Bileaflet heart valves are currently the most commonly implanted type of mechanical prosthetic valve, because of their low transvalvular pressure drop, centralised flow and durability. However, in common with all mechanical heart valves, implanted bileaflet valves show an inherent tendency for blood clot formation at the valve site. Fluid dynamical phenomena associated with blood clotting are elevated blood shear stresses and regions of persistent blood recirculation, particularly when both occur together. Using three-dimensional CFD modelling, combined with enlarged scale experimental modelling, we investigated the blood flow through the ATS bileaflet valve during forward flow, with particular attention to the leaflet pivot regions. Recirculating regions were found both within and downstream of the valve housing ring. Qualitative assessment of the entire cardiac cycle suggested that recirculating blood within the housing ring will be washed away whilst the valve is closed, but as with all bileaflet valve designs recirculating blood downstream of the valve may have a residence time much longer than one cardiac cycle.

Left-ventricular pressure gradients: a computer-model simulation.

Both invasive left-ventricular pressure measurements and non-invasive colour M-mode echographic measurements have shown the existence of intraventricular pressure gradients (IVPGs) during early filling. The mechanisms responsible for these IVPG cannot be completely explained by the experiments. Therefore a one-dimensional numerical model is developed and validated. The model describes filling (both velocities and pressures) along a left ventricular (LV) base-apex axis. Blood-wall interaction in the left ventricle with moving boundaries is taken into account. The computational results for a canine heart indicate that the observed IVPGs during filling are the consequence of a complex interaction between, on the one hand, pressure waves travelling in the LV and, on the other hand, LV geometry, relaxation and compliance. The computational results indicate the pressure dependency of wavespeed (0.77-1.90 m-1 s) for different mean intraventricular pressures (0.88-5.00 mmHg) and IVPGs up to 2 mmHg, independent of the ratio of end systolic volume and equilibrium volume. Increasing relaxation rate not only decreases minimum basal pressure (2.8 instead of 3.6 mmHg) but also has a strong influence on the time delay between the minimum basal and apical pressures (14 ms instead of 49 ms). The results sustain the hypothesis that pressure-wave propagation determines IVPGs and that IVPGs are no proof of elastic recoil.