Assessment of Reynolds stress components and turbulent pressure loss using 4D flow MRI with extended motion encoding.

PURPOSE: To measure the Reynolds stress tensor using 4D flow MRI, and to evaluate its contribution to computed pressure maps. METHODS: A method to assess both velocity and Reynolds stress using 4D flow MRI is presented and evaluated. The Reynolds stress is compared by cross-sectional integrals of the Reynolds stress invariants. Pressure maps are computed using the pressure Poisson equation-both including and neglecting the Reynolds stress. RESULT: Good agreement is seen for Reynolds stress between computational fluid dynamics, simulated MRI, and MRI experiment. The Reynolds stress can significantly influence the computed pressure loss for simulated (eg, -0.52% vs -15.34% error; P < 0.001) and experimental (eg, 306 ± 11 vs 203 ± 6 Pa; P < 0.001) data. A 54% greater pressure loss is seen at the highest experimental flow rate when accounting for Reynolds stress (P < 0.001). CONCLUSION: 4D flow MRI with extended motion-encoding enables quantification of both the velocity and the Reynolds stress tensor. The additional information provided by this method improves the assessment of pressure gradients across a stenosis in the presence of turbulence. Unlike conventional methods, which are only valid if the flow is laminar, the proposed method is valid for both laminar and disturbed flow, a common presentation in diseased vessels. Magn Reson Med 79:1962-1971, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Highly accelerated intracranial 4D flow MRI: evaluation of healthy volunteers and patients with intracranial aneurysms.

OBJECTIVES: To evaluate an accelerated 4D flow MRI method that provides high temporal resolution in a clinically feasible acquisition time for intracranial velocity imaging. MATERIALS AND METHODS: Accelerated 4D flow MRI was developed by using a pseudo-random variable-density Cartesian undersampling strategy (CIRCUS) with the combination of k-t, parallel imaging and compressed sensing image reconstruction techniques (k-t SPARSE-SENSE). Four-dimensional flow data were acquired on five healthy volunteers and eight patients with intracranial aneurysms using CIRCUS (acceleration factor of R = 4, termed CIRCUS4) and GRAPPA (R = 2, termed GRAPPA2) as the reference method. Images with three times higher temporal resolution (R = 12, CIRCUS12) were also reconstructed from the same acquisition as CIRCUS4. Qualitative and quantitative image assessment was performed on the images acquired with different methods, and complex flow patterns in the aneurysms were identified and compared. RESULTS: Four-dimensional flow MRI with CIRCUS was achieved in 5 min and allowed further improved temporal resolution of <30 ms. Volunteer studies showed similar qualitative and quantitative evaluation obtained with the proposed approach compared to the reference (overall image scores: GRAPPA2 3.2 ± 0.6; CIRCUS4 3.1 ± 0.7; CIRCUS12 3.3 ± 0.4; difference of the peak velocities: -3.83 ± 7.72 cm/s between CIRCUS4 and GRAPPA2, -1.72 ± 8.41 cm/s between CIRCUS12 and GRAPPA2). In patients with intracranial aneurysms, the higher temporal resolution improved capturing of the flow features in intracranial aneurysms (pathline visualization scores: GRAPPA2 2.2 ± 0.2; CIRCUS4 2.5 ± 0.5; CIRCUS12 2.7 ± 0.6). CONCLUSION: The proposed rapid 4D flow MRI with a high temporal resolution is a promising tool for evaluating intracranial aneurysms in a clinically feasible acquisition time.

Globally local: Hyper-local modeling for accurate forecast of COVID-19.

IMPORTANCE: Assumption of a well-mixed population during modeling is often erroneously made without due analysis of its validity. Ignoring the importance of the geo-spatial granularity at which the data is collected could have significant implications on the quality of forecasts and the actionable clinical recommendations that are based on it. OBJECTIVE: This paper's primary objective is to test the hypothesis that the characteristic dynamics defining the trajectory of the pandemic in a region is lost when the data is aggregated and modeled at higher geo-spatial levels. DESIGN: We use publicly available confirmed SARS-CoV-2 cases and deaths from January 1st, 2020 to August 3rd, 2020 in the United States at different geo-spatial granularities to conduct our experiments. To understand the impact of this hypothesis, the output of this study was implemented in Tampa General Hospital (TGH) to provide resource demand forecast. RESULTS: The Mean Absolute Percentage Error (MAPE) in the forecast confirmed cases can be 30% higher for modeling at the state-level than aggregating model results at the scale of counties or clusters of counties. Similarly, modeling at a state-level and crafting policy decisions based on them may not be effective - county-level forecasts made by partitioning state-level forecasts are 3x worse for confirmed cases and 20x worse for deaths relative to the same model at the county level. By leveraging these results, TGH was able to accurately allocate clinical resources to tackle COVID-19 cases, continue elective surgical procedures largely uninterrupted and avoid costly construction of overflow capacity in the first two epidemic waves. CONCLUSIONS AND RELEVANCE: Accurate forecasting at the county level requires hyper-local modeling with county resolution. State-level modeling does not accurately predict community spread in smaller sub-regions because state populations are not well mixed, resulting in large prediction errors. Actionable decisions such as deciding whether to cancel planned surgeries or construct overflow capacity require models with local specificity.

Predicting Nocturnal Hypoglycemia from Continuous Glucose Monitoring Data with Extended Prediction Horizon.

Nocturnal hypoglycemia is a serious complication of insulin-treated diabetes, which commonly goes undetected. Continuous glucose monitoring (CGM) devices have enabled prediction of impending nocturnal hypoglycemia, however, prior efforts have been limited to a short prediction horizon (~ 30 minutes). To this end, a nocturnal hypoglycemia prediction model with a 6-hour horizon (midnight-6 am) was developed using a random forest machine- learning model based on data from 10,000 users with more than 1 million nights of CGM data. The model demonstrated an overall nighttime hypoglycemia prediction performance of ROC AUC = 0.84, with AUC = 0.90 for early night (midnight-3 am) and AUC = 0.75 for late night (prediction at midnight, looking at 3-6 am window). While instabilities and the absence of late-night blood glucose patterns introduce predictability challenges, this 6-hour horizon model demonstrates good performance in predicting nocturnal hypoglycemia. Additional study and specific patient-specific features will provide refinements that further ensure safe overnight management of glycemia.

Flow patterns in the jugular veins of pulsatile tinnitus patients.

Pulsatile Tinnitus (PT) is a pulse-synchronous sound heard in the absence of an external source. PT is often related to abnormal flow in vascular structures near the cochlea. One vascular territory implicated in PT is the internal jugular vein (IJV). Using computational fluid dynamics (CFD) based on patient-specific Magnetic Resonance Imaging (MRI), we investigated the flow within the IJV of seven subjects, four symptomatic and three asymptomatic of PT. We found that there were two extreme anatomic types classified by the shape and position of the jugular bulbs: elevated and rounded. PT patients had elevated jugular bulbs that led to a distinctive helical flow pattern within the proximal internal jugular vein. Asymptomatic subjects generally had rounded jugular bulbs that neatly redirected flow from the sigmoid sinus directly into the jugular vein. These two flow patterns were quantified by calculating the length-averaged streamline curvature of the flow within the proximal jugular vein: 130.3±8.1m-1 for geometries with rounded bulbs, 260.7±29.4m-1 for those with elevated bulbs (P<0.005). Our results suggest that variations in the jugular bulb geometry lead to distinct flow patterns that are linked to PT, but further investigation is needed to determine if the vortex pattern is causal to sound generation.

The manifestation of vortical and secondary flow in the cerebral venous outflow tract: An in vivo MR velocimetry study.

Aberrations in flow in the cerebral venous outflow tract (CVOT) have been implicated as the cause of several pathologic conditions including idiopathic intracranial hypertension (IIH), multiple sclerosis (MS), and pulsatile tinnitus (PT). The advent of 4D flow magnetic resonance imaging (4D-flow MRI) has recently allowed researchers to evaluate blood flow patterns in the arterial structures with great success. We utilized similar imaging techniques and found several distinct flow characteristics in the CVOT of subjects with and without lumenal irregularities. We present the flow patterns of 8 out of 38 subjects who have varying heights of the internal jugular bulb and varying lumenal irregularities including stenosis and diverticulum. In the internal jugular vein (IJV) with an elevated jugular bulb (JB), 4Dflow MRI revealed a characteristic spiral flow that was dependent on the level of JB elevation. Vortical flow was also observed in the diverticula of the venous sinuses and IJV. The diversity of flow complexity in the CVOT illustrates the potential importance of hemodynamic investigations in elucidating venous pathologies.

Turbulence intensity measurements using particle image velocimetry in diseased carotid artery models: effect of stenosis severity, plaque eccentricity, and ulceration.

Clinical decision-making for the treatment of patients with diseased carotid artery is mainly based on the severity of the stenosis. However, stenosis severity alone is not a sensitive indicator, and other local factors for the assessment of stroke risk are required. Flow disturbance is of particular interest due to its proven association with increased thromboembolic activities. The objective of this study was to investigate the level of turbulence intensity (TI) with regards to certain geometrical features of the plaque - namely stenosis severity, eccentricity, and ulceration. A family of eight carotid-artery bifurcation models was examined using particle image velocimetry. Results showed a marked difference in turbulence intensity among these models; increasing degree of stenosis severity resulted in increased turbulence intensity, going from 0.12 m/s for mild stenosis to 0.37 m/s for severe stenosis (with concentric geometry). Moreover, independent of stenosis severity, eccentricity led to further elevations in turbulence intensity, increasing TI by 0.05-0.10 m/s over the counterpart concentric plaque. The presence of ulceration (in a 50% eccentric plaque) produced a larger portion of moderate turbulence intensity (~0.10 m/s) compared to the non-ulcerated model, more proximal to the bifurcation apex in the post-stenotic recirculation zone. The effect of plaque eccentricity and ulceration in enhancing the downstream turbulence has potential clinical implications for a more sensitive assessment of stroke risk beyond stenosis severity alone.

In vitro shear stress measurements using particle image velocimetry in a family of carotid artery models: effect of stenosis severity, plaque eccentricity, and ulceration.

Atherosclerotic disease, and the subsequent complications of thrombosis and plaque rupture, has been associated with local shear stress. In the diseased carotid artery, local variations in shear stress are induced by various geometrical features of the stenotic plaque. Greater stenosis severity, plaque eccentricity (symmetry) and plaque ulceration have been associated with increased risk of cerebrovascular events based on clinical trial studies. Using particle image velocimetry, the levels and patterns of shear stress (derived from both laminar and turbulent phases) were studied for a family of eight matched-geometry models incorporating independently varied plaque features - i.e. stenosis severity up to 70%, one of two forms of plaque eccentricity, and the presence of plaque ulceration). The level of laminar (ensemble-averaged) shear stress increased with increasing stenosis severity resulting in 2-16 Pa for free shear stress (FSS) and approximately double (4-36 Pa) for wall shear stress (WSS). Independent of stenosis severity, marked differences were found in the distribution and extent of shear stress between the concentric and eccentric plaque formations. The maximum WSS, found at the apex of the stenosis, decayed significantly steeper along the outer wall of an eccentric model compared to the concentric counterpart, with a 70% eccentric stenosis having 249% steeper decay coinciding with the large outer-wall recirculation zone. The presence of ulceration (in a 50% eccentric plaque) resulted in both elevated FSS and WSS levels that were sustained longer (∼20 ms) through the systolic phase compared to the non-ulcerated counterpart model, among other notable differences. Reynolds (turbulent) shear stress, elevated around the point of distal jet detachment, became prominent during the systolic deceleration phase and was widely distributed over the large recirculation zone in the eccentric stenoses.

Transitional flow analysis in the carotid artery bifurcation by proper orthogonal decomposition and particle image velocimetry.

Blood flow instabilities in the carotid artery bifurcation have been highly correlated to clot formation and mobilization resulting in ischemic stroke. In this work, PIV-measured flow velocities in normal and stenosed carotid artery bifurcation models were analyzed by means of proper orthogonal decomposition (POD). Through POD analysis, transition to more complex flow was visualized and quantified for increasing stenosis severity. While no evidence of transitional flow was seen in the normal model, the 50%-stenosed model started to show characteristics of transitional flow, which became highly evident in the 70% model, with greatest manifestation during the systolic phase of the cardiac cycle. By means of a model comparison, we demonstrate two quantitative measures of the flow complexity through the power-law decay slope of the energy spectrum and the global entropy. The more complex flow in the 70%-stenosed model showed a flatter slope of energy decay (-0.91 compared to -1.34 for 50% stenosis) and higher entropy values (0.26 compared to 0.17). Finally, the minimum temporal resolution required for POD analysis of carotid artery flow was found to be 100 Hz when determined through a more typical energy-mode convergence test, as compared to 400 Hz based on global entropy values.

3-D flow characterization and shear stress in a stenosed carotid artery bifurcation model using stereoscopic PIV technique.

The carotid artery bifurcation is a common site of atherosclerosis which is a major leading cause of ischemic stroke. The impact of stenosis in the atherosclerotic carotid artery is to disturb the flow pattern and produce regions with high shear rate, turbulence, and recirculation, which are key hemodynamic factors associated with plaque rupture, clot formation, and embolism. In order to characterize the disturbed flow in the stenosed carotid artery, stereoscopic PIV measurements were performed in a transparent model with 50% stenosis under pulsatile flow conditions. Simulated ECG gating of the flowrate waveform provides external triggering required for volumetric reconstruction of the complex flow patterns. Based on the three-component velocity data in the lumen region, volumetric shear-stress patterns were derived.