Physiological Pressure Waveforms Research Articles

A contact model based on the coefficient of restitution for simulations of bio-prosthetic heart valves.

A new general contact model is proposed for preventing inter-leaflet penetration of bio-prosthetic heart valves (BHV) at the end of the systole, which has the advantage of applying kinematic constraints directly and creating smooth free edges. At the end of each time step, the impenetrability constraints and momentum exchange between the impacting bodies are applied separately based on the coefficient of restitution. The contact method is implemented in a rotation-free, large deformation, and thin shell finite-element (FE) framework based on loop's subdivision surfaces. A nonlinear, anisotropic material model for a BHV is employed which uses Fung-elastic constitutive laws for in-plane and bending responses, respectively. The contact model is verified and validated against several benchmark problems. For a BHV-specific validation, the computed strains on different regions of a BHV under constant pressure are compared with experimentally measured data. Finally, dynamic simulations of BHV under physiological pressure waveform are performed for symmetrical and asymmetrical fiber orientations incorporating the new contact model and compared with the penalty contact method. The proposed contact model provides the coaptation area of a functioning BHV during the closing phase for both of the fiber orientations. Our results show that fiber orientation affects the dynamic of leaflets during the opening and closing phases. A swirling motion for the BHV with asymmetrical fiber orientation is observed, similar to experimental data. To include the fluid effects, fluid-structure interaction (FSI) simulation of the BHV is performed and compared to the dynamic results.

Computational fluid dynamics investigation on aortic hemodynamics in double aortic arch before and after ligation surgery

Double aortic arch (DAA) malformation is one of the reasons for symptomatic vascular rings, the hemodynamics of which is still poorly understood. This study aims to investigate the blood flow characteristics in patient-specific double aortic arches using computational fluid dynamics (CFD). Seven cases of infantile patients with DAA were collected and their computed tomography images were used to reconstruct 3D computational models. A modified Carreau model was used to consider the non-Newtonian effect of blood and a three-element Windkessel model taking the effect of the age of patients into account was applied to reproduce physiological pressure waveforms. Numerical results show that blood flow distribution and energy loss of DAA depends on relative sizes of the two aortic arches and their angles with the ascending aorta. Ligation of either aortic arch increases the energy loss of blood in the DAA, leading to the increase in cardiac workload. Generally, the rising rate of energy loss before and after the surgery is almost linear with the area ratio between the aortic arch without ligation and the ascending aorta.

Establishment of a Modular Hemodynamic Simulator for Accurate In Vitro Simulation of Physiological and Pathological Pressure Waveforms in Native and Bioartificial Blood Vessels

PurposeIn vitro stimulation of native and bioartificial vessels in perfusable systems simulating natural mechanical environments of the human vasculature represents an emerging approach in cardiovascular research. Promising results have been achieved for applications in both regenerative medicine and etiopathogenetic investigations. However, accurate and reliable simulation of the wide variety of physiological and pathological pressure environments observed in different vessels still remains an unmet challenge.MethodsWe established a modular hemodynamic simulator (MHS) with interchangeable and modifiable components suitable for the perfusion of native porcine—(i.e. the aorta, brachial and radial arteries and the inferior vena cava) and bioartificial fibrin-based vessels with anatomical site specific pressure curves. Additionally, different pathological pressure waveforms associated with cardiovascular diseases including hyper- and hypotension, tachy- and bradycardia, aortic valve stenosis and insufficiency, heart failure, obstructive cardiomyopathy and arterial stiffening were simulated. Pressure curves, cyclic distension and shear stress were measured for each vessel and compared to ideal clinical pressure waveforms.ResultsThe pressure waveforms obtained in the MHS showed high similarity to the ideal anatomical site specific pressure curves of different vessel types. Moreover, the system facilitated accurate emulation of physiological and different pathological pressure conditions in small diameter fibrin-based vessels.ConclusionThe MHS serves as a variable in vitro platform for accurate emulation of physiological and pathological pressure environments in biological probes. Potential applications of the system include bioartificial vessel maturation in cardiovascular tissue engineering approaches as well as etiopathogenetic investigations of various cardiovascular pathologies.

Biomechanical implications of the fenestration structure after thoracic endovascular aortic repair

The inadequate landing zone during thoracic endovascular aortic repair (TEVAR) could be resolved by fenestration technology. The fenestration structure consists of a main endograft and a left subclavian artery (LSA) stent-graft. The purpose of this study is to assess the biomechanical implications of the protruding segment (PS) of the LSA stent-graft after TEVAR with in situ fenestration (ISF-TEVAR). The PS is characterized by the protruding length and centerline angle between the LSA and PS. An idealized three-dimensional geometric model of the human aorta was constructed as the reference benchmark. We designed nine postoperative aortic geometries with different protruding lengths (5, 10, 15 mm) and centerline angles (-20°, 0°, +20°). The blood was assumed to be non-Newtonian and the three-element Windkessel model was applied to reproduce physiological pressure waveforms. The interaction between the blood and vessel wall was captured by a two-way fluid-structure method. We also considered the impact of the fenestration structure on the vessel wall. Long protruding length (15 mm) and the retrograde angle (+20°) result in a markedly reduced LSA flow ratio (1.3%). There is a pressure difference between the inner and outer walls of the PS. The region around the PS is more prone to thrombosis. The flow stability and turbulence intensity of downstream blood of the PS gradually deteriorate. The largest deformed region moves from the aortic arch to the ascending aorta after ISF-TEVAR. The postoperative hemodynamics largely depends on the protruding length and angle of the LSA stent-graft. The configurations should be carefully controlled during ISF-TEVAR.

Dynamic physiological temperature and pressure sensing with phase-resolved low-coherence interferometry

We report the development and characterisation of highly miniaturised fibre-optic sensors for simultaneous pressure and temperature measurement, and a compact interrogation system with a high sampling rate. The sensors, which have a maximum diameter of 250 µm, are based on multiple low-finesse optical cavities formed from polydimethylsiloxane (PDMS), positioned at the distal ends of optical fibres, and interrogated using phase-resolved low-coherence interferometry. At acquisition rates of 250 Hz, temperature and pressure changes of 0.0021 °C and 0.22 mmHg are detectable. An in vivo experiment demonstrated that the sensors had sufficient speed and sensitivity for monitoring dynamic physiological pressure waveforms. These sensors are ideally suited to various applications in minimally invasive surgery, where diminutive lateral dimensions, high sensitivity and low manufacturing complexities are particularly valuable.

Abstract 17080: A 3D Printed Ex Vivo Left Heart Simulator Quantifies and Validates Posterior Ventricular Anchoring Neochordoplasty

Introduction: The posterior ventricular anchoring neochordal (PVAN) repair is a nonresectional, single-suture technique for correcting posterior leaflet prolapse. While this technique has demonstrated clinical efficacy, a possible limitation is the stability of the suture anchored into myocardium as opposed to the fibrous portion of a papillary muscle. Hypothesis: We hypothesize that the PVAN suture serves only to position the leaflet for coaptation, after which systolic forces will be distributed throughout the valve, resulting in low peak forces on the suture. Methods: A left heart simulator was constructed using 3D printing, tuned to generate physiological pressure and flow waveforms, then validated. Porcine mitral valves (n=9) were dissected and mounted within the simulator. Chordal forces were measured using Fiber Bragg Grating (FBG) sensors, sewn in place using PTFE suture. FBG sensors are strain gauges made of 125 μ m optical fibers that use reflected peak wavelength changes to measure strain. Hemodynamic and echocardiographic data were also collected. Isolated severe mitral regurgitation (MR) was induced by cutting P2 primary chordae. The valve was repaired using the PVAN technique, anchoring the suture to a customized force-sensing post positioned to mimic in vivo placement. Results: Forces on 1° and 2° chordae of both anterior and posterior leaflets were significantly elevated in the prolapse condition ( P < 0.05). PVAN resulted in elimination of MR in all valves, as well as normalization of chordae forces to baseline levels for posterior primary ( P < 0.01 ) , posterior secondary ( P < 0.01 ) , and anterior primary chordae ( P < 0.05 ) , with reduction in anterior secondary chordal forces approaching significance ( P = 0.055 ) . Peak forces on the PVAN stitch were minimal, even compared to the forces experienced by primary chordae of normal, healthy valves ( P < 0.05). Conclusions: The PVAN technique eliminates MR by effectively positioning the posterior leaflet for optimal coaptation, distributing the forces amongst the subvalvular apparatus. Given the extremely low forces involved, the strength of the ventricular anchoring suture and myocardial anchoring point should not be a limiting factor.

On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics

Boundary conditions (BCs) are an essential part in computational fluid dynamics (CFD) simulations of blood flow in large arteries. Although several studies have investigated the influence of BCs on predicted flow patterns and hemodynamic wall parameters in various arterial models, there is a lack of comprehensive assessment of outlet BCs for patient-specific analysis of aortic flow. In this study, five different sets of outlet BCs were tested and compared using a subject-specific model of a normal aorta. Phase-contrast magnetic resonance imaging (PC-MRI) was performed on the same subject and velocity profiles extracted from the in vivo measurements were used as the inlet boundary condition. Computational results obtained with different outlet BCs were assessed in terms of their agreement with the PC-MRI velocity data and key hemodynamic parameters, such as pressure and flow waveforms and wall shear stress related indices. Our results showed that the best overall performance was achieved by using a well-tuned three-element Windkessel model at all model outlets, which not only gave a good agreement with in vivo flow data, but also produced physiological pressure waveforms and values. On the other hand, opening outlet BCs with zero pressure at multiple outlets failed to reproduce any physiologically relevant flow and pressure features.

Physiological outflow boundary conditions methodology for small arteries with multiple outlets: a patient-specific hepatic artery haemodynamics case study.

Physiological outflow boundary conditions are necessary to carry out computational fluid dynamics simulations that reliably represent the blood flow through arteries. When dealing with complex three-dimensional trees of small arteries, and therefore with multiple outlets, the robustness and speed of convergence are also important. This study derives physiological outflow boundary conditions for cases in which the physiological values at those outlets are not known (neither in vivo measurements nor literature-based values are available) and in which the tree exhibits symmetry to some extent. The inputs of the methodology are the three-dimensional domain and the flow rate waveform and the systolic and diastolic pressures at the inlet. The derived physiological outflow boundary conditions, which are a physiological pressure waveform for each outlet, are based on the results of a zero-dimensional model simulation. The methodology assumes symmetrical branching and is able to tackle the flow distribution problem when the domain outlets are at branches with a different number of upstream bifurcations. The methodology is applied to a group of patient-specific arteries in the liver. The methodology is considered to be valid because the pulsatile computational fluid dynamics simulation with the inflow flow rate waveform (input of the methodology) and the derived outflow boundary conditions lead to physiological results, that is, the resulting systolic and diastolic pressures at the inlet match the inputs of the methodology, and the flow split is also physiological.

Experimental study of hemodynamics in the Circle of Willis.

BackgroundThe Circle of Willis (CoW) is an important collateral pathway of the cerebral blood flow. An experimental study of the cerebral blood flow (CBF) distribution in different anatomical variations may help to a better understanding of the collateral mechanism of the CoW.MethodsAn in-vitro test rig was developed to simulate the physiological cerebral blood flow in the CoW. Ten anatomical variations were considered in this study, include a set of different degrees of stenosis in L-ICA and L-ICA occlusion coexist with common anatomical variations. Volume flow rates of efferent arteries and pressure signals at the end of communicating arteries of each case were recorded. Physiological pressure waveforms were applied as inlet boundary condition.ResultsIn the development of L-ICA stenosis, the total CBF decreases with the increase of stenosis degree. The blood supply of ipsilateral middle cerebral artery (MCA) was affected most by the stenosis of L-ICA. Anterior communicating artery (ACoA) and ipsilateral posterior communicating artery (PCoA) function as important collateral pathways of cerebral collateral circulation when unilateral stenosis occurred. The blood supply of anterior cerebral circulation was compensated by the posterior cerebral circulation through ipsilateral PCoA when L-ICA stenosis degree is greater than 40% and the affected side was compensated immediately by the unaffected side through ACoA. Blood flow of the anterior circulation and the total CBF reached the minimum among all cases studied when L-ICA occlusion coexist with the absence of PCoA.ConclusionThe results demonstrated the flow distribution patterns of the CoW under anatomical variations and clarified the collateral mechanism of the CoW. The flow ACoA is the most sensitive indexes to the morphology change of ipsilateral ICA. The relative independence of the circulation in anterior and posterior sections of the CoW is not broken and the function of ipsilateral PCoA is not activated until a severe stenosis of unilateral ICA occurs. PCoA is the most important collateral pathway of the collateral circulation and the missing of PCoA has the highest risk of stroke when the ipsilateral ICA has severe stenosis. These findings may provide the basis for future therapeutic and diagnosis applications.

Evaluation of the hemodynamics in straight 6-mm and tapered 6- to 8-mm grafts as upper arm hemodialysis vascular access.

The present study is intended to investigate and compare the hemodynamics in two different sizes of hemodialysis arteriovenous grafts for upper arm hemodialysis vascular access: 8-mm tapered to 6-mm at the arterial side and straight 6mm. A computational simulation approach is presented for this study, which is validated against the available experimental and numerical pressure measurements in the literature. The imposed boundary conditions at the arterial inlet and venous outlet boundaries of the models are physiological velocity and pressure waveforms, respectively. Blood flow fields and distribution patterns of the hemodynamic indices including wall shear stress (WSS) as one of the major hemodynamic parameters of the cardiovascular system and spatial wall shear stress gradient (SWSSG) as an indicator of disturbed flow patterns and hence susceptible sites of lesion developments are analyzed and compared between the two grafts. The tapered 6- to 8-mm graft seemingly is associated with less disturbed flow patterns within the venous anastomosis (VA) and the vein downstream while benefiting from higher blood flow rates within. Also, it shows a definitive advantage in terms of WSS and SWSSG distribution patterns around the VA and throughout the vein downstream with significantly lower values, which reduce the risk of thrombosis formation and stenotic lesion developments. The only disadvantage encountered in using 6- to 8-mm tapered graft is higher values of hemodynamic parameters at the arterial junction attributable to its significantly higher mean blood flow rate within. The results clearly indicate that the tapered 6- to 8-mm graft entirely outperforms straight 6-mm graft hemodynamically as an upper arm hemodialysis vascular access graft and confirms clinical data in the literature, which suggests advantageous use of tapered 6- to 8-mm grafts in the creation of upper arm brachioaxillary hemodialysis vascular access grafts in selected groups of patients with expectably higher patency rates and lower complications.

THE INFLUENCE OF USING THE NEEDLE ADAPTER TO REDUCE THE BIOMECHANICAL RISK FACTORS WITHIN HEMODIALYSIS ARTERIOVENOUS GRAFTS

Hemodialysis vascular access failure is related to increased morbidity and mortality in hemodialysis patients, representing a challenging clinical problem which results in a high percentage of hospital entrance and an important economic burden on government's disbursement. In this paper, the feasibility of using the needle adapter to reduce the biomechanical risk factors within arteriovenous grafts is considered. The three-dimensional (3D) tapered 6 to 8 mm loop graft in the presence of venous and arterial needles with and without adapter was numerically simulated. Navier–Stokes equations for incompressible Newtonian fluid are the governing equation of this problem. k – ω two equations turbulence modeling were applied to capture flow features of low Reynolds turbulent flow regions in this simulation. The physiological velocity waveform was used as an arterial inlet boundary condition. The venose outlet boundary condition was a time dependent physiological pressure waveform. The results for the dialysis without the adapter demonstrated that the graft wall experiences increased hemodynamic stresses as a result of the hitting needle jet flow. The dialysis with the adapter demonstrated that the venous anastomosis experiences lower biomechanical risk factors in comparison to the dialysis without the adapter and it reduced the vascular access failure. Using adapter caused less damage to endothelial cells during hemodialysis.

Unsteady blood flow and mass transfer of a human left coronary artery bifurcation: FSI vs. CFD

In this study, a Fluid Solid Interaction analysis (FSI) of a computerized tomography (CT) scan reconstructed left coronary artery was performed. The arterial wall was modeled as an isotropic hyperelastic material. The arterial wall shear stress (WSS) was computed in order to investigate a correlation between flow-induced wall shear stress and geometry of the artery. An unsteady state FSI analysis with a commercial finite element software was performed in order to evaluate the maximum and the minimum wall shear stress as a function of the flow regime and the arterial wall compliance in the left coronary. As boundary conditions, physiological pressure waveforms were applied. Comparison of the computational results between the FSI and rigid-wall models showed that the wall shear stress (WSS) distributions were substantially affected by the arterial wall compliance. In particular, the minimum and maximum WSS values significantly vary.

Immersed boundary model of aortic heart valve dynamics with physiological driving and loading conditions.

The immersed boundary (IB) method is a mathematical and numerical framework for problems of fluid–structure interaction, treating the particular case in which an elastic structure is immersed in a viscous incompressible fluid. The IB approach to such problems is to describe the elasticity of the immersed structure in Lagrangian form, and to describe the momentum, viscosity, and incompressibility of the coupled fluid–structure system in Eulerian form. Interaction between Lagrangian and Eulerian variables is mediated by integral equations with Dirac delta function kernels. The IB method provides a unified formulation for fluid–structure interaction models involving both thin elastic boundaries and also thick viscoelastic bodies. In this work, we describe the application of an adaptive, staggered-grid version of the IB method to the three-dimensional simulation of the fluid dynamics of the aortic heart valve. Our model describes the thin leaflets of the aortic valve as immersed elastic boundaries, and describes the wall of the aortic root as a thick, semi-rigid elastic structure. A physiological left-ventricular pressure waveform is used to drive flow through the model valve, and dynamic pressure loading conditions are provided by a reduced (zero-dimensional) circulation model that has been fit to clinical data. We use this model and method to simulate aortic valve dynamics over multiple cardiac cycles. The model is shown to approach rapidly a periodic steady state in which physiological cardiac output is obtained at physiological pressures. These realistic flow rates are not specified in the model, however. Instead, they emerge from the fluid–structure interaction simulation.

Aortic flow patterns resulting from right axillary artery cannulation

Right axillary artery (RAA) cannulation is increasingly used in cardiac surgery. Little is known about resulting flow patterns in the aorta. Therefore, flow was visualized and analyzed. A mock circulatory circuit was assembled based on a compliant transparent anatomical silicon aortic model. A RAA cannula was connected to a continuous flow rotary blood pump (RBP), pulsatile heart action was provided by a pneumatic ventricular assist device (PVAD). Peripheral vascular resistance, regional flow and vascular compliance were adjusted to obtain physiological flow and pressure waveforms. Colorants were injected automatically for flow visualization. Five flow distributions with a total flow of 4l/min were tested (%PVAD:%RBP): 100:0, 75:25, 50:50, 25:75, 0:100. Colorant distribution was assessed using quantitative 2D image processing. Continuous flow from the RAA divided in a retrograde and an antegrade portion. Retro- to antegrade flow ratio increased with increasing RAA-flow. At full RBP support flow was stagnant in the ascending aorta. There were distinct flow patterns between the right- and left-sided supra-aortic branches. At full RBP support retrograde flow was demonstrated in the right carotid and right vertebral arteries. Further studies are needed to confirm and evaluate the described flow patterns.

Theoretical modeling of airways pressure waveform for dual-controlled ventilation with physiological pattern and linear respiratory mechanics

The present paper describes the theoretical treatment performed for the geometrical optimization of advanced and improved-shape waveforms as airways pressure excitation for controlled breathings in dual-controlled ventilation applied to anaesthetized or severe brain injured patients, the respiratory mechanics of which can be assumed linear. Advanced means insensitive to patient breathing activity as well as to ventilator settings while improved-shape intends in comparison to conventional square waveform for a progressive approaching towards physiological transpulmonary pressure and respiratory airflow waveforms. Such functional features along with the best ventilation control for the specific therapeutic requirements of each patient can be achieved through the implementation of both diagnostic and compensation procedures effectively carried out by the Advance Lung Ventilation System (ALVS) already successfully tested for square waveform as airways pressure excitation. Triangular and trapezoidal waveforms have been considered as airways pressure excitation. The results shows that the latter fits completely the requirements for a physiological pattern of endoalveolar pressure and respiratory airflow waveforms, while the former exhibits a lower physiological behaviour but it is anyhow periodically recommended for performing adequately the powerful diagnostic procedure.

Mathematical modeling and experimental testing of three bioreactor configurations based on windkessel models

This paper presents an experimental study of three bioreactor configurations. The bioreactor is intended to be used for the development of tissue-engineered heart valve substitutes. Therefore it must be able to reproduce physiological flow and pressure waveforms accurately. A detailed analysis of three bioreactor arrangements is presented using mathematical models based on the windkessel (WK) approach. First, a review of the many applications of this approach in medical studies enhances its fundamental nature and its usefulness. Then the models are developed with reference to the actual components of the bioreactor. This study emphasizes different conflicting issues arising in the design process of a bioreactor for biomedical purposes, where an optimization process is essential to reach a compromise satisfying all conditions. Two important aspects are the need for a simple system providing ease of use and long-term sterility, opposed to the need for an advanced (thus more complex) architecture capable of a more accurate reproduction of the physiological environment. Three classic WK architectures are analyzed, and experimental results enhance the advantages and limitations of each one.

Fluid–structure interaction analysis of a patient‐specific right coronary artery with physiological velocity and pressure waveforms

AbstractCoupled fluid–structure interaction (FSI) analysis of the human right coronary artery (RCA) has been carried out to investigate the effects of wall compliance on coronary hemodynamics. A 3‐D model of a stenosed RCA was reconstructed based on multislice computerized tomography images. A velocity waveform in the proximal RCA and a pressure waveform in the distal RCA of a patient with a severe stenosis were acquired with a catheter delivered wire probe and applied as boundary conditions. The arterial wall was modeled as a Mooney–Rivlin hyperelastic material. The predicted maximum wall displacement (3.85 mm) was comparable with the vessel diameter (∼4 mm), but the diameter variation was much smaller, 0.134 mm at the stenosis and 0.486 mm in the distal region. Comparison of the computational results between the FSI and rigid‐wall models showed that the instantaneous wall shear stress (WSS) distributions were affected by diameter variation in the arterial wall; increasing systolic blood pressure dilated the vessel and consequently lowered WSS, whereas the opposite occurred when pressure started to decrease. However, the effects of wall compliance on time‐averaged WSS (TAWSS) and oscillatory shear index (OSI) were insignificant (4.5 and 2.7% difference in maximum TAWSS and OSI, respectively). Copyright © 2009 John Wiley & Sons, Ltd.

Neonatal Aortic Arch Hemodynamics and Perfusion During Cardiopulmonary Bypass

The objective of this study is to quantify the detailed three-dimensional (3D) pulsatile hemodynamics, mechanical loading, and perfusion characteristics of a patient-specific neonatal aortic arch during cardiopulmonary bypass (CPB). The 3D cardiac magnetic resonance imaging (MRI) reconstruction of a pediatric patient with a normal aortic arch is modified based on clinical literature to represent the neonatal morphology and flow conditions. The anatomical dimensions are verified from several literature sources. The CPB is created virtually in the computer by clamping the ascending aorta and inserting the computer-aided design model of the 10 Fr tapered generic cannula. Pulsatile (130 bpm) 3D blood flow velocities and pressures are computed using the commercial computational fluid dynamics (CFD) software. Second order accurate CFD settings are validated against particle image velocimetry experiments in an earlier study with a complex cardiovascular unsteady benchmark. CFD results in this manuscript are further compared with the in vivo physiological CPB pressure waveforms and demonstrated excellent agreement. Cannula inlet flow waveforms are measured from in vivo PC-MRI and 3 kg piglet neonatal animal model physiological experiments, distributed equally between the head-neck vessels and the descending aorta. Neonatal 3D aortic hemodynamics is also compared with that of the pediatric and fetal aortic stages. Detailed 3D flow fields, blood damage, wall shear stress (WSS), pressure drop, perfusion, and hemodynamic parameters describing the pulsatile energetics are calculated for both the physiological neonatal aorta and for the CPB aorta assembly. The primary flow structure is the high-speed canulla jet flow (approximately 3.0 m/s at peak flow), which eventually stagnates at the anterior aortic arch wall and low velocity flow in the cross-clamp pouch. These structures contributed to the reduced flow pulsatility (85%), increased WSS (50%), power loss (28%), and blood damage (288%), compared with normal neonatal aortic physiology. These drastic hemodynamic differences and associated intense biophysical loading of the pathological CPB configuration necessitate urgent bioengineering improvements--in hardware design, perfusion flow waveform, and configuration. This study serves to document the baseline condition, while the methodology presented can be utilized in preliminary CPB cannula design and in optimization studies reducing animal experiments. Coupled to a lumped-parameter model the 3D hemodynamic characteristics will aid the surgical decision making process of the perfusion strategies in complex congenital heart surgeries.

Dynamic Stress Analysis of the Arterial Wall Utilizing Physiological Pressure Waveforms

Arterial diseases, promoted by alterations in arterial wall properties, are among the main causes of mortality. Mechanical stresses of the arterial wall caused by pulsatile luminal pressure define arterial function in normal and pathological conditions. This study aims to determine dynamic stress distribution in the arterial wall subjected to physiological pressure waveforms. Finite element models of a typical artery are developed to evaluate Von Misses stress in the arterial wall due to physiological pressure waveforms and with differing mechanical properties. Mechanical parameters include Young's modulus of elasticity, non-linear stress-strain relationship and visco-elastic parameter. The appropriate boundary conditions are allocated to allow radial expansion. Application of physiological pulsatile pressure results in stress waves with the values and waveforms markedly influenced by mechanical properties of the arterial wall and blood pressure pulse. Elevated elastic modulus of the arterial wall results in significant increase in maximum stress. Viscoelastic property leads to reduction of the peak stress and smoothening of the stress waveform. The pressure waveform is also a major parameter affecting the stress pattern in arterial wall. Hypertensive arteries result in higher and sharper stress waves not only because of a higher systolic value but also because of the sharper waveform. The combination of these parameters produces the resultant stress pattern in the arterial wall.

Particle Image Velocimetry Study of Pulsatile Flow in Bi-leaflet Mechanical Heart Valves with Image Compensation Method

Particle Image Velocimetry (PIV) is an important technique in studying blood flow in heart valves. Previous PIV studies of flow around prosthetic heart valves had different research concentrations, and thus never provided the physical flow field pictures in a complete heart cycle, which compromised their pertinence for a better understanding of the valvular mechanism. In this study, a digital PIV (DPIV) investigation was carried out with improved accuracy, to analyse the pulsatile flow field around the bi-leaflet mechanical heart valve (MHV) in a complete heart cycle. For this purpose a pulsatile flow test rig was constructed to provide the necessary in vitro test environment, and the flow field around a St. Jude size 29 bi-leaflet MHV and a similar MHV model were studied under a simulated physiological pressure waveform with flow rate of 5.2 l/min and pulse rate at 72 beats/min. A phase-locking method was applied to gate the dynamic process of valve leaflet motions. A special image-processing program was applied to eliminate optical distortion caused by the difference in refractive indexes between the blood analogue fluid and the test section. Results clearly showed that, due to the presence of the two leaflets, the valvular flow conduit was partitioned into three flow channels. In the opening process, flow in the two side channels was first to develop under the presence of the forward pressure gradient. The flow in the central channel was developed much later at about the mid-stage of the opening process. Forward flows in all three channels were observed at the late stage of the opening process. At the early closing process, a backward flow developed first in the central channel. Under the influence of the reverse pressure gradient, the flow in the central channel first appeared to be disturbed, which was then transformed into backward flow. The backward flow in the central channel was found to be the main driving factor for the leaflet rotation in the valve closing process. After the valve was fully closed, local flow activities in the proximity of the valve region persisted for a certain time before slowly dying out. In both the valve opening and closing processes, maximum velocity always appeared near the leaflet trailing edges. The flow field features revealed in the present paper improved our understanding of valve motion mechanism under physiological conditions, and this knowledge is very helpful in designing the new generation of MHVs.