Modeling left ventricular dynamics with characteristic deformation modes

A computationally efficient method is described for simulating the dynamics of the left ventricle (LV) in three dimensions. LV motion is represented as a combination of a limited number of deformation modes, chosen to represent observed cardiac motions while conserving volume in the LV wall. The contribution of each mode to wall motion is determined by a corresponding time-dependent deformation variable. The principle of virtual work is applied to these deformation variables, yielding a system of ordinary differential equations for LV dynamics, including effects of muscle fiber orientations, active and passive stresses, and surface tractions. Passive stress is governed by a transversely isotropic elastic model. Active stress acts in the fiber direction and incorporates length–tension and force–velocity properties of cardiac muscle. Preload and afterload are represented by lumped vascular models. The variational equations and their numerical solutions are verified by comparison to analytic solutions of the strong form equations. Deformation modes are constructed using Fourier series with an arbitrary number of terms. Greater numbers of deformation modes increase deformable model resolution but at increased computational cost. Simulations of normal LV motion throughout the cardiac cycle are presented using models with 8, 23, or 46 deformation modes. Aggregate quantities that describe LV function vary little as the number of deformation modes is increased. Spatial distributions of stress and strain change as more deformation modes are included, but overall patterns are conserved. This approach yields three-dimensional simulations of the cardiac cycle on a clinically relevant time-scale.

     

Dependence of local left ventricular wall mechanics on myocardial fiber orientation: a model study.

The dependence of local left ventricular (LV) mechanics on myocardial muscle fiber orientation was investigated using a finite element model. In the model we have considered anisotropy of the active and passive components of myocardial tissue, dependence of active stress on time, strain and strain rate, activation sequence of the LV wall and aortic afterload. Muscle fiber orientation in the LV wall is quantified by the helix fiber angle, defined as the angle between the muscle fiber direction and the local circumferential direction. In a first simulation, a transmural variation of the helix fiber angle from +60 degrees at the endocardium through 0 degrees in the midwall layers to -60 degrees at the epicardium was assumed. In this simulation, at the equatorial level maximum active muscle fiber stress was found to vary from about 110 kPa in the subendocardial layers through about 30 kPa in the midwall layers to about 40 kPa in the subepicardial layers. Next, in a series of simulations, muscle fiber orientation was iteratively adapted until the spatial distribution of active muscle fiber stress was fairly homogeneous. Using a transmural course of the helix fiber angle of +60 degrees at the endocardium, +15 degrees in the midwall layers and -60 degrees at the epicardium, at the equatorial level maximum active muscle fiber stress varied from 52 kPa to 55 kPa, indicating a remarkable reduction of the stress range. Moreover, the change of muscle fiber strain with time was more similar in different parts of the LV wall than in the first simulation. It is concluded that (1) the distribution of active muscle fiber stress and muscle fiber strain across the LV wall is very sensitive to the transmural distribution of the helix fiber angle and (2) a physiological transmural distribution of the helix fiber angle can be found, at which active muscle fiber stress and muscle fiber strain are distributed approximately homogeneously across the LV wall.     

Model for left ventricular contraction combining the force length velocity relationship with the time varying elastance theory

A model for the contraction of the left ventricle (LV) is developed for a spheroidal geometry. The classical force-length-velocity relationship for a single muscle fiber is assumed. The linear maximum pressure volume relationship (maximum elastance), a measure of muscle contractility, is further extended into a time-varying function. This is achieved by utilizing a mechanical activation function, assumed as half a sinusoidal wave, to describe the time-dependent isometric stress for the activated cardiac muscle. This, in turn, results in the time-varying elastance function and represents the instantaneous activity of the muscle contractile proteins. The model is tested for a set of boundary conditions that determine preload, afterload, and the inherent properties of the muscle, i.e., the contractility. The computed results of the isovolumic contraction, auxotonic contraction, and isovolumic relaxation are in agreement with the expected behavior of the LV. The relations between the simulated variations on preload, afterload, and contractility, and the set of performance indexes of the LV, are presented and discussed.     

A model of the structural and functional development of the normal human fetal left ventricle based on a global growth law

The purpose of this research is to study the growth of the normal human left ventricle (LV) during the fetal period from 14 to 40 weeks of gestation. A new constitutive law for the active myocardium describing the mechanical properties of the active muscle during the whole cardiac cycle has been proposed. The LV model is a thick-walled, incompressible, hyperelastic cylinder, with families of helicoidal fibers running on cylindrical surfaces [1]. Based on the works of Lin and Taber [2] done on the embryonic chick heart, we use for the human fetal heart a growth law in which the growth rate depends on the wall stresses. The parameters of the growth law are adapted to agree with sizes and volumes inferred from two dimensional ultrasound measurements performed on 18 human fetuses.Then calculations are performed to extrapolate the cardiac performance during normal growth of the fetal LV. The results presented support the idea that a growth law in which the growth rate depends linearly on the mean wall stresses averaged through the space and during whole cardiac cycle, is adapted to the normal human fetal LV development.     

A Model of the Structural and Functional Development of the Normal Human Fetal Left Ventricle Based on a Global Growth Law

The purpose of this research is to study the growth of the normal human left ventricle (LV) during the fetal period from 14 to 40 weeks of gestation. A new constitutive law for the active myocardium describing the mechanical properties of the active muscle during the whole cardiac cycle has been proposed. The LV model is a thick-walled, incompressible, hyperelastic cylinder, with families of helicoidal fibers running on cylindrical surfaces [1] . Based on the works of Lin and Taber [2] done on the embryonic chick heart, we use for the human fetal heart a growth law in which the growth rate depends on the wall stresses. The parameters of the growth law are adapted to agree with sizes and volumes inferred from two dimensional ultrasound measurements performed on 18 human fetuses. Then calculations are performed to extrapolate the cardiac performance during normal growth of the fetal LV. The results presented support the idea that a growth law in which the growth rate depends linearly on the mean wall stresses averaged through the space and during whole cardiac cycle, is adapted to the normal human fetal LV development.     

Theoretical approach to blood ejection from the human left ventricle

An ejection dynamics mathematical model of human left ventricle (LV) based on physiological data of human heart is proposed in this study. The mathematical equations were expressed in terms of vorticity-stream function equations in a prolate spheroidal coordinate system. These equations combined with specified boundary conditions were numerically solved by using an alternating-direction-implicit (ADI) algorithm with second order accuracy. The unsteady aspects of the ejection process were subsequently introduced into the numerical simulation. The numerical results have shown that the present ellipsoidal model could be available to simulate the ejection process of the human LV. Such a model combined with cardiac muscle mechanics could be studied further to determine altered left ventricular function in cardiac diseases.     

A model study of intracellular oxygen gradients in a myoglobin-containing skeletal muscle fiber.

A theoretical two-dimensional model is used to investigate oxygen gradients in a red skeletal muscle fiber. The model describes the steady state, free and myoglobin-facilitated diffusion of oxygen into a respiring cylindrical muscle fiber cross section. The oxygen tension at the sarcolemma is assumed to vary along the sarcolemma as an approximation to the discrete capillary oxygen supply around the fiber. Maximal oxygen gradients are studied by considering parameters relevant to a maximally-respiring red muscle fiber. The model predicts that angular variations in the oxygen tension imposed at the sarcolemma due to the discrete capillary sources do not penetrate deeply into the fiber over a range of physiological values for myoglobin concentration, diffusion coefficients, number of surrounding capillaries, and oxygen tension level at the sarcolemma. Also, the oxygen tension in the core of the fiber is determined by the average oxygen tension at the sarcolemma. The drop in oxygen tension from fiber periphery to core, however, does depend significantly on the myoglobin concentration, the oxygen tension level at the sarcolemma, and the oxygen and myoglobin diffusivities. This dependence is summarized by calculating the minimum average sarcolemmal oxygen tension for maximal respiration without the development of an intracellular anoxic region. For a myoglobin-rich muscle fiber (0.5 mM myoglobin), the model predicts that maximal oxygen consumption can proceed with a relatively flat (less than 5 mm Hg) oxygen tension drop from fiber periphery to core over a large range for diffusion coefficients.     

Age-related changes in cardiac structure and function in Fischer 344 x Brown Norway hybrid rats

The effects of aging on cardiovascular function and cardiac structure were determined in a rat model recommended for gerontological studies. A cross-sectional analysis assessed cardiac changes in male Fischer 344 x Brown Norway F1 hybrid rats (FBN) from adulthood to the very aged (n = 6 per 12-, 18-, 21-, 24-, 27-, 30-, 33-, 36-, and 39-mo-old group). Rats underwent echocardiographic and hemodynamic analyses to determine standard values for left ventricular (LV) mass, LV wall thickness, LV chamber diameter, heart rate, LV fractional shortening, mitral inflow velocity, LV relaxation time, and aortic/LV pressures. Histological analyses were used to assess LV fibrotic infiltration and cardiomyocyte volume density over time. Aged rats had an increased LV mass-to-body weight ratio and deteriorated systolic function. LV systolic pressure declined with age. Histological analysis demonstrated a gradual increase in fibrosis and a decrease in cardiomyocyte volume density with age. We conclude that, although significant physiological and morphological changes occurred in heart function and structure between 12 and 39 mo of age, these changes did not likely contribute to mortality. We report reference values for cardiac function and structure in adult FBN male rats through very old age at 3-mo intervals.     

Dynamic myocardial contractile parameters from left ventricular pressure-volume measurements

A new dynamic model of left ventricular (LV) pressure-volume relationships in beating heart was developed by mathematically linking chamber pressure-volume dynamics with cardiac muscle force-length dynamics. The dynamic LV model accounted for >80% of the measured variation in pressure caused by small-amplitude volume perturbation in an otherwise isovolumically beating, isolated rat heart. The dynamic LV model produced good fits to pressure responses to volume perturbations, but there existed some systematic features in the residual errors of the fits. The issue was whether these residual errors would be damaging to an application where the dynamic LV model was used with LV pressure and volume measurements to estimate myocardial contractile parameters. Good agreement among myocardial parameters responsible for response magnitude was found between those derived by geometric transformations of parameters of the dynamic LV model estimated in beating heart and those found by direct measurement in constantly activated, isolated muscle fibers. Good agreement was also found among myocardial kinetic parameters estimated in each of the two preparations. Thus the small systematic residual errors from fitting the LV model to the dynamic pressure-volume measurements do not interfere with use of the dynamic LV model to estimate contractile parameters of myocardium. Dynamic contractile behavior of cardiac muscle can now be obtained from a beating heart by judicious application of the dynamic LV model to information-rich pressure and volume signals. This provides for the first time a bridge between the dynamics of cardiac muscle function and the dynamics of heart function and allows a beating heart to be used in studies where the relevance of myofilament contractile behavior to cardiovascular system function may be investigated.     

Low-intensity interval exercise training attenuates coronary vascular dysfunction and preserves Ca²⁺-sensitive K⁺ current in miniature swine with LV hypertrophy

Coronary vascular dysfunction has been observed in several models of heart failure (HF). Recent evidence indicates that exercise training is beneficial for patients with HF, but the precise intensity and underlying mechanisms are unknown. Left ventricular (LV) hypertrophy can play a significant role in the development of HF; therefore, the purpose of this study was to assess the effects of low-intensity interval exercise training on coronary vascular function in sedentary (HF) and exercise trained (HF-TR) aortic-banded miniature swine displaying LV hypertrophy. Six months postsurgery, in vivo coronary vascular responses to endothelin-1 (ET-1) and adenosine were measured in the left anterior descending coronary artery. Baseline and maximal coronary vascular conductance were similar between all groups. ET-1-induced reductions in coronary vascular conductance (P < 0.05) were greater in HF vs. sedentary control and HF-TR groups. Pretreatment with the ET type A (ET(A)) receptor blocker BQ-123 prevented ET-1 hypersensitivity in HF animals. Whole cell voltage clamp was used to characterize composite K(+) currents (I(K(+))) in coronary smooth muscle cells. Raising internal Ca(2+) from 200 to 500 nM increased Ca(2+)-sensitive K(+) current in HF-TR and control, but not HF animals. In conclusion, an ET(A)-receptor-mediated hypersensitivity to ET-1, elevated resting LV wall tension, and decreased coronary smooth muscle cell Ca(2+)-sensitive I(K(+)) was found in sedentary animals with LV hypertrophy. Low-intensity interval exercise training preserved normal coronary vascular function and smooth muscle cell Ca(2+)-sensitive I(K(+)), illustrating a potential mechanism underlying coronary vascular dysfunction in a large-animal model of LV hypertrophy. Our results demonstrate the potential clinical impact of exercise on coronary vascular function in HF patients displaying pathological LV hypertrophy.     

The effects of ADP and phosphate on the contraction of muscle fibers.

The products of MgATP hydrolysis bind to the nucleotide site of myosin and thus may be expected to inhibit the contraction of muscle fibers. We measured the effects of phosphate and MgADP on the isometric tensions and isotonic contraction velocities of glycerinated rabbit psoas muscle at 10 degrees C. Addition of phosphate decreased isometric force but did not affect the maximum velocity of shortening. To characterize the effects of ADP on fiber contractions, force-velocity curves were measured for fibers bathed in media containing various concentrations of MgATP (1.5-4 mM) and various concentrations of MgADP (1-4 mM). As the [MgADP]/[MgATP] ratio in the fiber increases, the maximum velocity achieved by the fiber decreases while the isometric tension increases. The inhibition of fiber velocities and the potentiation of fiber tension by MgADP is not altered by the presence of 12 mM phosphate. The concentration of both MgADP and MgATP within the fiber was calculated from the diffusion coefficient for nucleotides within the fiber, and the rate of MgADP production within the fiber. Using the calculated values for the nucleotide concentration inside the fiber, observed values of the maximum contraction velocity could be described, within experimental accuracy, by a model in which MgADP competed with MgATP and inhibited fiber velocity with an effective Ki of 0.2-0.3 mM. The average MgADP level generated by the fiber ATPase activity within the fiber was approximately 0.9 mM. In fatigued fibers MgADP and phosphate levels are known to be elevated, and tension and the maximum velocity of contraction are depressed. The results obtained here suggest that levels of MgADP in fatigued fibers play no role in these decreases in function, but the elevation of both phosphate and H+ is sufficient to account for much of the decrease in tension.     

Theoretical models in mechanics of the left ventricle

Different rheological concepts and theoretical studies have been recently presented using models of myocardial mechanics. Complex analysis of the mechanical behavior of the left ventricular wall have been developed in order to estimate the local stresses and deformations that occur during the heart cycle as well as the ventricular stroke volume and pressure. Theoretical models have taken into account non-linear and viscoelastic passive properties of the myocardium tissue, when subjected to large deformations, through given strain energy functions or stress-strain relations. Different prolate spheroid geometries have been considered for such thick shell cardiac structure. During the active state of the contraction, the rheological behavior of the fibers has been described using different muscle models and relationships between fiber tension and strain, and activation degree. A forthcoming approach for bridging the gap between the knowledge of the muscle fiber microrheological properties and the study of the mechanical behavior of the entire ventricle, consists in including anisotropic and inhomogeneous effects through fiber direction field.     

Effects of Long Term Supplementation of Anabolic Androgen Steroids on Human Skeletal Muscle

The effects of long-term (over several years) anabolic androgen steroids (AAS) administration on human skeletal muscle are still unclear. In this study, seventeen strength training athletes were recruited and individually interviewed regarding self-administration of banned substances. Ten subjects admitted having taken AAS or AAS derivatives for the past 5 to 15 years (Doped) and the dosage and type of banned substances were recorded. The remaining seven subjects testified to having never used any banned substances (Clean). For all subjects, maximal muscle strength and body composition were tested, and biopsies from the vastus lateralis muscle were obtained. Using histochemistry and immunohistochemistry (IHC), muscle biopsies were evaluated for morphology including fiber type composition, fiber size, capillary variables and myonuclei. Compared with the Clean athletes, the Doped athletes had significantly higher lean leg mass, capillary per fibre and myonuclei per fiber. In contrast, the Doped athletes had significantly lower absolute value in maximal squat force and relative values in maximal squat force (relative to lean body mass, to lean leg mass and to muscle fiber area). Using multivariate statistics, an orthogonal projection of latent structure discriminant analysis (OPLS-DA) model was established, in which the maximal squat force relative to muscle mass and the maximal squat force relative to fiber area, together with capillary density and nuclei density were the most important variables for separating Doped from the Clean athletes (regression  =  0.93 and prediction  =  0.92, p<0.0001). In Doped athletes, AAS dose-dependent increases were observed in lean body mass, muscle fiber area, capillary density and myonuclei density. In conclusion, long term AAS supplementation led to increases in lean leg mass, muscle fiber size and a parallel improvement in muscle strength, and all were dose-dependent. Administration of AAS may induce sustained morphological changes in human skeletal muscle, leading to physical performance enhancement.     

Stress Distribution in the Canine Left Ventricle during Diastole and Systole

A model is proposed for stress analysis of the left ventricular wall (LV wall) based on the realistic assumption that the myocardium is essentially composed of fiber elements which carry only axial tension and vary in orientation through the wall. Stress analysis based on such a model requires an extensive study of muscle fiber orientation and curvature through the myocardium. Accordingly, the principal curvatures were studied at a local site near the equator in ten dog hearts rapidly fixed in situ at end diastole and end systole; the fiber orientation for these hearts had already been established in a previous study. The principal radii of curvature were (a) measured by fitting templates to the endocardial and epicardial wall surfaces in the circumferential and longitudinal directions and (b) computed from measured lengths of semiaxes of ellipsoids of revolution representing the LV wall (“ellipsoid” data). The wall was regarded as a tethered set of nested shells, each having a unique fiber orientation. Results indicate the following. (a) Fiber curvature, k, is maximum at midwall at end systole; this peak shifts towards endocardium at end diastole. (b) The pressure or radial stress through the wall decreases more rapidly near the endocardium than near the epicardium at end diastole and at end systole when a constant tension is assumed for each fiber through the wall. (c) At end diastole the curve for the circumferential stress vs. wall thickness is convex with a maximum at midwall. In the longitudinal direction the stress distribution curve is concave with a minimum at midwall. Similar distributions are obtained at end systole when a constant tension is assumed for each fiber through the wall. (d) The curvature and stress distributions obtained by direct measurements at a selected local site agree well with those computed from “ellipsoid” data.     

Viscerosomatic interaction induced by myocardial ischemia in conscious dogs.

Studies tested the hypothesis that myocardial ischemia induces increased paraspinal muscular tone localized to the T(2)-T(5) region that can be detected by palpatory means. This is consistent with theories of manual medicine suggesting that disturbances in visceral organ physiology can cause increases in skeletal muscle tone in specific muscle groups. Clinical studies in manual and traditional medicine suggest this phenomenon occurs during episodes of myocardial ischemia and may have diagnostic potential. However, there is little direct evidence of a cardiac-somatic mechanism to explain these findings. Chronically instrumented dogs [12 neurally intact and 3 following selective left ventricular (LV) sympathectomy] were examined before, during, and after myocardial ischemia. Circumflex blood flow (CBF), left ventricular contractile function, electromyographic (EMG) analysis, and blinded manual palpatory assessments (MPA) of tissue over the transverse spinal processes at segments T(2)-T(5) and T(11)-T(12) (control) were performed. Myocardial ischemia was associated with a decrease in myocardial contractile function and an increase in heart rate. MPA revealed increases in muscle tension and texture/firmness during ischemia in the T(2)-T(5) segments on the left, but not on the right or in control segments. EMG demonstrated increased amplitude for the T(4)-T(5) segments. After LV sympathectomy, MPA and EMG evidence of increased muscle tone were absent. In conclusion, myocardial ischemia is associated with significant increased paraspinal muscle tone localized to the left side T(4)-T(5) myotomes in neurally intact dogs. LV sympathectomy eliminates the somatic response, suggesting that sympathetic neural traffic between the heart and somatic musculature may function as the mechanism for the interaction.     

Interaction between myosin heavy chain and troponin isoforms modulate cardiac myofiber contractile dynamics

Coordinated expression of species-specific myosin heavy chain (MHC) and troponin (Tn) isoforms may bring about a dynamic complementarity to match muscle contraction speed with species-specific heart rates. Contractile system function and dynamic force-length measurements were made in muscle fibers from mouse and rat hearts and in muscle fibers after reconstitution with either recombinant homologous Tn or orthologous Tn. The rate constants of length-mediated cross-bridge (XB) recruitment (b) and tension redevelopment (k(tr)) of mouse fibers were significantly faster than those of rat fibers. Both the tension cost (ATPase/tension) and rate constant of length-mediated XB distortion (c) were higher in the mouse than in the rat. Thus the mouse fiber was faster in all dynamic and functional aspects than the rat fiber. Mouse Tn significantly increased b and k(tr) in rat fibers; conversely, rat Tn significantly decreased b and k(tr) in mouse fibers. Thus the length-mediated recruitment of force-bearing XB occurs much more rapidly in the presence of mouse Tn than in the presence of rat Tn, demonstrating that the speed of XB recruitment is regulated by Tn. There was a significant interaction between Tn and MHC such that changes in either Tn or MHC affected the speed of XB recruitment. Our data demonstrate that the dynamics of myocardial contraction are different in the mouse and rat hearts because of sequence heterogeneity in MHC and Tn. At the myofilament level, coordinated expression of complementary regulatory contractile proteins produces a functional dynamic phenotype that allows the cardiovascular systems to function effectively at different heart rates.     

Optimizing ventricular fibers: uniform strain or stress,but not ATP consumption,leads to high efficiency

The aim of this study was to investigate the influence of fiber orientation in the left ventricular (LV) wall on the ejection fraction, efficiency, and heterogeneity of the distributions of developed fiber stress, strain and ATP consumption. A finite element model of LV mechanics was used with active properties of the cardiac muscle described by the Huxley-type cross-bridge model. The computed variances of sarcomere length (SL(var)), developed stress (DS(var)), and ATP consumption (ATP(var)) have several minima at different transmural courses of helix fiber angle. We identified only one region in the used design space with high ejection fraction, high efficiency of the LV and relatively small SL(var), DS(var), and ATP(var). This region corresponds to the physiological distribution of the helix fiber angle in the LV wall. Transmural fiber angle can be predicted by minimizing SL(var) and DS(var), but not ATP(var). If ATP(var) was minimized, then the transverse fiber angle was considerably underestimated. The results suggest that ATP consumption distribution is not regulating the fiber orientation in the heart.     

Modeling and simulation of musculoskeletal system of human lower limb based on tensegrity structure

Abstract

In this paper, a mechanical model of the skeletal muscle of human lower limb system is established by using the Hill muscle model and kinetic equation of the movement of lower extremities according to the attachment positions of skeletal muscle. State vector and neural control are delineated by the direct configuration method. Changes of gait and skeletal muscle stress during walking process are analyzed with energy consumption as objective function. Results illustrate that simulation data are in good agreement with actual walking gait data. Feasibility and correctness of the designed model and control behavior of skeletal muscle tension structure are also verified.     

A mathematical model of force transmission from intrafascicularly terminating muscle fibers

Many long skeletal muscles are comprised of fibers that terminate intrafascicularly. Force from terminating fibers can be transmitted through shear within the endomysium that surrounds fibers or through tension within the endomysium that extends from fibers to the tendon; however, it is unclear which pathway dominates in force transmission from terminating fibers. The purpose of this work was to develop mathematical models to (i) compare the efficacy of lateral (through shear) and longitudinal (through tension) force transmission in intrafascicularly terminating fibers, and (ii) determine how force transmission is affected by variations in the structure and properties of fibers and the endomysium. The models demonstrated that even though the amount of force that can be transmitted from an intrafascicularly terminating fiber is dependent on fiber resting length (the unstretched length at which passive stress is zero), endomysium shear modulus, and fiber volume fraction (the fraction of the muscle cross-sectional area that is occupied by fibers), fibers that have values of resting length, shear modulus, and volume fraction within physiologic ranges can transmit nearly all of their peak isometric force laterally through shearing of the endomysium. By contrast, the models predicted only limited force transmission ability through tension within the endomysium that extends from the fiber to the tendon. Moreover, when fiber volume fraction decreases to unhealthy ranges (less than 50%), the force-transmitting potential of terminating fibers through shearing of the endomysium decreases significantly. The models presented here support the hypothesis that lateral force transmission through shearing of the endomysium is an effective mode of force transmission in terminating fibers.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.