Computational simulations of the aortic valve validated by imaging data: evaluation of valve-sparing techniques☆

Multiscale structure and function of the aortic valve apparatus

Whereas studying the aortic valve in isolation has facilitated the development of life-saving procedures and technologies, the dynamic interplay of the aortic valve and its surrounding structures is vital to preserving their function across the wide range of conditions encountered in an active lifestyle. Our view is that these structures should be viewed as an integrated functional unit, here referred to as the aortic valve apparatus (AVA). The coupling of the aortic valve and root, left ventricular outflow tract, and blood circulation is crucial for AVA's functions: unidirectional flow out of the left ventricle, coronary perfusion, reservoir function, and support of left ventricular function. In this review, we explore the multiscale biological and physical phenomena that underlie the simultaneous fulfillment of these functions. A brief overview of the tools used to investigate the AVA, such as medical imaging modalities, experimental methods, and computational modeling, specifically fluid-structure interaction (FSI) simulations, is included. Some pathologies affecting the AVA are explored, and insights are provided on treatments and interventions that aim to maintain quality of life. The concepts explained in this article support the idea of AVA being an integrated functional unit and help identify unanswered research questions. Incorporating phenomena through the molecular, micro, meso, and whole tissue scales is crucial for understanding the sophisticated normal functions and diseases of the AVA.

Valve-Sparing Aortic Root Replacement State-of-the-Art Review, Part II: Surgical Techniques

Aortic valve disease is common, and valve-preserving operations are preferred whenever possible. Valve-sparing aortic root replacement (VSRR) has become an important tool for managing aortic root pathology in children and adults. The learning curve for this operation is challenging, but with increasing experience and technical modifications, early and late outcomes continue to improve. Durable long term results vary based on underlying anatomy, pathology, and patient selection, as well as surgeon expertise. Part II of this VSRR State-of-the-Art Review article provides technical pearls related to VSRR.

Aortic Root Vortex Formation During Left Ventricular Assist Device Support

The vortex that forms in the aortic sinus plays a vital role in optimizing blood flow. Disruption of the vortex can result in flow stagnation and activate thrombus formation in the aortic root, especially when aortic valve flow is reduced as during left ventricular assist device (LVAD) support. Our goal in this study was to visualize vortex formation in an experimental model of the aortic root as flow is progressively reduced. A mock circulatory loop that reproduces heart failure hemodynamics was combined with a HeartMate II LVAD and velocity measured in a transparent aortic root with a bioprosthetic valve. The aortic valve sinus vortices are clearly visible as counter-rotating structures in the velocity field at baseline and for all conditions with flow through the aortic valve. As LVAD speed increases, the central jet narrows but the vortices persist, disappearing only when the valve is completely closed. The vortices preserve fluid momentum and generate shear stress along the tissue surfaces which disrupts flow stasis. These features underscore the importance of maintaining "intermittent" aortic valve opening, as recommended for LVAD patients. This study is the first to report vortex formation in the aortic root during LVAD support, providing a motivation for further evaluation.

Heart Valve Biomechanics: The Frontiers of Modeling Modalities and the Expansive Capabilities of Ex Vivo Heart Simulation

The field of heart valve biomechanics is a rapidly expanding, highly clinically relevant area of research. While most valvular pathologies are rooted in biomechanical changes, the technologies for studying these pathologies and identifying treatments have largely been limited. Nonetheless, significant advancements are underway to better understand the biomechanics of heart valves, pathologies, and interventional therapeutics, and these advancements have largely been driven by crucial in silico, ex vivo, and in vivo modeling technologies. These modalities represent cutting-edge abilities for generating novel insights regarding native, disease, and repair physiologies, and each has unique advantages and limitations for advancing study in this field. In particular, novel ex vivo modeling technologies represent an especially promising class of translatable research that leverages the advantages from both in silico and in vivo modeling to provide deep quantitative and qualitative insights on valvular biomechanics. The frontiers of this work are being discovered by innovative research groups that have used creative, interdisciplinary approaches toward recapitulating in vivo physiology, changing the landscape of clinical understanding and practice for cardiovascular surgery and medicine.

Leaflet stress quantification of porcine vs bovine surgical bioprostheses: an in vitro study

Calcified aortic stenoses are among the most prevalent form of cardiovascular diseases in the industrialized countries. This progressive disease, with no effective medical therapy, ultimately requires aortic valve replacement - either a surgical or very recently transcatheter aortic valve implantation. Increase leaflet mechanical stress is one of the main determinants of the structural deterioration of bioprosthetic aortic valves. We applied a coupled in vitro/in silico method to compare the timing, magnitude, and regional distribution of leaflet mechanical stress in porcine versus pericardial bioprostheses (Mosaic and Trifecta). A double activation simulator was used for in vitro testing of a bioprosthesis with externally mounted pericardium (Abbott, Trifecta) and a bioprosthesis with internally mounted porcine valve (Medtronic, Mosaic). A non-contact system based on stereophotogammetry and digital image correlation (DIC) with high spatial and temporal resolution (2000 img/s) was used to visualize the valve leaflet motion and perform the three-dimensional analysis. A finite element model of the valve was developed, and the leaflet deformation obtained from the DIC analysis was applied to the finite element model calculate local leaflet mechanical stress throughout the cardiac cycle. The maximum leaflet stress was higher with the pericardial versus the porcine bioprosthesis (2.03 vs. 1.30 MPa) For both bioprostheses the highest values of leaflet stress occurred during diastole and were primarily observed in the upper leaflet edge near the commissures and to a lesser extent in the mid-portion of the leaflet body. In conclusion, the coupled in vitro/in silico method described in this study shows that the highest levels of leaflet stress occur in the regions of the commissures and mid-portion of the leaflet body. This method may have important insight with regard to bioprosthetic valve durability. Our results suggest that, compared to porcine bioprostheses, those with externally mounted pericardium have higher leaflet mechanical stress, which may translate into shorter durability.

Influence of aneurysmal aortic root geometry on mechanical stress to the aortic valve leaflet

AIMS

While mechanical stress caused by blood flow, e.g. wall shear stress (WSS), and related parameters, e.g. oscillatory shear index (OSI), are increasingly being recognized as key moderators of various cardiovascular diseases, studies on valves have been limited because of a lack of appropriate imaging modalities. We investigated the influence of aortic root geometry on WSS and OSI on the aortic valve (AV) leaflet.

METHODS AND RESULTS

We applied our novel approach of intraoperative epi-aortic echocardiogram to measure the haemodynamic parameters of WSS and OSI on the AV leaflet. Thirty-six patients were included, which included those who underwent valve-sparing aortic root replacement (VSARR) with no significant aortic regurgitation (n = 17) and coronary artery bypass graft (CABG) with normal AV (n = 19). At baseline, those who underwent VSARR had a higher systolic WSS (0.52 ± 0.12 vs. 0.32 ± 0.08 Pa, respectively, P < 0.001) and a higher OSI (0.37 ± 0.06 vs. 0.29 ± 0.04, respectively, P < 0.001) on the aortic side of the AV leaflet than those who underwent CABG. Multivariate regression analysis revealed that the size of the sinus of Valsalva had a significant association with WSS and OSI. Following VSARR, WSS and OSI values decreased significantly compared with the baseline values (WSS: 0.29 ± 0.12 Pa, P < 0.001; OSI: 0.26 ± 0.09, P < 0.001), and became comparable to the values in those who underwent CABG (WSS, P = 0.42; OSI, P = 0.15).

CONCLUSIONS

Mechanical stress on the AV gets altered in correlation with the size of the aortic root. An aneurysmal aortic root may expose the leaflet to abnormal fluid dynamics. The VSARR procedure appeared to reduce these abnormalities.

In Vitro Study of a Stentless Aortic Bioprosthesis Made of Bacterial Cellulose

PURPOSE

The paper present findings from an in vitro experimental study of a stentless human aortic bioprosthesis (HAB) made of bacterial cellulose (BC). Three variants of the basic model were designed and tested to identify the valve prosthesis with the best performance parameters. The modified models were made of BC, and the basic model of pericardium.

METHODS

Each model (named V1, V2 and V3) was implanted into a 90 mm porcine aorta. Effective Orifice Area (EOA), rapid valve opening time (RVOT) and rapid valve closing time (RVCT) were determined. The flow resistance of each bioprosthesis model during the simulated heart systole, i.e. for the mean differential pressure (ΔP) at the time of full valve opening was measured. All experimental specimens were exposed to a mean blood pressure (MBP) of 90.5 ± 2.3 mmHg.

RESULTS

The V3 model demonstrated the best performance. The index defining the maximum opening of the bioprosthesis during systole for models V1, V2 and V3 was 2.67 ± 0.59, 2.04 ± 0.23 and 2.85 ± 0.59 cm2, respectively. The mean flow rate through the V3 valve was 5.7 ± 1, 6.9 ± 0.7 and 8.9 ± 1.4 l/min for stroke volume (SV) of 65, 90 and 110 mL, respectively. The phase of immediate opening and closure for models V1, V2 and V3 was 8, 7 and 5% of the cycle duration, respectively. The mean flow resistance of the models was: 4.07 ± 2.1, 4.28 ± 2.51 and 5.6 ± 2.32 mmHg.

CONCLUSIONS

The V3 model of the aortic valve prosthesis is the most effective. In vivo tests using BC as a structural material for this model are recommended. The response time of the V3 model to changed work conditions is comparable to that of a healthy human heart. The model functions as an aortic valve prosthesis in in vitro conditions.

Fluid structure interaction modelling of aortic valve stenosis: Effects of valve calcification on coronary artery flow and aortic root hemodynamics

BACKGROUND AND OBJECTIVE

Coronary artery diseases and aortic valve stenosis are two of the main causes of mortality and morbidity worldwide. Stenosis of the aortic valve develops due to calcium deposition on the aortic valve leaflets during the cardiac cycle. Clinical investigations have demonstrated that aortic valve stenosis not only affects hemodynamic parameters inside the aortic root but also has a significant influence on the coronary artery hemodynamics and leads to the initiation of coronary artery disease. The aim of this study is to investigate the effect of calcification of the aortic valve on the variation of hemodynamic parameters in the aortic root and coronary arteries in order to find potential locations for initiation of the coronary stenoses.

METHODS

Fluid structure interaction modelling methodology was used to simulate aortic valve hemodynamics in the presence of coronary artery flow. A 2-D model of the aortic valve leaflets was developed in ANSYS Fluent based on the available echocardiography images in literature. The k-ω SST turbulence model was utilised to model the turbulent flow downstream of the leaflets.

RESULTS

The effects of calcification of the aortic valve on aortic root hemodynamics including transvalvular pressure gradient, valve orifice dimeter, vorticity magnitude in the sinuses and wall shear stress on the ventricularis and fibrosa layers of the leaflets were studied. Results revealed that the transvalvular pressure gradient increases from 792 Pa (∼ 6 mmHg) for a healthy aortic valve to 2885 Pa (∼ 22 mmHg) for a severely calcified one. Furthermore, the influence of the calcification of the aortic valve leaflets on the velocity profile and the wall shear stress in the coronary arteries was investigated and used for identification of potential locations of initiation of the coronary stenoses. Obtained results show that the maximum velocity inside the coronary arteries at early diastole decreases from 1 m/s for the healthy valve to 0.45 m/s for the severely calcified case.

CONCLUSIONS

Calcification significantly decreases the wall shear stress of the coronary arteries. This reduction in the wall shear stress can be a main reason for initiation of the coronary atherosclerosis process and eventually results in coronary stenoses.

Aortic annulus and the importance of annuloplasty

Dystrophic aortic insufficiency accounts for the majority of Western cases of aortic insufficiency and can be divided into the three phenotypes of isolated aortic insufficiency, dilated aortic root, and dilated ascending aorta. Each of these phenotypes is associated with a dilated annulus and/or sinotubular junction. Recent international guidelines recommend reimplantation or remodeling with aortic annuloplasty for valve-sparing root replacement, as well as consideration of aortic valve repair in cases of aortic insufficiency. A dilated aortic annulus is a major risk factor for failure of aortic valve repair procedures, indicating the need to address the annulus at the time of aortic valve or root repair. Calibrated annuloplasty should be performed at sub- and supravalular levels in order to restore the ratio of the sinotubular junction and annulus and be adapted according to the phenotype of the root and ascending aorta. Standardization of aortic valve repair techniques with use of a calibrated annuloplasty will improve dissemination of techniques and rate of aortic valve repair. Current medical evidence shows that aortic valve repair is safe, produces better quality of life, and reduces valve-related mortality compared to prosthetic valve replacement.

Aortic valve opening and closure: the clover dynamics

Background

Systolic aortic root expansion is reported to facilitate valve opening, but the precise dynamics remain unknown. A sonometric study with a high data sampling rate (200 to 800 Hz) was conducted in an acute ovine model to better understand the timing, mechanisms, and shape of aortic valve opening and closure.

Methods

Eighteen piezoelectric crystals were implanted in 8 sheep at each annular base, commissures, sinus of Valsalva, sinotubular junction, nodulus of Arantius, and ascending aorta (AA). Geometric changes were time related to pressures and flows.

Results

The aortic root was hemodynamically divided into left ventricular (LV) and aortic compartments situated, respectively, below and above the leaflets. During isovolumetric contraction (IVC), aortic root expansion started in the LV compartment, most likely due to volume redistribution in the LV outflow tract below the leaflets. This expansion initiated leaflet separation prior to ejection (2.1%±0.5% of total opening area). Aortic compartment expansion was delayed toward the end of IVC, likely related to volume redistribution above the leaflets due to accelerating aortic backflow toward the aortic valve and coronary flow reduction due to myocardial contraction. Maximum valve opening during the first third of ejection acquired a truncated cone shape [leaflet free edge area smaller than annular base area (-41.5%±5.5%)]. The distal orifice became clover shaped because the leaflet free edge area is larger than the commissural area by 16.3%±2.0%.

Conclusions

Aortic valve opening is initiated prior to ejection related to delicate balance between LV, aortic root, and coronary dynamics. It is clover shaped at maximum opening in systole. A better understanding of these mechanisms should stimulate more physiological surgical approaches of valve repair and replacement.

Rationale for aortic annuloplasty to standardise aortic valve repair

Available evidence shows that aortic valve repair reduces valve-related mortality and improves quality of life compared to prosthetic aortic valve replacement. One of the most important predictors of bicuspid and tricuspid aortic valve repair failure is the absence of treating a dilated aortic annulus greater than 25-28 mm. Competency of the aortic valve depends on multiple factors including the diameter of the annulus, sinotubular junction, valve cusps and commissures. Dystrophic aortic insufficiency (AI) is the commonest cause of AI in the Western world and is characterised by dilatation of the aortic annulus (≥25 mm), sinuses and/or sinotubular junction (≥30 mm). Depending on whether the sinuses of Valsalva and/or tubular ascending aorta are dilated, three phenotypes can be identified: dilated aortic root, dilated ascending aorta and isolated AI. All three phenotypes are associated with a dilated aortic annulus. Aortic annuloplasty reduces the dilated aortic annulus and improves the surface of coaptation, as in the case of mitral valve repair. In treating AI, it is also important to restore the physiological sinotubular junction/annulus ratio, which can be carried out with remodeling root repair + subvalvular annuloplasty (for dilated aortic root), tubular ascending aorta replacement + subvalvular annuloplasty (for dilated ascending aorta) and double sub- and supra-valvular annuloplasty (for isolated AI). Aortic annuloplasty is now considered an essential component of aortic valve repair and valve-sparing root surgery.

Flow Dynamics in the Aortic Arch and Its Effect on the Arterial Input Function in Cardiac Computed Tomography

The arterial input function (AIF)-time-density curve (TDC) of contrast at the coronary ostia-plays a central role in contrast enhanced computed tomography angiography (CTA). This study employs computational modeling in a patient-specific aorta to investigate mixing and dispersion of contrast in the aortic arch (AA) and to compare the TDCs in the coronary ostium and the descending aorta. Here, we examine the validity of the use of TDC in the descending aorta as a surrogate for the AIF. Computational fluid dynamics was used to study hemodynamics and contrast dispersion in a CTA-based patient model of the aorta. Variations in TDC between the aortic root, through the AA and at the descending aorta and the effect of flow patterns on contrast dispersion was studied via post-processing of the results. Simulations showed complex unsteady patterns of contrast mixing and dispersion in the AA that are driven by the pulsatile flow. However, despite the relatively long intra-aortic distance between the coronary ostia and the descending aorta, the TDCs at these two locations were similar in terms of rise-time and up-slope, and the time lag between the two TDCs was 0.19 seconds. TDC in the descending aorta is an accurate analog of the AIF. Methods that use quantitative metrics such as rise-time and slope of the AIF to estimate coronary flowrate and myocardial ischemia can continue with the current practice of using the TDC at the descending aorta as a surrogate for the AIF.

Validation and Extension of a Fluid-Structure Interaction Model of the Healthy Aortic Valve

Purpose

The understanding of the optimum function of the healthy aortic valve is essential in interpreting the effect of pathologies in the region, and in devising effective treatments to restore the physiological functions. Still, there is no consensus on the operating mechanism that regulates the valve opening and closing dynamics. The aim of this study is to develop a numerical model that can support a better comprehension of the valve function and serve as a reference to identify the changes produced by specific pathologies and treatments.

Methods

A numerical model was developed and adapted to accurately replicate the conditions of a previous in vitro investigation into aortic valve dynamics, performed by means of particle image velocimetry (PIV). The resulting velocity fields of the two analyses were qualitatively and quantitatively compared to validate the numerical model. In order to simulate more physiological operating conditions, this was then modified to overcome the main limitations of the experimental setup, such as the presence of a supporting stent and the non-physiological properties of the fluid and vessels.

Results

The velocity fields of the initial model resulted in good agreement with those obtained from the PIV, with similar flow structures and about 90% of the computed velocities after valve opening within the standard deviation of the equivalent velocity measurements of the in vitro model. Once the experimental limitations were removed from the model, the valve opening dynamics changed substantially, with the leaflets opening into the sinuses to a much greater extent, enlarging the effective orifice area by 11%, and reducing greatly the vortical structures previously observed in proximity of the Valsalva sinuses wall.

Conclusions

The study suggests a new operating mechanism for the healthy aortic valve leaflets considerably different from what reported in the literature to date and largely more efficient in terms of hydrodynamic performance. This work also confirms the crucial role that numerical approaches, complemented with experimental findings, can play in overcoming some of the limitations inherent in experimental techniques, supporting the full understanding of complex physiological phenomena.

Electronic supplementary material

The online version of this article (doi:10.1007/s13239-018-00391-1) contains supplementary material, which is available to authorized users.

Fluid-structure interaction of a pulsatile flow with an aortic valve model: A combined experimental and numerical study

The complex fluid-structure interaction problem associated with the flow of blood through a heart valve with flexible leaflets is investigated both experimentally and numerically. In the experimental test rig, a pulse duplicator generates a pulsatile flow through a biomimetic rigid aortic root where a model of aortic valve with polymer flexible leaflets is implanted. High-speed recordings of the leaflets motion and particle image velocimetry measurements were performed together to investigate the valve kinematics and the dynamics of the flow. Large eddy simulations of the same configuration, based on a variant of the immersed boundary method, are also presented. A massively parallel unstructured finite-volume flow solver is coupled with a finite-element solid mechanics solver to predict the fluid-structure interaction between the unsteady flow and the valve. Detailed analysis of the dynamics of opening and closure of the valve are conducted, showing a good quantitative agreement between the experiment and the simulation regarding the global behavior, in spite of some differences regarding the individual dynamics of the valve leaflets. A multicycle analysis (over more than 20 cycles) enables to characterize the generation of turbulence downstream of the valve, showing similar flow features between the experiment and the simulation. The flow transitions to turbulence after peak systole, when the flow starts to decelerate. Fluctuations are observed in the wake of the valve, with maximum amplitude observed at the commissure side of the aorta. Overall, a very promising experiment-vs-simulation comparison is shown, demonstrating the potential of the numerical method.

Development of an Off-the-Shelf Tissue-Engineered Sinus Valve for Transcatheter Pulmonary Valve Replacement: a Proof-of-Concept Study

Tissue-engineered heart valves with self-repair and regeneration properties may overcome the problem of long-term degeneration of currently used artificial prostheses. The aim of this study was the development and in vivo proof-of-concept of next-generation off-the-shelf tissue-engineered sinus valve (TESV) for transcatheter pulmonary valve replacement (TPVR). Transcatheter implantation of off-the-shelf TESVs was performed in a translational sheep model for up to 16 weeks. Transapical delivery of TESVs was successful and showed good acute and short-term performance (up to 8 weeks), which then worsened over time most likely due to a non-optimized in vitro valve design. Post-mortem analyses confirmed the remodelling potential of the TESVs, with host cell infiltration, polymer degradation, and collagen and elastin deposition. TESVs proved to be suitable as TPVR in a preclinical model, with encouraging short-term performance and remodelling potential. Future studies will enhance the clinical translation of such approach by improving the valve design to ensure long-term functionality.

INFLUENCE OF LEAFLET’S MATRIX STIFFNESS AND FIBER ORIENTATION ON THE OPENING DYNAMICS OF A PROSTHETIC TRILEAFLET HEART VALVE

Biological valves are employed for aortic valve substitution since a long time but there is a growing effort toward the development of new engineered tissues, in which the complex mechanical response of native leaflets is replicated using composite materials consisting of a soft matrix with embedded reinforcing fibers. The main goal of the present study is to investigate the influence that variations on fiber orientation and matrix stiffness may have on valve dynamics. To this aim, a fluid–structure interaction (FSI) model of a trileaflet valve was implemented in which the opening phase was simulated and leaflet matrix stiffness and fiber orientation were varied in the framework of an anisotropic hyperelastic strain energy function. Results show that both parameters may affect significantly transvalvular pressure gradient and effective orifice area (EOA). For the opening phase of the valve examined, less favorable flow conditions were found when preferred fiber orientation is circumferential, due to lower maximum EOA achievable. Such configuration in combination with stiffer matrix may result in significant degradation of valve performances. Overall fiber orientation can potentially be taylored to optimize valve dynamics, provided also structural aspects that may be prominent in the closure phase, are considered.

Fiber-reinforced computational model of the aortic root incorporating thoracic aorta and coronary structures

Cardiovascular diseases are still the leading causes of death in the developed world. The decline in the mortality associated with circulatory system diseases is accredited to development of new diagnostic and prognostic tools. It is well known that there is an inter relationship between the aortic valve impairment and pathologies of the aorta and coronary vessels. However, due to the limitations of the current tools, the possible link is not fully elucidated. Following our previous model of the aortic root including the coronaries, in this study, we have further developed the global aspect of the model by incorporating the anatomical structure of the thoracic aorta. This model is different from all the previous studies in the sense that inclusion of the coronary structures and thoracic aorta into the natural aortic valve introduces the notion of globality into the model enabling us to explore the possible link between the regional pathologies. The developed model was first validated using the available data in the literature under physiological conditions. Then, to provide a support for the possible association between the localized cardiovascular pathologies and global variations in hemodynamic conditions, we simulated the model for two pathological conditions including moderate and severe aortic valve stenoses. The findings revealed that malformations of the aortic valve are associated with development of low wall shear stress regions and helical blood flow in thoracic aorta that are considered major contributors to aortic pathologies.

A patient-specific aortic valve model based on moving resistive immersed implicit surfaces

In this paper, we propose a full computational framework to simulate the hemodynamics in the aorta including the valve. Closed and open valve surfaces, as well as the lumen aorta, are reconstructed directly from medical images using new ad hoc algorithms, allowing a patient-specific simulation. The fluid dynamics problem that accounts from the movement of the valve is solved by a new 3D-0D fluid-structure interaction model in which the valve surface is implicitly represented through level set functions, yielding, in the Navier-Stokes equations, a resistive penalization term enforcing the blood to adhere to the valve leaflets. The dynamics of the valve between its closed and open position is modeled using a reduced geometric 0D model. At the discrete level, a finite element formulation is used and the SUPG stabilization is extended to include the resistive term in the Navier-Stokes equations. Then, after time discretization, the 3D fluid and 0D valve models are coupled through a staggered approach. This computational framework, applied to a patient-specific geometry and data, allows to simulate the movement of the valve, the sharp pressure jump occurring across the leaflets, and the blood flow pattern inside the aorta.

3D physiological model of the aortic valve incorporating small coronary arteries

The diseases of the coronary arteries and the aortic root are still the leading causes of mortality and morbidity worldwide. In this study, a 3D global fluid-structure interaction of the aortic root with inclusion of anatomically inspired small coronary arteries using the finite element method is presented. This innovative model allows to study the impact and interaction of root biomechanics on coronary hemodynamics and brings a new understanding to small coronary vessels hemodynamics. For the first time, the velocity profiles and shear stresses are reported in distal coronary arteries as a result of the aortic flow conditions in a global fluid-structure interaction model.

Evaluation of an aortic valve prosthesis: Fluid-structure interaction or structural simulation?

Bio-inspired polymeric heart valves (PHVs) are excellent candidates to mimic the structural and the fluid dynamic features of the native valve. PHVs can be implanted as prosthetic alternative to currently clinically used mechanical and biological valves or as potential candidate for a minimally invasive treatment, like the transcatheter aortic valve implantation. Nevertheless, PHVs are not currently used for clinical applications due to their lack of reliability. In order to investigate the main features of this new class of prostheses, pulsatile tests in an in-house pulse duplicator were carried out and reproduced in silico with both structural Finite-Element (FE) and Fluid-Structure interaction (FSI) analyses. Valve kinematics and geometric orifice area (GOA) were evaluated to compare the in vitro and the in silico tests. Numerical results showed better similarity with experiments for the FSI than for the FE simulations. The maximum difference between experimental and FSI GOA at maximum opening time was only 5%, as compared to the 46.5% between experimental and structural FE GOA. The stress distribution on the valve leaflets clearly reflected the difference in valve kinematics. Higher stress values were found in the FSI simulations with respect to those obtained in the FE simulation. This study demonstrates that FSI simulations are more appropriate than FE simulations to describe the actual behaviour of PHVs as they can replicate the valve-fluid interaction while providing realistic fluid dynamic results.

The impact of the aortic valve impairment on the distant coronary arteries hemodynamics: a fluid-structure interaction study

Atherosclerosis is still the leading cause of death in the developed world. Although its initiation and progression is a complex multifactorial process, it is well known that blood flow-induced wall shear stress (WSS) is an important factor involved in early atherosclerotic plaque initiation. In recent clinical studies, it was established that the regional pathologies of the aortic valve can be involved in the formation of atherosclerotic plaques. However, the impact of hemodynamic effects is not yet fully elucidated for disease initiation and progression. In this study, our developed 3D global fluid-structure interaction model of the aortic root incorporating coronary arteries is used to investigate the possible interaction between coronary arteries and aortic valve pathologies. The coronary hemodynamics was examined and quantified for different degrees of aortic stenosis varying from nonexistent to severe. For the simulated healthy model, the calculated WSS varied between 0.41 and 1.34 Pa which is in the atheroprotective range. However, for moderate and severe aortic stenoses, wide regions of the coronary structures, especially the proximal sections around the first bifurcation, were exposed to lower values of WSS and therefore they were prone to atherosclerosis even in the case of healthy coronary arteries.

Fluid-Structure Interaction Study of Transcatheter Aortic Valve Dynamics Using Smoothed Particle Hydrodynamics

Computational modeling of heart valve dynamics incorporating both fluid dynamics and valve structural responses has been challenging. In this study, we developed a novel fully-coupled fluid-structure interaction (FSI) model using smoothed particle hydrodynamics (SPH). A previously developed nonlinear finite element (FE) model of transcatheter aortic valves (TAV) was utilized to couple with SPH to simulate valve leaflet dynamics throughout the entire cardiac cycle. Comparative simulations were performed to investigate the impact of using FE-only models vs. FSI models, as well as an isotropic vs. an anisotropic leaflet material model in TAV simulations. From the results, substantial differences in leaflet kinematics between FE-only and FSI models were observed, and the FSI model could capture the realistic leaflet dynamic deformation due to its more accurate spatial and temporal loading conditions imposed on the leaflets. The stress and the strain distributions were similar between the FE and FSI simulations. However, the peak stresses were different due to the water hammer effect induced by the fluid inertia in the FSI model during the closing phase, which led to 13-28% lower peak stresses in the FE-only model compared to that of the FSI model. The simulation results also indicated that tissue anisotropy had a minor impact on hemodynamics of the valve. However, a lower tissue stiffness in the radial direction of the leaflets could reduce the leaflet peak stress caused by the water hammer effect. It is hoped that the developed FSI models can serve as an effective tool to better assess valve dynamics and optimize next generation TAV designs.

Physiological vortices in the sinuses of Valsalva: An in vitro approach for bio-prosthetic valves

PURPOSE

The physiological flow dynamics within the Valsalva sinuses, in terms of global and local parameters, are still not fully understood. This study attempts to identify the physiological conditions as closely as possible, and to give an explanation of the different and sometime contradictory results in literature.

METHODS

An in vitro approach was implemented for testing porcine bio-prosthetic valves operating within different aortic root configurations. All tests were performed on a pulse duplicator, under physiological pressure and flow conditions. The fluid dynamics established in the various cases were analysed by means of 2D Particle Image Velocimetry, and related with the achieved hydrodynamic performance.

RESULTS

Each configuration is associated with substantially different flow dynamics, which significantly affects the valve performance. The configuration most closely replicating healthy native anatomy was characterised by the best hemodynamic performance, and any mismatch in size and position between the valve and the root produced substantial modification of the fluid dynamics downstream of the valve, hindering the hydrodynamic performance of the system. The worst conditions were observed for a configuration characterised by the total absence of the Valsalva sinuses.

CONCLUSION

This study provides an explanation for the different vortical structures described in the literature downstream of bioprosthetic valves, enlightening the experimental complications in valve testing. Most importantly, the results clearly identify the fluid mechanisms promoted by the Valsalva sinuses to enhance the ejection and closing phases, and this study exposes the importance of an optimal integration of the valve and root, to operate as a single system.

Multiphysics simulation of the effect of leaflet thickness inhomogeneity and material anisotropy on the stress-strain distribution on the aortic valve

This study developed a realistic 3D FSI computational model of the aortic valve using the fixed-grid method, which was eventually employed to investigate the effect of the leaflet thickness inhomogeneity and leaflet mechanical nonlinearity and anisotropy on the simulation results. The leaflet anisotropy and thickness inhomogeneity were found to significantly affect the valve stress-strain distribution. However, their effect on valve dynamics and fluid flow through the valve were minor. Comparison of the simulation results against in-vivo and in-vitro data indicated good agreement between the computational models and experimental data. The study highlighted the importance of simulating multi-physics phenomena (such as fluid flow and structural deformation), regional leaflet thickness inhomogeneity and anisotropic nonlinear mechanical properties, to accurately predict the stress-strain distribution on the natural aortic valve.

A Novel Method for Optical High Spatiotemporal Strain Analysis for Transcatheter Aortic Valves In Vitro

The transcatheter aortic valve implantation (TAVI) valve is a bioprosthetic valve within a metal stent frame. Like traditional surgical bioprosthetic valves, the TAVI valve leaflet tissue is expected to calcify and degrade over time. However, clinical studies of TAVI valve longevity are still limited. In order to indirectly assess the longevity of TAVI valves, an estimate of the mechanical wear and tear in terms of valvular deformation and strain of the leaflets under various conditions is warranted. The aim of this study was, therefore, to develop a platform for noncontact TAVI valve deformation analysis with both high temporal and spatial resolutions based on stereophotogrammetry and digital image correlation (DIC). A left-heart pulsatile in vitro flow loop system for mounting of TAVI valves was designed. The system enabled high-resolution imaging of all three TAVI valve leaflets simultaneously for up to 2000 frames per second through two high-speed cameras allowing three-dimensional analyses. A coating technique for applying a stochastic pattern on the leaflets of the TAVI valve was developed. The technique allowed a pattern recognition software to apply frame-by-frame cross correlation based deformation measurements from which the leaflet motions and the strain fields were derived. The spatiotemporal development of a very detailed strain field was obtained with a 0.5 ms time resolution and a spatial resolution of 72 μm/pixel. Hence, a platform offering a new and enhanced supplementary experimental evaluation of tissue valves during various conditions in vitro is presented.

Hemodynamic performance of Edwards Intuity valve in a compliant aortic root model

Numerous designs of bioprosthetic valves exist. The sutureless surgical valve is a newer design concept which combines elements of the transcatheter valve technology with surgical valves. This design aims at shorter and easier implantation. It was the aim of this study to perform hemodynamic and kinematic measurements for this type of valves to serve as a baseline for following studies which investigate the effect of the aortic root on the valve performance. To this end, the Edwards Intuity aortic valve was investigated in a new in vitro flow loop mimicking the left heart. The valve was implanted in a transparent, compliant aortic root model, and the valve kinematics was investigated using a high speed camera together with synchronized hemodynamic measurements of pressures and flows. The valve closure was asynchronous (one by one leaflet), and the valve started to close before the deceleration of the fluid. The aortic root model showed a dilation of the sinuses which was different to the ascending aorta, and the annulus was found to move towards the left ventricle during diastole and towards the aorta during systole.

Early results of extra-aortic annuloplasty ring implantation on aortic annular dimensions

OBJECTIVES

Dilatation of the aortic annulus is a cause of recurrent aortic regurgitation following the Ross or Yacoub procedures. Use of an extra-aortic annuloplasty ring is a potentially useful adjunct procedure. The aim of this study was to analyze the early effectiveness and mid-term stability of this surgical approach.

METHODS

From 2011 to 2015, 50 patients (mean age, 43 ± 14 years) underwent adjunct extra-aortic annuloplasty ring implantation (n = 39 Dacron rings and n = 11 ExAo rings [CORONEO Inc, Montreal, Canada]). Median ring size was 28 mm (range, 27-32 mm). All patients had aortic regurgitation or a dilated aortic annulus. Concomitant surgical procedure was a valve-sparing remodeling procedure (n = 32) or a Ross procedure (n = 18). Baseline and follow-up echocardiographic systolic and diastolic aortic annular dimensions were prospectively collected. Longitudinal analyses were performed using mixed-effect models. Median follow-up was 12 months (98% complete).

RESULTS

Use of an extra-aortic annuloplasty ring resulted in a significant decrease in both systolic (27.9 ± 0.5 mm preoperatively vs 23.6 ± 0.3 mm at discharge, P < .001) and diastolic (24.8 ± 0.4 mm preoperatively vs 20.3 ± 0.3 mm at discharge, P < .001) dimensions. Mean systolic and diastolic dimensions remained statistically unchanged up to 2 years postoperatively, compared with their predischarge values. Systolic expansion of the annulus was conserved early after surgery (16% systolic expansion) and preserved up to 2 years after ring implantation.

CONCLUSIONS

Use of an extra-aortic annuloplasty ring is effective at reducing annular diameters. This remains stable at mid-term follow-up, with preservation of aortic annular dynamics. Longer-term studies are required to determine the continued stability and impact on long-term clinical outcomes.

Aortic valve dynamics using a fluid structure interaction model--The physiology of opening and closing

Comparative study among aortic valves requires the use of an unbiased and relevant boundary condition. Pressure and flow boundary conditions used in literature are not sufficient for an unbiased analysis. We need a different boundary condition to analyze the valves in an unbiased, relevant environment. The proposed boundary condition is a combination of the pressure and flow boundary condition methods, which is chosen considering the demerits of the pressure and flow boundary conditions. In order to study the valve in its natural environment and to give a comparative analysis between different boundary conditions, a fluid-structure interaction analysis is made using the pressure and the proposed boundary conditions for a normal aortic valve. Commercial software LS-DYNA is used in all our analysis. The proposed boundary condition ensures a full opening of the valve with reduced valve regurgitation. It is found that for a very marginal raise in the ventricular pressure caused by pumping a fixed stroke volume, the cardiac output is considerably raised. The mechanics of the valve is similar between these two boundary conditions, however we observe that the importance of the root to raise the cardiac output may be overstated, considering the importance of the fully open nodule of arantius. Our proposed boundary condition delivers all the insights offered by the pressure and flow boundary conditions, along with providing an unbiased framework for the analysis of different valves and hence, more suitable for comparative analysis.

Computational modeling of cardiac valve function and intervention

In the past two decades, major advances have been made in the clinical evaluation and treatment of valvular heart disease owing to the advent of noninvasive cardiac imaging modalities. In clinical practice, valvular disease evaluation is typically performed on two-dimensional (2D) images, even though most imaging modalities offer three-dimensional (3D) volumetric, time-resolved data. Such 3D data offer researchers the possibility to reconstruct the 3D geometry of heart valves at a patient-specific level. When these data are integrated with computational models, native heart valve biomechanical function can be investigated, and preoperative planning tools can be developed. In this review, we outline the advances in valve geometry reconstruction, tissue property modeling, and loading and boundary definitions for the purpose of realistic computational structural analysis of cardiac valve function and intervention.

Experimental and numerical analysis of soft tissue stiffness measurement using manual indentation device--significance of indentation geometry and soft tissue thickness

BACKGROUND

Indentation techniques haves been applied to measure stiffness of human soft tissues. Tissue properties and geometry of the indentation instrument control the measured response.

METHODS

Mechanical roles of different soft tissues were characterized to understand the performance of the indentation instrument. An optimal instrument design was investigated. Experimental indentations in forearm of human subjects (N = 11) were conducted. Based on peripheral quantitative computed tomography imaging, a finite element (FE) model for indentation was created. The model response was matched with the experimental data.

RESULTS

Optimized values for the elastic modulus of skin and adipose tissue were 130.2 and 2.5 kPa, respectively. The simulated indentation response was 3.9 ± 1.2 (mean ± SD) and 4.9 ± 2.0 times more sensitive to changes in the elastic modulus of the skin than to changes in the elastic modulus of adipose tissue and muscle, respectively. Skin thickness affected sensitivity of the instrument to detect changes in stiffness of the underlying tissues.

CONCLUSION

Finite element modeling provides a feasible method to quantitatively evaluate the geometrical aspects and the sensitivity of an indentation measurement device. Systematically, the skin predominantly controlled the indentation response regardless of the indenter geometry or variations in the volume of different soft tissues.

In vitro characterization of an aortic bioprosthetic valve using Doppler echocardiography and qualitative flow visualization

A 19 mm diameter prototype bioprosthetic valve mounted in a cardiac pulse duplicator was characterized using Doppler echocardiography and qualitative flow visualization at a heart rate of 72 bpm. Analysis of the flow visualization images revealed that the prototype and control valve leaflets open symmetrically but close asymmetrically. The asymmetry in the closing of the valves is likely due to the large pressure gradients across the valves and may have implications for the long term mechanical failure of the valves. The relatively high peak systolic velocity of 309.9 cm/s, which was measured in the prototype 19 mm valve, can be attributed to the small valve diameter and the high cardiac output used in the current study.

Toward patient-specific simulations of cardiac valves: state-of-the-art and future directions

Recent computational methods enabling patient-specific simulations of native and prosthetic heart valves are reviewed. Emphasis is placed on two critical components of such methods: (1) anatomically realistic finite element models for simulating the structural dynamics of heart valves; and (2) fluid structure interaction methods for simulating the performance of heart valves in a patient-specific beating left ventricle. It is shown that the significant progress achieved in both fronts paves the way toward clinically relevant computational models that can simulate the performance of a range of heart valves, native and prosthetic, in a patient-specific left heart environment. The significant algorithmic and model validation challenges that need to be tackled in the future to realize this goal are also discussed.

INVESTIGATION OF LEAFLET GEOMETRY IN A PERCUTANEOUS AORTIC VALVE WITH THE USE OF FLUID-STRUCTURE INTERACTION SIMULATION

Prosthetic aortic valves have been used for the replacement of dysfunctional native aortic valves in humans for more than fifty years. Current prosthetic valves have significant limitations and the development of improved aortic valve prostheses remains an important research focus area. This paper investigates one of the newer additions to the family of replacement valves, namely the stented percutaneous valve. An important design aspect of stented percutaneous valves, is the configuration of the leaflet's attachment to the surrounding stent. There are essentially two possible configurations: The first method is attaching the leaflets in a straight configuration, and the second method is to attach the leaflets in a curved configuration. Finite element models of both configurations were created, and the behavior of these configurations was then studied using a fluid-structure interaction (FSI) simulation. The FSI simulation was validated by means of comparing simulation results to actual measurements from a pulse duplicator using prototype valves of both configurations. The FSI results showed no significant difference between the valves' opening and closing behaviors. The von Mises stress distributions proved to be the largest differentiating and decisive factor between the two valves. The FSI simulations did however show that the leaflets that are attached in the straight configuration form folds that resembles that of the curved configuration as well as the native valve, but to a larger scale. The effect that these folds might have on valve tissue fatigue leaves room for future investigation.

Therapeutic vascular compliance change may cause significant variation in coronary perfusion: a numerical study

In some pathological conditions like aortic stiffening and calcific aortic stenosis (CAS), the microstructure of the aortic root and the aortic valve leaflets are altered in response to stress resulting in changes in tissue thickness, stiffness, or both. This aortic stiffening and CAS are thought to affect coronary blood flow. The goal of the present paper was to include the flow in the coronary ostia in the previous fluid structure interaction model we have developed and to analyze the effect of diseased tissues (aortic root stiffening and CAS) on coronary perfusion. Results revealed a significant impact on the coronary perfusion due to a moderate increase in the aortic wall stiffness and CAS (increase of the aortic valve leaflets thickness). A marked drop of coronary peak velocity occurred when the values of leaflet thickness and aortic wall stiffness were above a certain threshold, corresponding to a threefold of their normal value. Consequently, mild and prophylactic treatments such as smoking cessation, exercise, or diet, which have been proven to increase the aortic compliance, may significantly improve the coronary perfusion.

Aortic valve leaflet mechanical properties facilitate diastolic valve function

This work was concerned with the numerical simulation of the behaviour of aortic valves whose material can be modelled as non-linear elastic anisotropic. Linear elastic models for the valve leaflets with parameters used in previous studies were compared with hyperelastic models, incorporating leaflet anisotropy with pronounced stiffness in the circumferential direction through a transverse isotropic model. The parameters for the hyperelastic models were obtained from fits to results of orthogonal uniaxial tensile tests on porcine aortic valve leaflets. The computational results indicated the significant impact of transverse isotropy and hyperelastic effects on leaflet mechanics; in particular, increased coaptation with peak values of stress and strain in the elastic limit. The alignment of maximum principal stresses in all models follows approximately the coarse collagen fibre distribution found in aortic valve leaflets. The non-linear elastic leaflets also demonstrated more evenly distributed stress and strain which appears relevant to long-term scaffold stability and mechanotransduction.

Structural analysis of the natural aortic valve in dynamics: from unpressurized to physiologically loaded

A novel finite element model of the natural aortic valve was developed implementing anisotropic hyperelastic material properties for the leaflets and aortic tissues, and starting from the unpressurized geometry. Static pressurization of the aortic root, silicone rubber moulds and published data helped to establish the model parameters, while high-speed video recording of the leaflet motion in a left-heart simulator allowed for comparisons with simulations. The model was discretized with brick elements and loaded with time-varying pressure using an explicit commercial solver. The aortic valve model produced a competent valve whose dynamic behavior (geometric orifice area vs. time) closely matched that observed in the experiment. In both cases, the aortic valve took approximately 30 ms to open to an 800 mm(2) orifice and remained completely or more than half open for almost 200 ms, after which it closed within 30-50 ms. The highest values of stress were along the leaflet attachment line and near the commissure during diastole. Von Mises stress in the leaflet belly reached 600-750 kPa from early to mid-diastole. While the model using the unpressurized geometry as initial configuration was specially designed to satisfy the requirements of continuum mechanics for large deformations of hyperelastic materials, it also clearly demonstrated that dry models can be adequate to analyze valve dynamics. Although improvements are still needed, the advanced modeling and validation techniques used herein contribute toward improved and quantified accuracy over earlier simplified models.

Computational analysis of an aortic valve jet with Lagrangian coherent structures

Important progress has been achieved in recent years in simulating the fluid-structure interaction around cardiac valves. An important step in making these computational tools useful to clinical practice is the development of postprocessing techniques to extract clinically relevant information from these simulations. This work focuses on flow through the aortic valve and illustrates how the computation of Lagrangian coherent structures can be used to improve insight into the transport mechanics of the flow downstream of the valve, toward the goal of aiding clinical decision making and the understanding of pathophysiology.