Assessment of calcified aortic valve leaflet deformations and blood flow dynamics using fluid-structure interaction modeling

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.

Multiscale computational modeling of aortic valve calcification

Calcific aortic valve disease (CAVD) is a common cardiovascular disease that affects millions of people worldwide. The disease is characterized by the formation of calcium nodules on the aortic valve leaflets, which can lead to stenosis and heart failure if left untreated. The pathogenesis of CAVD is still not well understood, but involves several signaling pathways, including the transforming growth factor beta (TGF β ) pathway. In this study, we developed a multiscale computational model for TGF β -stimulated CAVD. The model framework comprises cellular behavior dynamics, subcellular signaling pathways, and tissue-level diffusion fields of pertinent chemical species, where information is shared among different scales. Processes such as endothelial to mesenchymal transition (EndMT), fibrosis, and calcification are incorporated. The results indicate that the majority of myofibroblasts and osteoblast-like cells ultimately die due to lack of nutrients as they become trapped in areas with higher levels of fibrosis or calcification, and they subsequently act as sources for calcium nodules, which contribute to a polydispersed nodule size distribution. Additionally, fibrosis and calcification processes occur more frequently in regions closer to the endothelial layer where the cell activity is higher. Our results provide insights into the mechanisms of CAVD and TGF β signaling and could aid in the development of novel therapeutic approaches for CAVD and other related diseases such as cancer. More broadly, this type of modeling framework can pave the way for unraveling the complexity of biological systems by incorporating several signaling pathways in subcellular models to simulate tissue remodeling in diseases involving cellular mechanobiology.

Patient-specific analysis of bicuspid aortic valve hemodynamics using a fully coupled fluid-structure interaction model

Bicuspid aortic valve (BAV), the most common congenital heart disease, is prone to develop significant valvular dysfunction and aortic wall abnormalities such as ascending aortic aneurysm. Growing evidence has suggested that abnormal BAV hemodynamics could contribute to disease progression. In order to investigate BAV hemodynamics, we performed 3D patient-specific fluid-structure interaction (FSI) simulations with fully coupled blood flow dynamics and valve motion throughout the cardiac cycle. Results showed that the hemodynamics during systole can be characterized by a systolic jet and two counter-rotating recirculation vortices. At peak systole, the jet was usually eccentric, with asymmetric recirculation vortices and helical flow motion in the ascending aorta. The flow structure at peak systole was quantified using the vorticity, flow rate reversal ratio and local normalized helicity (LNH) at four locations from the aortic root to the ascending aorta. The systolic jet was evaluated with the peak velocity, normalized flow displacement, and jet angle. It was found that peak velocity and normalized flow displacement (rather than jet angle) gave a strong correlation with the vorticity and LNH in the ascending aorta, which suggests that these two metrics could be used for clinical noninvasive evaluation of abnormal blood flow patterns in BAV patients.

Validation of fluid-structure interaction simulations of the opening phase of phantom mitral heart valves under physiologically inspired conditions

BACKGROUND AND OBJECTIVE

Atrioventricular valve disease is a common cause of heart failure, and successful surgical or interventional outcomes are crucial. Patient-specific fluid-structure interaction (FSI) modeling may provide valuable insights into valve dynamics and guidance of valve repair strategies. However, lack of validation has kept FSI modeling from clinical implementation. Therefore, this study aims to validate FSI simulations against in vitro benchmarking data, based on clinically relevant parameters for evaluating heart valve disease.

METHODS

An FSI model that mimics the left heart was developed. The domain included a deformable mitral valve of different stiffnesses run with different inlet velocities. Five different cases were simulated and compared to in vitro data based on the pressure difference across the valve, the valve opening, and the velocity in the flow domain.

RESULTS

The simulations underestimate the pressure difference across the valve by 6.8-14 % compared to catheter measurements. Evaluation of the valve opening showed an underprediction of 5.4-7.3 % when compared to cine MRI, 2D Echo, and 3D Echo data. Additionally, the simulated velocity through the valve showed a 7.9-8.4 % underprediction in relation to Doppler Echo measurements. Qualitative assessment of the velocity profile in the ventricle and the streamlines of the flow in the domain showed good agreement of the flow behavior.

CONCLUSIONS

Parameters relevant to the diagnosis of heart valve disease estimated by FSI simulations showed good agreement when compared to in vitro benchmarking data, with differences small enough not to affect the grading of heart valve disease. The FSI model is thus deemed good enough for further development toward patient-specific cases.

Biomechanical Evaluation of Aortic Valve Stenosis by Means of a Virtual Stress Test: A Fluid-Structure Interaction Study

The impact of aortic valve stenosis (AS) extends beyond the vicinity of the narrowed leaflets into the left ventricle (LV) and into the systemic vasculature because of highly unpredictable valve behavior and complex blood flow in the ascending aorta that can be attributed to the strong interaction between the narrowed cusps and the ejected blood. These effects can become exacerbated during exercise and may have implications for disease progression, accurate diagnosis, and timing of intervention. In this 3-D patient-specific study, we employ strongly coupled fluid-structure interaction (FSI) modeling to perform a comprehensive biomechanical evaluation of systolic ejection dynamics in a stenosed aortic valve (AV) during increasing LV contraction. Our model predictions reveal that the heterogeneous ∆P vs. Q relationship that was observed in our previous clinical study can be attributed to a non-linear increase (by ~ 1.5-fold) in aortic valve area as LV heart rate increases from 70 to 115 bpm. Furthermore, our results show that even for a moderately stenotic valve, increased LV contraction during exercise can lead to high-velocity flow turbulence (Re = 11,700) in the aorta similar to that encountered with a severely stenotic valve (Re ~ 10,000), with concomitant greater viscous loss (~3-fold increase) and elevated wall stress in the ascending aorta. Our FSI predictions also reveal that individual valve cusps undergo distinct and highly non-linear increases (>100%) in stress during exercise, potentially contributing to progressive calcification. Such quantitative biomechanical evaluations from realistic FSI workflows provide insights into disease progression and can be integrated with current stress testing for AS patients to comprehensively predict hemodynamics and valve function under both baseline and exercise conditions.

A Doppler-exclusive non-invasive computational diagnostic framework for personalized transcatheter aortic valve replacement

Given the associated risks with transcatheter aortic valve replacement (TAVR), it is crucial to determine how the implant will affect the valve dynamics and cardiac function, and if TAVR will improve or worsen the outcome of the patient. Effective treatment strategies, indeed, rely heavily on the complete understanding of the valve dynamics. We developed an innovative Doppler-exclusive non-invasive computational framework that can function as a diagnostic tool to assess valve dynamics in patients with aortic stenosis in both pre- and post-TAVR status. Clinical Doppler pressure was reduced by TAVR (52.2 ± 20.4 vs. 17.3 ± 13.8 [mmHg], p < 0.001), but it was not always accompanied by improvements in valve dynamics and left ventricle (LV) hemodynamics metrics. TAVR had no effect on LV workload in 4 patients, and LV workload post-TAVR significantly rose in 4 other patients. Despite the group level improvements in maximum LV pressure (166.4 ± 32.2 vs 131.4 ± 16.9 [mmHg], p < 0.05), only 5 of the 12 patients (41%) had a decrease in LV pressure. Moreover, TAVR did not always improve valve dynamics. TAVR did not necessarily result in a decrease (in 9 out of 12 patients investigated in this study) in major principal stress on the aortic valve leaflets which is one of the main contributors in valve degeneration and, consequently, failure of heart valves. Diastolic stresses increased significantly post-TAVR (34%, 109% and 81%, p < 0.001) for each left, right and non-coronary leaflets respectively. Moreover, we quantified the stiffness and material properties of aortic valve leaflets which correspond with the reduced calcified region average stiffness among leaflets (66%, 74% and 62%; p < 0.001; N = 12). Valve dynamics post-intervention should be quantified and monitored to ensure the improvement of patient conditions and prevent any further complications. Improper evaluation of biomechanical valve features pre-intervention as well as post-intervention may result in harmful effects post-TAVR in patients including paravalvular leaks, valve degeneration, failure of TAVR and heart failure.

Multiscale model for blood flow after a bileaflet artificial aortic valve implantation

Cardiovascular diseases are the leading cause of mortality in the world, mainly due to atherosclerosis and its consequences. The article presents the numerical model of the blood flow through artificial aortic valve. The overset mesh approach was applied to simulate the valve leaflets motion and to realize the moving mesh, in the aortic arch and the main branches of cardiovascular system. To capture the cardiac system's response and the effect of vessel compliance on the outlet pressure, the lumped parameter model has been also included within the solution procedure. Three different turbulence modeling approaches were used and compared - the laminar, k-ϵ and k-ω model. The simulation results were also compared with the model excluding the moving valve geometry and the importance of the lumped parameter model for the outlet boundary condition was analyzed. Proposed numerical model and protocol was found as suitable for performing the virtual operations on the real patient vasculature geometry. The time-efficient turbulence model and overall solving procedure allows to support the clinicians in making decisions about the patient treatment and to predict the results of the future surgery.

An ultrasound-exclusive non-invasive computational diagnostic framework for personalized cardiology of aortic valve stenosis

Aortic stenosis (AS) is an acute and chronic cardiovascular disease and If left untreated, 50% of these patients will die within two years of developing symptoms. AS is characterized as the stiffening of the aortic valve leaflets which restricts their motion and prevents the proper opening under transvalvular pressure. Assessments of the valve dynamics, if available, would provide valuable information about the patient's state of cardiac deterioration as well as heart recovery and can have incredible impacts on patient care, planning interventions and making critical clinical decisions with life-threatening risks. Despite remarkable advancements in medical imaging, there are no clinical tools available to quantify valve dynamics invasively or noninvasively. In this study, we developed a highly innovative ultrasound-based non-invasive computational framework that can function as a diagnostic tool to assess valve dynamics (e.g. transient 3-D distribution of stress and displacement, 3-D deformed shape of leaflets, geometric orifice area and angular positions of leaflets) for patients with AS at no risk to the patients. Such a diagnostic tool considers the local valve dynamics and the global circulatory system to provide a platform for testing the intervention scenarios and evaluating their effects. We used clinical data of 12 patients with AS not only to validate the proposed framework but also to demonstrate its diagnostic abilities by providing novel analyses and interpretations of clinical data in both pre and post intervention states. We used transthoracic echocardiogram (TTE) data for the developments and transesophageal echocardiography (TEE) data for validation.

Importance of Non-Newtonian Computational Fluid Modeling on Severely Calcified Aortic Valve Geometries - Insights from Quasi-Steady State Simulations

The Newtonian model has commonly been used to represent the viscosity of blood in the aorta, despite blood itself being a non-Newtonian fluid. This is justified where shear rates tend to be large. However, we hypothesized that using the Newtonian model to predict the hemodynamics on the aortic valve, particularly in those with severe calcifications, is inaccurate owing to valve leaflet geometry irregularities inducing multiple regions of low shear rates, &lt;100 s-1, where a Newtonian model is invalid. We investigated the use of three fluid models: Newtonian, Carreau, and Quemada on a healthy and calcified valve and compared their ability to capture important hemodynamic parameters of time averaged wall shear stress (TAWSS) and the oscillatory shear index (OSI). Our findings indicate that when the shear rates were large enough, &gt;100 s-1, use of a Newtonian model was justified for the ventricularis side of the leaflet. However, in spatial regions of low shear rates, &lt;100 s-1, specifically on the inner cusps of the fibrosa side of the valve, TAWSS calculations under a Newtonian model were found to be much lower when compared with their non-Newtonian Carreau and Quemada counterparts. We hereby conclude that in order to facilitate more accurate computational flow simulations in severe aortic valve calcification which is subjected to relatively large spatial regions of low shear (&lt;100 s-1), a non-Newtonian model should be applied.

Improving transcatheter aortic valve interventional predictability via fluid-structure interaction modelling using patient-specific anatomy

Transcatheter aortic valve replacement (TAVR) is now a standard treatment for high-surgical-risk patients with severe aortic valve stenosis. TAVR is being explored for broader indications including degenerated bioprosthetic valves, bicuspid valves and for aortic valve (AV) insufficiency. It is, however, challenging to predict whether the chosen valve size, design or its orientation would produce the most-optimal haemodynamics in the patient. Here, we present a novel patient-specific evaluation framework to realistically predict the patient's AV performance with a high-fidelity fluid-structure interaction analysis that included the patient's left ventricle and ascending aorta (AAo). We retrospectively evaluated the pre- and post-TAVR dynamics of a patient who underwent a 23 mm TAVR and evaluated against the patient's virtually de-calcified AV serving as a hypothetical benchmark. Our model predictions were consistent with clinical data. Stenosed AV produced a turbulent flow during peak-systole, while aortic flow with TAVR and de-calcified AV were both in the laminar-to-turbulent transitional regime with an estimated fivefold reduction in viscous dissipation. For TAVR, dissipation was highest during early systole when valve deformation was the greatest, suggesting that an efficient valve opening may reduce energy loss. Our study demonstrates that such patient-specific modelling frameworks can be used to improve predictability and in the planning of AV interventions.

Impact of calcific aortic valve disease on valve mechanics

The aortic valve is a highly dynamic structure characterized by a transvalvular flow that is unsteady, pulsatile, and characterized by episodes of forward and reverse flow patterns. Calcific aortic valve disease (CAVD) resulting in compromised valve function and increased pressure overload on the ventricle potentially leading to heart failure if untreated, is the most predominant valve disease. CAVD is a multi-factorial disease involving molecular, tissue and mechanical interactions. In this review, we aim at recapitulating the biomechanical loads on the aortic valve, summarizing the current and most recent research in the field in vitro, in-silico, and in vivo, and offering a clinical perspective on current strategies adopted to mitigate or approach CAVD.

Numerical Investigation of the Fetal Left Heart Hemodynamics During Gestational Stages

Flow-driven hemodynamic forces on the cardiac tissues have critical importance, and have a significant role in the proper development of the heart. These mechanobiological mechanisms govern the cellular responses for the growth and remodeling of the heart, where the altered hemodynamic environment is believed to be a major factor that is leading to congenital heart defects (CHDs). In order to investigate the mechanobiological development of the normal and diseased hearts, identification of the blood flow patterns and wall shear stresses (WSS) on these tissues are required for an accurate hemodynamic assessment. In this study, we focus on the left heart hemodynamics of the human fetuses throughout the gestational stages. Computational fetal left heart models are created for the healthy fetuses using the ultrasound images at various gestational weeks. Realistic inflow boundary conditions are implemented in the models using the Doppler ultrasound measurements for resolving the specific blood flow waveforms in the mitral valve. Obtained results indicate that WSS and vorticity levels in the fetal left heart decrease with the development of the fetus. The maximum WSS around the mitral valve is determined around 36 Pa at the gestational week of 16. This maximum WSS decreases to 11 Pa at the gestational week of 27, indicating nearly three-times reduction in the peak shear stress. These findings reveal the highly dynamic nature of the left heart hemodynamics throughout the development of the human fetus and shed light into the relevance of hemodynamic environment and development of CHDs.

Computational Analysis of Wall Shear Stress Patterns on Calcified and Bicuspid Aortic Valves: Focus on Radial and Coaptation Patterns

Calcification and bicuspid valve formation are important aortic valve disorders that disturb the hemodynamics and the valve function. The detailed analysis of aortic valve hemodynamics would lead to a better understanding of the disease’s etiology. We computationally modeled the aortic valve using simplified three-dimensional geometry and inlet velocity conditions obtained via echocardiography. We examined various calcification severities and bicuspid valve formation. Fluid-structure interaction (FSI) analyses were adapted using ANSYS Workbench to incorporate both flow dynamics and leaflet deformation accurately. Simulation results were validated by comparing leaflet movements in B-mode echo recordings. Results indicate that the biomechanical environment is significantly changed for calcified and bicuspid valves. High flow jet velocities are observed in the calcified valves which results in high transvalvular pressure difference (TPG). Wall shear stresses (WSS) increased with the calcification on both fibrosa (aorta side) and ventricularis (left ventricle side) surfaces of the leaflet. The WSS distribution is regular on the ventricularis, as the WSS values proportionally increase from the base to the tip of the leaflet. However, WSS patterns are spatially complex on the fibrosa side. Low WSS levels and spatially complex WSS patterns on the fibrosa side are considered as promoting factors for further calcification and valvular diseases.

Fluid Flow Characteristics of Healthy and Calcified Aortic Valves Using Three-Dimensional Lagrangian Coherent Structures Analysis

Aortic valve calcification is an important cardiovascular disorder that deteriorates the accurate functioning of the valve leaflets. The increasing stiffness due to the calcification prevents the complete closure of the valve and therefore leads to significant hemodynamic alterations. Computational fluid dynamics (CFD) modeling enables the investigation of the entire flow domain by processing medical images from aortic valve patients. In this study, we computationally modeled and simulated a 3D aortic valve using patient-specific dimensions of the aortic root and aortic sinus. Leaflet stiffness is deteriorated in aortic valve disease due to calcification. In order to investigate the influence of leaflet calcification on flow dynamics, three different leaflet-stiffness values were considered for healthy, mildly calcified, and severely calcified leaflets. Time-dependent CFD results were used for applying the Lagrangian coherent structures (LCS) technique by performing finite-time Lyapunov exponent (FTLE) computations along with Lagrangian particle residence time (PRT) analysis to identify unique vortex structures at the front and backside of the leaflets. Obtained results indicated that the peak flow velocity at the valve orifice increased with the calcification rate. For the healthy aortic valve, a low-pressure field was observed at the leaflet tips. This low-pressure field gradually expanded through the entire aortic sinus as the calcification level increased. FTLE field plots of the healthy and calcified valves showed a variety of differences in terms of flow structures. When the number of fluid particles in the healthy valve model was taken as reference, 1.59 and 1.74 times more particles accumulated in the mildly and severely calcified valves, respectively, indicating that the calcified valves were not sufficiently opened to allow normal mass flow rates.

Modeling of aortic valve stenosis using fluid-structure interaction method

BACKGROUND AND OBJECTIVES

Fluid structure interaction (FSI) is defined as interaction of the structures with contacting fluids. The aortic valve experiences the interaction with blood flow in systolic phase. In this study, we have tried to predict the hemodynamics of blood flow through a normal and stenotic aortic valve in two relaxation and exercise conditions using a three-dimensional FSI method.

METHODS

The aorta valve was modeled as a three-dimensional geometry including a normal model and two others with 25% and 50% stenosis. The geometry of the aortic valve was extracted from CT images and the models were generated by MMIMCS software and then they were implemented in ANSYS software. The pulsatile flow rate was used for all cases and the numerical simulations were conducted based on a time-dependent domain.

RESULTS

The obtained results including the velocity, pressure, and shear stress contours in different systolic time sequences were explained and discussed. The maximum blood flow velocity in relaxation phase was obtained 1.62 m/s (normal valve), 3.78 m/s (25% stenosed valve), and 4.73 m/s (50% stenosed valve). In exercise condition, the maximum velocities are 2.86, 4.32, and 5.42 m/s respectively. The maximum blood pressure in relaxation phase was calculated 111.45 mmHg (normal), 148.66 mmHg (25% stenosed), and 164.21 mmHg (50% stenosed). However, the calculated values in exercise situation were 129.57, 163.58, and 191.26 mmHg.The validation of the predicted results was also conducted using existing literature.

CONCLUSIONS

We believe that such model are useful tools for biomechanical experts. The further studies should be done using experimental data and the data are implemented on the boundary conditions for better comparison of the results.

Computational Modeling of Blood Flow Hemodynamics for Biomechanical Investigation of Cardiac Development and Disease

The heart is the first functional organ in a developing embryo. Cardiac development continues throughout developmental stages while the heart goes through a serious of drastic morphological changes. Previous animal experiments as well as clinical observations showed that disturbed hemodynamics interfere with the development of the heart and leads to the formation of a variety of defects in heart valves, heart chambers, and blood vessels, suggesting that hemodynamics is a governing factor for cardiogenesis, and disturbed hemodynamics is an important source of congenital heart defects. Therefore, there is an interest to image and quantify the flowing blood through a developing heart. Flow measurement in embryonic fetal heart can be performed using advanced techniques such as magnetic resonance imaging (MRI) or echocardiography. Computational fluid dynamics (CFD) modeling is another approach especially useful when the other imaging modalities are not available and in-depth flow assessment is needed. The approach is based on numerically solving relevant physical equations to approximate the flow hemodynamics and tissue behavior. This approach is becoming widely adapted to simulate cardiac flows during the embryonic development. While there are few studies for human fetal cardiac flows, many groups used zebrafish and chicken embryos as useful models for elucidating normal and diseased cardiogenesis. In this paper, we explain the major steps to generate CFD models for simulating cardiac hemodynamics in vivo and summarize the latest findings on chicken and zebrafish embryos as well as human fetal hearts.

Computational Assessment of Valvular Dysfunction in Discrete Subaortic Stenosis: A Parametric Study

PURPOSE

Discrete subaortic stenosis (DSS) is a left-ventricular outflow tract (LVOT) obstruction caused by a membranous lesion. DSS is associated with steep aortoseptal angles (AoSAs) and is a risk factor for aortic regurgitation (AR). However, the etiology of AR secondary to DSS remains unknown. This study aimed at quantifying computationally the impact of AoSA steepening and DSS on aortic valve (AV) hemodynamics and AR.

METHODS

An LV geometry reconstructed from cine-MRI data was connected to an AV geometry to generate a unified 2D LV-AV model. Six geometrical variants were considered: unobstructed (CTRL) and DSS-obstructed LVOT (DSS), each reflecting three AoSA variations (110°, 120°, 130°). Fluid-structure interaction simulations were run to compute LVOT flow, AV leaflet dynamics, and regurgitant fraction (RF).

RESULTS

AoSA steepening and DSS generated vortex dynamics alterations and stenotic flow conditions. While the CTRL-110° model generated the highest degree of leaflet opening asymmetry, DSS preferentially altered superior leaflet kinematics, and caused leaflet-dependent alterations in systolic fluttering. LVOT steepening and DSS subjected the leaflets to increasing WSS overloads (up to 94% increase in temporal shear magnitude), while DSS also increased WSS bidirectionality on the inferior leaflet belly (+ 0.30-point in oscillatory shear index). Although AoSA steepening and DSS increased diastolic transvalvular backflow, regurgitant fractions (RF < 7%) remained below the threshold defining clinical mild AR.

CONCLUSIONS

The mechanical interactions between AV leaflets and LVOT steepening/DSS hemodynamic derangements do not cause AR. However, the leaflet WSS abnormalities predicted in those anatomies provide new support to a mechanobiological etiology of AR secondary to DSS.

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.

Using Discrete Multiphysics Modelling to Assess the Effect of Calcification on Hemodynamic and Mechanical Deformation of Aortic Valve

This study proposes a 3D particle-based (discrete) multiphysics approach for modelling calcification in the aortic valve. Different stages of calcification (from mild to severe) were simulated, and their effects on the cardiac output were assessed. The cardiac flow rate decreases with the level of calcification. In particular, there is a critical level of calcification below which the flow rate decreases dramatically. Mechanical stress on the membrane is also calculated. The results show that, as calcification progresses, spots of high mechanical stress appear. Firstly, they concentrate in the regions connecting two leaflets; when severe calcification is reached, then they extend to the area at the basis of the valve.

Pathologic Shear and Elongation Rates Do Not Cause Cleavage of Von Willebrand Factor by ADAMTS13 in a Purified System

INTRODUCTION

Pathological flows in patients with severe aortic stenosis are associated with acquired von Willebrand syndrome. This syndrome is characterized by excessive cleavage of von Willebrand factor by its main protease, A Disintegrin and Metalloproteinase with a Thrombospondin Type 1 Motif, Member 13 (ADAMTS13) leading to decreased VWF function and mucocutaneous bleeding. Aortic valve replacement and correction of the flow behavior to physiological levels reverses the syndrome, supporting the association between pathological flow and acquired von Willebrand syndrome. We investigated the effects of shear and elongational rates on von Willebrand factor cleavage in the presence of ADAMTS13.

METHODS

We identified acquired von Willebrand syndrome in five patients with severe aortic stenosis. Doppler echography values from these patients were used to develop three computational fluid dynamic (CFD) aortic valve models (normal, mild and severe stenosis). Shear, elongational rates and exposure times identified in the CFD simulations were used as parameters for the design of microfluidic devices to test the effects of pathologic shear and elongational rates on the structure and function of von Willebrand factor.

RESULTS

The shear rates (0-10,000s-1), elongational rates (0-1000 s-1) and exposure times (1-180 ms) tested in our microfluidic designs mimicked the flow features identified in patients with aortic stenosis. The shear and elongational rates tested in vitro did not lead to excessive cleavage or decreased function of von Willebrand factor in the presence of the protease.

CONCLUSIONS

High shear and elongational rates in the presence of ADAMTS13 are not sufficient for excessive cleavage of von Willebrand Factor.

Foam-in-vein: A review of rheological properties and characterization methods for optimization of sclerosing foams

Varicose veins are chronic venous defects that affect >20% of the population in developed countries. Among potential treatments, sclerotherapy is one of the most commonly used. It involves endovenous injection of a surfactant solution (or foam) in varicose veins, inducing damage to the endothelial layer and subsequent vessel sclerosis. Treatments have proven to be effective in the short-term, however recurrence is reported at rates of up to 64% 5-year post-treatment. Thus, once diagnosed with varicosities there is a high probability of a permanently reduced quality of life. Recently, foam sclerotherapy has become increasingly popular over its liquid counterpart, since foams can treat larger and longer varicosities more effectively, they can be imaged using ultrasound, and require lower amounts of sclerosing agent. In order to minimize recurrence rates however, an investigation of current treatment methods should lead to more effective and long-lasting effects. The literature is populated with studies aimed at characterizing the fundamental physics of aqueous foams; nevertheless, there is a significant need for appropriate product development platforms. Despite successfully capturing the microstructural evolution of aqueous foams, the complexity of current models renders them inadequate for pharmaceutical development. This review article will focus on the physics of foams and the attempts at optimizing them for sclerotherapy. This takes the form of a discussion of the most recent numerical and experimental models, as well as an overview of clinically relevant parameters. This holistic approach could contribute to better foam characterization methods that patients may eventually derive long term benefit from.

Effect of Aortic Wall Deformation with Healthy and Calcified Annulus on Hemodynamic Performance of Implanted On-X Valve

INTRODUCTION

In this research, the hemodynamic performance of a 23-mm On-X bileaflet mechanical heart valve (BMHV) was investigated with the realistic geometry model of the valve and the deformable aorta in accelerating systole. In addition, the effect of ascending aorta flexibility and aortic annulus calcification on the complex blood flow characteristics were investigated.

METHODS

The geometry of the aorta is derived from the medical images, and the Ogden model has been utilized for the mechanical behavior of the ascending aorta. The 3D numerical simulation by a two-way Fluid-Structure Interaction (FSI) analysis using the Arbitrary Lagrangian-Eulerian (ALE) method was performed throughout the accelerating systolic phase.

RESULTS

The dynamics of the leaflets are investigated, and blood flow characteristics such as velocities, vorticities as well as viscous and turbulent shear stress were precisely captured in the flow domain specifically in the hinge region. Streamline results are in accordance with the previously reported data, which show that the flared On-X valves inlet yields a more uniform flow in accelerating systole. Simulations show that aorta flexibility or valve annulus calcification causes variations up to 7% in maximum fluid velocity and 20% in Turbulence Kinetic Energy (TKE).

CONCLUSIONS

In this study, the complex flow field characteristics in the new generation of BMHVs considering aorta flexibility with healthy and calcified annulus were investigated. It was found that the blood flow around the hinges region is in the danger of hemolysis and platelet activation and subsequently thromboembolism. Furthermore, the results show that similar to vessel wall deformation, considering the probable annulus calcification after valve replacement is also essential.

Biomechanical Investigation of Disturbed Hemodynamics-Induced Tissue Degeneration in Abdominal Aortic Aneurysms Using Computational and Experimental Techniques

Abdominal aortic aneurysm (AAA) is the dilatation of the aorta beyond 50% of the normal vessel diameter. It is reported that 4-8% of men and 0.5-1% of women above 50 years of age bear an AAA and it accounts for ~15,000 deaths per year in the United States alone. If left untreated, AAA might gradually expand until rupture; the most catastrophic complication of the aneurysmal disease that is accompanied by a striking overall mortality of 80%. The precise mechanisms leading to AAA rupture remains unclear. Therefore, characterization of disturbed hemodynamics within AAAs will help to understand the mechanobiological development of the condition which will contribute to novel therapies for the condition. Due to geometrical complexities, it is challenging to directly quantify disturbed flows for AAAs clinically. Two other approaches for this investigation are computational modeling and experimental flow measurement. In computational modeling, the problem is first defined mathematically, and the solution is approximated with numerical techniques to get characteristics of flow. In experimental flow measurement, once the setup providing physiological flow pattern in a phantom geometry is constructed, velocity measurement system such as particle image velocimetry (PIV) enables characterization of the flow. We witness increasing number of applications of these complimentary approaches for AAA investigations in recent years. In this paper, we outline the details of computational modeling procedures and experimental settings and summarize important findings from recent studies, which will help researchers for AAA investigations and rupture mechanics.

Adaptation of a Mice Doppler Echocardiography Platform to Measure Cardiac Flow Velocities for Embryonic Chicken and Adult Zebrafish

Ultrasonography is the most widely used imaging technique in cardiovascular medicine. In this technique, a piezoelectric crystal produces, sends, and receives high frequency ultrasound waves to the body to create an image of internal organs. It enables practical real time visualization in a non-invasive manner, making the modality especially useful to image dynamic cardiac structures. In the last few decades, echocardiography has been applied to in vivo cardiac disease models, mainly to rodents. While clinical echocardiography platforms can be used for relatively large animals such as pigs and rats, specialized systems are needed for smaller species. Theoretically, as the size of the imaged sample decreases, the frequency of the ultrasound transducer needed to image the sample increases. There are multiple modes of echocardiography imaging. In Doppler mode, erythrocytes blood flow velocities are measured from the frequency shift of the sent ultrasound waves compared to received echoes. Recorded data are then used to calculate cardiac function parameters such as cardiac output, as well as the hemodynamic shear stress levels in the heart and blood vessels. The multi-mode (i.e., b-mode, m-mode, Pulsed Doppler, Tissue Doppler, etc.) small animal ultrasound systems in the market can be used for most in vivo cardiac disease models including mice, embryonic chick and zebrafish. These systems are also associated with significant costs. Alternatively, there are more economical single-mode echocardiography platforms. However, these are originally built for mice studies and they need to be tested and evaluated for smaller experimental models. We recently adapted a mice Doppler echocardiography system to measure cardiac flow velocities for adult zebrafish and embryonic chicken. We successfully assessed cardiac function and hemodynamic shear stress for normal as well as for diseased embryonic chicken and zebrafish. In this paper, we will present our detailed protocols for Doppler flow measurements and further cardiac function analysis on these models using the setup. The protocols will involve detailed steps for animal stabilization, probe orientation for specific measurements, data acquisition, and data analysis. We believe this information will help cardiac researchers to establish similar echocardiography platforms in their labs in a practical and economical manner.

Modeling mitral valve stenosis: A parametric study on the stenosis severity level

New computational techniques providing more accurate representation of human heart pathologies could help uncovering relevant physical phenomena and improve the outcome of medical therapies. In this framework, the present work describes an efficient computational model for the evaluation of the ventricular flow alteration in presence of mitral valve stenosis. The model is based on the direct numerical simulation of the Navier-Stokes equations two-way coupled with a structural solver for the left ventricle and mitral valve dynamics. The presence of mitral valve stenosis is mimicked by a single-parameter constraint acting on the kinematics of the mitral leaflets. Four different degrees of mitral valve stenosis are considered focusing on the hemodynamic alterations occurring in pathologic conditions. The mitral jet, generated during diastole, is seen to shrink and strengthen when the stenosis gets more severe. As a consequence, the kinetic energy of the flow, the tissues shear stresses, the transvalvular pressure drop and mitral regurgitation increase. It results that, as the stenosis severity level increases, the geometric and effective orifice areas decrease up to 50% with respect the normal case due to the reduced leaflets mobility and stronger blood acceleration during the diastolic phase. The modified intraventricular hemodynamics is also related to a stronger pressure gradient that, for severe stenosis, can be more than ten times larger than the healthy valve case. These computational results are fully consistent with the available clinical literature and open the way to the virtual assessment of surgical procedures and to the evaluation of prosthetic devices.