Assessing the impact of turbulent kinetic energy boundary conditions on turbulent flow simulations using computational fluid dynamics

Computational fluid dynamics has been widely used to study hemodynamics, but accurately determining boundary conditions for turbulent blood flow remains challenging. This study aims to investigate the effect of patient-specific turbulence boundary conditions on the accuracy of turbulent flow simulation. Using a stenosis model with 50% severity in diameter, the post-stenosis turbulence flow region was simulated with different planes to obtain inlet boundary conditions and simulate downstream flows. The errors of simulated flow fields obtained with turbulence kinetic energy (TKE) boundary data and arbitrary turbulence intensity were compared. Additionally, the study tested various TKE data resolutions and noise levels to simulate experimental environments. The mean absolute error of velocity and TKE was investigated with various turbulence intensities and TKE mapping. While voxel size and signal-to-noise ratio of the TKE data affected the results, simulation with SNR > 5 and voxel size < 10% resulted in better accuracy than simulations with turbulence intensities. The simulation with appropriate TKE boundary data resulted in a more accurate velocity and turbulence field than those with arbitrary turbulence intensity boundary conditions. The study demonstrated the potential improvement of turbulent blood flow simulation with patient-specific turbulence boundary conditions, which can be obtained from recent measurement techniques.

4D Flow cardiovascular magnetic resonance consensus statement: 2023 update

Hemodynamic assessment is an integral part of the diagnosis and management of cardiovascular disease. Four-dimensional cardiovascular magnetic resonance flow imaging (4D Flow CMR) allows comprehensive and accurate assessment of flow in a single acquisition. This consensus paper is an update from the 2015 '4D Flow CMR Consensus Statement'. We elaborate on 4D Flow CMR sequence options and imaging considerations. The document aims to assist centers starting out with 4D Flow CMR of the heart and great vessels with advice on acquisition parameters, post-processing workflows and integration into clinical practice. Furthermore, we define minimum quality assurance and validation standards for clinical centers. We also address the challenges faced in quality assurance and validation in the research setting. We also include a checklist for recommended publication standards, specifically for 4D Flow CMR. Finally, we discuss the current limitations and the future of 4D Flow CMR. This updated consensus paper will further facilitate widespread adoption of 4D Flow CMR in the clinical workflow across the globe and aid consistently high-quality publication standards.

Characterization of baseline hemodynamics after the Fontan procedure: a retrospective cohort study on the comparison of 4D Flow MRI and computational fluid dynamics

Introduction: The aim of this study was to characterize the hemodynamics of Fontan patients using both four-dimensional flow magnetic resonance imaging (4D Flow MRI) and computational fluid dynamics (CFD). Methods: Twenty-nine patients (3.5 ± 0.5 years) who had undergone the Fontan procedure were enrolled, and the superior vena cava (SVC), left pulmonary artery (LPA), right pulmonary artery (RPA), and conduit were segmented based on 4D Flow MRI images. Velocity fields from 4D Flow MRI were used as boundary conditions for CFD simulations. Hemodynamic parameters such as peak velocity (Vmax), pulmonary flow distribution (PFD), kinetic energy (KE), and viscous dissipation (VD) were estimated and compared between the two modalities. Results and discussion: The Vmax, KE, VD, PFD_{Total to LPA}, and PFD_{Total to RPA} of the Fontan circulation were 0.61 ± 0.18 m/s, 0.15 ± 0.04 mJ, 0.14 ± 0.04 mW, 41.3 ± 15.7%, and 58.7 ± 15.7% from 4D Flow MRI; and 0.42 ± 0.20 m/s, 0.12 ± 0.05 mJ, 0.59 ± 0.30 mW, 40.2 ± 16.4%, and 59.8 ± 16.4% from CFD, respectively. The overall velocity field, KE, and PFD from the SVC were in agreement between modalities. However, PFD from the conduit and VD showed a large discrepancy between 4D Flow MRI and CFD, most likely due to spatial resolution and data noise. This study highlights the necessity for careful consideration when analyzing hemodynamic data from different modalities in Fontan patients.

Non-invasive cardiovascular magnetic resonance assessment of pressure recovery distance after aortic valve stenosis

BACKGROUND

Decisions in the management of aortic stenosis are based on the peak pressure drop, captured by Doppler echocardiography, whereas gold standard catheterization measurements assess the net pressure drop but are limited by associated risks. The relationship between these two measurements, peak and net pressure drop, is dictated by the pressure recovery along the ascending aorta which is mainly caused by turbulence energy dissipation. Currently, pressure recovery is considered to occur within the first 40-50 mm distally from the aortic valve, albeit there is inconsistency across interventionist centers on where/how to position the catheter to capture the net pressure drop.

METHODS

We developed a non-invasive method to assess the pressure recovery distance based on blood flow momentum via 4D Flow cardiovascular magnetic resonance (CMR). Multi-center acquisitions included physical flow phantoms with different stenotic valve configurations to validate this method, first against reference measurements and then against turbulent energy dissipation (respectively n = 8 and n = 28 acquisitions) and to investigate the relationship between peak and net pressure drops. Finally, we explored the potential errors of cardiac catheterisation pressure recordings as a result of neglecting the pressure recovery distance in a clinical bicuspid aortic valve (BAV) cohort of n = 32 patients.

RESULTS

In-vitro assessment of pressure recovery distance based on flow momentum achieved an average error of 1.8 ± 8.4 mm when compared to reference pressure sensors in the first phantom workbench. The momentum pressure recovery distance and the turbulent energy dissipation distance showed no statistical difference (mean difference of 2.8 ± 5.4 mm, R2 = 0.93) in the second phantom workbench. A linear correlation was observed between peak and net pressure drops, however, with strong dependences on the valvular morphology. Finally, in the BAV cohort the pressure recovery distance was 78.8 ± 34.3 mm from vena contracta, which is significantly longer than currently accepted in clinical practise (40-50 mm), and 37.5% of patients displayed a pressure recovery distance beyond the end of the ascending aorta.

CONCLUSION

The non-invasive assessment of the distance to pressure recovery is possible by tracking momentum via 4D Flow CMR. Recovery is not always complete at the ascending aorta, and catheterised recordings will overestimate the net pressure drop in those situations. There is a need to re-evaluate the methods that characterise the haemodynamic burden caused by aortic stenosis as currently clinically accepted pressure recovery distance is an underestimation.

Fluid-structure interaction simulation of visceral perfusion and impact of different cannulation methods on aortic dissection

Hemodynamics in aortic dissection (AD) is closely associated with the risk of aortic aneurysm, rupture, and malperfusion. Altered blood flow in patients with AD can lead to severe complications such as visceral malperfusion. In this study, we aimed to investigate the effect of cannulation flow on hemodynamics in AD using a fluid-structure interaction simulation. We developed a specific-idealized AD model that included an intimal tear in the descending thoracic aorta, a re-entry tear in the left iliac artery, and nine branches. Two different cannulation methods were tested: (1) axillary cannulation (AC) only through the brachiocephalic trunk and (2) combined axillary and femoral cannulation (AFC) through the brachiocephalic trunk and the right common iliac artery. AC was found to result in the development of a pressure difference between the true lumen and false lumen, owing to the difference in the flow rate through each lumen. This pressure difference collapsed the true lumen, disturbing blood flow to the celiac and superior mesenteric arteries. However, in AFC, the pressure levels between the two lumens were similar, and no collapse occurred. Moreover, the visceral flow was higher than that in AC. Lastly, the stiffness of the intimal flap affected the true lumen's collapse.

Modelling and simulation of fluid flow through stenosis and aneurysm blood vessel: a computational hemodynamic analysis

In this article, the hemodynamics of nanofluid flow through the modelled stenosis-aneurysm models in the presence of the catheter has been studied. The eight stenosis-aneurysm models are developed to mimic biological observations and thus make the model more realistic. The mathematical understanding helps in treating the stenosis in the blood vessel by targeting the unhealthy region to the drug, which is coated on nanoparticles. The catheter achieves the active drug release to the aimed organs by coating on the catheter surface, which adds additional benefits. In the present hemodynamic study, the blood is modeled as a couple stress fluid; as a result, the highly non-linear momentum, temperature, and concentration equations were obtained. The fluid flow equations' complexity is further increased by incorporating the variable viscosity effects that arose due to the suspension of nanoparticles. The resultant mathematical model is solved by using the homotopy perturbation method. The convergence of the perturbed solutions is studied and depicted the degree of deformation in the case of temperature and concentration. The effects of the porous nature of the stenosis, no-slip at the catheter surface, and the free slip at the blood vessel boundary in the non-stenotic region are also considered in the model. The essential physiological property like surface shear stress is computed, and various parameters' influence on shear stress is analyzed. The present analysis can be helpful in understanding the enhancement in mass dispersion and heat transfer in unhealthy blood vessels, which could be used for drug delivery in the treatment of stenotic conditions.

In-vitro and In-Vivo Assessment of 4D Flow MRI Reynolds Stress Mapping for Pulsatile Blood Flow

Imaging hemodynamics play an important role in the diagnosis of abnormal blood flow due to vascular and valvular diseases as well as in monitoring the recovery of normal blood flow after surgical or interventional treatment. Recently, characterization of turbulent blood flow using 4D flow magnetic resonance imaging (MRI) has been demonstrated by utilizing the changes in signal magnitude depending on intravoxel spin distribution. The imaging sequence was extended with a six-directional icosahedral (ICOSA6) flow-encoding to characterize all elements of the Reynolds stress tensor (RST) in turbulent blood flow. In the present study, we aimed to demonstrate the feasibility of full RST analysis using ICOSA6 4D flow MRI under physiological conditions. First, the turbulence analysis was performed through in vitro experiments with a physiological pulsatile flow condition. Second, a total of 12 normal subjects and one patient with severe aortic stenosis were analyzed using the same sequence. The in-vitro study showed that total turbulent kinetic energy (TKE) was less affected by the signal-to-noise ratio (SNR), however, maximum principal turbulence shear stress (MPTSS) and total turbulence production (TP) had a noise-induced bias. Smaller degree of the bias was observed for TP compared to MPTSS. In-vivo study showed that the subject-variability on turbulence quantification was relatively low for the consistent scan protocol. The in vivo demonstration of the stenosis patient showed that the turbulence analysis could clearly distinguish the difference in all turbulence parameters as they were at least an order of magnitude larger than those from the normal subjects.

Quantification of visceral perfusion and impact of femoral cannulation: in vitro model of aortic dissection

OBJECTIVES

We aimed to simulate blood flow at an aortic dissection in an in vitro vascular model and assess the impact of the cannulation method on visceral perfusion.

METHODS

An aortic-dissection model with an acrylic aortic wall and silicone intimal flap was developed to study visceral perfusion under various cannulation conditions. The primary tear was placed in the proximal descending aorta and the re-entry site in the left common iliac artery. A cardiovascular pump was used to reproduce a normal pulsatile aortic flow and a steady cannulation flow. Axillary and axillary plus femoral cannulation were compared at flow rates of 3-7 l/min. Haemodynamics were analysed by using four-dimensional flow magnetic resonance imaging.

RESULTS

Axillary cannulation (AC) was found to collapse the true lumen at the coeliac and superior mesentery arteries, while combined axillary and femoral cannulation did not change the size of the true lumen. Combined axillary and femoral cannulation resulted in a larger visceral flow than did AC alone. When axillary plus femoral cannulation was used, the visceral flow increased by 125% at 3 l/min, by 89% at 4 l/min, by 67% at 5 L/min, by 98% at 6 l/min and by 101% at 7 l/min, respectively, compared to those with the AC only.

CONCLUSIONS

Our model was useful to understanding the haemodynamics in aortic dissection. In this specific condition, we confirmed that the intimal flap motion can partially block blood flow to the coeliac and superior mesenteric arteries and that additional femoral cannulation can increase visceral perfusion.

Improving Blood Flow Visualization of Recirculation Regions at Carotid Bulb in 4D Flow MRI Using Semi-Automatic Segmentation with ITK-SNAP

Assessment of carotid bulb hemodynamics using four-dimensional (4D) flow magnetic resonance imaging (MRI) requires accurate segmentation of recirculation regions that is frequently hampered by limited resolution. This study aims to improve the accuracy of 4D flow MRI carotid bulb segmentation and subsequent recirculation regions analysis. Time-of-flight (TOF) MRI and 4D flow MRI were performed on bilateral carotid artery bifurcations in seven healthy volunteers. TOF-MRI data was segmented into 3D geometry for computational fluid dynamics (CFD) simulations. ITK-SNAP segmentation software was included in the workflow for the semi-automatic generation of 4D flow MRI angiographic data. This study compared the velocities calculated at the carotid bifurcations and the 3D blood flow visualization at the carotid bulbs obtained by 4D flow MRI and CFD. By applying ITK-SNAP segmentation software, an obvious improvement in the 4D flow MRI visualization of the recirculation regions was observed. The 4D flow MRI images of the recirculation flow characteristics of the carotid artery bulbs coincided with the CFD. A reasonable agreement was found in terms of velocity calculated at the carotid bifurcation between CFD and 4D flow MRI. However, the dispersion of velocity data points relative to the local errors of measurement in 4D flow MRI remains. Our proposed strategy showed the feasibility of improving recirculation regions segmentation and the potential for reliable blood flow visualization in 4D flow MRI. However, quantitative analysis of recirculation regions in 4D flow MRI with ITK-SNAP should be enhanced for use in clinical situations.

Comparison of Hemodynamic Visualization in Cerebral Arteries: Can Magnetic Resonance Imaging Replace Computational Fluid Dynamics?

A multimodality approach was applied using four-dimensional flow magnetic resonance imaging (4D flow MRI), time-of-flight magnetic resonance angiography (TOF-MRA) signal intensity gradient (SIG), and computational fluid dynamics (CFD) to investigate the 3D blood flow characteristics and wall shear stress (WSS) of the cerebral arteries. TOF-MRA and 4D flow MRI were performed on the major cerebral arteries in 16 healthy volunteers (mean age 34.7 ± 7.6 years). The flow rate measured with 4D flow MRI in the internal carotid artery, middle cerebral artery, and anterior cerebral artery were 3.8, 2.5, and 1.2 mL/s, respectively. The 3D blood flow pattern obtained through CFD and 4D flow MRI on the cerebral arteries showed reasonable consensus. CFD delivered much greater resolution than 4D flow MRI. TOF-MRA SIG and CFD WSS of the major cerebral arteries showed reasonable consensus with the locations where the WSS was relatively high. However, the visualizations were very different between TOF-MRA SIG and CFD WSS at the internal carotid artery bifurcations, the anterior cerebral arteries, and the anterior communicating arteries. 4D flow MRI, TOF-MRA SIG, and CFD are complementary methods that can provide additional insight into the hemodynamics of the human cerebral artery.

In vitro experiments on ICOSA6 4D flow MRI measurement for the quantification of velocity and turbulence parameters

PURPOSE

To perform comprehensive in vitro experiments using six-directional icosahedral flow encoding (ICOSA6) 4D flow magnetic resonance imaging (MRI) under various scan conditions to analyze the robustness of velocity and turbulence quantification.

MATERIALS AND METHODS

In vitro flow phantoms with steady flow rates of 10 and 20 L/min were scanned using both conventional 4D flow MRI and ICOSA6. Experiments focused on comparisons between ICOSA6 and conventional four point (4P) methods, and the effects of contrast agents, velocity encoding range (Venc), and scan direction on velocity and turbulence quantification.

RESULTS

The results demonstrated that 1) ICOSA6 improves the velocity-to-noise ratio (VNR) of velocity estimation by 33% (on average) and results in similar turbulent kinetic energy (TKE) estimation as the 4P method. 2) Measurements with a contrast agent resulted in more than a 2.5 fold increase in average VNR. However, the improvement of total TKE quantification was not obvious. 3) TKE estimation was less affected by Venc and the scan direction, whereas turbulence production (TP) estimation was largely affected by these measurement conditions. The effects of Venc and scan direction accounted for less than 11.63% of TKE estimation, but up to 33.89% of TP estimation.

CONCLUSION

The ICOSA6 scheme is compatible with conventional 4D flow MRI for velocity and TKE measurement. Contrast agents are effective at increasing VNR, but not signal-to-noise ratio for TKE quantification. The effects of Venc and scan direction influence total TP more than total TKE.

First report of the use of mucosal swabs of the palatine tonsillar crypt for detection of Mycoplasma bovis in naturally infected calves

Aims: To compare detection by real-time PCR of DNA from Mycoplasma bovis on mucosal swabs taken from the palatine tonsillar crypt and the mainstem bronchi of clinically asymptomatic calves after slaughter. Methods: We compared the sensitivity of mucosal swabs taken from two sites: the palatine tonsillar crypt and the mainstem bronchi. Paired samples were taken post-mortem at slaughter from 55 clinically well calves from an infected herd and were tested by real-time PCR for the presence of M. bovis-specific DNA. Results: Mycoplasma bovis DNA was detected in 51 palatine tonsillar crypt swabs (92.7 (95% CI = 82.4-98.0)%) and seven mainstem bronchial swabs (12.7 (95% CI = 5.3-24.5)%). All seven calves with positive mainstem bronchial swabs also had positive palatine tonsillar crypt swabs. Conclusions: When compared to mucosal swabs of the mainstem bronchi, mucosal swabs of the palatine tonsillar crypt were seven times more sensitive for the post-mortem detection of M. bovis DNA. The viability of detected M. bovis was not assessed, because any cattle carrying viable or non-viable M. bovis DNA were determined to be a potential risk to eradication. Palatine tonsillar crypt mucosa may be a useful anatomical site for real-time PCR detection of M. bovis DNA in naturally infected calves. More work is needed to define the persistence and viability of M. bovis at this anatomical site. Clinical relevance: The results of this study helped form the basis of surveillance tools used in M. bovis control and eradication efforts. Familiarity with these results may help veterinarians better communicate with their clients about the science behind the eradication efforts.

Non-invasive estimation of relative pressure in turbulent flow using virtual work-energy

Vascular pressure differences are established risk markers for a number of cardiovascular diseases. Relative pressures are, however, often driven by turbulence-induced flow fluctuations, where conventional non-invasive methods may yield inaccurate results. Recently, we proposed a novel method for non-turbulent flows, νWERP, utilizing the concept of virtual work-energy to accurately probe relative pressure through complex branching vasculature. Here, we present an extension of this approach for turbulent flows: νWERP-t. We present a theoretical method derivation based on flow covariance, quantifying the impact of flow fluctuations on relative pressure. νWERP-t is tested on a set of in-vitro stenotic flow phantoms with data acquired by 4D flow MRI with six-directional flow encoding, as well as on a patient-specific in-silico model of an acute aortic dissection. Over all tests νWERP-t shows improved accuracy over alternative energy-based approaches, with excellent recovery of estimated relative pressures. In particular, the use of a guaranteed divergence-free virtual field improves accuracy in cases where turbulent flows skew the apparent divergence of the acquired field. With the original νWERP allowing for assessment of relative pressure into previously inaccessible vasculatures, the extended νWERP-t further enlarges the method's clinical scope, underlining its potential as a novel tool for assessing relative pressure in-vivo.

Extent of Subprosthetic Pannus after Aortic Valve Replacement: Changes Over Time and Relationship with Echocardiographic Findings

Purpose

This study aimed to evaluate changes of subprosthetic pannus on cardiac CT and determine its relationship to echocardiographic findings in patients with mechanical aortic valve replacement (AVR).

Materials and Methods

Between April 2011 and November 2017, 17 AVR patients (56.8 ± 8.9 years, 12% male) who showed pannus formation on CT and had undergone both follow-up CT and echocardiography were included. The mean interval from AVR to the date of pannus detection was 10.5 ± 7.1 years. In the initial and follow-up CT and echocardiography, the pannus extent and echocardiographic parameters were compared using paired t-tests. The relationship between the opening angle of the prosthetic valve and the pannus extent was evaluated using Pearson correlation analysis.

Results

The pannus extent was significantly increased on CT (p < 0.05). The peak velocity (3.9 ± 0.8 m/s vs. 4.2 ± 0.8 m/s, p = 0.03) and mean pressure gradient (36.4 ± 15.5 mm Hg vs. 42.1 ± 15.8 mm Hg, p = 0.03) were significantly increased. The mean opening angles of the mechanical aortic leaflets were slightly decreased, but there was no statistical significance (73.1 ± 8.3° vs. 69.4 ± 12.1°, p = 0.12). The opening angle of the prosthetic leaflets was inversely correlated with the pannus extent (r = −0.57, p < 0.001).

Conclusion

The pannus extent increases over time, increasing transvalvular peak velocity and the pressure gradient. CT can be used to evaluate the pannus extent associated with hemodynamic changes that need to be managed by surgical intervention.

Impact of coronary lumen reconstruction on the estimation of endothelial shear stress: in vivo comparison of three-dimensional quantitative coronary angiography and three-dimensional fusion combining optical coherent tomography

Aims

It is not clearly elucidated how the fusion technique improves the accuracy of endothelial shear stress (ESS) prediction, in comparison with that of three-dimensional (3D) quantitative coronary angiography (QCA) alone. We aimed to evaluate the difference in geometric measurements and haemodynamic estimation between 3D QCA and a 3D fusion model combining 3D QCA and optical coherence tomography (OCT).

Methods and results

Computational fluid dynamics was assessed in the coronary models of 20 patients. In the plane-per-plane comparison, the difference and agreement were assessed using a generalized linear mixed model and concordance correlation coefficient (CCC), respectively. The haemodynamic feature around minimum-lumen-diameter (MLD) was characterized using CCC values calculated for 1-mm segments. In comparison with the 3D fusion model, 3D QCA showed a shorter maximum lumen diameter (2.54 ± 0.67 mm vs. 2.78 ± 0.73 mm, P < 0.001) and smaller lumen area (4.81 ± 2.56 mm2 vs. 5.66 ± 2.97 mm2, P < 0.001), resulting in a significantly higher ESS (4.64 Pa vs. 3.78 Pa, p = 0.029). A more asymmetric lumen shape of the 3D fusion model was more likely associated with under- and over-estimation of the maximum and minimum lumen diameters in the 3D QCA model, respectively. The circumferential ESS variations, which were blunted by 3D QCA, showed the worst concordance near the MLD site (CCC = 0.370) on segment-based comparison.

Conclusion

The 3D fusion technique may be a more relevant tool for the haemodynamic simulation of coronary arteries through providing more accurate lumen characterization than 3D QCA.

Validation of pressure drop assessment using 4D flow MRI-based turbulence production in various shapes of aortic stenoses

PURPOSE

To validate pressure drop measurements using 4D flow MRI-based turbulence production in various shapes of stenotic stenoses.

METHODS

In vitro flow phantoms with seven different 3D-printed aortic valve geometries were constructed and scanned with 4D flow MRI with six-directional flow encoding (ICOSA6). The pressure drop through the valve was non-invasively predicted based on the simplified Bernoulli, the extended Bernoulli, the turbulence production, and the shear-scaling methods. Linear regression and agreement of the predictions with invasively measured pressure drop were analyzed.

RESULTS

All pressure drop predictions using 4D Flow MRI were linearly correlated to the true pressure drop but resulted in different regression slopes. The regression slope and 95% limits of agreement for the simplified Bernoulli method were 1.35 and 11.99 ± 21.72 mm Hg. The regression slope and 95% limits of agreement for the extended Bernoulli method were 1.02 and 0.74 ± 8.48 mm Hg. The regression slope and 95% limits of agreement for the turbulence production method were 0.89 and 0.96 ± 8.01 mm Hg. The shear-scaling method presented good correlation with an invasively measured pressure drop, but the regression slope varied between 0.36 and 1.00 depending on the shear-scaling coefficient.

CONCLUSION

The pressure drop assessment based on the turbulence production method agrees well with the extended Bernoulli method and invasively measured pressure drop in various shapes of the aortic valve. Turbulence-based pressure drop estimation can, as a complement to the conventional Bernoulli method, play a role in the assessment of valve diseases.

Effect of pannus formation on the prosthetic heart valve: In vitro demonstration using particle image velocimetry

Although hemodynamic influence of the subprosthetic tissue, termed as pannus, may contribute to prosthetic aortic valve dysfunction, the relationship between pannus extent and hemodynamics in the prosthetic valve has rarely been reported. We investigated the fluid dynamics of pannus formation using in vitro experiments with particle image velocimetry. Subvalvular pannus formation caused substantial changes in prosthetic valve transvalvular peak velocity, transvalvular pressure gradient (TPG) and opening angle. Maximum flow velocity and corresponding TPG were mostly affected by pannus width. When the pannus width was 25% of the valve diameter, pannus formation elevated TPG to >2.5 times higher than that without pannus formation. Opening dysfunction was observed only for a pannus involvement angle of 360°. Although circumferential pannus with an involvement angle of 360° decreased the opening angle of the valve from approximately 82° to 58°, eccentric pannus with an involvement angle of 180° did not induce valve opening dysfunction. The pannus involvement angle largely influenced the velocity flow field at the aortic sinus and corresponding hemodynamic indices, including wall shear stress, principal shear stress and viscous energy loss distributions. Substantial discrepancy between the velocity-based TPG estimation and direct pressure measurements was observed for prosthetic valve flow with pannus formation.

Age-Related Vascular Changes Affect Turbulence in Aortic Blood Flow

Turbulent blood flow is implicated in the pathogenesis of several aortic diseases but the extent and degree of turbulent blood flow in the normal aorta is unknown. We aimed to quantify the extent and degree of turbulece in the normal aorta and to assess whether age impacts the degree of turbulence. 22 young normal males (23.7 ± 3.0 y.o.) and 20 old normal males (70.9 ± 3.5 y.o.) were examined using four dimensional flow magnetic resonance imaging (4D Flow MRI) to quantify the turbulent kinetic energy (TKE), a measure of the intensity of turbulence, in the aorta. All healthy subjects developed turbulent flow in the aorta, with total TKE of 3–19 mJ. The overall degree of turbulence in the entire aorta was similar between the groups, although the old subjects had about 73% more total TKE in the ascending aorta compared to the young subjects (young = 3.7 ± 1.8 mJ, old = 6.4 ± 2.4 mJ, p < 0.001). This increase in ascending aorta TKE in old subjects was associated with age-related dilation of the ascending aorta which increases the volume available for turbulence development. Conversely, age-related dilation of the descending and abdominal aorta decreased the average flow velocity and suppressed the development of turbulence. In conclusion, turbulent blood flow develops in the aorta of normal subjects and is impacted by age-related geometric changes. Non-invasive assessment enables the determination of normal levels of turbulent flow in the aorta which is a prerequisite for understanding the role of turbulence in the pathophysiology of cardiovascular disease.

Estimating the irreversible pressure drop across a stenosis by quantifying turbulence production using 4D Flow MRI

The pressure drop across a stenotic vessel is an important parameter in medicine, providing a commonly used and intuitive metric for evaluating the severity of the stenosis. However, non-invasive estimation of the pressure drop under pathological conditions has remained difficult. This study demonstrates a novel method to quantify the irreversible pressure drop across a stenosis using 4D Flow MRI by calculating the total turbulence production of the flow. Simulation MRI acquisitions showed that the energy lost to turbulence production can be accurately quantified with 4D Flow MRI within a range of practical spatial resolutions (1-3 mm; regression slope = 0.91, R² = 0.96). The quantification of the turbulence production was not substantially influenced by the signal-to-noise ratio (SNR), resulting in less than 2% mean bias at SNR > 10. Pressure drop estimation based on turbulence production robustly predicted the irreversible pressure drop, regardless of the stenosis severity and post-stenosis dilatation (regression slope = 0.956, R² = 0.96). In vitro validation of the technique in a 75% stenosis channel confirmed that pressure drop prediction based on the turbulence production agreed with the measured pressure drop (regression slope = 1.15, R² = 0.999, Bland-Altman agreement = 0.75 ± 3.93 mmHg).

Flow characteristics around proximal and distal stenoses in a series of tandem stenosed vessels

The flow characteristics around the proximal and distal stenoses in tandem vessel models are experimentally investigated with varying flow rates (Q=0.25, 0.5, 1.0L/min), interspacing distances (L=3, 6, 10 of diameter D) and severities (S=50%, 75% reduction in diameter). When the interspacing L is larger than 10 D, no fluid-dynamic interaction is observed. The flow between the proximal and distal stenoses becomes stabilized (turbulence intensity of <3%) as the interspacing distance decreases. When the severity S is 75%, the transition from laminar to turbulent flow occurs at a flow rate higher than 0.5L/min, although the interspacing distance L is 3 D. Formation of recirculation flow is restricted by the presence of distal stenosis as the interspacing distance decreases. In this case, the flow between the stenoses is focused on the central region. The center-line velocity at the neck of the distal stenosis is approximately 10-15% higher than that of the proximal stenosis with equal severity of S=50%. When the inlet flow is center-focused, the lengths of the recirculation and the jet core behind the distal stenosis increase with decrease in interspacing distance L. When the inlet flow is turbulent, the transition from laminar to turbulent flow occurs early as the interspacing distance L is reduced. When the upstream proximal stenosis exhibits increased severity, the pressure drop is measured to be 20% compared with that when the severity of the downstream distal stenosis is increased at the flow rate of Q=1.0L/min.

The influence of the aortic valve angle on the hemodynamic features of the thoracic aorta

Since the first observation of a helical flow pattern in aortic blood flow, the existence of helical blood flow has been found to be associated with various pathological conditions such as bicuspid aortic valve, aortic stenosis, and aortic dilatation. However, an understanding of the development of helical blood flow and its clinical implications are still lacking. In our present study, we hypothesized that the direction and angle of aortic inflow can influence helical flow patterns and related hemodynamic features in the thoracic aorta. Therefore, we investigated the hemodynamic features in the thoracic aorta and various aortic inflow angles using patient-specific vascular phantoms that were generated using a 3D printer and time-resolved, 3D, phase-contrast magnetic resonance imaging (PC-MRI). The results show that the rotational direction and strength of helical blood flow in the thoracic aorta largely vary according to the inflow direction of the aorta, and a higher helical velocity results in higher wall shear stress distributions. In addition, right-handed rotational flow conditions with higher rotational velocities imply a larger total kinetic energy than left-handed rotational flow conditions with lower rotational velocities.

Hemodynamic Measurement Using Four-Dimensional Phase-Contrast MRI: Quantification of Hemodynamic Parameters and Clinical Applications

Recent improvements have been made to the use of time-resolved, three-dimensional phase-contrast (PC) magnetic resonance imaging (MRI), which is also named four-dimensional (4D) PC-MRI or 4D flow MRI, in the investigation of spatial and temporal variations in hemodynamic features in cardiovascular blood flow. The present article reviews the principle and analytical procedures of 4D PC-MRI. Various fluid dynamic biomarkers for possible clinical usage are also described, including wall shear stress, turbulent kinetic energy, and relative pressure. Lastly, this article provides an overview of the clinical applications of 4D PC-MRI in various cardiovascular regions.

Three-Dimensional Printing: Basic Principles and Applications in Medicine and Radiology

The advent of three-dimensional printing (3DP) technology has enabled the creation of a tangible and complex 3D object that goes beyond a simple 3D-shaded visualization on a flat monitor. Since the early 2000s, 3DP machines have been used only in hard tissue applications. Recently developed multi-materials for 3DP have been used extensively for a variety of medical applications, such as personalized surgical planning and guidance, customized implants, biomedical research, and preclinical education. In this review article, we discuss the 3D reconstruction process, touching on medical imaging, and various 3DP systems applicable to medicine. In addition, the 3DP medical applications using multi-materials are introduced, as well as our recent results.

Post-stenotic Recirculating Flow May Cause Hemodynamic Perforator Infarction

BACKGROUND AND PURPOSE

The primary mechanism underlying paramedian pontine infarction (PPI) is atheroma obliterating the perforators. Here, we encountered a patient with PPI in the post-stenotic area of basilar artery (BA) without a plaque, shown by high-resolution magnetic resonance imaging (HR-MRI). We performed an experiment using a 3D-printed BA model and a particle image velocimetry (PIV) to explore the hemodynamic property of the post-stenotic area and the mechanism of PPI.

METHODS

3D-model of a BA stenosis was reconstructed with silicone compound using a 3D-printer based on the source image of HR-MRI. Working fluid seeded with fluorescence particles was used and the velocity of those particles was measured horizontally and vertically. Furthermore, microtubules were inserted into the posterior aspect of the model to measure the flow rates of perforators (pre-and post-stenotic areas). The flow rates were compared between the microtubules.

RESULTS

A recirculating flow was observed from the post-stenotic area in both directions forming a spiral shape. The velocity of the flow in these regions of recirculation was about one-tenth that of the flow in other regions. The location of recirculating flow well corresponded with the area with low-signal intensity at the time-of-flight magnetic resonance angiography and the location of PPI. Finally, the flow rate through the microtubule inserted into the post-stenotic area was significantly decreased comparing to others (P<0.001).

CONCLUSIONS

Perforator infarction may be caused by a hemodynamic mechanism altered by stenosis that induces a recirculation flow. 3D-printed modeling and PIV are helpful understanding the hemodynamics of intracranial stenosis.

Multi-VENC acquisition of four-dimensional phase-contrast MRI to improve precision of velocity field measurement

PURPOSE

The present study aims to improve precision of four-dimensional (4D) phase-contrast (PC) MRI technique by using multiple velocity encoding (VENC) parameters.

THEORY AND METHODS

The 3D flow fields in an in vitro stenosis phantom and an in vivo ascending aorta were determined using a 4D PC-MRI sequence with multiple VENC values. The velocity field obtained for large VENC was combined with that from small VENC, unless velocity data were lost by phase aliasing and phase dispersion. Noise levels of the combined velocity fields were compared with the increasing overlapping number of VENC parameters.

RESULTS

The phantom measurement showed that the multi-VENC acquisition reduced the noise levels in radial and axial velocities (> 24 cm/s at VENC = 300 cm/s) down to 0.80 ± 0.45 cm/s and 5.60 ± 2.63 cm/s, respectively. This increased the velocity-to-noise ratio (VNR) by approximately two-fold to six-fold depending on the locations. As a result, the multi-VENC measurement could visualize the low-velocity recirculating flows more clearly.

CONCLUSION

The multi-VENC measurement of 4D PC-MRI sequence increased the VNR distribution by reducing velocity noise. The improved VNR can be beneficial for investigating blood flow structures in a flow field with a high velocity dynamic range.

Effect of swirling blood flow on vortex formation at post-stenosis

Various clinical observations reported that swirling blood flow is a normal physiological flow pattern in various vasculatures. The swirling flow has beneficial effects on blood circulation through the blood vessels. It enhances oxygen transfer and reduces low-density lipoprotein concentration in the blood vessel by enhancing cross-plane mixing of the blood. However, the fluid-dynamic roles of the swirling flow are not yet fully understood. In this study, inhibition of material deposition at the post-stenosis region by the swirling flow was observed. To reveal the underlying fluid-dynamic characteristics, pathline flow visualization and time-resolved particle image velocimetry measurements were conducted. Results showed that the swirling inlet flow increased the development of vortices at near wall region of the post-stenosis, which can suppress further development of stenosis by enhancing transport and mixing of the blood flow. The fluid-dynamic characteristics obtained in this study would be useful for improving hemodynamic characteristics of vascular grafts and stents in which the stenosis frequently occurred. Moreover, the time-resolved particle image velocimetry measurement technique and vortex identification method employed in this study would be useful for investigating the fluid-dynamic effects of the swirling flow on various vascular environments.

Google unveils a glimpse of allergic rhinitis in the real world

Google Trends (GT) is a Web-based surveillance tool used to explore the searching trends of specific queries on Google. Recent studies have suggested the utility of GT in predicting outbreaks of influenza and other diseases. However, this utility has not been thoroughly evaluated for allergic diseases. Therefore, we investigated the utility of GT for predicting the epidemiology of allergic rhinitis. In the USA, GT for allergic rhinitis showed repetitive seasonality that peaked in late April and early May and then rapidly decreased, and a second small peak occurred in September. These trends are highly correlated with the searching trends for other queries such as 'pollen count', antihistamines such as loratadine and cetirizine (all r > 0.88 and all P < 0.001), and even the total pollen count collected from 21 pollen counters across the USA (r = 0.928, P < 0.001). Google Trends for allergic rhinitis was similar to the monthly changes in rhinitis symptoms according to the US National Health and Nutrition Examination Survey III, sales for Claritin(®) and all over-the-counter antihistamines, and the number of monthly page views of 'claritin.com'. In conclusion, GT closely reflects the real-world epidemiology of allergic rhinitis in the USA and could potentially be used as a monitoring tool for allergic rhinitis.

α-secretase cleaved amyloid precursor protein (APP) accumulates in cholinergic dystrophic neurites in normal, aged hippocampus

AIMS

Dystrophic neurites are associated with β-amyloid (Aβ) plaques in the brains of Alzheimer's disease (AD) patients and are also found in some specific areas of normal, aged brains. This study assessed the molecular characteristics of dystrophic neurites in normal ageing and its difference from AD.

METHODS

We compared the dystrophic neurites in normal aged human brains (age 20-70 years) and AD brains (Braak stage 4-6) by immunostaining against ChAT, synaptophysin, γ-tubulin, cathepsin-D, Aβ1-16, Aβ17-24, amyloid precursor protein (APP)-CT695 and APP-NT. We then tested the reproducibility in C57BL/6 mice neurone cultures.

RESULTS

In normal, aged mice and humans, we found an increase in clustered dystrophic neurites of cholinergic neurones in CA1 regions of the hippocampus and layer II and III regions of the entorhinal cortex, which are the major and earliest affected areas in AD. These dystrophic neurites showed accumulation of sAPPα peptides cleaved from the amyloid precursor protein by α-secretase rather than Aβ or C-terminal fragments. In contrast, Aβ and APP-CTFs accumulated in the dystrophic neurites in and around Aβ plaques of AD patients. Several experiments suggested that the accumulation of sAPPα resulted from ageing-related proteasomal dysfunction.

CONCLUSIONS

Ageing-associated impairment of the proteasomal system and accumulation of sAPPα at cholinergic neurites in specific areas of brain regions associated with memory could be associated with the normal decline of memory in aged individuals. In addition, these age-related changes might be the most vulnerable targets of pathological insults that result in pathological accumulation of Aβ and/or APP-CTFs and lead to neurodegenerative conditions such as AD.

Fluid-dynamic optimal design of helical vascular graft for stenotic disturbed flow

Although a helical configuration of a prosthetic vascular graft appears to be clinically beneficial in suppressing thrombosis and intimal hyperplasia, an optimization of a helical design has yet to be achieved because of the lack of a detailed understanding on hemodynamic features in helical grafts and their fluid dynamic influences. In the present study, the swirling flow in a helical graft was hypothesized to have beneficial influences on a disturbed flow structure such as stenotic flow. The characteristics of swirling flows generated by helical tubes with various helical pitches and curvatures were investigated to prove the hypothesis. The fluid dynamic influences of these helical tubes on stenotic flow were quantitatively analysed by using a particle image velocimetry technique. Results showed that the swirling intensity and helicity of the swirling flow have a linear relation with a modified Germano number (Gn*) of the helical pipe. In addition, the swirling flow generated a beneficial flow structure at the stenosis by reducing the size of the recirculation flow under steady and pulsatile flow conditions. Therefore, the beneficial effects of a helical graft on the flow field can be estimated by using the magnitude of Gn*. Finally, an optimized helical design with a maximum Gn* was suggested for the future design of a vascular graft.