Phase-contrast velocimetry with hyperpolarized3He for in vitro and in vivo characterization of airflow

Human upper-airway respiratory airflow: In vivo comparison of computational fluid dynamics simulations and hyperpolarized 129Xe phase contrast MRI velocimetry

Computational fluid dynamics (CFD) simulations of respiratory airflow have the potential to change the clinical assessment of regional airway function in health and disease, in pulmonary medicine and otolaryngology. For example, in diseases where multiple sites of airway obstruction occur, such as obstructive sleep apnea (OSA), CFD simulations can identify which sites of obstruction contribute most to airway resistance and may therefore be candidate sites for airway surgery. The main barrier to clinical uptake of respiratory CFD to date has been the difficulty in validating CFD results against a clinical gold standard. Invasive instrumentation of the upper airway to measure respiratory airflow velocity or pressure can disrupt the airflow and alter the subject's natural breathing patterns. Therefore, in this study, we instead propose phase contrast (PC) velocimetry magnetic resonance imaging (MRI) of inhaled hyperpolarized 129Xe gas as a non-invasive reference to which airflow velocities calculated via CFD can be compared. To that end, we performed subject-specific CFD simulations in airway models derived from 1H MRI, and using respiratory flowrate measurements acquired synchronously with MRI. Airflow velocity vectors calculated by CFD simulations were then qualitatively and quantitatively compared to velocity maps derived from PC velocimetry MRI of inhaled hyperpolarized 129Xe gas. The results show both techniques produce similar spatial distributions of high velocity regions in the anterior-posterior and foot-head directions, indicating good qualitative agreement. Statistically significant correlations and low Bland-Altman bias between the local velocity values produced by the two techniques indicates quantitative agreement. This preliminary in vivo comparison of respiratory airway CFD and PC MRI of hyperpolarized 129Xe gas demonstrates the feasibility of PC MRI as a technique to validate respiratory CFD and forms the basis for further comprehensive validation studies. This study is therefore a first step in the pathway towards clinical adoption of respiratory CFD.

Improved pulmonary 129 Xe ventilation imaging via 3D-spiral UTE MRI

PURPOSE

Hyperpolarized ¹²⁹ Xe MRI characterizes regional lung ventilation in a variety of disease populations, with high sensitivity to airway obstruction in early disease. However, ventilation images are usually limited to a single breath-hold and most-often acquired using gradient-recalled echo sequences with thick slices (~10-15 mm), which increases partial-volume effects, limits ability to observe small defects, and suffers from imperfect slice selection. We demonstrate higher-resolution ventilation images, in shorter breath-holds, using FLORET (Fermat Looped ORthogonally Encoded Trajectories), a center-out 3D-spiral UTE sequence.

METHODS

In vivo human adult (N = 4; 2 healthy, 2 with cystic fibrosis) ¹²⁹ Xe images were acquired using 2D gradient-recalled echo, 3D radial, and FLORET. Each sequence was acquired at its highest possible resolution within a 16-second breath-hold with a minimum voxel dimension of 3 mm. Images were compared using ¹²⁹ Xe ventilation defect percentage, SNR, similarity coefficients, and vasculature cross-sections.

RESULTS

The FLORET sequence obtained relative normalized SNR, 40% greater than 2D gradient-recalled echo (P = .012) and 26% greater than 3D radial (P = .067). Moreover, the FLORET images were acquired with 3-fold-higher nominal resolution in a 15% shorter breath-hold. Finally, vasculature was less prominent in FLORET, likely due to diminished susceptibility-induced dephasing at shorter TEs afforded by UTE sequences.

CONCLUSION

The FLORET sequence yields higher SNR for a given resolution with a shorter breath-hold than traditional ventilation imaging techniques. This sequence more accurately measures ventilation abnormalities and enables reduced scan times in patients with poor compliance and severe lung disease.

Effect of upper airway on tracheobronchial fluid dynamics

The upper airways play a significant role in the tracheal flow dynamics. Despite many previous studies, however, the effect of the upper airways on the ventilation distribution in distal airways has remained a challenge. The aim of this study is to experimentally and computationally investigate the dynamic behaviour in the intratracheal flow induced by the upper respiratory tract and to assess its influence on the subsequent tributaries. Patient-specific images from 2 different modalities (magnetic resonance imaging of the upper airways and computed tomography of the lower airways) were segmented and combined. An experimental phantom of patient-specific airways (including the oral cavity, larynx, trachea, down to generations 6-8) was generated using 3D printing. The flow velocities in this phantom model were measured by the flow-sensitised phase contrast magnetic resonance imaging technique and compared with the computational fluid dynamics simulations. Both experimental and computational results show a good agreement in the time-averaged velocity fields as well as fluctuating velocity. The flows in the proximal trachea were complex and unsteady under both lower- and higher-flow rate conditions. Computational fluid dynamics simulations were also performed with an airways model without the upper airways. Although the flow near the carina remained unstable only when the inflow rate was high, the influence of the upper airways caused notable changes in distal flow distributions when the 2 airways models were compared with and without the upper airways. The results suggest that the influence of the upper airways should be included in the respiratory flow assessment as the upper airways extensively affect the flows in distal airways and consequent ventilation distribution in the lungs.

Fast dynamic ventilation MRI of hyperpolarized 129 Xe using spiral imaging

Purpose

To develop and optimize a rapid dynamic hyperpolarized ¹²⁹Xe ventilation (DXeV) MRI protocol and investigate the feasibility of capturing pulmonary signal‐time curves in human lungs.

Theory and Methods

Spiral k‐space trajectories were designed with the number of interleaves N_int = 1, 2, 4, and 8 corresponding to voxel sizes of 8 mm, 5 mm, 4 mm, and 2.5 mm, respectively, for field of view = 15 cm. DXeV images were acquired from a gas‐flow phantom to investigate the ability of N_int = 1, 2, 4, and 8 to capture signal‐time curves. A finite element model was constructed to investigate gas‐flow dynamics corroborating the experimental signal‐time curves. DXeV images were also carried out in six subjects (three healthy and three chronic obstructive pulmonary disease subjects).

Results

DXeV images and numerical modelling of signal‐time curves permitted the quantification of temporal and spatial resolutions for different numbers of spiral interleaves. The two‐interleaved spiral (N_int = 2) was found to be the most time‐efficient to obtain DXeV images and signal‐time curves of whole lungs with a temporal resolution of 624 ms for 13 slices. Signal‐time curves were well matched in three healthy volunteers. The Spearman's correlations of chronic obstructive pulmonary disease subjects were statistically different from three healthy subjects (P < 0.05).

Conclusion

The N_int = 2 spiral demonstrates the successful acquisition of DXeV images and signal‐time curves in healthy subjects and chronic obstructive pulmonary disease patients. Magn Reson Med 79:2597–2606, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Observation of cardiogenic flow oscillations in healthy subjects with hyperpolarized 3He MRI

Recently, dynamic MRI of hyperpolarized (3)He during inhalation revealed an alternation of the image intensity between left and right lungs with a cardiac origin (Sun Y, Butler JP, Ferrigno M, Albert MS, Loring SH. Respir Physiol Neurobiol 185: 468-471, 2013). This effect is investigated further using dynamic and phase-contrast flow MRI with inhaled (3)He during slow inhalations (flow rate ∼100 ml/s) to elucidate airflow dynamics in the main lobes in six healthy subjects. The ventilation MR signal and gas inflow in the left lower lobe (LLL) of the lungs were found to oscillate clearly at the cardiac frequency in all subjects, whereas the MR signals in the other parts of the lungs had a similar oscillatory behavior but were smaller in magnitude and in anti-phase to the signal in the left lower lung. The airflow in the main bronchi showed periodic oscillations at the frequency of the cardiac cycle. In four of the subjects, backflows were observed for a short period of time of the cardiac cycle, demonstrating a pendelluft effect at the carina bifurcation between the left and right lungs. Additional (1)H structural MR images of the lung volume and synchronized ECG recording revealed that maximum inspiratory flow rates in the LLL of the lungs occurred during systole when the corresponding left lung volume increased, whereas the opposite effect was observed during diastole, with gas flow redirected to the other parts of the lung. In conclusion, cardiogenic flow oscillations have a significant effect on regional gas flow and distribution within the lungs.

Phase-contrast helium-3 MRI of aerosol deposition in human airways

One of the key challenges in the study of health-related aerosols is predicting and monitoring sites of particle deposition in the respiratory tract. The potential health risks of ambient exposure to environmental or workplace aerosols and the beneficial effects of medical aerosols are strongly influenced by the site of aerosol deposition along the respiratory tract. Nuclear medicine is the only current modality that combines quantification and regional localization of aerosol deposition, and this technique remains limited by its spatial and temporal resolutions and by patient exposure to radiation. Recent work in MRI has shed light on techniques to quantify micro-sized magnetic particles in living bodies by the measurement of associated static magnetic field variations. With regard to lung MRI, hyperpolarized helium-3 may be used as a tracer gas to compensate for the lack of MR signal in the airways, so as to allow assessment of pulmonary function and morphology. The extrathoracic region of the human respiratory system plays a critical role in determining aerosol deposition patterns, as it acts as a filter upstream from the lungs. In the present work, aerosol deposition in a mouth-throat phantom was measured using helium-3 MRI and compared with single-photon emission computed tomography. By providing high sensitivity with high spatial and temporal resolutions, phase-contrast helium-3 MRI offers new insights for the study of particle transport and deposition.

Phase-contrast velocity mapping for highly diffusive fluids: optimal bipolar gradient pulse parameters for hyperpolarized helium-3

PURPOSE

In MR-velocity phase-contrast measurements, increasing the encoding bipolar gradient, i.e., decreasing the field of speed, usually improves measurement precision. However, in gases, fast diffusion during the bipolar gradient pulses may dramatically decrease the signal-to-noise ratio, thus degrading measurement precision. These two effects are contradictory. This work aims at determining the optimal sequence parameters to improve the velocity measurement precision.

THEORY AND METHODS

This work presents the theoretical optimization of bipolar gradient parameters (duration and amplitude) to improve velocity measurement precision. An analytical approximation is given as well as a numerical optimization. It is shown that the solution depends on the diffusion coefficient and T2 *. Experimental validation using hyperpolarized (3) He diluted in various buffer gases ((4) He, N2 , and SF6 ) is presented at 1.5 Tesla (T) in a straight pipe.

RESULTS

Excellent agreement was found with the theoretical results for prediction of optimal field of speed and good agreement was found for the precision in measured velocity, but for SF6 buffered gas.

CONCLUSION

The theoretical predictions were validated, providing a way to optimize velocity mapping in gases.

Phase-contrast MRI and CFD modeling of apparent ³He gas flow in rat pulmonary airways

Phase-contrast (PC) magnetic resonance imaging (MRI) with hyperpolarized ³He is potentially useful for developing and testing patient-specific models of pulmonary airflow. One challenge, however, is that PC-MRI provides apparent values of local ³He velocity that not only depend on actual airflow but also on gas diffusion. This not only blurs laminar flow patterns in narrow airways but also introduces anomalous airflow structure that reflects gas-wall interactions. Here, both effects are predicted in a live rat using computational fluid dynamics (CFD), and for the first time, simulated patterns of apparent ³He gas velocity are compared with in vivo PC-MRI. Results show (1) that correlations (R²) between measured and simulated airflow patterns increase from 0.23 to 0.79 simply by accounting for apparent ³He transport, and (2) that remaining differences are mainly due to uncertain airway segmentation and partial volume effects stemming from relatively coarse MRI resolution. Higher-fidelity testing of pulmonary airflow predictions should therefore be possible with future imaging improvements.

Hyperpolarized gas MR Imaging of the lung: current status as a research tool

Hyperpolarized gas magnetic resonance imaging has been explored extensively as a promising tool for the quantitative evaluation of regional pulmonary pathophysiology. This noninvasive technique is capable of providing both structural information down to the level of the alveolar microstructure and functional information, such as dynamic ventilation, intrapulmonary partial pressure of oxygen, and alveolar surface area. This study reviews the role of hyperpolarized 3-helium and 129-xenon magnetic resonance imaging in this research.

MR imaging of apparent 3He gas transport in narrow pipes and rodent airways

High sensitivity makes hyperpolarized (3)He an attractive signal source for visualizing gas flow with magnetic resonance (MR) imaging. Its rapid Brownian motion, however, can blur observed flow lamina and alter measured diffusion rates when excited nuclei traverse shear-induced velocity gradients during data acquisition. Here, both effects are described analytically, and predicted values for measured transport during laminar flow through a straight, 3.2-mm diameter pipe are validated using two-dimensional (2D) constant-time images of different binary gas mixtures. Results show explicitly how measured transport in narrow conduits is characterized by apparent values that depend on underlying gas dynamics and imaging time. In ventilated rats, this is found to obscure acquired airflow images. Nevertheless, flow splitting at airway branches is still evident and use of 3D vector flow mapping is shown to reveal surprising detail that highlights the correlation between gas dynamics and lung structure.

In vitro validation of computational fluid dynamic simulation in human proximal airways with hyperpolarized 3He magnetic resonance phase-contrast velocimetry

Computational fluid dynamics (CFD) and magnetic resonance (MR) gas velocimetry were concurrently performed to study airflow in the same model of human proximal airways. Realistic in vivo-based human airway geometry was segmented from thoracic computed tomography. The three-dimensional numerical description of the airways was used for both generation of a physical airway model using rapid prototyping and mesh generation for CFD simulations. Steady laminar inspiratory experiments (Reynolds number Re = 770) were performed and velocity maps down to the fourth airway generation were extracted from a new velocity mapping technique based on MR velocimetry using hyperpolarized (3)He gas. Full two-dimensional maps of the velocity vector were measured within a few seconds. Numerical simulations were carried out with the experimental flow conditions, and the two sets of data were compared between the two modalities. Flow distributions agreed within 3%. Main and secondary flow velocity intensities were similar, as were velocity convective patterns. This work demonstrates that experimental and numerical gas velocity data can be obtained and compared in the same complex airway geometry. Experiments validated the simulation platform that integrates patient-specific airway reconstruction process from in vivo thoracic scans and velocity field calculation with CFD, hence allowing the results of this numerical tool to be used with confidence in potential clinical applications for lung characterization. Finally, this combined numerical and experimental approach of flow assessment in realistic in vivo-based human airway geometries confirmed the strong dependence of airway flow patterns on local and global geometrical factors, which could contribute to gas mixing.