3D Fluid-Structure Interaction Simulation of Aortic Valves Using a Unified Continuum ALE FEM Model

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.

Nature of aortic annulus: Influence of annulus dynamic on the aortic valve hemodynamics

Accurate imaging reports of the aortic valve indicate that the diameter of the aortic annulus changes regularly during a cardiac cycle. Most of these studies aim to demonstrate the proper method for estimating the aortic annulus diameter before performing TAVR surgery, revealing that the aortic annulus is dynamic and not constant throughout the cardiac cycle. This raises the question of how fixing the aortic annulus might affect valve function, which is a question that still needs to be addressed. Therefore, the present study seeks to address this question and elucidate the dynamic impact of the aortic annulus on aortic valve hemodynamics. Two computational models based on this hypothesis were created and solved, and then their results were compared. Both models are identical, except for the intrinsic dynamic nature of the aortic annulus. One model consists of the dynamic behavior, and the other simulates a fixed annulus, resembling the effect of a TAVR operation, SAVR, or any phenomenon that eliminates the dynamic nature of the annulus. Our research findings indicate that the dynamic nature of the annulus enhances blood flow (+2.7 %), increases mean velocity (+11.9) and kinetic energy density (+34 %), prolongs momentum retention during systole, stabilizes the flow jet at the end of systole, reduces the required pressure to keep the leaflets open (-40.9 % at 0.3s), and sustains ventricular pressure superiority (+9.4 %) over the aorta for a longer duration (+17.7 % of systole), a crucial factor in preventing backflow during aortic valve closure. Based on these results, more attention should be paid to the dynamic nature of the annulus.

Image-Based Computational Fluid Dynamics to Compare Two Repair Techniques for Mitral Valve Prolapse

OBJECTIVE : The treatment of mitral valve prolapse involves two distinct repair techniques: chordal replacement (Neochordae technique) and leaflet resection (Resection technique). However, there is still a debate in the literature about which is the optimal one. In this context, we performed an image-based computational fluid dynamic study to evaluate blood dynamics in the two surgical techniques. METHODS : We considered a healthy subject (H) and two patients (N and R) who underwent surgery for prolapse of the posterior leaflet and were operated with the Neochordae and Resection technique, respectively. Computational Fluid Dynamics (CFD) was employed with prescribed motion of the entire left heart coming from cine-MRI images, with a Large Eddy Simulation model to describe the transition to turbulence and a resistive method for managing valve dynamics. We created three different virtual scenarios where the operated mitral valves were inserted in the same left heart geometry of the healthy subject to study the differences attributed only to the two techniques. RESULTS : We compared the three scenarios by quantitatively analyzing ventricular velocity patterns and pressures, transition to turbulence, and the ventricle ability to prevent thrombi formation. From these results, we found that the operative techniques affected the ventricular blood dynamics in different ways, with variations attributed to the reduced mobility of the Resection posterior leaflet. Specifically, the Resection technique resulted in turbulent forces, related with the risk of hemolysis formation, up to 640 Pa, while the other two scenarios exhibited a maximum of 240 Pa. Moreover, in correspondence of the ventricular apex, the Resection technique reduced the areas with low velocity to 15%, whereas the healthy case and the Neochordae case maintained these areas at 30 and 48%, respectively. Our findings suggest that the Neochordae technique developed a more physiological flow with respect to the Resection technique. CONCLUSION: Resection technique gives rise to a different direction of the mitral jet during diastole increasing the ability to washout the ventricular apex preventing from thrombi formation, but at the same time it promotes turbulence formation that is associated with ventricular effort and risk of hemolysis.

Fluid-structure interaction simulation of mechanical aortic valves: a narrative review exploring its role in total product life cycle

Over the last years computer modelling and simulation has emerged as an effective tool to support the total product life cycle of cardiovascular devices, particularly in the device preclinical evaluation and post-market assessment. Computational modelling is particularly relevant for heart valve prostheses, which require an extensive assessment of their hydrodynamic performance and of risks of hemolysis and thromboembolic complications associated with mechanically-induced blood damage. These biomechanical aspects are typically evaluated through a fluid-structure interaction (FSI) approach, which enables valve fluid dynamics evaluation accounting for leaflets movement. In this context, the present narrative review focuses on the computational modelling of bileaflet mechanical aortic valves through FSI approach, aiming to foster and guide the use of simulations in device total product life cycle. The state of the art of FSI simulation of heart valve prostheses is reviewed to highlight the variety of modelling strategies adopted in the literature. Furthermore, the integration of FSI simulations in the total product life cycle of bileaflet aortic valves is discussed, with particular emphasis on the role of simulations in complementing and potentially replacing the experimental tests suggested by international standards. Simulations credibility assessment is also discussed in the light of recently published guidelines, thus paving the way for a broader inclusion of in silico evidence in regulatory submissions. The present narrative review highlights that FSI simulations can be successfully framed within the total product life cycle of bileaflet mechanical aortic valves, emphasizing that credible in silico models evaluating the performance of implantable devices can (at least) partially replace preclinical in vitro experimentation and support post-market biomechanical evaluation, leading to a reduction in both time and cost required for device development.

A mechanically consistent unified formulation for fluid-porous-structure-contact interaction

Fluid-structure interaction with contact poses profound mathematical and numerical challenges, particularly when considering realistic contact scenarios and the influence of surface roughness. Computationally, contact introduces challenges in altering the fluid domain topology and preserving stress balance. This work introduces a new mathematical framework for a unified continuum description of fluid-porous-structure-contact interaction (FPSCI), leveraging the Navier-Stokes-Brinkman (NSB) equations to incorporate porous effects within the surface asperities in the contact region. Our approach maintains mechanical consistency during contact, circumventing issues associated with contact models and complex interface coupling conditions, allowing for the modeling of tangential creeping flows due to surface roughness. The unified continuum and variational multiscale formulation ensure robustness by enabling stable and unified integration of fluid, porous, and solid sub-problems. Computational efficiency and ease of implementation - key advantages of our approach - are demonstrated by solving two benchmark problems of a falling ball and an idealized heart valve. This research has broad implications for fields reliant on accurate fluid-structure interactions and promising advancements in modeling and numerical simulation techniques.

Blood flow assessment technology in aortic surgery: a narrative review

Background and Objective

Blood flow assessment is an emerging technique that allows for assessment of hemodynamics in the heart and blood vessels. Recent advances in cardiovascular imaging technologies have made it possible for this technique to be more accessible to clinicians and researchers. Blood flow assessment typically refers to two techniques: measurement-based flow visualization using echocardiography or four-dimensional flow magnetic resonance imaging (4D flow MRI), and computer-based flow simulation based on computational fluid dynamics modeling. Using these methods, blood flow patterns can be visualized and quantitative measurements of mechanical stress on the walls of the ventricles and blood vessels, most notably the aorta, can be made. Thus, blood flow assessment has been enhancing the understanding of cardiac and aortic diseases; however, its introduction to clinical practice has been negligible yet. In this article, we aim to discuss the clinical applications and future directions of blood flow assessment in aortic surgery. We then provide our unique perspective on the technique's translational impact on the surgical management of aortic disease.

Methods

Articles from the PubMed database and Google Scholar regarding blood flow assessment in aortic surgery were reviewed. For the initial search, articles published between 2013 and 2023 were prioritized, including original articles, clinical trials, case reports, and reviews. Following the initial search, additional articles were considered based on manual searches of the references from the retrieved literature.

Key Content and Findings

In aortic root pathology and ascending aortic aneurysms, blood flow assessment can elucidate postoperative hemodynamic changes after surgical reconfiguration of the aortic valve complex or ascending aorta. In cases of aortic dissection, analysis of blood flow can predict future aortic dilatation. For complicated congenital aortic anomalies, surgeons may use preoperative imaging to perform "virtual surgery", in which blood flow assessment can predict postoperative hemodynamics for different surgical reconstructions and assist in procedural planning even before entering the operating room.

Conclusions

Blood flow assessment and computational modeling can evaluate hemodynamics and flow patterns by visualizing blood flow and calculating biomechanical forces in patients with aortic disease. We anticipate that blood flow assessment will become an essential tool in the treatment planning and understanding of the progression of aortic disease.

Design, manufacturing and testing of a green non-isocyanate polyurethane prosthetic heart valve

The sole effective treatment for most patients with heart valve disease is valve replacement by implantation of mechanical or biological prostheses. However, mechanical valves represent high risk of thromboembolism, and biological prostheses are prone to early degeneration. In this work, we aim to determine the potential of novel environmentally-friendly non-isocyanate polyurethanes (NIPUs) for manufacturing synthetic prosthetic heart valves. Polyhydroxyurethane (PHU) NIPUs are synthesized via an isocyanate-free route, tested in vitro, and used to produce aortic valves. PHU elastomers reinforced with a polyester mesh show mechanical properties similar to native valve leaflets. These NIPUs do not cause hemolysis. Interestingly, both platelet adhesion and contact activation-induced coagulation are strongly reduced on NIPU surfaces, indicating low thrombogenicity. Fibroblasts and endothelial cells maintain normal growth and shape after indirect contact with NIPUs. Fluid-structure interaction (FSI) allows modeling of the ideal valve design, with minimal shear stress on the leaflets. Injection-molded valves are tested in a pulse duplicator and show ISO-compliant hydrodynamic performance, comparable to clinically-used bioprostheses. Poly(tetrahydrofuran) (PTHF)-NIPU patches do not show any evidence of calcification over a period of 8 weeks. NIPUs are promising sustainable biomaterials for the manufacturing of improved prosthetic valves with low thrombogenicity.

Smoothed particle hydrodynamics based FSI simulation of the native and mechanical heart valves in a patient-specific aortic model

The failure of the aortic heart valve is common, resulting in deterioration of the pumping function of the heart. For the end stage valve failure, bi-leaflet mechanical valve (most popular artificial valve) is implanted. However, due to its non-physiological behaviour, a significant alteration is observed in the normal haemodynamics of the aorta. While in-vivo experimentation of a human heart valve (native and artificial) is a formidable task, in-silico study using computational fluid dynamics (CFD) with fluid structure interaction (FSI) is an effective and economic tool for investigating the haemodynamics of natural and artificial heart valves. In the present work, a haemodynamic model of a natural and mechanical heart valve has been developed using meshless particle-based smoothed particle hydrodynamics (SPH). In order to further enhance its clinical relevance, this study employs a patient-specific vascular geometry and presents a successful validation against traditional finite volume method and 4D magnetic resonance imaging (MRI) data. The results have demonstrated that SPH is ideally suited to simulate the heart valve function due to its Lagrangian description of motion, which is a favourable feature for FSI. In addition, a novel methodology for the estimation of the wall shear stress (WSS) and other related haemodynamic parameters have been proposed from the SPH perspective. Finally, a detailed comparison of the haemodynamic parameters has been carried out for both native and mechanical aortic valve, with a particular emphasis on the clinical risks associated with the mechanical valve.

Turbulent blood dynamics in the left heart in the presence of mitral regurgitation: a computational study based on multi-series cine-MRI

In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution.

Modeling isovolumetric phases in cardiac flows by an Augmented Resistive Immersed Implicit Surface method

A major challenge in the computational fluid dynamics modeling of the heart function is the simulation of isovolumetric phases when the hemodynamics problem is driven by a prescribed boundary displacement. During such phases, both atrioventricular and semilunar valves are closed: consequently, the ventricular pressure may not be uniquely defined, and spurious oscillations may arise in numerical simulations. These oscillations can strongly affect valve dynamics models driven by the blood flow, making unlikely to recovering physiological dynamics. Hence, prescribed opening and closing times are usually employed, or the isovolumetric phases are neglected altogether. In this article, we propose a suitable modification of the Resistive Immersed Implicit Surface (RIIS) method (Fedele et al., Biomech Model Mechanobiol 2017, 16, 1779-1803) by introducing a reaction term to correctly capture the pressure transients during isovolumetric phases. The method, that we call Augmented RIIS (ARIIS) method, extends the previously proposed ARIS method (This et al., Int J Numer Methods Biomed Eng 2020, 36, e3223) to the case of a mesh which is not body-fitted to the valves. We test the proposed method on two different benchmark problems, including a new simplified problem that retains all the characteristics of a heart cycle. We apply the ARIIS method to a fluid dynamics simulation of a realistic left heart geometry, and we show that ARIIS allows to correctly simulate isovolumetric phases, differently from standard RIIS method. Finally, we demonstrate that by the new method the cardiac valves can open and close without prescribing any opening/closing times.

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.

A mathematical model that integrates cardiac electrophysiology, mechanics, and fluid dynamics: Application to the human left heart

We propose a mathematical and numerical model for the simulation of the heart function that couples cardiac electrophysiology, active and passive mechanics and hemodynamics, and includes reduced models for cardiac valves and the circulatory system. Our model accounts for the major feedback effects among the different processes that characterize the heart function, including electro-mechanical and mechano-electrical feedback as well as force-strain and force-velocity relationships. Moreover, it provides a three-dimensional representation of both the cardiac muscle and the hemodynamics, coupled in a fluid-structure interaction (FSI) model. By leveraging the multiphysics nature of the problem, we discretize it in time with a segregated electrophysiology-force generation-FSI approach, allowing for efficiency and flexibility in the numerical solution. We employ a monolithic approach for the numerical discretization of the FSI problem. We use finite elements for the spatial discretization of partial differential equations. We carry out a numerical simulation on a realistic human left heart model, obtaining results that are qualitatively and quantitatively in agreement with physiological ranges and medical images.

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.

Assessment of valve implantation in the descending aorta as an alternative for aortic regurgitation patients not treatable with conventional procedures

BACKGROUND

Aortic Regurgitation (AR) produces the entrance of an abnormal amount of blood in the left ventricle. This disease is responsible for high morbidity and mortality worldwide and may be caused by an aortic valve dysfunction. Surgical and transcatheter aortic valve replacement (TAVR) are the current options for treating AR. They have replaced older procedures such as Hufnagel's one. However, some physicians have reconsidered this procedure as a less aggressive alternative for patients not eligible for surgical or TAVR. Although Hufnagel suggested a 75% regurgitation reduction when a valve is placed in the descending aorta, a quantification of this value has not been reported.

METHODS

In this paper, CFD/FSI numerical simulation is conducted on an idealized geometry. We quantify the effect of placing a bileaflet mechanical heart valve in the descending aorta on a moderate-severe AR case. A three-element Windkessel model is employed to prescribe pressure outlet boundary conditions. We calculate the resulting flow rates and pressures at the aorta and first-generation vessels. Moreover, we evaluate several indices to assess the improvement due to the valve introduction.

RESULTS AND CONCLUSIONS

Regurgitation fraction (RF) is reduced from 37.5% (without valve) to 18.0% (with valve) in a single cardiac cycle. This reduction clearly shows the remarkable efficacy of the rescued technique. It will further ameliorate the left ventricle function in the long-term. Moreover, the calculations show that the implantation in that location introduces fewer incompatibilities' risks than a conventional one. The proposed methodology can be extended to any particular conditions (pressure waveforms/geometry) and is designed to assess usual clinical parameters employed by physicians.

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

Computational Methods for Fluid-Structure Interaction Simulation of Heart Valves in Patient-Specific Left Heart Anatomies

Given the complexity of human left heart anatomy and valvular structures, the fluid–structure interaction (FSI) simulation of native and prosthetic valves poses a significant challenge for numerical methods. In this review, recent numerical advancements for both fluid and structural solvers for heart valves in patient-specific left hearts are systematically considered, emphasizing the numerical treatments of blood flow and valve surfaces, which are the most critical aspects for accurate simulations. Numerical methods for hemodynamics are considered under both the continuum and discrete (particle) approaches. The numerical treatments for the structural dynamics of aortic/mitral valves and FSI coupling methods between the solid Ωs and fluid domain Ωf are also reviewed. Future work toward more advanced patient-specific simulations is also discussed, including the fusion of high-fidelity simulation within vivo measurements and physics-based digital twining based on data analytics and machine learning techniques.

Predictive Model for Thrombus Formation After Transcatheter Valve Replacement

PURPOSE

Leaflet thrombosis is a significant adverse event after transcatheter aortic valve (TAV) replacement (TAVR). The purpose of our study was to present a semi-empirical, mathematical model that links patient-specific anatomic, valve, and flow parameters to predict likelihood of leaflet thrombosis.

METHODS

The two main energy sources of neo-sinus (NS) washout after TAVR include the jet flow downstream of the TAV and NS geometric change in volume due to the leaflets opening and closing. Both are highly dependent on patient anatomic and hemodynamic factors. As rotation of blood flow is prevalent in both the sinus of Valsalva and then the NS, we adopted the vorticity flux or circulation (Г) as a metric quantifying overall washout. Leaflet thrombus volumes were segmented based on hypo-attenuating leaflet thickening (HALT) in post-TAVR patient's gated computed tomography. Г was assessed using dimensional scaling as well as computational fluid dynamics (CFD) respectively and correlated to the thrombosis volumes using sensitivity and specificity analysis.

RESULTS

Г in the NS, that accounted for patient flow and anatomic conditions derived from scaling arguments significantly better predicted the occurrence of leaflet thrombus than CFD derived measures such as stasis volumes or wall shear stress. Given results from the six patient datasets considered herein, a threshold Г value of 28.0 yielded a sensitivity and specificity of 100% where patients with Gamma < 28 developed valve thrombosis. A 10% error in measurements of all variables can bring the sensitivity specificity down to 87%.

CONCLUSION

A predictive model relating likelihood of valve thrombosis using Г in the NS was developed with promising sensitivity and specificity. With further studies and improvements, this predictive technology may lead to alerting physicians on the risk for thrombus formation following TAVR.

A design-based model of the aortic valve for fluid-structure interaction

This paper presents a new method for modeling the mechanics of the aortic valve and simulates its interaction with blood. As much as possible, the model construction is based on first principles, but such that the model is consistent with experimental observations. We require that tension in the leaflets must support a pressure, then derive a system of partial differential equations governing its mechanical equilibrium. The solution to these differential equations is referred to as the predicted loaded configuration; it includes the loaded leaflet geometry, fiber orientations and tensions needed to support the prescribed load. From this configuration, we derive a reference configuration and constitutive law. In fluid-structure interaction simulations with the immersed boundary method, the model seals reliably under physiological pressures and opens freely over multiple cardiac cycles. Further, model closure is robust to extreme hypo- and hypertensive pressures. Then, exploiting the unique features of this model construction, we conduct experiments on reference configurations, constitutive laws and gross morphology. These experiments suggest the following conclusions: (1) The loaded geometry, tensions and tangent moduli primarily determine model function. (2) Alterations to the reference configuration have little effect if the predicted loaded configuration is identical. (3) The leaflets must have sufficiently nonlinear material response to function over a variety of pressures. (4) Valve performance is highly sensitive to free edge length and leaflet height. These conclusions suggest appropriate gross morphology and material properties for the design of prosthetic aortic valves. In future studies, our aortic valve modeling framework can be used with patient-specific models of vascular or cardiac flow.

Local and global growth and remodeling in calcific aortic valve disease and aging

Aging and calcific aortic valve disease (CAVD) are the main factors leading to aortic stenosis. Both processes are accompanied by growth and remodeling pathways that play a crucial role in aortic valve pathophysiology. Herein, a computational growth and remodeling (G&R) framework was developed to investigate the effects of aging and calcification on aortic valve dynamics. Particularly, an algorithm was developed to couple the global growth and stiffening of the aortic valve due to aging and the local growth and stiffening due to calcification with the aortic valve transient dynamics. The aortic valve dynamics during baseline were validated with available data in the literature. Subsequently, the changes in aortic valve dynamic patterns during aging and CAVD progression were studied. The results revealed the patterns in geometric orifice area reduction and an increase in the valve stress during local and global growth and remodeling of the aortic valve. The proposed algorithm provides a framework to couple mechanobiology models of disease growth with tissue-scale transient structural mechanics models to study the biomechanical changes during cardiovascular disease growth and aging.

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.

Biomechanical Assessment of Bicuspid Aortic Valve Phenotypes: A Fluid-Structure Interaction Modelling Approach

PURPOSE

Bicuspid aortic valve (BAV) is a congenital heart malformation with phenotypic heterogeneity. There is no prior computational study that assesses the haemodynamic and valve mechanics associated with BAV type 2 against a healthy tricuspid aortic valve (TAV) and other BAV categories.

METHODS

A proof-of-concept study incorporating three-dimensional fluid-structure interaction (FSI) models with idealised geometries (one TAV and six BAVs, namely type 0 with lateral and anterior-posterior orientations, type 1 with R-L, N-R and N-L leaflet fusion and type 2) has been developed. Transient physiological boundary conditions have been applied and simulations were run using an Arbitrary Lagrangian-Eulerian formulation.

RESULTS

Our results showed the presence of abnormal haemodynamics in the aorta and abnormal valve mechanics: type 0 BAVs yielded the best haemodynamical and mechanical outcomes, but cusp stress distribution varied with valve orifice orientation, which can be linked to different cusp calcification location onset; type 1 BAVs gave rise to similar haemodynamics and valve mechanics, regardless of raphe position, but this position altered the location of abnormal haemodynamic features; finally, type 2 BAV constricted the majority of blood flow, exhibiting the most damaging haemodynamic and mechanical repercussions when compared to other BAV phenotypes.

CONCLUSION

The findings of this proof-of-concept work suggest that there are specific differences across haemodynamics and valve mechanics associated with BAV phenotypes, which may be critical to subsequent processes associated with their pathophysiology processes.

Fluid-Structure Interaction Models of Bioprosthetic Heart Valve Dynamics in an Experimental Pulse Duplicator

Computer modeling and simulation is a powerful tool for assessing the performance of medical devices such as bioprosthetic heart valves (BHVs) that promises to accelerate device design and regulation. This study describes work to develop dynamic computer models of BHVs in the aortic test section of an experimental pulse-duplicator platform that is used in academia, industry, and regulatory agencies to assess BHV performance. These computational models are based on a hyperelastic finite element extension of the immersed boundary method for fluid-structure interaction (FSI). We focus on porcine tissue and bovine pericardial BHVs, which are commonly used in surgical valve replacement. We compare our numerical simulations to experimental data from two similar pulse duplicators, including a commercial ViVitro system and a custom platform related to the ViVitro pulse duplicator. Excellent agreement is demonstrated between the computational and experimental results for bulk flow rates, pressures, valve open areas, and the timing of valve opening and closure in conditions commonly used to assess BHV performance. In addition, reasonable agreement is demonstrated for quantitative measures of leaflet kinematics under these same conditions. This work represents a step towards the experimental validation of this FSI modeling platform for evaluating BHVs.

New Insights into Valve Hemodynamics

Heart valve diseases are common disorders with five million annual diagnoses being made in the United States alone. All heart valve disorders alter cardiac hemodynamic performance; therefore, treatments aim to restore normal flow. This paper reviews the state-of-the-art clinical and engineering advancements in heart valve treatments with a focus on hemodynamics. We review engineering studies and clinical literature on the experience with devices for aortic valve treatment, as well as the latest advancements in mitral valve treatments and the pulmonic and tricuspid valves on the right side of the heart. Upcoming innovations will potentially revolutionize treatment of heart valve disorders. These advancements, and more gradual enhancements in the procedural techniques and imaging modalities, could improve the quality of life of patients suffering from valvular disease who currently cannot be treated.

A partition of unity approach to fluid mechanics and fluid-structure interaction

For problems involving large deformations of thin structures, simulating fluid-structure interaction (FSI) remains a computationally expensive endeavour which continues to drive interest in the development of novel approaches. Overlapping domain techniques have been introduced as a way to combine the fluid-solid mesh conformity, seen in moving-mesh methods, without the need for mesh smoothing or re-meshing, which is a core characteristic of fixed mesh approaches. In this work, we introduce a novel overlapping domain method based on a partition of unity approach. Unified function spaces are defined as a weighted sum of fields given on two overlapping meshes. The method is shown to achieve optimal convergence rates and to be stable for steady-state Stokes, Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for FSI in the case of 2D flow past an elastic beam simulation. These initial results point to the potential applicability of the method to a wide range of FSI applications, enabling boundary layer refinement and large deformations without the need for re-meshing or user-defined stabilization.

Augmented resistive immersed surfaces valve model for the simulation of cardiac hemodynamics with isovolumetric phases

In order to reduce the complexity of heart hemodynamics simulations, uncoupling approaches are often considered for the modeling of the immersed valves as an alternative to complex fluid-structure interaction (FSI) models. A possible shortcoming of these simplified approaches is the difficulty to correctly capture the pressure dynamics during the isovolumetric phases. In this work, we propose an enhanced resistive immersed surfaces (RIS) model of cardiac valves, which overcomes this issue. The benefits of the model are investigated and tested in blood flow simulations of the left heart where the physiological behavior of the intracavity pressure during the isovolumetric phases is recovered without using fully coupled fluid-structure models and without important alteration of the associated velocity field.

Simulation of semilunar valve function: computer-aided design, 3D printing and flow assessment with MR

BACKGROUND

The structure of the valve leaflets and sinuses are crucial in supporting the proper function of the semilunar valve and ensuring leaflet durability. Therefore, an enhanced understanding of the structural characteristics of the semilunar valves is fundamental to the evaluation and staging of semilunar valve pathology, as well as the development of prosthetic or bioprosthetic valves. This paper illustrates the process of combining computer-aided design (CAD), 3D printing and flow assessment with 4-dimensional flow magnetic resonance imaging (MRI) to provide detailed assessment of the structural and hemodynamic characteristics of the normal semilunar valve.

METHODS

Previously published geometric data on the aortic valve was used to model the 'normal' tricuspid aortic valve with a CAD software package and 3D printed. An MRI compatible flow pump with the capacity to mimic physiological flows was connected to the phantom. A peak flow rate of 100 mL/s and heart rate of 60 beats per minute were used. MRI measurements included cine imaging, 2D and 4D phase-contrast imaging to assess valve motion, flow velocity and complex flow patterns.

RESULTS

Cine MRI data showed normal valve function and competency throughout the cardiac cycle in the 3D-printed phantom. Quantitative analysis of 4D Flow data showed net flow through 2D planes proximal and distal to the valve were very consistent (26.03 mL/s and 26.09 mL/s, respectively). Measurements of net flow value agreed closely with the flow waveform provided to the pump (27.74 mL/s), confirming 4D flow acquisition in relation to the pump output. Peak flow values proximal and distal to the valve were 78.4 mL/s and 63.3 mL/s, respectively. Particle traces of flow from 4D-phase contrast MRI data demonstrated flow through the valve into the ascending aorta and vortices within the aortic sinuses, which are expected during ventricular diastole.

CONCLUSION

In this proof of concept study, we have demonstrated the ability to generate physiological 3D-printed aortic valve phantoms and evaluate their function with cine- and 4D Flow MRI. This technology can work synergistically with promising tissue engineering research to develop optimal aortic valve replacements, which closely reproduces the complex function of the normal aortic valve.

Computational investigation of hemodynamics in hardshell venous reservoirs: A comparative study

Extracorporeal circulation using heart-lung-machines is associated with a profound activation of corpuscular and plasmatic components of circulating blood, which can also lead to deleterious events such as systemic inflammatory response and hemolysis. Individual components used to install the extracorporeal circulation have an impact on the level of activation, most predominantly membrane oxygenators and hardshell venous reservoirs as used in extracorporeal systems. The blood flows in two different hardshell reservoirs are computationally investigated. A special emphasis is placed on the prediction of an onset of transition and turbulence generation. Reynolds-averaged numerical simulations (RANS) based on a transitional turbulence model, as well as large eddy simulations (LES) are applied to achieve an accurate prediction. In the LES analysis, the non-Newtonian behavior of the blood is considered via the Carreau model. Blood damage potential is quantified applying the Modified Index of Hemolysis (MIH) based on the predicted flow fields. The results indicate that the flows in both reservoirs remain predominantly laminar. For one of the reservoirs, considerable turbulence generation is observed near the exit site, caused by the specific design for the connection with the drainage tube. This difference causes the MIH of this reservoir to be nearly twice as large as compared to the alternative design. However, a substantial improvement of these performance criteria can be expected by a local geometry modification.