Numerical simulation of the dynamics of a bileaflet prosthetic heart valve using a fluid-structure interaction approach

The main purpose of this study is to reproduce in silico the dynamics of a bileaflet mechanical heart valve (MHV; St Jude Hemodynamic Plus, 27mm characteristic size) by means of a fully implicit fluid-structure interaction (FSI) method, and experimentally validate the results using an ultrafast cinematographic technique. The computational model was constructed to realistically reproduce the boundary condition (72 beats per minute (bpm), cardiac output 4.5l/min) and the geometry of the experimental setup, including the valve housing and the hinge configuration. The simulation was carried out coupling a commercial computational fluid dynamics (CFD) package based on finite-volume method with user-defined code for solving the structural domain, and exploiting the parallel performance of the whole numerical setup. Outputs are leaflets excursion from opening to closure and the fluid dynamics through the valve. Results put in evidence a favorable comparison between the computed and the experimental data: the model captures the main features of the leaflet motion during the systole. The use of parallel computing drastically limited the computational costs, showing a linear scaling on 16 processors (despite the massive use of user-defined subroutines to manage the FSI process). The favorable agreement obtained between in vitro and in silico results of the leaflet displacements confirms the consistency of the numerical method used, and candidates the application of FSI models to become a major tool to optimize the MHV design and eventually provides useful information to surgeons.     

Validation of a Fluid–Structure Interaction Model of a Heart Valve using the Dynamic Mesh Method in Fluent

Simulations of coupled problems such as fluid–structure interaction (FSI) are becoming more and more important for engineering purposes. This is particularly true when modeling the aortic valve, where the FSI between the blood and the valve determines the valve movement and the valvular hemodynamics. Nevertheless only a few studies are focusing on the opening and closing behavior during the ejection phase (systole). In this paper, we present the validation of a FSI model using the dynamic mesh method of Fluent for the two-dimensional (2D) simulation of mechanical heart valves during the ejection phase of the cardiac cycle. The FSI model is successfully validated by comparing simulation results to experimental data obtained from in vitro studies using a CCD camera.     

Evaluation of a transient,simultaneous, arbitrary Lagrange–Euler based multi-physics method for simulating the mitral heart valve

A transient multi-physics model of the mitral heart valve has been developed, which allows simultaneous calculation of fluid flow and structural deformation. A recently developed contact method has been applied to enable simulation of systole (the stage when blood pressure is elevated within the heart to pump blood to the body). The geometry was simplified to represent the mitral valve within the heart walls in two dimensions. Only the mitral valve undergoes deformation. A moving arbitrary Lagrange–Euler mesh is used to allow true fluid–structure interaction (FSI). The FSI model requires blood flow to induce valve closure by inducing strains in the region of 10–20%. Model predictions were found to be consistent with existing literature and will undergo further development.     

Fluid-Structure Interaction Simulation of Prosthetic Aortic Valves: Comparison between Immersed Boundary and Arbitrary Lagrangian-Eulerian Techniques for the Mesh Representation

In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations’ outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results.     

Systolic fluid–structure interaction model of the congenitally bicuspid aortic valve: assessment of modelling requirements

A transient fluid–structure interaction (FSI) model of a congenitally bicuspid aortic valve has been developed which allows simultaneous calculation of fluid flow and structural deformation. The valve is modelled during the systolic phase (the stage when blood pressure is elevated within the heart to pump blood to the body). The geometry was simplified to represent the bicuspid aortic valve in two dimensions. A congenital bicuspid valve is compared within the aortic root only and within the aortic arch. Symmetric and asymmetric cusps were simulated, along with differences in mechanical properties. A moving arbitrary Lagrange–Euler mesh was used to allow FSI. The FSI model requires blood flow to induce valve opening and induced strains in the region of 10%. It was determined that bicuspid aortic valve simulations required the inclusion of the ascending aorta and aortic arch. The flow patterns developed were sensitive to cusp asymmetry and differences in mechanical properties. Stiffening of the valve amplified peak velocities, and recirculation which developed in the ascending aorta. Model predictions demonstrate the need to take into account the category, including any existing cusp asymmetry, of a congenital bicuspid aortic valve when simulating its fluid flow and mechanics.     

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

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

Validation of a fluid-structure interaction model of a heart valve using the dynamic mesh method in fluent

Simulations of coupled problems such as fluid-structure interaction (FSI) are becoming more and more important for engineering purposes. This is particularly true when modeling the aortic valve, where the FSI between the blood and the valve determines the valve movement and the valvular hemodynamics. Nevertheless only a few studies are focusing on the opening and closing behavior during the ejection phase (systole). In this paper, we present the validation of a FSI model using the dynamic mesh method of Fluent for the two-dimensional (2D) simulation of mechanical heart valves during the ejection phase of the cardiac cycle. The FSI model is successfully validated by comparing simulation results to experimental data obtained from in vitro studies using a CCD camera.     

Effects of fluid structure interaction in a three dimensional model of the spinal subarachnoid space

It is unknown whether spinal cord motion has a significant effect on cerebrospinal fluid (CSF) pressure and therefore the importance of including fluid structure interaction (FSI) in computational fluid dynamics models (CFD) of the spinal subarachnoid space (SAS) is unclear. This study aims to determine the effects of FSI on CSF pressure and spinal cord motion in a normal and in a stenosis model of the SAS. A three-dimensional patient specific model of the SAS and spinal cord were constructed from MR anatomical images and CSF flow rate measurements obtained from a healthy human being. The area of SAS at spinal level T4 was constricted by 20% to represent the stenosis model. FSI simulations in both models were performed by running ANSYS CFX and ANSYS Mechanical in tandem. Results from this study show that the effect of FSI on CSF pressure is only about 1% in both the normal and stenosis models and therefore show that FSI has a negligible effect on CSF pressure.     

Does clinical data quality affect fluid-structure interaction simulations of patient-specific stenotic aortic valve models?

Numerical models are increasingly used in the cardiovascular field to reproduce, study and improve devices and clinical treatments. The recent literature involves a number of patient-specific models replicating the transcatheter aortic valve implantation procedure, a minimally invasive treatment for high-risk patients with aortic diseases. The representation of the actual patient’s condition with truthful anatomy, materials and working conditions is the first step toward the simulation of the clinical procedure.The aim of this work is to quantify how the quality of routine clinical data, from which the patient-specific models are built, affects the outputs of the numerical models representing the pathological condition of stenotic aortic valve.Seven fluid–structure interaction (FSI) simulations were performed, completed with a sensitivity analysis on patient-specific reconstructed geometries and boundary conditions. The structural parts of the models consisted of the aortic root, native tri-leaflets valve and calcifications. Ventricular and aortic pressure curves were applied to the fluid domain.The differences between clinical data and numerical results for the aortic valve area were less than 2% but reached 12% when boundary conditions and geometries were changed. The difference in the aortic stenosis jet velocity between measured and simulated values was less than 11% reaching 27% when the geometry was changed. The CT slice thickness was found to be the most sensitive parameter on the presented FSI numerical model.In conclusion, the results showed that the segmentation and reconstruction phases need to be carefully performed to obtain a truthful patient-specific domain to be used in FSI analyses.     

A mathematical model of aortic valve vibration

It has been shown experimentally that the second heart sound is produced by diastolic vibrations of the closed aortic valve. In the present paper a mathematical model of this vibration is developed from first principles. The model assumes a one dimensional but non-linear fluid behavior. The problem is coupled through a non-linear, planar valve. A solution is obtained using the method of characteristics developed in finite difference form. The resulting valve frequency and amplitude are in good agreement with patient data. The model predicts a strong dependency of response on the valve forcing function and valve stiffness; and a weaker dependency of response on valve mass.     

An integrated fluid–structure interaction and thrombosis model for type B aortic dissection

False lumen thrombosis (FLT) in type B aortic dissection has been associated with the progression of dissection and treatment outcome. Existing computational models mostly assume rigid wall behavior which ignores the effect of flap motion on flow and thrombus formation within the FL. In this study, we have combined a fully coupled fluid–structure interaction (FSI) approach with a shear-driven thrombosis model described by a series of convection–diffusion reaction equations. The integrated FSI-thrombosis model has been applied to an idealized dissection geometry to investigate the interaction between vessel wall motion and growing thrombus. Our simulation results show that wall compliance and flap motion can influence the progression of FLT. The main difference between the rigid and FSI models is the continuous development of vortices near the tears caused by drastic flap motion up to 4.45 mm. Flap-induced high shear stress and shear rates around tears help to transport activated platelets further to the neighboring region, thus speeding up thrombus formation during the accelerated phase in the FSI models. Reducing flap mobility by increasing the Young’s modulus of the flap slows down the thrombus growth. Compared to the rigid model, the predicted thrombus volume is 25% larger using the FSI-thrombosis model with a relatively mobile flap. Furthermore, our FSI-thrombosis model can capture the gradual effect of thrombus growth on the flow field, leading to flow obstruction in the FL, increased blood viscosity and reduced flap motion. This model is a step closer toward simulating realistic thrombus growth in aortic dissection, by taking into account the effect of intimal flap and vessel wall motion.

     

Evaluation of dentinal fluid flow behaviours: a fluid-structure interaction simulation

This study uses the fluid-structure interaction (FSI) method to investigate the fluid flow in dental pulp. First, the FSI method is used for the biomechanical simulation of dental intrapulpal responses during force loading (50, 100 and 150 N) on a tooth. The results are validated by comparison with experimental outcomes. Second, the FSI method is used to investigate an intact tooth subjected to a mechanical stimulus during loading at various loading rates. Force loading (0–100 N) is applied gradually to an intact tooth surface with loading rates of 125, 62.5, 25 and 12.5 N/s, respectively, and the fluid flow changes in the pulp are evaluated. FSI analysis is found to be suitable for examining intrapulpal biomechanics. An external force applied to a tooth with a low loading rate leads to a low fluid flow velocity in the pulp chamber, thus avoiding tooth pain.     

The effect of implantation orientation of a bileaflet mechanical heart valve on kinematics and hemodynamics in an anatomic aorta

We carry out three-dimensional high-resolution numerical simulations of a bileaflet mechanical heart valve under physiologic pulsatile flow conditions implanted at different orientations in an anatomic aorta obtained from magnetic resonance imaging (MRI) of a volunteer. We use the extensively validated for heart valve flow curvilinear-immersed boundary (CURVIB) fluid-structure interaction (FSI) solver in which the empty aorta is discretized with a curvilinear, aorta-conforming grid while the valve is handled as an immersed boundary. The motion of the valve leaflets are calculated through a strongly coupled FSI algorithm implemented in conjunction with the Aitken convergence acceleration technique. We perform simulations for three valve orientations, which differ from each other by 45 deg and compare the results in terms of leaflet motion and flow field. We show that the valve implanted symmetrically relative to the symmetry plane of the ascending aorta curvature exhibits the smallest overall asymmetry in the motion of its two leaflets and lowest rebound during closure. Consequently, we hypothesize that this orientation is beneficial to reduce the chance of intermittent regurgitation. Furthermore, we find that the valve orientation does not significantly affect the shear stress distribution in the aortic lumen, which is in agreement with previous studies.     

Numerical Simulation of Callus Healing for Optimization of Fracture Fixation Stiffness

The stiffness of fracture fixation devices together with musculoskeletal loading defines the mechanical environment within a long bone fracture, and can be quantified by the interfragmentary movement. In vivo results suggested that this can have acceleratory or inhibitory influences, depending on direction and magnitude of motion, indicating that some complications in fracture treatment could be avoided by optimizing the fixation stiffness. However, general statements are difficult to make due to the limited number of experimental findings. The aim of this study was therefore to numerically investigate healing outcomes under various combinations of shear and axial fixation stiffness, and to detect the optimal configuration. A calibrated and established numerical model was used to predict fracture healing for numerous combinations of axial and shear fixation stiffness under physiological, superimposed, axial compressive and translational shear loading in sheep. Characteristic maps of healing outcome versus fixation stiffness (axial and shear) were created. The results suggest that delayed healing of 3 mm transversal fracture gaps will occur for highly flexible or very rigid axial fixation, which was corroborated by in vivo findings. The optimal fixation stiffness for ovine long bone fractures was predicted to be 1000–2500 N/mm in the axial and >300 N/mm in the shear direction. In summary, an optimized, moderate axial stiffness together with certain shear stiffness enhances fracture healing processes. The negative influence of one improper stiffness can be compensated by adjustment of the stiffness in the other direction.     

Probing the effect of morphology on lymphatic valve dynamic function

The lymphatic system is vital to the circulatory and immune systems, performing a range of important functions such as transport of interstitial fluid, fatty acid, and immune cells. Lymphatic vessels are composed of contractile walls and lymphatic valves, allowing them to pump lymph against adverse pressure gradients and to prevent backflow. Despite the importance of the lymphatic system, the contribution of mechanical and geometric changes of lymphatic valves and vessels in pathologies of lymphatic dysfunction, such as lymphedema, is not well understood. We develop a fully coupled fluid–solid, three-dimensional computational model to interrogate the various parameters thought to influence valve behavior and the consequences of these changes to overall lymphatic function. A lattice Boltzmann model is used to simulate the lymph, while a lattice spring model is used to model the mechanics of lymphatic valves. Lymphatic valve functions such as enabling lymph flow and preventing backflow under varied lymphatic valve geometries and mechanical properties are investigated to provide an understanding of the function of lymphatic vessels and valves. The simulations indicate that lymphatic valve function is optimized when valves are of low aspect ratio and bending stiffness, so long as these parameters are maintained at high enough values to allow for proper valve closing. This suggests that valve stiffening could have a profound effect on overall lymphatic pumping performance. Furthermore, dynamic valve simulations showed that this model captures the delayed response of lymphatic valves to dynamic flow conditions, which is an essential feature of valve operation. Thus, our model enhances our understanding of how lymphatic pathologies, specifically those exhibiting abnormal valve morphologies such as has been suggested to occur in cases of primary lymphedema, can lead to lymphatic dysfunctions.     

Anatomy of the oral valve in nymphalid butterflies and a functional model for fluid uptake in Lepidoptera

The food canal of the proboscis of Lepidoptera serves for the uptake of nutrient fluids and the discharge of saliva. A valve was discovered at the entrance to the sucking pump in the head that separates these countercurrent flows in nymphalid butterflies. Three species of Nymphalidae were examined by dissections and light microscopic serial semithin sections. The sucking pump is a unit composed of three structures: (1) the oral valve, which is a projection of the epipharynx extending into the anterior cibarial lumen, (2) the expandable lumen, and (3) the posterior sphincter valve which controls influx into the oesophagus. Based on the microanatomical results, a functional model is presented to account for the uptake and swallowing of fluids and for the control of the salivary flow into the food canal of the proboscis. Dilator muscles of the sucking pump expand the lumen by pulling on the muscular dorso-anterior side. This opens the oral valve and fluid can be drawn into the lumen from the food canal of the proboscis. Circular compressor muscles which attach to both sides of the sclerotized ventro-posterior wall of the sucking pump reduce the size of the lumen; passively they close the oral valve and press fluid through the relaxed posterior sphincter opening into the oesophagus. According to this model saliva can be discharged into the food canal during the swallowing phase. The oral valve and pumping unit are similar in all studied species despite the fact that saliva presumably plays a special role in the derived pollen-feeding behaviour of one of them, viz. Heliconius melpomene.     

Measurement of a spinal motion segment stiffness matrix

The six-degrees-of-freedom elastic behavior of spinal motion segments can be approximated by a stiffness matrix. A method is described to measure this stiffness matrix directly with the motion segment held under physiological conditions of axial preload and in an isotonic fluid bath by measuring the forces and moments associated with each of the six orthogonal translations and rotations. The stiffness matrix was obtained from the load-displacement measurements by linear least squares assuming a symmetric matrix. Results from a pig lumbar spinal motion segment in an isotonic bath, with and without a 500 N axial preload, showed a large stiffening effect with axial preload.     

Velocity profiles in the human ductus venosus: a numerical fluid structure interaction study

The veins distributing oxygenated blood from the placenta to the fetal body have been given much attention in clinical Doppler velocimetry studies, in particular the ductus venosus. The ductus venosus is embedded in the left liver lobe and connects the intra-abdominal portion of the umbilical vein (IUV) directly to the inferior vena cava, such that oxygenated blood can bypass the liver and flow directly to the fetal heart. In the current work, we have developed a mathematical model to assist the clinical assessment of volumetric flow rate at the inlet of the ductus venosus. With a robust estimate of the velocity profile shape coefficient (VC), the volumetric flow rate may be estimated as the product of the time-averaged cross-sectional area, the time-averaged cross-sectional maximum velocity and the VC. The time average quantities may be obtained from Doppler ultrasound measurements, whereas the VC may be estimated from numerical simulations. The mathematical model employs a 3D fluid structure interaction model of the bifurcation formed by the IUV, the ductus venosus and the left portal vein. Furthermore, the amniotic portion of the umbilical vein, the right liver lobe and the inferior vena cava were incorporated as lumped model boundary conditions for the fluid structure interaction model. A hyperelastic material is used to model the structural response of the vessel walls, based on recently available experimental data for the human IUV and ductus venous. A parametric study was constructed to investigate the VC at the ductus venosus inlet, based on a reference case for a human fetus at 36 weeks of gestation. The VC was found to be \(0.687\,\pm \,0.023\) (Mean \(\pm \) SD of parametric case study), which confirms previous studies in the literature on the VC at the ductus venosus inlet. Additionally, CFD simulations with rigid walls were performed on a subsection of the parametric case study, and only minor changes in the predicted VCs were observed compared to the FSI cases. In conclusion, the presented mathematical model is a promising tool for the assessment of ductus venosus Doppler velocimetry.     

Leaching losses of nitrogen from a clay soil under grass and cereal crops in Finland

Summary A 16-plot experimental field was established in 1975 on a clay soil in Jokioinen, Finland. The water discharge through tile drains was measured and its ammonium and nitrate N contents determined for each plot separately. The surface runoff was also measured and analysed. The annual runoff and the N leached from the surface of moderately fertilized (100 kg/ha/y N) cereal plots varied during 1976–1982 from 21 to 301 mm and from 2 to 7 kg/ha, respectively. The discharge of water and leaching of N through subdrains varied from 65 to 225 mm and from 1 to 38 kg/ha, respectively. The highest leaching was probably caused by a previous fallow. The annual N uptake by the crop varied between 41 and 122 kg/ha.Of the fertilizer-N used for cereals, 20% of that applied in the autumn was lost, but only 1 to 4 per cent was lost when applied in the spring. There was much less N leaching from ley than from barley plots, although the former was given twice as much N. The rate of N fertilization had only a very slight effect on N leaching from both ley and barley plots.The results were compared with those obtained in lysimeters filled with clay, silt, sand and peat soils. No definite conclusions can be drawn because the lysimeter experimental data are only for the first year.