Fluid-structure interaction model of aortic valve with porcine-specific collagen fiber alignment in the cusps

On the role of aortic valve architecture for physiological hemodynamics and valve replacement, Part I: Flow configuration and vortex dynamics

Aortic valve replacement has become an increasing concern due to the rising prevalence of aortic stenosis in an ageing population. Existing replacement options have limitations, necessitating the development of improved prosthetic aortic valves. In this study, flow characteristics during systole in a stenotic aortic valve case are compared with those downstream of two newly designed surgical bioprosthetic aortic valves (BioAVs). To do so, advanced three-dimensional fluid-structure interaction simulations are conducted and dedicated analysis methods to investigate jet flow configuration and vortex dynamics are developed. Our findings reveal that the stenotic case maintains a high jet flow eccentricity due to a fixed orifice geometry, resulting in flow separation and increased vortex stretching and tilting in the commissural low-flow regions. One BioAV design introduces non-axisymmetric leaflet motion, which reduces the maximum jet velocity and forms more vortical structures. The other BioAV design produces a fixed symmetric triangular jet shape due to non-moving leaflets and exhibits favourable vorticity attenuation, revealed by negative temporally and spatially averaged projected vortex stretching values, and significantly reduced drag. Therefore, this study highlights the benefits of custom-designed aortic valves in the context of their replacement through comprehensive and novel flow analyses. The results emphasise the importance of analysing jet flow, vortical structures, momentum balance and vorticity transport for thoroughly evaluating aortic valve performance.

Systematic adjustment of root dimensions to cusp size in aortic valve repair: a computer simulation

OBJECTIVES

Aortic valve repair requires the creation of a normal geometry of cusps and aortic root. Of the different dimensions, geometric cusp height is the most difficult to change while annular and sinotubular dimensions can be easily modified. The objective of this study was to investigate, by computer simulation, ideal combinations of annular and sinotubular junction size for a given geometric height.

METHODS

Based on a literature review of anatomical data, a computational biomechanics model was generated for a tricuspid aortic valve. We aimed to determine the ideal relationships for the root dimensions, keeping geometric height constant and creating different combinations of the annular and sinotubular junction dimensions. Using this model, 125 virtual anatomies were created, with 25 different combinations of annulus and sinotubular junction. Effective height, coaptation height and mechanical cusp stress were calculated with the valves in closed configuration.

RESULTS

Generally, within the analysed range of geometric heights, changes to the annular diameter yielded a stronger impact than sinotubular junction diameter changes for optimal valve configuration. The best results were obtained with the sinotubular junction being 2-4 mm larger than the annulus, leading to higher effective height, normal coaptation height and lower stress. Within the range tested, stenosis did not occur due to annular reduction.

CONCLUSIONS

In tricuspid aortic valves, the geometric height can be used to predict ideal post-repair annular and sinotubular junction dimensions for optimal valve configuration. Such an ideal configuration is associated with reduced cusp stress.

Fluid-structure interaction modeling of compliant aortic valves using the lattice Boltzmann CFD and FEM methods

The lattice Boltzmann method (LBM) has been increasingly used as a stand-alone CFD solver in various biomechanical applications. This study proposes a new fluid-structure interaction (FSI) co-modeling framework for the hemodynamic-structural analysis of compliant aortic valves. Toward that goal, two commercial software packages are integrated using the lattice Boltzmann (LBM) and finite element (FE) methods. The suitability of the LBM-FE hemodynamic FSI is examined in modeling healthy tricuspid and bicuspid aortic valves (TAV and BAV), respectively. In addition, a multi-scale structural approach that has been employed explicitly recognizes the heterogeneous leaflet tissues and differentiates between the collagen fiber network (CFN) embedded within the elastin matrix of the leaflets. The CFN multi-scale tissue model is inspired by monitoring the distribution of the collagen in 15 porcine leaflets. Different simulations have been examined, and structural stresses and resulting hemodynamics are analyzed. We found that LBM-FE FSI approach can produce good predictions for the flow and structural behaviors of TAV and BAV and correlates well with those reported in the literature. The multi-scale heterogeneous CFN tissue structural model enhances our understanding of the mechanical roles of the CFN and the elastin matrix behaviors. The importance of LBM-FE FSI also emerges in its ability to resolve local hemodynamic and structural behaviors. In particular, the diastolic fluctuating velocity phenomenon near the leaflets is explicitly predicted, providing vital information on the flow transient nature. The full closure of the contacting leaflets in BAV is also demonstrated. Accordingly, good structural kinematics and deformations are captured for the entire cardiac cycle.

Aortic Leaflet Stresses Are Substantially Lower Using Pulmonary Visceral Pleura Than Pericardial Tissue

Background: Porcine heart and bovine pericardium valves, which are collagen-based with relatively little elastin, have been broadly utilized to construct bioprosthetic heart valves (BHVs). With a larger proportion of elastin, the pulmonary visceral pleura (PVP) has greater elasticity and could potentially serve as an advantageous biomaterial for the construction/repair of BHVs. The question of how the aortic valve’s performance is affected by its bending rigidity has not been well studied.

Methods: Based on the stress–strain relationships of the pericardium and PVP determined by planar uni-axial tests, a three-dimensional (3D) computational fluid–structure interaction (FSI) framework is employed to numerically investigate the aortic valve’s performance by considering three different cases with Young’s modulus as follows: E=375, 750, and 1500 kPa, respectively.

Results: The stroke volumes are 112, 99.6, and 91.4 ml as Young’s modulus increases from 375 to 750 and 1500 kPa, respectively. Peak geometric opening area (GOA) values are 2.3, 2.2, and 2.0 cm² for E=375, 750, and 1500 kPa, respectively. The maximum value of the aortic leaflet stress is about 271 kPa for E=375 kPa, and it increases to about 383 and 540 kPa for E=750 and 1500 kPa in the belly region at the peak systole, while it reduces from 550 kPa to 450 and 400 kPa for E=375, 750, and 1500 kPa, respectively, at the instant of peak “water-hammer”.

Conclusion: A more compliant PVP aortic leaflet valve with a smaller Young’s modulus, E, has a higher cardiac output, larger GOA, and lower hemodynamic resistance. Most importantly, the aortic leaflet stresses are substantially lower in the belly region within the higher compliance PVP aortic valve tissue during the systole phase, even though some stress increase is also found during the fast-closing phase due to the “water-hammer” effect similar to that in the pericardial tissue. Future clinical studies will be conducted to test the hypothesis that the PVP-based valve leaflets with higher compliance will have lower fatigue or calcification rates due to the overall lower stress.

The effect of the fibrocalcific pathological process on aortic valve stenosis in female patients: a finite element study

Calcific aortic valve disease (CAVD) is the most common heart valvular disease in the developed world. Most of the relevant research has been sex-blind, ignoring sex-related biological variables and thus under-appreciate sex differences. However, females present pronounced fibrosis for the same aortic stenosis (AS) severity compared with males, who exhibit more calcification. Herein, we present a computational model of fibrocalcific AV, aiming to investigate its effect on AS development. A parametric study was conducted to explore the influence of the total collagen fiber volume and its architecture on the aortic valve area (AVA). Towards that goal, computational models were generated for three females with stenotic AVs and different volumes of calcium. We have tested the influence of fibrosis on various parameters as fiber architecture, fibrosis location, and transvalvular pressure. We found that increased fiber volume with a low calcium volume could actively contribute to AS and reduce the AVA similarly to high calcium volume. Thus, the computed AVAs for our fibrocalcific models were 0.94 and 0.84 cm²and the clinical (Echo) AVAs were 0.82 and 0.8 cm². For the heavily calcified model, the computed AVA was 0.8 cm²and the clinical AVA was 0.73 cm². The proposed models demonstrated how collagen thickening influence the fibrocalcific-AS process in female patients. These models can assist in the clinical decision-making process and treatment development in valve therapy for female patients.

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.

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.

Structural Responses of Integrated Parametric Aortic Valve in an Electro-Mechanical Full Heart Model

The aortic valve (AV) is located between the left ventricle and the aorta and responsible for maintaining an outward unidirectional flow. Many AV hemodynamic and structural aspects of have been extensively studied, however, more sophisticated models are needed to better understand the AV biomechanical behavior. This study deals with integrating a new parametric AV structural model with the electro-mechanical Living Heart Human Model® (LHHM). The LHHM is a finite element model simulating human heart capable of realistic electro-mechanical simulations. Different geometric metrics of AV have been examined. New integrated structural AV model within the LHHM better predict local stresses during the cardiac cycle due to the realistic boundary condition derived from the LHHM. It was found that ellipticity index (EI), calculated as the ratio between the maximal (Max) and minimal (Min) aortic annulus (AA) diameters, well correlates with measured clinical data obtained from patients undergoing computed tomography (CT) while the annular perimeter (Perim) matches the same trend. This increases the confidence in the predicted kinematic behavior, leaflets coaptation, and the overall stresses. From the clinical aspect, the new proposed coupled and integrated AV modeling can serve as a platform for design and implementation of pre-transcatheter aortic valve replacement (TAVR) procedures.

Fluid-Structure Interaction Models of Bicuspid Aortic Valves: The Effects of Nonfused Cusp Angles

Bicuspid aortic valve (BAV) is the most common type of congenital heart disease, occurring in 0.5-2% of the population, where the valve has only two rather than the three normal cusps. Valvular pathologies, such as aortic regurgitation and aortic stenosis, are associated with BAVs, thereby increasing the need for a better understanding of BAV kinematics and geometrical characteristics. The aim of this study is to investigate the influence of the nonfused cusp (NFC) angle in BAV type-1 configuration on the valve's structural and hemodynamic performance. Toward that goal, a parametric fluid-structure interaction (FSI) modeling approach of BAVs is presented. Four FSI models were generated with varying NFC angles between 120 deg and 180 deg. The FSI simulations were based on fully coupled structural and fluid dynamic solvers and corresponded to physiologic values, including the anisotropic hyper-elastic behavior of the tissue. The simulated angles led to different mechanical behavior, such as eccentric jet flow direction with a wider opening shape that was found for the smaller NFC angles, while a narrower opening orifice followed by increased jet flow velocity was observed for the larger NFC angles. Smaller NFC angles led to higher concentrated flow shear stress (FSS) on the NFC during peak systole, while higher maximal principal stresses were found in the raphe region during diastole. The proposed biomechanical models could explain the early failure of BAVs with decreased NFC angles, and suggests that a larger NFC angle is preferable in suture annuloplasty BAV repair surgery.

The Genetic Regulation of Aortic Valve Development and Calcific Disease

Heart valves are dynamic, highly organized structures required for unidirectional blood flow through the heart. Over an average lifetime, the valve leaflets or cusps open and close over a billion times, however in over 5 million Americans, leaflet function fails due to biomechanical insufficiency in response to wear-and-tear or pathological stimulus. Calcific aortic valve disease (CAVD) is the most common valve pathology and leads to stiffening of the cusp and narrowing of the aortic orifice leading to stenosis and insufficiency. At the cellular level, CAVD is characterized by valve endothelial cell dysfunction and osteoblast-like differentiation of valve interstitial cells. These processes are associated with dysregulation of several molecular pathways important for valve development including Notch, Sox9, Tgfβ, Bmp, Wnt, as well as additional epigenetic regulators. In this review, we discuss the multifactorial mechanisms that contribute to CAVD pathogenesis and the potential of targeting these for the development of novel, alternative therapeutics beyond surgical intervention.

A New Growth Model for Aortic Valve Calcification

Calcific aortic valve disease (CAVD) is a progressive disease in which minerals accumulate in the tissue of the aortic valve cusps, stiffening them and preventing valve opening and closing. The process of valve calcification was found to be similar to that of bone formation including cell differentiation to osteoblast-like cells. Studies have shown the contribution of high strains to calcification initiation and growth process acceleration. In this paper, a new strain-based calcification growth model is proposed. The model aims to explain the unique shape of the calcification and other disease characteristics. The calcification process was divided into two stages: Calcification initiation and calcification growth. The initiation locations were based on previously published findings and a reverse calcification technique (RCT), which uses computed tomography (CT) scans of patients to reveal the calcification initiation point. The calcification growth process was simulated by a finite element model of one aortic valve cusp loaded with cyclic loading. Similar to Wolff's law, describing bone response to stress, our model uses strains to drive calcification formation. The simulation grows calcification from its initiation point to its full typical stenotic shape. Study results showed that the model was able to reproduce the typical calcification growth pattern and shape, suggesting that strain is the main driving force behind calcification progression. The simulation also sheds light on other disease characteristics, such as calcification growth acceleration as the disease progresses, as well as sensitivity to hypertension.

Image-based immersed boundary model of the aortic root

Each year, approximately 300,000 heart valve repair or replacement procedures are performed worldwide, including approximately 70,000 aortic valve replacement surgeries in the United States alone. Computational platforms for simulating cardiovascular devices such as prosthetic heart valves promise to improve device design and assist in treatment planning, including patient-specific device selection. This paper describes progress in constructing anatomically and physiologically realistic immersed boundary (IB) models of the dynamics of the aortic root and ascending aorta. This work builds on earlier IB models of fluid-structure interaction (FSI) in the aortic root, which previously achieved realistic hemodynamics over multiple cardiac cycles, but which also were limited to simplified aortic geometries and idealized descriptions of the biomechanics of the aortic valve cusps. By contrast, the model described herein uses an anatomical geometry reconstructed from patient-specific computed tomography angiography (CTA) data, and employs a description of the elasticity of the aortic valve leaflets based on a fiber-reinforced constitutive model fit to experimental tensile test data. The resulting model generates physiological pressures in both systole and diastole, and yields realistic cardiac output and stroke volume at physiological Reynolds numbers. Contact between the valve leaflets during diastole is handled automatically by the IB method, yielding a fully competent valve model that supports a physiological diastolic pressure load without regurgitation. Numerical tests show that the model is able to resolve the leaflet biomechanics in diastole and early systole at practical grid spacings. The model is also used to examine differences in the mechanics and fluid dynamics yielded by fresh valve leaflets and glutaraldehyde-fixed leaflets similar to those used in bioprosthetic heart valves. Although there are large differences in the leaflet deformations during diastole, the differences in the open configurations of the valve models are relatively small, and nearly identical hemodynamics are obtained in all cases considered.

A study of extracellular matrix remodeling in aortic heart valves using a novel biaxial stretch bioreactor

In aortic valves, biaxial cyclic stretch is known to modulate cell differentiation, extracellular matrix (ECM) synthesis and organization. We designed a novel bioreactor that can apply independent and precise stretch along radial and circumferential directions in a tissue culture environment. While this bioreactor can be used for either native or engineered tissues, this study determined matrix remodeling and strain distribution of aortic cusps after culturing under biaxial stretch for 14 days. The contents of collagen and glycosaminoglycans were determined using standard biochemical assays and compared with fresh controls. Strain fields in static cusps were more uniform than those in stretched cusps, which indicated degradation of the ECM fibers. The glycosaminoglycan content was significantly elevated in the static control as compared to fresh or stretched cusps, but no difference was observed in collagen content among the groups. The strain profile of freshly isolated fibrosa vs. ventricularis and left, right, and noncoronary cusps were also determined by Digital Image Correlation technique. Distinct strain patterns were observed under stretch on fibrosa and ventricularis sides and among the three cusps. This work highlights the critical role of the anisotropic ECM structure for proper functions of native aortic valves and the beneficial effects of biaxial stretch for maintenance of the native ECM structure.

Inverted orientation improves decellularization of whole porcine hearts

In structurally heterogeneous organs, such as heart, it is challenging to retain extracellular matrix integrity in the thinnest regions (eg, valves) during perfusion decellularization and completely remove cellular debris from thicker areas. The high inflow rates necessary to maintain physiologic pressure can distend or damage thin tissues, but lower pressures prolong the process and increase the likelihood of contamination. We examined two novel retrograde decellularization methods for porcine hearts: inverting the heart or venting the apex to decrease inflow rate. We measured flow dynamics through the aorta (Ao) and pulmonary artery (PA) at different Ao pressures and assessed the heart's appearance, turbidity of the outflow solutions, and coronary perfusion efficiency. We used rectangle image fitting of decellularized heart images to obtain a heart shape index. Using nonlinear optical microscopy, we determined the microstructure of collagen and elastin fibers of the aortic valve cusps. DNA, glycosaminoglycan, and residual detergent levels were compared. The inverted method was superior to the vented method, as shown by a higher coronary perfusion efficiency, more cell debris outflow, higher collagen and elastin content inside the aortic valve, lower DNA content, and better retention of the heart shape after decellularization. To our knowledge, this is the first study to use flow dynamics in a whole heart throughout the decellularization procedure to provide real-time information about the success of the process and the integrity of the vulnerable regions of the matrix. Heart orientation was important in optimizing decellularization efficiency and maintaining extracellular matrix integrity.

STATEMENT OF SIGNIFICANCE

The use of decellularized tissue as a suitable scaffold for engineered tissue has emerged over the past decade as one of the most promising biofabrication platforms. The decellularization process removes all native cells, leaving the natural biopolymers, extracellular matrix materials and native architecture intact. This manuscript describes heart orientation as important in optimizing decellularization efficiency and maintaining extracellular matrix integrity. To our knowledge, this is the first study to assess flow dynamics in a whole heart throughout the decellularization procedure. Our findings compared to currently published methods demonstrate that continuous complex real-time measurements and analyses are required to produce an optimal scaffold for cardiac regeneration.

Numerical simulation of closure performance for neo-aortic valve for arterial switch operation

Background

Modeling neo-aortic valve for arterial switch surgical planning to simulate the neo-aortic valve closure performance.

Methods

We created five geometrical models of neo-aortic valve, namely model A, model B, model C, model D and model E with different size of sinotubular junction or sinus. The nodes at the ends of aorta and left ventricle duct fixed all the degrees of freedom. Transvalvular pressure of normal diastolic blood pressure of 54 mmHg was applied on the neo-aortic valve cusps. The neo-aortic valve closure performance was investigated by the parameters, such as stress of neo-aortic root, variation of neo-aortic valve ring as well as aortic valve cusps contact force in the cardiac diastole.

Results

The maximum stress of the five neo-aortic valves were 96.29, 98.34, 96.28, 98.26, and 90.60 kPa, respectively. Compared among five neo-aortic valve, aortic valve cusps contact forces were changed by 43.33, −10.00% enlarging or narrowing the sinotubular junction by 20% respectively based on the reference model A. The cusps contact forces were changed by 6.67, −23.33% with sinus diameter varying 1.2 times and 0.8 times respectively.

Conclusions

Comparing with stress of healthy adult subjects, the neo-aortic valve of infant creates lower stress. It is evident that enlarging or narrowing the sinotubular junction within a range of 20% can increase or decrease the maximum stress and aortic valve cusps contact force of neo-aortic valve.

Fluid-structure interaction modeling of calcific aortic valve disease using patient-specific three-dimensional calcification scans

Calcific aortic valve disease (CAVD) is characterized by calcification accumulation and thickening of the aortic valve cusps, leading to stenosis. The importance of fluid flow shear stress in the initiation and regulation of CAVD progression is well known and has been studied recently using fluid-structure interaction (FSI) models. While cusp calcifications are three-dimensional (3D) masses, previously published FSI models have represented them as either stiffened or thickened two-dimensional (2D) cusps. This study investigates the hemodynamic effect of these calcifications employing FSI models using 3D patient-specific calcification masses. A new reverse calcification technique (RCT) is used for modeling different stages of calcification growth based on the spatial distribution of calcification density. The RCT is applied to generate the 3D calcification deposits reconstructed from a patient-specific CT scans. Our results showed that consideration of 3D calcification deposits led to both higher fluid shear stresses and unique fluid shear stress distribution on the aortic side of the cusps that may have an impact on the calcification growth rate. However, the flow did not seem to affect the geometry of the calcification during the growth phase.

Fibrous scaffolds for building hearts and heart parts

Extracellular matrix (ECM) structure and biochemistry provide cell-instructive cues that promote and regulate tissue growth, function, and repair. From a structural perspective, the ECM is a scaffold that guides the self-assembly of cells into distinct functional tissues. The ECM promotes the interaction between individual cells and between different cell types, and increases the strength and resilience of the tissue in mechanically dynamic environments. From a biochemical perspective, factors regulating cell-ECM adhesion have been described and diverse aspects of cell-ECM interactions in health and disease continue to be clarified. Natural ECMs therefore provide excellent design rules for tissue engineering scaffolds. The design of regenerative three-dimensional (3D) engineered scaffolds is informed by the target ECM structure, chemistry, and mechanics, to encourage cell infiltration and tissue genesis. This can be achieved using nanofibrous scaffolds composed of polymers that simultaneously recapitulate 3D ECM architecture, high-fidelity nanoscale topography, and bio-activity. Their high porosity, structural anisotropy, and bio-activity present unique advantages for engineering 3D anisotropic tissues. Here, we use the heart as a case study and examine the potential of ECM-inspired nanofibrous scaffolds for cardiac tissue engineering. We asked: Do we know enough to build a heart? To answer this question, we tabulated structural and functional properties of myocardial and valvular tissues for use as design criteria, reviewed nanofiber manufacturing platforms and assessed their capabilities to produce scaffolds that meet our design criteria. Our knowledge of the anatomy and physiology of the heart, as well as our ability to create synthetic ECM scaffolds have advanced to the point that valve replacement with nanofibrous scaffolds may be achieved in the short term, while myocardial repair requires further study in vitro and in vivo.

Imaging analysis of collagen fiber networks in cusps of porcine aortic valves: effect of their local distribution and alignment on valve functionality

The cusps of native aortic valve (AV) are composed of collagen bundles embedded in soft tissue, creating a heterogenic tissue with asymmetric alignment in each cusp. This study compares native collagen fiber networks (CFNs) with a goal to better understand their influence on stress distribution and valve kinematics. Images of CFNs from five porcine tricuspid AVs are analyzed and fluid-structure interaction models are generated based on them. Although the valves had similar overall kinematics, the CFNs had distinctive influence on local mechanics. The regions with dilute CFN are more prone to damage since they are subjected to higher stress magnitudes.

Three-dimensional macro-scale assessment of regional and temporal wall shear stress characteristics on aortic valve leaflets

The aortic valve (AV) achieves unidirectional blood flow between the left ventricle and the aorta. Although hemodynamic stresses have been shown to regulate valvular biology, the native wall shear stress (WSS) experienced by AV leaflets remains largely unknown. The objective of this study was to quantify computationally the macro-scale leaflet WSS environment using fluid-structure interaction modeling. An arbitrary Lagrangian-Eulerian approach was implemented to predict valvular flow and leaflet dynamics in a three-dimensional AV geometry subjected to physiologic transvalvular pressure. Local WSS characteristics were quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI) and temporal shear gradient (TSG). The dominant radial WSS predicted on the leaflets exhibited high amplitude and unidirectionality on the ventricularis (TSM>7.50 dyn/cm(2), OSI < 0.17, TSG>325.54 dyn/cm(2) s) but low amplitude and bidirectionality on the fibrosa (TSM < 2.73 dyn/cm(2), OSI>0.38, TSG < 191.17 dyn/cm(2) s). The radial WSS component computed in the leaflet base, belly and tip demonstrated strong regional variability (ventricularis TSM: 7.50-22.32 dyn/cm(2), fibrosa TSM: 1.26-2.73 dyn/cm(2)). While the circumferential WSS exhibited similar spatially dependent magnitude (ventricularis TSM: 1.41-3.40 dyn/cm(2), fibrosa TSM: 0.42-0.76 dyn/cm(2)) and side-specific amplitude (ventricularis TSG: 101.73-184.43 dyn/cm(2) s, fibrosa TSG: 41.92-54.10 dyn/cm(2) s), its temporal variations were consistently bidirectional (OSI>0.25). This study provides new insights into the role played by leaflet-blood flow interactions in valvular function and critical hemodynamic stress data for the assessment of the hemodynamic theory of AV disease.

Biomechanics and mechanobiology in functional tissue engineering

The field of tissue engineering continues to expand and mature, and several products are now in clinical use, with numerous other preclinical and clinical studies underway. However, specific challenges still remain in the repair or regeneration of tissues that serve a predominantly biomechanical function. Furthermore, it is now clear that mechanobiological interactions between cells and scaffolds can critically influence cell behavior, even in tissues and organs that do not serve an overt biomechanical role. Over the past decade, the field of "functional tissue engineering" has grown as a subfield of tissue engineering to address the challenges and questions on the role of biomechanics and mechanobiology in tissue engineering. Originally posed as a set of principles and guidelines for engineering of load-bearing tissues, functional tissue engineering has grown to encompass several related areas that have proven to have important implications for tissue repair and regeneration. These topics include measurement and modeling of the in vivo biomechanical environment; quantitative analysis of the mechanical properties of native tissues, scaffolds, and repair tissues; development of rationale criteria for the design and assessment of engineered tissues; investigation of the effects biomechanical factors on native and repair tissues, in vivo and in vitro; and development and application of computational models of tissue growth and remodeling. Here we further expand this paradigm and provide examples of the numerous advances in the field over the past decade. Consideration of these principles in the design process will hopefully improve the safety, efficacy, and overall success of engineered tissue replacements.

Computational modeling of cardiac valve function and intervention

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

Fluid structure interaction (FSI) simulation of the left ventricle (LV) during the early filling wave (E-wave), diastasis and atrial contraction wave (A-wave)

In this paper, the hemodynamic characteristics inside a physiologically correct three-dimensional LV model using fluid structure interaction scheme are examined under various heart beat conditions during early filling wave (E-wave), diastasis and atrial contraction wave (A-wave). The time dependent and incompressible viscous fluid, nonlinear viscous fluid and the stress tensor equations are coupled with the full Navier-Stoke's equations together with the Arbitrary Lagrangian-Eulerian and elasticity in the solid domain are used in the analysis. The results are discussed in terms of the variation in the intraventricular pressure, wall shear stress (WSS) and the fluid flow patterns inside the LV model. Moreover, changes in the magnitude of displacements on the LV are also observed during diastole period. The results obtained demonstrate that the magnitude of the intraventricle pressure is found higher in the basal region of the LV during the beginning of the E-wave and A-wave, whereas the Ip is found much higher in the apical region when the flow propagation is in peak E-wave, peak A-wave and diastasis. The magnitude of the pressure is found to be 5.4E2 Pa during the peak E-wave. Additionally, WSS elevates with the rise in the E-wave and A-wave but the magnitude decreases during the diastasis phase. During the peak E-wave, maximum WSS is found to be 5.7 Pa. Subsequent developments, merging and shifting of the vortices are observed throughout the filling wave. Formations of clockwise vortices are evident during the peak E-wave and at the onset of the A-wave, but counter clockwise vortices are found at the end of the diastasis and at the beginning of the A-wave. Moreover, the maximum magnitude of the structural displacement is seen in the ventricle apex with the value of 3.7E-5 m.