The effects of supravalvular aortic stenosis on realistic three-dimensional left ventricular blood ejection

Intraventricular pressure drop and aortic blood acceleration as indices of cardiac inotropy: a comparison with the first derivative of aortic pressure based on computer fluid dynamics

This paper presents a computational approach to ventricular fluid mechanics to evaluate three inotropic indices of early ejection: the intraventricular pressure drop (deltap). the first derivative of aortic flow rate (df/dt) and the first derivative of aortic pressure dp/dt. dp/dt is one of the most frequently used indices for assessing myocardial inotropy. Deltap and df/dt are characteristic of inertia driven flows and reflect the impulsive nature of the flow inside the ventricle during the ejection phase. The study is based on an axisymmetric fluid dynamics model of the left ventricle, developed according to the finite element approach. The fluid cavity is bounded by a shell containing two sets of counter-rotating contractile fibres. Two simulation sets were performed: the former to investigate the sensitivity of deltap and df/dt peaks (deltap(max) and df/dt(max)) with respect to changes in the inotropic state of the fibre. The latter allows the evaluation of the dependency of deltap(max) and df/dt(max) on afterload by means of two supravalvular stenoses of 50% and 70%. The model simulates the inertial features of ventricle behaviour. The calculated values of the indices investigated are in close agreement with those reported in the literature. The sensitivities of deltap(max) df/dt(max) and dp/dt(max) are calculated for the two simulation sets. Data are normalised with respect to the maximum values reached in the simulation set. The comparison indicates that deltap(max) has a greater sensitivity (3.4 vs. 3.1 ) and a more linear pattern than dp/dt(max) for changes in the inotropic state of the fibre. df/dt(max), shows a sensitivity close to dp/dt(max). Results confirm that the afterload does not affect dp/dt(max), in accordance with experimental observations, while deltap(max) and, to a major degree, df/dt(max) decrease when the afterload is increased.

Computational evaluation of intraventricular pressure gradients based on a fluid-structure approach

The dynamics of intraventricular blood flow, i.e. its rapid evolution, implies the rise of intraventricular pressure gradients (IPGs) characteristic of the inertia-driven events as experimentally observed by Pasipoularides (1987, 1990) and by Falsetti et al. (1986). The IPG time course is determined by the wall contraction which, in turn, depends on the load applied, namely the intraventricular pressure which is the sum of the aortic pressure (i.e., the systemic net response) and the IPG. Hence the IPGs account, at least in part, for the wall movement. These considerations suggest the necessity of a comprehensive analysis of the ventricular mechanics involving both ventricular wall mechanics and intraventricular fluid dynamics as each domain determines the boundary conditions of the other. This paper presents a computational approach to ventricular ejection mechanics based on a fluid-structure interaction calculation for the evaluation of the IPG time course. An axisymmetric model of the left ventricle is utilized. The intraventricular fluid is assumed to be Newtonian. The ventricle wall is thin and is composed of two sets of counter-rotating fibres which behave according to the modified version of Wong's sarcomere model proposed by Montevecchi and Pietrabissa and Pietrabissa et al. (1987, 1991). The full Navier-Stokes equations describing the fluid domain are solved using Galerkin's weighted residual approach in conjunction with finite element approximation (FIDAP). The wall displacement is solved using the multiplane quasi-Newton method proposed by Buzzi Ferraris and Tronconi (1985). The interaction procedure is performed by means of an external macro which compares the flow fields and the wall displacement and appropriately modifies the boundary conditions to reach the simultaneous and congruous convergence of the two problems. The results refer to a simulation of the ventricular ejection with a heart rate of 72 bpm. In this phase the ventricle ejects 61 cm3 (ejection fraction equal to 54 percent) and the ventricular pressure varies from 78 mmHg to 140 mmHg. The IPG show an oscillating behaviour with two major peaks at the beginning (11.09 mmHg) and at the end (4.32 mmHg) of the ejection phase, when the flow rate hardly changes, according to the experimental data. Furthermore the wall displacement, the wall stress and strain, the pressure and velocity fields are calculated and reported.

Realistic three-dimensional left ventricular ejection determined from computational fluid dynamics

A realistic model of the left ventricle of the human heart was constructed using a cast from a dog heart which was in diastole. A coordinate measuring machine was used to measure and digitize the coordinates of the left ventricle. From the complex measured left ventricle shape values, a three-dimensional finite volume representation was constructed using a simulation package. The left ventricular walls moved towards the centre of the aortic outlet in order to study the effects of time-varying left ventricular ejection. The left ventricular wall motion was assumed to follow the blood flow and the wall grid was reformed 25 times during the calculation. The 25.8 cm3 ventricular volume was reduced by 75% in 0.25 s. Centreline and cross-sectional velocity vectors greatly increased in magnitude at the aortic outlet, and most of the pressure occurred in the top 15% of the heart. The computational method should make it possible to compare simulation results with important measurement techniques such as ultrasound and magnetic resonance imaging, and this should allow a finer detail of flow understanding than is presently available using either a modelling or imaging method alone.