Numerical investigation of oxygen mass transfer in a helical-type artery bypass graft

On the identification of hypoxic regions in subject-specific cerebral vasculature by combined CFD/MRI

A long-time exposure to lack of oxygen (hypoxia) in some regions of the cerebrovascular system is believed to be one of the causes of cerebral neurological diseases. In the present study, we show how a combination of magnetic resonance imaging (MRI) and computational fluid dynamics (CFD) can provide a non-invasive alternative for studying blood flow and transport of oxygen within the cerebral vasculature. We perform computer simulations of oxygen mass transfer in the subject-specific geometry of the circle of Willis. The computational domain and boundary conditions are based on four-dimensional (4D)-flow MRI measurements. Two different oxygen mass transfer models are considered: passive (where oxygen is treated as a dilute chemical species in plasma) and active (where oxygen is bonded to haemoglobin) models. We show that neglecting haemoglobin transport results in a significant underestimation of the arterial wall mass transfer of oxygen. We identified the hypoxic regions along the arterial walls by introducing the critical thresholds that are obtained by comparison of the estimated range of Damköhler number (Da ⊂ 〈9; 57〉) with the local Sherwood number. Finally, we recommend additional validations of the combined MRI/CFD approach proposed here for larger groups of subject- or patient-specific brain vasculature systems.

Hemodynamics and Oxygen Transport through Pararenal Aortic Aneurysm Treated with Multilayer Stent: A Numerical Study

BACKGROUND

As opposed to an endoluminal stent graft, a multilayer stent (MS) consists of a porous mesh, which allows for the possibility of treating pararenal aortic aneurysms (PRAAs) that involve a significant branch vessel. However, the choice of the density of the MS plays a vital role in isolating the aneurysm and allowing unobstructed blood flow in the branch vessel.

METHOD

In the present study, we examined 3 cases (without a stent and with single-layer and double-layer stents) via numerical simulations to explore the feasibility of the MSs used in the treatment of such aneurysms and estimate whether there is a more appropriate or optimal stent density.

RESULTS

With stent intervention, the velocity of blood flow in the sac decreased, but the pressure on the surface of the aneurysm did not decrease although it became more uniform. In addition, the "region of double low" (with low wall shear stress and a low Sherwood number) enlarged after stent implantation. Even with the double-layer stent, however, the flux of the branch vessel was still above normal, and we could predict that the optimal stent porosity was approximately 49.9%.

CONCLUSIONS

Unlike in previous studies, an MS could not be feasibly applied to high-risk PRAAs. However, an MS can induce sac thrombosis in the later stages while maintaining visceral vessel patency, and our results suggest that the optimal stent may be a 4-layer stent.

Study of helical flow inducers with different thread pitches and diameters in vena cava

Pulmonary embolism is a severe, potentially life-threatening condition. Inferior vena cava filters have been used to prevent recurrent pulmonary embolisms. However, the build-up of thrombosis in vena cava filters after deployment presents a severe problem to patients. Previous studies proposed that filters with helical flow are beneficial and capable of alleviating this problem. In this study, the hemodynamic performances of four typical helical flow inducers in the vena cava are determined using computational fluid dynamics simulations (steady-state and pulsatile flow) and compared. Pilot in vitro experiments were also conducted. The simulation results demonstrate that large-diameter inducers produce helical flow. Among inducers with identical diameter, those with a smaller thread pitch are more likely to induce increased helical flow. We also observed that the small thread pitch inducers can yield higher shear rates. Furthermore, a large diameter, small thread pitch helical flow inducer increases the time-averaged wall shear stress and reduces the oscillating shear index and relative residence time on the vessel wall in the vicinity of the helical flow inducer. In vitro experiments also verify that large diameter inducers generate a helical flow. A notable observation of this study is that the diameter is the key parameter that affects the induction of a helical flow. This study will likely provide important guidance for the design of interventional treatments and the deployment of filters associated with helical flow in the vena cava.

Numerical identification of the rupture locations in patient-specific abdominal aortic aneurysmsusing hemodynamic parameters

The rupture of an abdominal aortic aneurysm (AAA) is generally an unexpected event. Up to now, there is no agreement on an accurate criteria to predict the rupture risk of AAAs. This paper aims to numerically investigate the hemodynamics of three ruptured and one non-ruptured patient-specific AAA models to correlate local hemodynamic parameters with the rupture sites, and for the first time, this study introduced helicity as a potential index for the rupture potential of AAAs.3D reconstructions from CT scans were done. The simulation revealed that all the rupture sites were in regions of stagnation with near zero wall shear stress (WSS) but large WSS gradient (WSSG), which may explain the observation by the former researchers that the rupture site in the ruptured AAA has the lowest recorded wall thickness compared to other non-ruptured regions. Moreover, all the ruptures occurred at regions of zero helicity which represents a purely axial or circumferential flow. In addition, this study revealed that the double low region for the non-ruptured AAA was present with a thick layer of plaques, it suggests that the AAA rupture and the formation of atherosclerotic plaques may share a lot common physiological features. However, the fact that there are no plaques present in the walls of three RAAAs also indicates that AAA is not always a result of atherosclerosis. The current computational study may complement the maximum diameter, peak wall stress and other clinically relevant factors in AAA ruptures to identify the rupture sites of AAAs.

Analysis of Systolic Backflow and Secondary Helical Blood Flow in the Ascending Aorta Using Vector Flow Imaging

Secondary rotational flow and systolic backflow are seen in the ascending aorta and, in this study, were analyzed with the vector velocity method transverse oscillation. Twenty-five patients were scanned intra-operatively, and the vector velocities were related to estimates of transesophageal echocardiography and pulmonary artery catheter thermodilution, and associated with gender, age, aortic diameter, atherosclerotic plaques, left ventricular ejection fraction and previous myocardial infarctions. Secondary flow was present for all patients. The duration and rotational frequency (p < 0.001) and the duration and flow direction of the secondary flow (p < 0.002) were associated. Systolic backflow was present in 40% of the patients and associated with systolic velocities (p < 0.002) and the presence of atherosclerotic plaques (p < 0.001). No other significant associations were observed. The study indicates that backflow is injurious and that secondary flow is a normal flow phenomenon. The study also shows that transverse oscillation can provide new information on blood flow in the ascending aorta.

Finite element analysis of helical flows in human aortic arch: a novel index

This study investigates the helical secondary flows in the aortic arch using finite element analysis. The relationship between helical flow and the configuration of the aorta in patients of whose three-dimensional images constructed from computed tomography scans was examined. A finite element model of the pressurized root, arch, and supra-aortic vessels was developed to simulate the pattern of helical secondary flows. Calculations indicate that most of the helical secondary flow was formed in the ascending aorta. Angle α between the zero reference point and the aortic ostium (correlation coefficient (r) = -0.851, P = 0.001), the dispersion index of the cross section of the ascending (r = 0.683, P = 0.021) and descending aorta (r = 0.732, P = 0.010), all correlated closely with the presence of helical flow (P < 0.05). Stepwise multiple linear regression analysis confirmed angel α to be independently associated with the helical flow pattern in therein (standardized coefficients = -0.721, P = 0.023). The presence of helical fluid motion based on the atherosclerotic risks of patients, including those associated with diabetes, hypertension, hyperlipidemia, or renal insufficiency, was also evaluated. Numerical simulation of the flow patterns in aortas incorporating the atherosclerotic risks may better explain the mechanism of formation of helical flows and provide insight into causative factors that underlie them.

Computational fluid dynamic simulations of image-based stented coronary bifurcation models

One of the relevant phenomenon associated with in-stent restenosis in coronary arteries is an altered haemodynamics in the stented region. Computational fluid dynamics (CFD) offers the possibility to investigate the haemodynamics at a level of detail not always accessible within experimental techniques. CFD can quantify and correlate the local haemodynamics structures which might lead to in-stent restenosis. The aim of this work is to study the fluid dynamics of realistic stented coronary artery models which replicate the complete clinical procedure of stent implantation. Two cases of pathologic left anterior descending coronary arteries with their bifurcations are reconstructed from computed tomography angiography and conventional coronary angiography images. Results of wall shear stress and relative residence time show that the wall regions more prone to the risk of restenosis are located next to stent struts, to the bifurcations and to the stent overlapping zone for both investigated cases. Considering a bulk flow analysis, helical flow structures are generated by the curvature of the zone upstream from the stent and by the bifurcation regions. Helical recirculating microstructures are also visible downstream from the stent struts. This study demonstrates the feasibility to virtually investigate the haemodynamics of patient-specific coronary bifurcation geometries.