Computing the ankle-brachial index with parallel computational fluid dynamics

The ankle-brachial index (ABI), a ratio of arterial blood pressure in the ankles and upper arms, is used to diagnose and monitor circulatory conditions such as coarctation of the aorta and peripheral artery disease. Computational simulations of the ABI can potentially determine the parameters that produce an ABI indicative of ischemia or other abnormalities in blood flow. However, 0- and 1-D computational methods are limited in describing a 3-D patient-derived geometry. Thus, we present a massively parallel framework for computational fluid dynamics (CFD) simulations in the full arterial system. Using the lattice Boltzmann method to solve the Navier–Stokes equations, we employ highly parallelized and scalable methods to generate the simulation domain and efficiently distribute the computational load among processors. For the first time, we compute an ABI with 3-D CFD. In this proof-of-concept study, we investigate the dependence of ABI on the presence of stenoses, or narrowed regions of the arteries, by directly modifying the arterial geometry. As a result, our framework enables the computation a hemodynamic factor characterizing flow at the scale of the full arterial system, in a manner that is extensible to patient-specific imaging data and holds potential for treatment planning.     

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

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

Numerical modeling of hemodynamics scenarios of patient-specific coronary artery bypass grafts

A fast computational framework is devised to the study of several configurations of patient-specific coronary artery bypass grafts. This is especially useful to perform a sensitivity analysis of the hemodynamics for different flow conditions occurring in native coronary arteries and bypass grafts, the investigation of the progression of the coronary artery disease and the choice of the most appropriate surgical procedure. A complete pipeline, from the acquisition of patient-specific medical images to fast parameterized computational simulations, is proposed. Complex surgical configurations employed in the clinical practice, such as Y-grafts and sequential grafts, are studied. A virtual surgery platform based on model reduction of unsteady Navier–Stokes equations for blood dynamics is proposed to carry out sensitivity analyses in a very rapid and reliable way. A specialized geometrical parameterization is employed to compare the effect of stenosis and anastomosis variation on the outcome of the surgery in several relevant cases.     

Meniscal Tear Film Fluid Dynamics Near Marx’s Line

Extensive studies have explored the dynamics of the ocular surface fluid, though theoretical investigations are typically limited to the use of the lubrication approximation, which is not guaranteed to be uniformly valid a-priori throughout the tear meniscus. However, resolving tear film behaviour within the meniscus and especially its apices is required to characterise the flow dynamics where the tear film is especially thin, and thus most susceptible to evaporatively induced hyperosmolarity and subsequent epithelial damage. Hence, we have explored the accuracy of the standard lubrication approximation for the tear film by explicit comparisons with the 2D Navier–Stokes model, considering both stationary and moving eyelids. Our results demonstrate that the lubrication model is qualitatively accurate except in the vicinity of the eyelids. In particular, and in contrast to lubrication theory, the solution of the full Navier–Stokes equations predict a distinct absence of fluid flow, and thus convective mixing in the region adjacent to the tear film contact line. These observations not only support emergent hypotheses concerning the formation of Marx’s line, a region of epithelial cell staining adjacent to the contact line on the eyelid, but also enhance our understanding of the pathophysiological consequences of the flow profile near the tear film contact line.     

Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.     

Investigation of left heart flow using a numerical correlation to model heart wall motion

A three-dimensional computational fluid dynamics (CFD) method has been developed to model the flow in the left heart including atrium and ventricle. Since time resolution of the medical scans does not fit the requirements of the CFD calculations, the main challenge in a numerical simulation of heart chambers is wall motion modeling. This study employs a novel three-dimensional approximation scheme to correlate the wall boundary and grid movement in systole and diastole. It uses a geometry extracted from medical images in the literature and deformed based on the reported flow rates. The opening and closing of the mitral (MV) and the aortic valve (AV) considered as simultaneous events. Unstructured tetragonal grids were used for the meshing of the domain. The calculation was performed by a Navier–Stokes solver using the arbitrary Lagrange–Euler (ALE) formulation. Results show that the proposed correlation for the wall motion could predict the main features of heart flows.     

A computational fluid dynamics (CFD) based method for assessing the hydrodynamic impact of animal borne data loggers on host marine mammals

Animal‐borne data loggers (ABDLs) or “tags” are regularly used to elucidate animal ecology and physiology, but current literature highlights the need to assess associated deleterious impacts including increased resistive force to motion. Previous studies have used computational fluid dynamics (CFD) to estimate this impact, but many suffer limitations (e.g., inaccurate turbulence modeling, neglecting boundary layer transition, neglecting added mass effects, and analyzing the ABDL in isolation from the animal). A novel CFD‐based method is presented in which a “tag impact envelope” is defined utilizing simulations with and without transition modeling to define upper and lower drag limits, respectively, and added mass coefficients are found via simulations with sinusoidally varying inlet velocity, with modified Navier‐Stokes conservation of momentum equations enforcing a shift to the animal's noninertial reference frame. The method generates coefficients for calculating total resistive force for any velocity and acceleration combination, and is validated against theory for a prolate spheroid. An example case shows ABDL drag impact on a harp seal of 11.21%–16.24%, with negligible influence on added mass. By considering the effects of added mass and boundary layer transition, the approach presented is an enhancement to the CFD‐based ABDL impact assessment methods previously applied by researchers.     

Viscous flow simulation in a stenosis model using discrete particle dynamics: a comparison between DPD and CFD

Flow and stresses induced by blood flow acting on the blood cellular constituents can be represented to a certain extent by a continuum mechanics approach down to the order of the?μm level. However, the molecular effects of, e.g., adhesion/aggregation bonds of blood clotting can be on the order of nm. The coupling of the disparate length and timescales between such molecular levels and macroscopic transport represents a major computational challenge. To address this challenge, a multiscale numerical approach based on discrete particle dynamics (DPD) methodology derived from molecular dynamics (MD) principles is proposed. The feasibility of the approach was firstly tested for its ability to simulate viscous flow conditions. Simulations were conducted in low Reynolds numbers flows (Re = 25–33) through constricted tubes representing blood vessels with various degrees of stenosis. Multiple discrete particles interacting with each other were simulated, with 1.24–1.36 million particles representing the flow domain and 0.4 million particles representing the vessel wall. The computation was carried out on the massive parallel supercomputer NY BlueGene/L employing NAMD-a parallel MD package for high performance computing (HPC). Typical recirculation zones were formed distal to the stenoses. The velocity profiles and recirculation zones were in excellent agreement with computational fluid dynamics (CFD) 3D Navier–Stokes viscous fluid flow simulations and with classic numerical and experimental results by YC Fung in constricted tubes. This feasibility analysis demonstrates the potential of a methodology that widely departs from a continuum approach to simulate multiscale phenomena such as flow induced blood clotting.     

Multiscale bio-chemo-mechanical model of intimal hyperplasia

We consider a computational multiscale framework of a bio-chemo-mechanical model for intimal hyperplasia. With respect to existing models, we investigate the interactions between hemodynamics, cellular dynamics and biochemistry on the development of the pathology. Within the arterial wall, we propose a mathematical model consisting of kinetic differential equations for key vascular cell types, collagen and growth factors. The luminal hemodynamics is modeled with the Navier–Stokes equations. Coupling hypothesis among time and space scales are proposed to build a tractable modeling of such a complex multifactorial and multiscale pathology. A one-dimensional numerical test-case is presented for validation by comparing the results of the framework with experiments at short and long timescales. Our model permits to capture many cellular phenomena which have a central role in the physiopathology of intimal hyperplasia. Results are quantitatively and qualitatively consistent with experimental findings at both short and long timescales.

     

Patient-specific simulation of a stentless aortic valve implant: the impact of fibres on leaflet performance

In some cases of aortic valve leaflet disease, the implant of a stentless biological prosthesis represents an excellent option for aortic valve replacement (AVR). In particular, if compared with the implant of mechanical valves, it provides a more physiological haemodynamic performance and a reduced thrombogeneticity, avoiding the use of anticoagulants. The clinical outcomes of AVR are strongly dependent on an appropriate choice of both prosthesis size and replacement technique, which is, at present, strictly related to surgeon's experience and skill. This represents the motivation for patient-specific finite element analysis able to virtually reproduce stentless valve implantation. With the aim of performing reliable patient-specific simulations, we remark that, on the one hand, it is not well established in the literature whether bioprosthetic leaflet tissue is isotropic or anisotropic; on the other hand, it is of fundamental importance to incorporate an accurate material model to realistically predict post-operative performance. Within this framework, using a novel computational methodology to simulate stentless valve implantation, we test the impact of using different material models on both the stress pattern and post-operative coaptation parameters (i.e. coaptation area, length and height). As expected, the simulation results suggest that the material properties of the valve leaflets affect significantly the post-operative prosthesis performance.     

Effects of vessel compliance on flow pattern in porcine epicardial right coronary arterial tree

The compliance of the vessel wall affects hemodynamic parameters which may alter the permeability of the vessel wall. Based on experimental measurements, the present study established a finite element (FE) model in the proximal elastic vessel segments of epicardial right coronary arterial (RCA) tree obtained from computed tomography. The motion of elastic vessel wall was measured by an impedance catheter and the inlet boundary condition was measured by an ultrasound flow probe. The Galerkin FE method was used to solve the Navier–Stokes and Continuity equations, where the convective term in the Navier–Stokes equation was changed in the arbitrary Lagrangian–Eulerian (ALE) framework to incorporate the motion due to vessel compliance. Various hemodynamic parameters (e.g., wall shear stress—WSS, WSS spatial gradient—WSSG, oscillatory shear index—OSI) were analyzed in the model. The motion due to vessel compliance affects the time-averaged WSSG more strongly than WSS at bifurcations. The decrease of WSSG at flow divider in elastic bifurcations, as compared to rigid bifurcations, implies that the vessel compliance decreases the permeability of vessel wall and may be atheroprotective. The model can be used to predict coronary flow pattern in subject-specific anatomy as determined by noninvasive imaging.     

Modeling aqueous humor collection from the human eye

Glaucoma is a common cause of blindness. Studies of this disease can involve collection of aqueous humor (AH) fluid from eyes of patients undergoing surgery, which involves aspirating a small fluid volume from the anterior region of the eye through a fine-bore needle. Unfortunately, the composition of the AH is spatially non-uniform in the eye, and thus the composition of the aspirated fluid is uncertain. Our goal was to numerically simulate the AH aspiration process to determine where the aspirated fluid was being collected from and thus gain insight into the composition of the collected fluid.A 3D computational model of the anterior region of the human eye was created and the Navier–Stokes equations were numerically solved during the aspiration process for a set of typical (baseline) conditions: 40 μl aspirated volume and needle placement in the central anterior chamber. We also ran variations of this baseline simulation.The main finding was that the aspirated fluid comes from a very localized region around the needle tip, so that for typical conditions, almost no aspirated fluid is withdrawn from the angle region of the anterior chamber. This is important because the AH in this angle region is protein-rich and directly interacts with the tissues that control fluid drainage from the eye. Recommendations for standardizing aspiration conditions are given.     

Assessing the impact of including leaflets in the simulation of TAVI deployment into a patient-specific aortic root

Computational simulation of transcatheter aortic valve implantation (TAVI) device deployment presents a significant challenge over and above similar simulations for percutaneous coronary intervention due to the presence of prosthetic leaflets. In light of the complexity of these leaflets, simulations have been performed to assess the effect of including the leaflets in a complete model of a balloon-expandable TAVI device when deployed in a patient-specific aortic root. Using an average model discrepancy metric, the average frame positions (with and without the leaflets) are shown to vary by 0.236% of the expanded frame diameter (26 mm). This relatively small discrepancy leads to the conclusion that for a broad range of replacement valve studies, including new frame configurations and designs, patient-specific assessment of apposition, paravalvular leakage and tissue stress, modelling of the prosthetic leaflets is likely to have a marginal effect on the results     

In vitro assessment of a combined radiofrequency ablation and cryo-anchoring catheter for treatment of mitral valve prolapse

Percutaneous approaches to mitral valve repair are an attractive alternative to surgical repair or replacement. Radiofrequency ablation has the potential to approximate surgical leaflet resection by using resistive heating to reduce leaflet size, and cryogenic temperatures on a percutaneous catheter can potentially be used to reversibly adhere to moving mitral valve leaflets for reliable application of radiofrequency energy. We tested a combined cryo-anchoring and radiofrequency ablation catheter using excised porcine mitral valves placed in a left heart flow loop capable of reproducing physiologic pressure and flow waveforms. Transmitral flow and pressure were monitored during the cryo-anchoring procedure and compared to baseline flow conditions, and the extent of radiofrequency energy delivery to the mitral valve was assessed post-treatment. Long term durability of radiofrequency ablation treatment was assessed using statically treated leaflets placed in a stretch bioreactor for four weeks. Transmitral flow and pressure waveforms were largely unaltered during cryo-anchoring. Parameter fitting to mechanical data from leaflets treated with radiofrequency ablation and cryo-anchoring revealed significant mechanical differences from untreated leaflets, demonstrating successful ablation of mitral valves in a hemodynamic environment. Picrosirius red staining showed clear differences in morphology and collagen birefringence between treated and untreated leaflets. The durability study indicated that statically treated leaflets did not significantly change size or mechanics over four weeks. A cryo-anchoring and radiofrequency ablation catheter can adhere to and ablate mitral valve leaflets in a physiologic hemodynamic environment, providing a possible percutaneous alternative to surgical leaflet resection of mitral valve tissue.     

Helical flows and asymmetry of blood jet in dilated ascending aorta with normally functioning bicuspid valve

Bicuspid aortic valve (BAV) is associated with aortic dilatation and aneurysm. Several studies evidenced an eccentric systolic flow in ascending aorta associated with increased wall shear stresses (WSS) and the occurrence of an helical systolic flow. This study seeks to elucidate the connections between jet asymmetry and helical flow in patients with normally functioning BAV and dilated ascending aorta. We performed a computational parametric study by varying, for a patient-specific geometry, the valve area and the flow rate entering the aorta and drawing also a tricuspid valve (TAV). We considered also phase-contrast magnetic resonance imaging of four BAV and TAV patients. Measurement of normalized flow asymmetry index, systolic WSS and of a new index (positive helix fraction, PHF) quantifying the presence of a single a single helical flow were performed. In our computation, BAV cases featured higher values of all indices with respect to TAV in both numerical and imaged-based results. Moreover, all indices increased with decreasing valve area and/or with increasing flow rate. This allowed to separate the BAV and TAV cases with respect to the jet asymmetry, WSS localization and helical flow. Interestingly, these results were obtained without modeling the leaflets.     

Fundamental mechanics of aortic heart valve closure

Stresses in a prosthetic heart valve at closure are determined by its geometrical and structural characteristics, by the mechanical support environment, and by the momentum of the valve leaflets or occluder and of the blood at the instant of closure. The mass of blood to be arrested is significantly greater than that of the leaflets or occluder, and is therefore likely to dominate the closure impulse. The kinetic energy of the blood must be transduced into potential energy in the structural components (valve leaflets, aortic root and aorta). This paper presents a methodology for computation and parameterisation of the blood momentum associated with a valve in the aortic position. It is suggested that the influence of physiological parameters, such as systolic waveform and systemic impedance, on the closure characteristics can be investigated based on the fluid dynamic implications. Detailed results are presented for a single leaflet mechanical valve (Bjork-Shiley 60 degrees Convexo-Concave). It is demonstrated that a simple analytical method can yield results that might be adequate for the purposes of valve design.     

An anisotropic constitutive model for immersogeometric fluid–structure interaction analysis of bioprosthetic heart valves

This paper considers an anisotropic hyperelastic soft tissue model, originally proposed for native valve tissue and referred to herein as the Lee–Sacks model, in an isogeometric thin shell analysis framework that can be readily combined with immersogeometric fluid–structure interaction (FSI) analysis for high-fidelity simulations of bioprosthetic heart valves (BHVs) interacting with blood flow. We find that the Lee–Sacks model is well-suited to reproduce the anisotropic stress–strain behavior of the cross-linked bovine pericardial tissues that are commonly used in BHVs. An automated procedure for parameter selection leads to an instance of the Lee–Sacks model that matches biaxial stress–strain data from the literature more closely, over a wider range of strains, than other soft tissue models. The relative simplicity of the Lee–Sacks model is attractive for computationally-demanding applications such as FSI analysis and we use the model to demonstrate how the presence and direction of material anisotropy affect the FSI dynamics of BHV leaflets.     

Automatic noninvasive measurement of systolic blood pressure using photoplethysmography

Background

Systolic blood flow has been simulated in the abdominal aorta and the superior mesenteric artery. The simulations were carried out using two different computational hemodynamic methods: the finite element method to solve the Navier Stokes equations and the lattice Boltzmann method.

Results

We have validated the lattice Boltzmann method for systolic flows by comparing the veloCity and pressure profiles of simulated blood flow between methods. We have also analyzed flow-specific characteristics such as the formation of a vortex at curvatures and traces of flow.

Conclusion

The lattice Boltzmann Method is as accurate as a Navier Stokes solver for computing complex blood flows. As such it is a good alternative for computational hemodynamics, certainly in situation where coupling to other models is required.