A Newton-Krylov method with an approximate analytical Jacobian for implicit solution of Navier-Stokes equations on staggered overset-curvilinear grids with immersed boundaries

A contact model based on the coefficient of restitution for simulations of bio-prosthetic heart valves

A new general contact model is proposed for preventing inter-leaflet penetration of bio-prosthetic heart valves (BHV) at the end of the systole, which has the advantage of applying kinematic constraints directly and creating smooth free edges. At the end of each time step, the impenetrability constraints and momentum exchange between the impacting bodies are applied separately based on the coefficient of restitution. The contact method is implemented in a rotation-free, large deformation, and thin shell finite-element (FE) framework based on loop's subdivision surfaces. A nonlinear, anisotropic material model for a BHV is employed which uses Fung-elastic constitutive laws for in-plane and bending responses, respectively. The contact model is verified and validated against several benchmark problems. For a BHV-specific validation, the computed strains on different regions of a BHV under constant pressure are compared with experimentally measured data. Finally, dynamic simulations of BHV under physiological pressure waveform are performed for symmetrical and asymmetrical fiber orientations incorporating the new contact model and compared with the penalty contact method. The proposed contact model provides the coaptation area of a functioning BHV during the closing phase for both of the fiber orientations. Our results show that fiber orientation affects the dynamic of leaflets during the opening and closing phases. A swirling motion for the BHV with asymmetrical fiber orientation is observed, similar to experimental data. To include the fluid effects, fluid-structure interaction (FSI) simulation of the BHV is performed and compared to the dynamic results.

The Effects of Implantation Orientation of a Bileaflet Mechanical Heart Valve in an Anatomic Left Ventricle-Aorta Configuration

We have performed three-dimensional high-resolution numerical simulations of a bi-leaflet mechanical heart valve (BMHV) implanted at different orientations in an anatomic left ventricle-aorta obtained from magnetic resonance imaging (MRI) of a volunteer. The thoroughly validated overset curvilinear-immersed boundary (overset-CURVIB) fluid-structure interaction (FSI) flow solver is used in which the aorta and LV are discretized with boundary-conforming and non-conforming curvilinear grids, respectively. The motion of the LV wall is prescribed based on a lumped parameter model while the motion of the leaflets are calculated using a strong coupled FSI algorithm enhanced with Aitken convergence technique. We carried out simulations for three valve orientations, which differ from each other by 45 degrees and compared the leaflet motion and flow field for multiple cycles. Our results show reproducible and relatively symmetrical opening for all valve orientations. The presence of small-scale vortical structures after peak systole, cause significant cycle-to-cycle variations in valve kinematics during the closing phase for all valve orientations. Furthermore, our results show that valve orientation does not have a significant effect on the distribution of viscous shear stress in the ascending aorta. Additionally, two different mathematical activation models including linear level of activation and Soares model are used to quantify the platelet activation in the ascending aorta. The results show that the valve orientation does not significantly affect (less than 8%) the total platelet activation in the ascending aorta.

Fluid-structure coupled biotransport processes in aortic valve disease

Biological transport processes near the aortic valve play a crucial role in calcific aortic valve disease initiation and bioprosthetic aortic valve thrombosis. Hemodynamics coupled with the dynamics of the leaflets regulate these transport patterns. Herein, two-way coupled fluid-structure interaction (FSI) simulations of a 2D bicuspid aortic valve and a 3D mechanical heart valve were performed and coupled with various convective mass transport models that represent some of the transport processes in calcification and thrombosis. Namely, five different continuum transport models were developed to study biochemicals that originate from the blood and the leaflets, as well as residence-time and flow stagnation. Low-density lipoprotein (LDL) and platelet activation were studied for their role in calcification and thrombosis, respectively. Coherent structures were identified using vorticity and Lagrangian coherent structures (LCS) for the 2D and 3D models, respectively. A very close connection between vortex structures and biochemical concentration patterns was shown where different vortices controlled the concentration patterns depending on the transport mechanism. Additionally, the relationship between leaflet concentration and wall shear stress was revealed. Our work shows that blood flow physics and coherent structures regulate the flow-mediated biological processes that are involved in aortic valve calcification and thrombosis, and therefore could be used in the design process to optimize heart valve replacement durability.

A hybrid echocardiography-CFD framework for ventricular flow simulations

Image-based CFD is a powerful tool to study cardiovascular flows while 2D echocardiography (echo) is the most widely used noninvasive imaging modality for the diagnosis of heart disease. Here, echo is combined with CFD, that is, an echo-CFD framework, to study ventricular flows. To achieve this, the previous 3D reconstruction from multiple 2D echo at standard cross sections is extended by: (a) reconstructing aortic and mitral valves from 2D echo and closing the left-ventricle (LV) geometry by approximating a superior wall; (b) incorporating the physiological assumption of the fixed apex as a reference (fixed) point in the 3D reconstruction; and (c) incorporating several smoothing algorithms to remove the nonphysical oscillations (ringing) near the basal section. The method is applied to echo from a baseline LV and one after inducing acute myocardial ischemia (AMI). The 3D reconstruction is validated by comparing it against a reference reconstruction from many echo sections while flow simulations are validated against the Doppler ultrasound velocity measurements. The sensitivity study shows that the choice of the smoothing algorithm does not change the flow pattern inside the LV. However, the presence of the mitral valve can significantly change the flow pattern during the diastole phase. In addition, the abnormal shape of a LV with AMI can drastically change the flow during diastole. Furthermore, the hemodynamic energy loss, as an indicator of the LV pumping performance, for different test cases is calculated, which shows a larger energy loss for a LV with AMI compared to the baseline one.

Numerical Simulations of Flow around Copepods: Challenges and Future Directions

Copepods are small aquatic creatures which are abundant in oceans as a major food source for fish, thereby playing a vital role in marine ecology. Because of their role in the food chain, copepods have been subject to intense research through different perspectives from anatomy, form-function biology, to ecology. Numerical simulations can uniquely support such investigations by quantifying: (i) the force and flow generated by different parts of the body, thereby clarify the form-function relation of each part; (ii) the relation between the small-scale flow around animal and the large-scale (e.g., oceanic) flow of its surroundings; and (iii) the flow and its energetics, thereby answering ecological questions, particularly, the three major survival tasks, i.e., feeding, predator avoidance, and mate-finding. Nevertheless, such numerical simulations need to overcome challenges involving complex anatomic shape of copepods, multiple moving appendages, resolving different scales (appendage-, animal- to large-scale). The numerical methods capable of handling such problems and some recent simulations are reviewed. At the end, future developments necessary to simulate copepods from animal- to surrounding-scale are discussed.

A non-dimensional parameter for classification of the flow in intracranial aneurysms. I. Simplified geometries

Non-dimensional parameters are routinely used to classify different flow regimes. We propose a non-dimensional parameter, called Aneurysm number (An), which depends on both geometric and flow characteristics, to classify the flow inside aneurysm-like geometries (sidewalls and bifurcations). The flow inside aneurysm-like geometries can be widely classified into (i) the vortex mode in which a vortex ring is formed and (ii) the cavity mode in which a stationary shear layer acts similar to a moving lid of a lid-driven cavity. In these modes, two competing time scales exist: (a) a transport time scale, T t , which is the time scale to develop a shear layer by transporting a fluid particle across the expansion region, and (b) the vortex formation time scale,T v . Consequently, a relevant non-dimensional parameter is the ratio of these two time scales, which is called Aneurysm number: An = T t /T v . It is hypothesized, based on this definition, that the flow is in the vortex mode if the time required for vortex ring formationT v is less than the transport time T t (An ≳ 1). Otherwise, the flow is in the cavity mode (An ≲ 1). This hypothesis is systematically tested through numerical simulations on simplified geometries and shown to be true through flow visualizations and identification of the main vortex and shear layer. The main vortex is shown to evolve when An ≳ 1 but stationary when An ≲ 1. In fact, it is shown that the flows with An ≲ 1 (cavity mode) are characterized by much smaller fluctuations of wall shear stress and oscillatory shear index relative to flows with An ≳ 1 (vortex mode) because of their quasi-stationary flow pattern (cavity mode) compared to the evolution and breakdown of the formed vortex ring (vortex mode).

A non-dimensional parameter for classification of the flow in intracranial aneurysms. II. Patient-specific geometries

A simple parameter, called the Aneurysm number (An) which is defined as the ratio of transport to vortex time scales, has been shown to classify the flow mode in simplified aneurysm geometries. Our objective is to test the hypothesis that An can classify the flow in patient-specific intracranial aneurysms (IA). Therefore, the definition of this parameter is extended to anatomic geometries by using hydraulic diameter and the length of expansion area in the approximate direction of the flow. The hypothesis is tested using image-based flow simulations in five sidewall and four bifurcation geometries, i.e., if An ≲ 1 (shorter transport time scale), then the fluid is transported across the neck before the vortex could be formed, creating a quasi-stationary shear layer (cavity mode). By contrast, if An ≳ 1 (shorter vortex time scale), a vortex is formed. The results show that if An switches from An ≲ 1 to An ≳ 1, then the flow mode switches from the cavity mode to the vortex mode. However, if An does not switch, then the IAs stay in the same mode. It is also shown that IAs in the cavity mode have significantly lower An, temporal fluctuations of wall shear stress and oscillatory shear index (OSI) compared to the vortex mode (p < 0.01). In addition, OSI correlates with An in each flow mode and with pulsatility index in each IA. This suggests An to be a viable hemodynamic parameter which can be easily calculated without the need for detailed flow measurements/ simulations.