Numerical Modeling of Intraventricular Flow during Diastole after Implantation of BMHV

This work presents a numerical simulation of intraventricular flow after the implantation of a bileaflet mechanical heart valve at the mitral position. The left ventricle was simplified conceptually as a truncated prolate spheroid and its motion was prescribed based on that of a healthy subject. The rigid leaflet rotation was driven by the transmitral flow and hence the leaflet dynamics were solved using fluid-structure interaction approach. The simulation results showed that the bileaflet mechanical heart valve at the mitral position behaved similarly to that at the aortic position. Sudden area expansion near the aortic root initiated a clockwise anterior vortex, and the continuous injection of flow through the orifice resulted in further growth of the anterior vortex during diastole, which dominated the intraventricular flow. This flow feature is beneficial to preserving the flow momentum and redirecting the blood flow towards the aortic valve. To the best of our knowledge, this is the first attempt to numerically model intraventricular flow with the mechanical heart valve incorporated at the mitral position using a fluid-structure interaction approach. This study facilitates future patient-specific studies.     

Modeling mitral valve stenosis: A parametric study on the stenosis severity level

New computational techniques providing more accurate representation of human heart pathologies could help uncovering relevant physical phenomena and improve the outcome of medical therapies. In this framework, the present work describes an efficient computational model for the evaluation of the ventricular flow alteration in presence of mitral valve stenosis. The model is based on the direct numerical simulation of the Navier–Stokes equations two-way coupled with a structural solver for the left ventricle and mitral valve dynamics. The presence of mitral valve stenosis is mimicked by a single-parameter constraint acting on the kinematics of the mitral leaflets.Four different degrees of mitral valve stenosis are considered focusing on the hemodynamic alterations occurring in pathologic conditions. The mitral jet, generated during diastole, is seen to shrink and strengthen when the stenosis gets more severe. As a consequence, the kinetic energy of the flow, the tissues shear stresses, the transvalvular pressure drop and mitral regurgitation increase. It results that, as the stenosis severity level increases, the geometric and effective orifice areas decrease up to 50% with respect the normal case due to the reduced leaflets mobility and stronger blood acceleration during the diastolic phase. The modified intraventricular hemodynamics is also related to a stronger pressure gradient that, for severe stenosis, can be more than ten times larger than the healthy valve case. These computational results are fully consistent with the available clinical literature and open the way to the virtual assessment of surgical procedures and to the evaluation of prosthetic devices.     

Role of diastolic vortices in flow and energy dynamics during systolic ejection

MRI-based computational fluid dynamics simulations were performed in the left ventricles of two adult porcine subjects with varying physiological states (before and after an induced infarction). The hypothesis that diastolic vortices store kinetic energy and assist systolic ejection was tested, by performing systolic simulations in the presence and absence of diastolic vortices. The latter was achieved by reinitializing the entire velocity field to be zero at the beginning of systole. A rudimentary prescribed motion model of a mitral valve was included in the simulations to direct the incoming mitral jet towards the apex. Results showed that the presence or absence of diastolic vortex rings had insignificant impact on the energy expended by walls of the left ventricles for systolic ejection for both the porcine subjects, under all physiological conditions. Although substantial kinetic energy was stored in diastolic vortices by end diastole, it provided no appreciable savings during systolic ejection, and most likely continued to complete dissipation during systole. The role of diastolic vortices in apical washout was investigated by studying the cumulative mass fraction of passive dye that was ejected during systole in the presence and absence of vortices. Results indicated that the diastolic vortices play a crucial role in ensuring efficient washout of apical blood during systolic ejection.     

A theoretical model for the noninvasive assessment of the transmitral pressure-flow relation.

The purpose of this paper is to formulate from the equations of fluid mechanics an equation which describes the transmitral pressure-flow relationship. According to the linear momentum equation applied to the atrioventricular coupling, the left-atrium-left-ventricle pressure difference (Pa-Pv) can be written as Pa-P v = A delta v/delta t + B v 2 + C v, where v is the transmitral blood velocity and A, B, and C are variables related to the geometry of the atrium, ventricle and mitral orifice, respectively. Based on this theory, Pa-Pv is calculated noninvasively in a patient with a nonobstructive mitral valve. Mitral flow and cardiac dimensions recorded by Doppler echocardiography are digitized and analyzed. Calculation shows that Pa-Pv reaches its peak value at the time of flow peak acceleration and has already considerably decreased at the time of peak velocity. The time course of calculated Pa-Pv is in close agreement with the published experimental catherization data. Numerical computation of early diastolic left atrium and left ventricle pressure curves based on the experimental data of others for the time constant of left ventricular relaxation, left atrial and ventricular chambers stiffness constants, combined with sine-waveform-simulated mitral flow, verifies the time course and the magnitude of Pa-Pv as predicted from flow equations. This paper provides a theoretical method for the noninvasive assessment of the transmitral pressure-flow relationship using ultrasound technique and might help to achieve a better understanding of the diastolic function as assessed by Doppler echocardiography.     

A computational procedure for prediction of structural effects of edge-to-edge repair on mitral valve

Edge-to-edge technique is a surgical procedure for the correction of mitral valve leaflets prolapse by suturing the edge of the prolapsed leaflet to the free edge of the opposing one. Suture presence modifies valve mechanical behavior and orifice flow area in the diastolic phase, when the valve opens and blood flows into the ventricle. In the present work, in order to support identification of potentially critical conditions, a computational procedure is described to evaluate the effects of changing suture length and position in combination with valve size and shape. The procedure is based on finite element method analyses applied to a range of different mitral valves, investigating for each configuration the influence of repair on functional parameters, such as mitral valve orifice area and transvalvular pressure gradient, and on structural parameters, such as stress in the leaflets and stitch tension. This kind of prediction would ideally require a coupled fluid-structural analysis, where the interactions between blood flows and mitral apparatus deformation are simultaneously considered. In the present study, however, an alternative approach is proposed, in which results obtained by purely structural finite element analyses are elaborated and interpreted taking into account the Bernoulli type equations available in literature to describe blood flow through mitral orifice. In this way, the effects of each parameter in terms of orifice flow area, suture loads, and leaflets stresses can be expressed as functions of atrioventricular pressure gradient and then correlated to blood flow rate. Results obtained by using this procedure for different configurations are finally discussed.     

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).     

Model and influence of mitral valve opening during the left ventricular filling

The flow inside a model left ventricle during filling (diastole) is simulated by the numerical solution of the equations of motion under the axisymmetric approximation. The left ventricle is taken with a truncated ellipsoid geometry, and a simple conceptual model is introduced to simulate the presence of the moving mitral valve. A relevant role during the left ventricle diastolic flow, as already discussed by other authors, is played by the travelling vortex wake that is formed from the transmitral jet during the early filling acceleration phase. The presence of a moving valve is found to produce a non-simultaneous spatial development of the entering bulk flow and a slightly more complex vortex wake structure; the results are discussed in comparison with fixed valve ones. They are analysed also in terms of M-mode representation suggesting a physical interpretation of the pattern detected in the clinical measurements that extends the one given previously on the basis of fixed valve models.     

Validation of a 3D computational fluid-structure interaction model simulating flow through an elastic aperture

This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two.     

Jet collisions and vortex reversal in the human left ventricle

Unnatural dynamics of the notorious vortex in the left ventricle is often associated with cardiac disease. Understanding how different cardiac diseases alter the flow physics in the left ventricle may therefore provide a powerful tool for disease detection. In this work, the fluid dynamics in the left ventricle subject to different severities of aortic regurgitation is experimentally investigated by performing time-resolved particle image velocimetry in a left heart duplicator. Diastolic vortex reversal was observed in the left ventricle accompanied by an increase in viscous energy dissipation. Vortex dynamics and energy dissipation may provide useful insights on sub-optimal flow patterns in the left ventricle.     

Diastolic right ventricular filling vortex in normal and volume overload states

Functional imaging computational fluid dynamics simulations of right ventricular (RV) inflow fields were obtained by comprehensive software using individual animal-specific dynamic imaging data input from three-dimensional (3-D) real-time echocardiography (RT3D) on a CRAY T-90 supercomputer. Chronically instrumented, lightly sedated awake dogs (n = 7) with normal wall motion (NWM) at control and normal or diastolic paradoxical septal motion (PSM) during RV volume overload were investigated. Up to the E-wave peak, instantaneous inflow streamlines extended from the tricuspid orifice to the RV endocardial surface in an expanding fanlike pattern. During the descending limb of the E-wave, large-scale (macroscopic or global) vortical motions ensued within the filling RV chamber. Both at control and during RV volume overload (with or without PSM), blood streams rolled up from regions near the walls toward the base. The extent and strength of the ring vortex surrounding the main stream were reduced with chamber dilatation. A hypothesis is proposed for a facilitatory role of the diastolic vortex for ventricular filling. The filling vortex supports filling by shunting inflow kinetic energy, which would otherwise contribute to an inflow-impeding convective pressure rise between inflow orifice and the large endocardial surface of the expanding chamber, into the rotational kinetic energy of the vortical motion that is destined to be dissipated as heat. The basic information presented should improve application and interpretation of noninvasive (Doppler color flow mapping, velocity-encoded cine magnetic resonance imaging, etc.) diastolic diagnostic studies and lead to improved understanding and recognition of subtle, flow-associated abnormalities in ventricular dilatation and remodeling.     

Determinants of the occurrence of vortex rings in the left ventricle during diastole

This study employs classical inviscid fluid dynamics theory to investigate whether LV diastolic inflow volume and the size of the LV play a role in vortex ring formation. Fluid injection across an orifice into a large container results in the generation of a vortex ring having a constant size and speed. Relations between the vortex size and speed and the injection were obtained by applying conservation laws regarding kinetic energy, impulse and vorticity; the initial state was computed using a bolus injection model, and the final state by using the Kelvin vortex model. An important parameter in the equations is the relative injection length, i.e., the ratio of the length of the injected bolus and the radius of the orifice (L/R). Its estimated highest value in man, L/R = 15, produces a rather thick vortex ring (relative thickness 0.77). Comparable results following from the Hill vortex model convinced us that the Kelvin vortex model can be applied in the whole range of injection lengths in the human left ventricle. In an in vitro model it is shown experimentally that vortex rings can be generated for L/R in the range from 2 to 16. The measured traveling speed of the vortex ring is in fair agreement with the theory, as well as the ring radius for large injections. A vortex ring located in a narrow channel cannot reach its proper traveling speed. The method of images is used to estimate the speed reduction of vortex rings within a cylinder. It turns out that propagation of vortex rings is possible when the ratio of orifice to cylinder radius is less than about 0.5.(ABSTRACT TRUNCATED AT 250 WORDS)     

Intraventricular flow patterns and stasis in the LVAD-assisted heart

Left ventricular assist device (LVAD) support disrupts the natural blood flow path through the heart, introducing flow patterns associated with thrombosis, especially in the presence of medical devices. The aim of this study was to quantitatively evaluate the flow patterns in the left ventricle (LV) of the LVAD-assisted heart, with a focus on alterations in vortex development and stasis. Particle image velocimetry of a LVAD-supported LV model was performed in a mock circulatory loop. In the Pre-LVAD flow condition, a vortex ring initiating from the LV base migrated toward the apex during diastole and remained in the LV by the end of ejection. During LVAD support, vortex formation was relatively unchanged although vortex circulation and kinetic energy increased with LVAD speed, particularly in systole. However, as pulsatility decreased and aortic valve opening ceased, a region of fluid stasis formed near the left ventricular outflow tract. These findings suggest that LVAD support does not substantially alter vortex dynamics unless cardiac function is minimal. The altered blood flow introduced by the LVAD results in stasis adjacent to the LV outflow tract, which increases the risk of thrombus formation in the heart.     

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

BackgroundAs 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.MethodsIn 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.ResultsThe 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 2 mmHg is consistent with reported data.ConclusionWe 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.     

Wave intensity analysis of left atrial mechanics and energetics in anesthetized dogs

The left atrium (LA) acts as a booster pump during late diastole, generating the Doppler transmitral A wave and contributing incrementally to left ventricular (LV) filling. However, after volume loading and in certain disease states, LA contraction fills the LV less effectively, and retrograde flow (i.e., the Doppler Ar wave) into the pulmonary veins increases. The purpose of this study was to provide an energetic analysis of LA contraction to clarify the mechanisms responsible for changes in forward and backward flow. Wave intensity analysis was performed at the mitral valve and a pulmonary vein orifice. As operative LV stiffness increased with progressive volume loading, the reflection coefficient (i.e., energy of reflected wave/energy of incident wave) also increased. This reflected wave decelerated the forward movement of blood through the mitral valve and was transmitted through the LA, accelerating retrograde blood flow in the pulmonary veins. Although total LA work increased with volume loading, the forward hydraulic work decreased and backward hydraulic work increased. Thus wave reflection due to increased LV stiffness accounts for the decrease in the A wave and the increase in the Ar wave measured by Doppler.     

Computational study on the evolution of an intraventricular vortical flow during early diastole for the interpretation of color M-mode Doppler echocardiograms

A computational fluid dynamics study of intraventricular flow during early diastole was carried out using a 3D model of the human left ventricle (LV). It was found that a vortical flow formed under the aortic orifice and then grew in size and extended laterally along the ventricular wall towards the posterior side. With further expansion of the LV, it developed into an annular vortex asymmetrically enlarged on the side of the aortic orifice, narrowing the passage of blood inflow and thus causing a shift of the high-velocity portion of inflow towards the apex. This appeared as an elongation of the aliasing area when the velocity of the inflow was expressed as a spatiotemporal map in the same manner as a color M-mode Doppler (CMD) echocardiogram. Based on these findings, it was concluded that the shape of the aliasing area in a CMD echocardiogram shows the change in the velocity of blood inflow affected by the development of an annular vortex formed in the LV.     

Fluid–structure interaction in a fully coupled three-dimensional mitral–atrium–pulmonary model

This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.

     

Fluid dynamics of the left ventricular filling in dilated cardiomyopathy

Modifications in diastolic function occur in a broad range of cardiovascular diseases and there is an increasing evidence that abnormalities in left ventricular function may contribute significantly to the symptomatology. The flow inside the left ventricle during the diastole is here investigated by numerical solution of the Navier-Stokes equations under the axisymmetric assumption. The equation are written in a body-fitted, moving prolate spheroid, system of coordinates and solved using a fractional step method. The system is forced by a given volume time-law derived from clinical data, and varying the two-degrees-of-freedom ventricle geometry on the basis of a simple model. The solution under healthy conditions is analysed in terms of vorticity dynamics, showing that the flow field is characterised by the presence of a vortex wake; it is attached to the mitral valve during the accelerating phase of the E-wave, and it detaches and translate towards the ventricle apex afterwards. The flow evolution is discussed, results are also reported as an M-mode representation of colour-coded Doppler velocity maps. In the presence of ventricle dilatation the mitral jet extends farther inside the ventricle, propagation velocity decreases, and the fluid stagnates longer at the apex.     

Three-Directional Evaluation of Mitral Flow in the Rat Heart by Phase-Contrast Cardiovascular Magnetic Resonance

Purpose

Determination of mitral flow is an important aspect in assessment of cardiac function. Traditionally, mitral flow is measured by Doppler echocardiography which suffers from several challenges, particularly related to the direction and the spatial inhomogeneity of flow. These challenges are especially prominent in rodents. The purpose of this study was to establish a cardiovascular magnetic resonance (CMR) protocol for evaluation of three-directional mitral flow in a rodent model of cardiac disease.

Materials and Methods

Three-directional mitral flow were evaluated by phase contrast CMR (PC-CMR) in rats with aortic banding (AB) (N = 7) and sham-operated controls (N = 7). Peak mitral flow and deceleration rate from PC-CMR was compared to conventional Doppler echocardiography. The accuracy of PC-CMR was investigated by comparison of spatiotemporally integrated mitral flow with left ventricular stroke volume assessed by cine CMR.

Results

PC-CMR portrayed the spatial distribution of mitral flow and flow direction in the atrioventricular plane throughout diastole. Both PC-CMR and echocardiography demonstrated increased peak mitral flow velocity and higher deceleration rate in AB compared to sham. Comparison with cine CMR revealed that PC-CMR measured mitral flow with excellent accuracy. Echocardiography presented significantly lower values of flow compared to PC-CMR.

Conclusions

For the first time, we show that PC-CMR offers accurate evaluation of three-directional mitral blood flow in rodents. The method successfully detects alterations in the mitral flow pattern in response to cardiac disease and provides novel insight into the characteristics of mitral flow.     

Unsteady flow field around a human hand and propulsive force in swimming

Much effort has been undertaken for the estimation of propulsive force of swimmers in the front crawl. Estimation is typically based on steady flow theory: the so-called quasi-steady analysis. Flow fields around a swimmer, however, are extremely unsteady because the change direction of hand produces unsteady vortex motions. To evaluate the force correctly, it is necessary to know the unsteady properties determined from the vortex dynamics because that unsteadiness is known to make the force greater. Unsteady flow measurements were made for this study using a sophisticated technique called particle image velocimetry (PIV) in several horizontal planes for subjects swimming in a flume. Using that method, a 100 time-sequential flow fields are obtainable simultaneously. Each flow field was calculated from two particle images using the cross-correlation method. The intensity of vortices and their locations were identified. A strong vortex was generated near the hand and then shed by directional change of the hand in the transition phase from in-sweep to out-sweep. When the vortex was shed, a new vortex rotating in the opposite direction around the hand was created. The pair of vortices induced the velocity component in the direction opposite to the swimming. Results of this study show that the momentum change attributable to the increase in this velocity component is the origin of thrust force by the hand.     

RV instantaneous intraventricular diastolic pressure and velocity distributions in normal and volume overload awake dog disease models

Intraventricular diastolic right ventricular (RV) flow field dynamics were studied by functional imaging using three-dimensional (3D) real-time echocardiography with sonomicrometry and computational fluid dynamics in seven awake dogs at control with normal wall motion (NWM) and RV volume overload with diastolic paradoxical septal motion. Burgeoning flow cross section between inflow anulus and chamber walls induces a convective pressure rise, which represents a "convective deceleration load" (CDL). High spatiotemporal resolution dynamic pressure and velocity distributions of the intraventricular RV flow field revealed time-dependent, subtle interactions between intraventricular local acceleration and convective pressure gradients. During the E-wave upstroke, the total pressure gradient along intraventricular flow is the algebraic sum of a pressure decrease contributed by local acceleration and a pressure rise contributed by a convective deceleration that partially counterbalances the local acceleration gradient. This underlies the smallness of early diastolic intraventricular gradients. At peak volumetric inflow, local acceleration vanishes and the total adverse intraventricular gradient is convective. During the E-wave downstroke, the strongly adverse gradient embodies the streamwise pressure augmentations from both local and convective decelerations. It induces flow separation and large-scale vortical motions, stronger in NWM. Their dynamic corollaries on intraventricular pressure and velocity distributions were ascertained. In the NWM pattern, the strong ring-like vortex surrounding the central core encroaches on the area available for flow toward the apex. This results in higher linear velocities later in the downstroke of the E wave than at peak inflow rate. The augmentation of CDL by ventriculoannular disproportion may contribute to E wave and E-to-A ratio depression with chamber dilatation.