An LDA and flow visualization study of pulsatile flow in an aortic bifurcation model

The structure of pulsatile flow in a rigid aortic bifurcation model was studied by means of a flow visualization technique and three-dimensional laser-Doppler anemometry. The model was made of glass, having the same shape as that of the average human aortic bifurcation. It was installed into a mock circulatory loop that generated physiological pulsatile flow. Flow separation was observed during accelerated and decelerated flow periods. Double helical flow existed inside the flow separation in the early accelerated flow period. In the decelerated flow period, disturbed flow appeared behind the separation zone. Flow was strongly disturbed during the back flow period, and then was gradually stabilized in the forward flow period. The flow separation and the disturbances released from the flow separation zone greatly influenced near-wall velocities along the lateral wall. The wave form of the near-wall velocity in the flow separation zone was much different from that observed in the aortic portion and behind the separation zone; for example, the magnitude of the negative peak velocity in the direction of the tube axis was larger than that of the positive one, and mean velocity in a cycle was very low. This abnormal phasic change of the near-wall velocity may be associated with atherogenesis. The three-dimensional velocity measurement is very useful for the detailed analysis of near-wall velocity patterns.     

Turbulence in the canine ascending aorta and the blood pressure.

Turbulent velocity fluctuations were measured and analyzed in the canine ascending aorta using a hot-film anemometer. Blood flow rate and temperature were stabilized using a special bypass technique. Blood pressure was elevated by Methoxamine infusion. Turbulence components were extracted from measured data using an ensemble averaging technique. Turbulence intensity correlated best with blood flow rate although the variance was relatively large, especially when the blood flow velocity was high. When pooled data were grouped into subclasses using peak aortic flow velocity as the criteria, turbulence intensity correlated well with aortic systolic blood pressure in each of the subclasses. Spectral bandwidth correlated with aortic pressure in the same manner. In summary, turbulence in the aorta developed when blood pressure was high. Both an increase of turbulence intensity and an widening of turbulence spectra may be ascribed to a stiffening of the aortic wall due to an elevation of blood pressure.     

Instantaneous blood flow, impedance and elastic properties computed in man from aortic pulse waves

A mathematical model of the pressure-flow relationship in the arterial circulation and its possible use in routine hemodynamics in man are described. The instantaneous blood flow velocity in the ascending aorta can be calculated from two pressure curves simultaneously recorded 5 cm apart. The mechanical aortic input impedance is computed from the recorded pressure and the calculated blood flow velocity curves. Projection of the pulse waves on a time-length plane leads to the determination of the pulse wave velocity and then an estimation of the elastic modulus of the aortic wall.     

Continuous measurement of aortic blood velocity,after cardiac surgery,by means of an extractable doppler ultrasound probe

A new method has been developed for the continuous measurement of aortic blood velocity in patients following cardiac surgery. Using an extractable Doppler ultrasound probe placed on the ascending aorta, the changes in aortic velocity were recorded up to 24 h postoperatively, in 14 patients undergoing coronary bypass surgery. Volume flow rate is calculated from the mean velocity, the diameter of the aorta and the angle between the ultrasound beam and the direction of the blood flow, by means of an analogue flow calculator. Estimation of aortic flow showed a correlation of r = 0.79 with cardiac output measured by a thermodilution technique. The main advantage of the system is that it allows continuous monitoring of cardiac output, as well as short and long-term trend analyses, during the early postoperative period.     

Spatial velocity profile in mouse embryonic aorta and Doppler-derived volumetric flow: a preliminary model

Characterizing embryonic circulatory physiology requires accurate cardiac output and flow data. Despite recent applications of high-frequency ultrasound Doppler to the study of embryonic circulation, current Doppler analysis of volumetric flow is relatively crude. To improve Doppler derivation of volumetric flow, we sought a preliminary model of the spatial velocity profile in the mouse embryonic dorsal aorta using ultrasound biomicroscopy (UBM)-Doppler data. Embryonic hematocrit is 0.05-0.10 so rheologic properties must be insignificant. Low Reynolds numbers (<500) and Womersley parameters (<0.76) suggest laminar flow. UBM demonstrated a circular dorsal aortic cross section with no significant tapering. Low Dean numbers (<100) suggest the presence of minimal skewing of the spatial velocity profile. The inlet length allows for fully developed flow. There is no apparent aortic wall pulsatility. Extrapolation of prior studies to these vessel diameters (300-350 microm) and flow velocities (~50-200 mm/s) suggests parabolic spatial velocity profiles. Therefore, mouse embryonic dorsal aortic blood flow may correspond to Poiseuille flow in a straight rigid tube with parabolic spatial velocity profiles. As a first approximation, these results are an important step toward precise in utero ultrasound characterization of blood flow within the developing mammalian circulation.     

Velocity distribution along an elastic model of human arterial tree

An experimental investigation of an elastic model of the human arterial tree, has been performed for physiological type flow by pulsed Doppler ultrasonic velocimetry. The arterial tree model, fabricated in clear polyurethane, includes the aortic arch, with a Starr-Edwards ball valve mounted in the root of the aorta, the descending aorta and the iliac bifurcation. Our study showed that the velocity profile, a few centimeters beyond the valve, is skewed, with higher velocities towards the top and the inner wall (anatomically the posterior and left lateral wall). An inward shift of the maximum velocity and reverse flow are denoted along the inner wall of the aortic arch. The velocity profiles in the descending aorta are blunted. Downstream from the vertex of the iliac bifurcation, there is vorticity creation, but the branching effect is quickly damped by the pulsatility of the flow and the elasticity of the wall.     

Computational analysis of blood flow in an integrated model of the left ventricle and the aorta

To study the effects of intraventricular flow dynamics on the aortic flow, we created an integrated model of the left ventricle and aorta and conducted a computer simulation of diastolic and systolic blood flow within this model. The results demonstrated that the velocity profile at the aortic annulus changed dynamically, and was influenced by the intraventricular flow dynamics. The profile was almost flat in early systole but became nonuniform as systole progressed, and was skewed toward the posterior side in midsystole and toward the anterior side in later systole. At a distance from the aortic annulus, a different velocity profile was induced by the twisting and torsion of the aorta. In the ascending aorta, the fastest flow was initially located in the posteromedial sector, and it moved to the posterior section along the circumference as systole progressed. The nonuniformity of the aortic inflow gave rise to a complex wall shear stress (WSS) distribution in the aorta. A comparison of the WSS distribution obtained in this integrated analysis with that obtained in flow calculations using an isolated aorta model with Poiseuille and flat inlet conditions showed that intraventricular flow affected the WSS distribution in the ascending aorta. These results address the importance of an integrated analysis of flow in the left ventricle and aorta.     

A new approach for the evaluation of the severity of coarctation of the aorta using Doppler velocity index and effective orifice area: in vitro validation and clinical implications

Early detection and accurate estimation of COA severity are the most important predictors of successful long-term outcome. However, current clinical parameters used for the evaluation of the severity of COA have several limitations and are flow dependent. The objectives of this study are to evaluate the limitations of current existing parameters for the evaluation of the severity of coarctation of the aorta (COA) and suggest two new parameters: COA Doppler velocity index and COA effective orifice area. Three different severities of COAs were tested in a mock flow circulation model under various flow conditions and in the presence of normal and stenotic aortic valves. Catheter trans-COA pressure gradients and Doppler echocardiographic trans-COA pressure gradients were evaluated. COA Doppler velocity index was defined as the ratio of pre-COA to post-COA peak velocities measured by Doppler echocardiography. COA Doppler effective orifice area was determined using continuity equation. The results show that peak-to-peak trans-COA pressure gradient significantly increased with flow rate (from 83% to 85%). Peak Doppler pressure gradient also significantly increased with flow rate (80-85%). A stenotic or bicuspid aortic valve increased peak Doppler pressure gradient by 20-50% for a COA severity of 75%. Both COA Doppler velocity index and COA effective orifice area did not demonstrate significant flow dependence or dependence upon aortic valve condition. As a conclusion, COA Doppler velocity index and COA effective orifice area are flow independent and do not depend on aortic valve conditions. They can, then, more accurately predict the severity of COA.     

Experimental investigation of steady flow through a model of the human aortic arch

Steady flow through a model of the human aortic arch has been studied with hot-film anemometry. A three sensor hot-film velocity probe was inserted into an acrylic flow chamber fabricated from the in situ casting of a human aorta, and the axial, radial and tangential velocity profiles were determined for steady flows in the region of the aortic arch. These studies demonstrated the presence of a potential core throughout the arch region, with a concomitant boundary layer adjacent to the inner wall of curvature of the arch. Trapped secondary flows in this fluid layer along the inner wall were quantitatively determined. Our steady flow studies in the model human aortic arch suggests that a shear-dependent mass transfer mechanism may play a significant role in the development and propagation of atherosclerotic lesions in this segment of the human cardiovascular system.     

Kinematic and Dynamic Characteristics of Pulsating Flow in 180^o Tube

Kinematic and dynamic characteristics of pulsating flow in a model of human aortic arch are obtained by a computational analysis. Three-dimensional flow processes are summarized by pressure distributions on the symmetric plane together with velocity and pressure contours on a few cross sections for systolic acceleration and deceleration. Without considering the effects of aortic tapering and the carotid arteries, the development of tubular boundary layer with centrifugal forces and pulsation are also analyzed for flow separation and backflow during systolic deceleration.     

Numerical simulations of steady flow inside a three dimensional aortic bifurcation model

A finite element formulation of the Navier-Stokes equations for three dimensional flow is presented. The equations are solved using the finite element method. The model is constructed from a cast of a human aortic bifurcation. The numerical problems introduced by solving the equation system are discussed and special attention is paid to the selection of the linear equation solver. The simulations of the steady blood flow patterns in an aortic bifurcation is shown. The results of the numerical analysis are presented as three dimensional plots of velocity vectors, wall shear vectors, streamlines and pressure isobars. The flow simulations are done for Reynolds number 10. The flow patterns found in the bifurcation model are discussed in connection with proposed theories to explain the event of early atherosclerosis.     

稳恒流下腹主动脉旁瘤超声多普勒血流信号的仿真分析

腹主动脉旁瘤超声多普勒血流信号的仿真研究,可以为采用超声多普勒技术检测腹主动脉旁瘤的形成、生长过程和估计动脉旁瘤瘤体大小提供指导。先通过有限元数值计算方法得到稳恒流下腹主动脉旁瘤区域内的血液流场分布,然后采用余弦叠加的合成方法仿真出相应的超声多普勒血流信号,最后对仿真信号进行频谱分析,计算其平均频率,并研究它与腹主动脉旁瘤瘤体大小的关系。结果表明：当动脉旁瘤较小时,平均频率的幅度变化较小;当动脉旁瘤较大时,平均频率的幅度变化较大。因此,采用平均频率的幅度变化可以在一定程度上估计动脉旁瘤的瘤体大小。     

Estimation of valvular resistance of segmented aortic valves using computational fluid dynamics

Aortic valve stenosis is associated with an elevated left ventricular pressure and transaortic pressure drop. Clinicians routinely use Doppler ultrasound to quantify aortic valve stenosis severity by estimating this pressure drop from blood velocity. However, this method approximates the peak pressure drop, and is unable to quantify the partial pressure recovery distal to the valve. As pressure drops are flow dependent, it remains difficult to assess the true significance of a stenosis for low-flow low-gradient patients. Recent advances in segmentation techniques enable patient-specific Computational Fluid Dynamics (CFD) simulations of flow through the aortic valve. In this work a simulation framework is presented and used to analyze data of 18 patients. The ventricle and valve are reconstructed from 4D Computed Tomography imaging data. Ventricular motion is extracted from the medical images and used to model ventricular contraction and corresponding blood flow through the valve. Simplifications of the framework are assessed by introducing two simplified CFD models: a truncated time-dependent and a steady-state model. Model simplifications are justified for cases where the simulated pressure drop is above 10 mmHg. Furthermore, we propose a valve resistance index to quantify stenosis severity from simulation results. This index is compared to established metrics for clinical decision making, i.e. blood velocity and valve area. It is found that velocity measurements alone do not adequately reflect stenosis severity. This work demonstrates that combining 4D imaging data and CFD has the potential to provide a physiologically relevant diagnostic metric to quantify aortic valve stenosis severity.     

A model investigation of the velocity and pressure spectra in vascular murmurs

A physical model consisting of an axisymmetrical jet in a rigid plexiglass pipe was used to study the flow and pressure fluctuations downstream from an aortic stenosis. The fluctuating velocity components, u and v, at several locations in the steady liquid jet were measured using a laser Doppler anemometer system. Simultaneous wall pressure fluctuations were monitored by an array of nine miniature pressure transducers wall mounted in the axial direction. This paper presents the detailed measurements of mean velocity profiles, turbulent intensity distributions and RMS pressure fluctuations. The energy spectra obtained for the pressure fluctuations and the u and v velocity components are compared. Contrary to earlier works, we found that the differences between peak frequencies of the pressure spectra and the characteristic frequencies of the velocity spectra vary with positions downstream from the nozzle. These differences are discussed in light of pseudosound generation by the eddy structures in the stenotic flow field.     

Some flow visualization and laser-Doppler-velocity measurements in a true-to-scale elastic model of a human aortic arch--a new model technique.

Flow studies were done in an elastic true-to-scale silicone rubber model of an aortic arch to study further hemodynamic influences on atherosclerosis. The model was prepared from a cast of a young woman. A revised model technique was used. The model had a compliance similar to that of the human aortic arch. Velocity measurements were done in the model with a two component laser-Doppler-anemometer in steady and pulsatile flow using a calcium chloride solution with a viscosity of eta = 3.18 mPas and density of rho = 1.28 kg/m3 at 20 degrees C. The time average Reynolds numbers over a whole cycle in the ascending aorta was Re = 1350. The Womersley parameter for pulsatile flow was a = 20. The pulse wave velocity in the ascending aorta was about c = 5.4 m/sec. The secondary flow behavior was discussed for steady and pulsatile flow. Reverse flows were found, especially along the inner radius of the aortic arch in the descending aorta in steady and pulsatile flow and also in small areas of the ascending aorta and at the branches of the aortic arch. The formation of atherosclerotic plaques at preferred local flow regions is discussed.     

Two-dimensional color-mapping of turbulent shear stress distribution downstream of two aortic bioprosthetic valves in vitro.

Since artificial heart valve related complications such as thrombus formation, hemolysis and calcification are considered related to flow disturbances caused by the inserted valve, a thorough hemodynamic characterization of heart valve prostheses is essential. In a pulsatile flow model, fluid velocities were measured one diameter downstream of a Hancock Porcine (HAPO) and a Ionescu-Shiley Pericardial Standard (ISPS) aortic valve. Hot-film anemometry (HFA) was used for velocity measurements at 41 points in the cross-sectional area of the ascending aorta. Three-dimensional visualization of the velocity profiles, at 100 different instants during one mean pump cycle, was performed. Turbulence analysis was performed as a function of time by calculating the axial turbulence energy within 50 ms overlapping time windows during the systole. The turbulent shear stresses were estimated by using the correlation equation between Reynolds normal stress and turbulent (Reynolds) shear stress. The turbulent shear stress distribution was visualized by two-dimensional color-mapping at different instants during one mean pump cycle. Based on the velocity profiles and the turbulent shear stress distribution, a relative blood damage index (RBDI) was calculated. It has the feature of combining the magnitude and exposure time of the estimated shear stresses in one index, covering the entire cross-sectional area. The HAPO valve showed a skewed jet-type velocity profile with the highest velocities towards the left posterior aortic wall. The ISPS valve revealed a more parabolic-shaped velocity profile during systole. The turbulent shear stresses were highest in areas of high or rapidly changing velocity gradients. For the HAPO valve the maximum estimated turbulent shear stress was 194 N m-2 and for the ISPS valve 154 Nm-2. The RBDI was the same for the two valves. The turbulent shear stresses had magnitudes and exposure times that might cause endothelial damage and sublethal or lethal damage to blood corpuscules. The RBDI makes comparison between different heart valves easier and may prove important when making correlation with clinical observations.     

Comparison of blood flow parameters measured by Doppler ultrasound and radionuclide techniques in the baboon (Papio ursinus)

Ascending aortic blood velocity was measured in the baboon (Papio ursinus) by using continuous wave Doppler ultrasound. The blood flow parameters thus obtained were compared to those by the standardized radionuclide technique. It appears that, due to the anatomical position of the ascending aorta and brachiocephalic trunks in relation to the ultrasound beam, Doppler ultrasound does not provide an accurate method of measuring aortic blood velocity in the baboon, which could be the reason for the poor correlation of the results from the two techniques.     

Comparison of hinge microflow fields of bileaflet mechanical heart valves implanted in different sinus shape and downstream geometry

The characterization of the bileaflet mechanical heart valves (BMHVs) hinge microflow fields is a crucial step in heart valve engineering. Earlier in vitro studies of BMHV hinge flow at the aorta position in idealized straight pipes have shown that the aortic sinus shapes and sizes may have a direct impact on hinge microflow fields. In this paper, we used a numerical study to look at how different aortic sinus shapes, the downstream aortic arch geometry, and the location of the hinge recess can influence the flow fields in the hinge regions. Two geometric models for sinus were investigated: a simplified axisymmetric sinus and an idealized three-sinus aortic root model, with two different downstream geometries: a straight pipe and a simplified curved aortic arch. The flow fields of a 29-mm St Jude Medical BMHV with its four hinges were investigated. The simulations were performed throughout the entire cardiac cycle. At peak systole, recirculating flows were observed in curved downsteam aortic arch unlike in straight downstream pipe. Highly complex three-dimensional leakage flow through the hinge gap was observed in the simulation results during early diastole with the highest velocity at 4.7 m/s, whose intensity decreased toward late diastole. Also, elevated wall shear stresses were observed in the ventricular regions of the hinge recess with the highest recorded at 1.65 kPa. Different flow patterns were observed between the hinge regions in straight pipe and curved aortic arch models. We compared the four hinge regions at peak systole in an aortic arch downstream model and found that each individual hinge did not vary much in terms of the leakage flow rate through the valves.     

Numerical simulation of fluid–structure interaction in bypassed DeBakey III aortic dissection

It was found that bypass graft alone could achieve great effects in treating aortic dissection. In order to investigate the mechanical mechanism and the haemodynamic validity of the bypassing treatment for DeBakey III aortic dissection, patient-specific models of DeBakey III aortic dissection treated with different bypassing strategies were constructed. One of the bypassing strategies is bypassing between ascending aorta and abdominal aorta, and the other is bypassing between left subclavian artery and abdominal aorta. Numerical simulations under physiological flow conditions based on fluid–structure interaction were performed using finite element method. The results show that blood flow velocity, pressure and vessel wall displacement of false lumen are all reduced after bypassing. This phenomenon indicates that bypassing is an effective surgery for the treatment of DeBakey III aortic dissection. The effectiveness to cure through lumen is better when bypassing between left subclavian artery and abdominal aorta, while the effectiveness to cure blind lumen is better when bypassing between ascending aorta and abdominal aorta.