Recent advances in flow MRI

Performance of respiratory gated 4D flow MRI with adaptive k-space reordering in healthy controls and aortic dissection: reproducibility and agreement with 2D phase contrast MRI

A four-dimensional phase-contrast magnetic resonance imaging sequence with respiratory-controlled adaptive k-space reordering (ReCAR-4DPC) offers potential benefits of improved scan efficiency and motion robustness. The purpose of this study was to evaluate the reproducibility of flow measurement using this technique and to compare hemodynamic metrics obtained to two-dimensional phase contrast MRI (2DPC)-derived metrics of the thoracic aorta. ReCAR-4DPC was performed with identical scan parameters in 15 healthy volunteers (6M,9F, mean [range] 37 [23-47] years) and 11 patients with thoracic aortic dissection (6M,5F, 56 [31-81] years) and acquisition time was recorded. Peak systolic velocity (PSV), average flow (AF) and net forward volume (NFV) were quantified by two readers for ReCAR-4DPC at ascending, descending and diaphragmatic aorta levels. Reference standard 2DPC measurements at the same levels were performed by a separate experienced cardiovascular radiologist. ReCAR-4DPC intra-reader agreement, inter-reader agreement, inter-scan repeatability and concordance with 2DPC-derived metrics (all segments combined) were evaluated with Lin's concordance correlation coefficient (LCCC) and reduced major axis regression. The overall average ± SD MRI acquisition time of all subjects was 11:59 ± 3:57 min, with shorter average times (9:37 ± 1:57 min) in healthy volunteers compared to patients (15:13 ± 3:44 min). There was near-perfect intra-reader, inter-reader and inter-scan concordance (LCCC for all metrics > 0.97, > 0.98 and > 0.92 respectively) for ReCAR-4DPC. Concordance with 2DPC was also high (LCCC all > 0.89), with overall minimally lower PSV, AF and NFV values derived from ReCAR-4DPC compared to reference 2DPC derived metrics. ReCAR-4DPC is a reproducible and relatively fast approach for comprehensive measurement of thoracic aortic flow metrics, with robust correlation to conventional 2DPC.

A low-field ceramic magnet design for magnetic resonance

We describe the design of a low-field portable magnet, based on two ceramic magnets, separated by a distance, with their magnetic poles aligned to create a large homogeneous region with a field strength of 425 gauss. Ceramic magnets are an uncommon choice compared to Neodymium Iron Boron magnets for low-field magnetic resonance but are preferable for our purposes to create a homogeneous region at lower field strength. The low cost of large ceramic magnets results in an inexpensive design with a large measurement volume. The magnets rest in a 3D-printed structure, which allows for the magnets to be moved by hand so the experimentalist has control over the field topology. To test the utility of the design, we explored an Overhauser dynamic nuclear polarization experiment with an aqueous solution of 4-Hydroxy-TEMPO. We also explored a simple flow measurement employing the ceramic magnets at a 6-degree pitch, creating a 14.6 gauss/cm constant gradient.

Optimal control flow encoding for time-efficient magnetic resonance velocimetry

Phase contrast velocimetry relies on bipolar gradients to establish a direct and linear relationship between the phase of the magnetic resonance signal, and the corresponding fluid motion. Despite its utility, several limitations and drawbacks have been reported, the most important being the extended echo time due to the encoding after the excitation. In this study, we elucidate a new approach based on optimal control theory that circumvents some of these disadvantages. An excitation pulse, termed FAUCET (flow analysis under controlled encoding transients), is designed to encode velocity into phase already during the radiofrequency excitation. As a result of concurrent excitation and flow encoding, and hence elimination of post-excitation flow encoding, FAUCET achieves a shorter echo time than the conventional method. This achievement is a matter of significance not only because it decreases the loss of signal due to spin-spin relaxation and B₀ inhomogeneity, but also because a shorter echo time is always preferred in order to reduce the dimensionless dephasing parameter and the required residence time of the flowing sample in the detection coil. The method is able to establish a non-linear bijective relationship between phase and velocity, which can be employed to enhance the resolution over a specific range of velocities, for example along flow boundaries. A computational comparison between the phase contrast and optimal control methods reveals that the latter's encoding is more robust against remnant higher-order-moment terms of the Taylor expansion for faster voxels, such as acceleration, jerk, and snap.

Net-phase flow NMR for compact applications

The net phase of the NMR signal is proposed as a robust mechanism for the encoding of fluid flow velocity into phase, showing local bijectivity. While magnitude-based or imaging-based methods suffer from loss of signal, by increasing the flow rate, the present method enables us to maintain the high SNR even for the case of fast flow. In addition, it is shown that a well-engineered flow channel is also necessary, which is not the case for traditional cylindrical flow channels. In this contribution, we report on implementing this approach in a low-cost NMR-based flowmeter for use in a low field (1 T) setting, for example, for monitoring reaction flow industrial processes.

Towards Tomography-Based Real-Time Control of Multiphase Flows: A Proof of Concept in Inline Fluid Separation

The performance of multiphase flow processes is often determined by the distribution of phases inside the equipment. However, controllers in the field are typically implemented based on flow variables, which are simpler to measure, but indirectly connected to performance (e.g., pressure). Tomography has been used in the study of the distribution of phases of multiphase flows for decades, but only recently, the temporal resolution of the technique was sufficient for real-time reconstructions of the flow. Due to the strong connection between the performance and distribution of phases, it is expected that the introduction of tomography to the real-time control of multiphase flows will lead to substantial improvements in the system performance in relation to the current controllers in the field. This paper uses a gas–liquid inline swirl separator to analyze the possibilities and limitations of tomography-based real-time control of multiphase flow processes. Experiments were performed in the separator using a wire-mesh sensor (WMS) and a high-speed camera to show that multiphase flows have two components in their dynamics: one intrinsic to its nonlinear physics, occurring independent of external process disturbances, and one due to process disturbances (e.g., changes in the flow rates of the installation). Moreover, it is shown that the intrinsic dynamics propagate from upstream to inside the separator and can be used in predictive and feedforward control strategies. In addition to the WMS experiments, a proportional–integral feedback controller based on electrical resistance tomography (ERT) was implemented in the separator, with successful results in relation to the control of the distribution of phases and impact on the performance of the process: the capture of gas was increased from 76% to 93% of the total gas with the tomography-based controller. The results obtained with the inline swirl separator are extended in the perspective of the tomography-based control of quasi-1D multiphase flows.

Integrating Pore-Scale Flow MRI and X-ray μCT for Validation of Numerical Flow Simulations in Porous Sedimentary Rocks

Single-phase fluid flow velocity maps in Ketton and Estaillades carbonate rock core plugs are computed at a pore scale, using the lattice Boltzmann method (LBM) simulations performed directly on three-dimensional (3D) X-ray micro-computed tomography (µCT) images (≤ 7 µm spatial resolution) of the core plugs. The simulations are then benchmarked on a voxel-by-voxel and pore-by-pore basis to quantitative, 3D spatially resolved magnetic resonance imaging (MRI) flow velocity maps, acquired at 35 µm isotropic spatial resolution for flow of water through the same rock samples. Co-registration of the 3D experimental and simulated velocity maps and coarse-graining of the simulation to the same resolution as the experimental data allowed the data to be directly compared. First, the results are demonstrated for Ketton limestone rock, for which good qualitative and quantitative agreement was found between the simulated and experimental velocity maps. The flow-carrying microstructural features in Ketton rock are mostly larger than the spatial resolution of the µCT images, so that the segmented images are an adequate representation of the pore space. Second, the flow data are presented for Estaillades limestone, which presents a more heterogeneous case with microstructural features below the spatial resolution of the µCT images. Still, many of the complex flow patterns were qualitatively reproduced by the LBM simulation in this rock, although in some pores, noticeable differences between the LBM and MRI velocity maps were observed. It was shown that 80% of the flow (fractional summed z-velocities within pores) in the Estaillades rock sample is carried by just 10% of the number of macropores, which is an indication of the high structural heterogeneity of the rock; in the more homogeneous Ketton rock, 50% of the flow is carried by 10% of the macropores. By analysing the 3D MRI velocity map, it was found that approximately one-third of the total flow rate through the Estaillades rock is carried by microporosity—a porosity that is not captured at the spatial resolution of the µCT image.

Multiphase Flow Regime Characterization and Liquid Flow Measurement Using Low-Field Magnetic Resonance Imaging

Multiphase flow metering with operationally robust, low-cost real-time systems that provide accuracy across a broad range of produced volumes and fluid properties, is a requirement across a range of process industries, particularly those concerning petroleum. Especially the wide variety of multiphase flow profiles that can be encountered in the field provides challenges in terms of metering accuracy. Recently, low-field magnetic resonance (MR) measurement technology has been introduced as a feasible solution for the petroleum industry. In this work, we study two phase air-water horizontal flows using MR technology. We show that low-field MR technology applied to multiphase flow has the capability to measure the instantaneous liquid holdup and liquid flow velocity using a constant gradient low flip angle CPMG (LFA-CPMG) pulse sequence. LFA-CPMG allows representative sampling of the correlations between liquid holdup and liquid flow velocity, which allows multiphase flow profiles to be characterized. Flow measurements based on this method allow liquid flow rate determination with an accuracy that is independent of the multiphase flow profile observed in horizontal pipe flow for a wide dynamic range in terms of the average gas and liquid flow rates.

Full-Field Comparison of MRV and CFD of Gas Flow through Regular Catalytic Monolithic Structures

Understanding the influence of gas flow maldistribution in honeycombs can be beneficial for the process design in various technical applications. Although recent studies have investigated the effect of maldistribution by comparing the results of numerical simulations with experimental measurements, an exhaustive 3D full-field comparison is still lacking. Such full-field comparisons are required to identify and eliminate possible limitations of numerical and experimental tools. For that purpose, spatially resolved flow patterns were simulated by computational fluid dynamics (CFD) and measured experimentally by non-invasive NMR velocimetry (MRV). While the latter might suffer from a misinterpretation of artefacts, the reliability of CFD is linked to correctly chosen boundary conditions. Here, a full-field numerical and experimental analysis of the gas flow within catalytic honeycombs is presented. The velocity field of thermally polarized methane gas was measured in a regular 3D-printed honeycomb and a commercial monolith using an optimized MRV pulse sequence to enhance the obtained signal-to-noise ratio. A second pulse sequence was used to show local flow propagators along the axial and radial direction of the honeycomb to quantify the contribution of diffusion to mass transport. A quantitative comparison of the axially averaged convective flow as determined by MRV and CFD shows a very good matching with an agreement of ±5% and 10% for printed and commercial samples, respectively. The impact of maldistribution on the gas flow pattern can be observed in both simulation and experiments, confirming the existence of an entrance effect. Gas displacement measurements, however, revealed that diffusive interchannel transport can also contribute to maldistribution, as was shown for the commercial sample. The good agreement between the simulation and experiments underpins the reliability of both methods for studying gas hydrodynamics within opaque monolith structures.

Characterizing pore-scale structure-flow correlations in sedimentary rocks using magnetic resonance imaging

Quantitative, three-dimensional (3D) spatially resolved magnetic resonance flow imaging (flow MRI) methods are presented to characterize structure-flow correlations in a 4-mm-diameter plug of Ketton limestone rock using undersampled k- and q-space data acquisition methods combined with compressed sensing (CS) data reconstruction techniques. The acquired MRI data are coregistered with an X-ray microcomputed tomography (μCT) image of the same rock sample, allowing direct correlation of the structural features of the rock with local fluid transport characteristics. First, 3D velocity maps acquired at 35 μm isotropic spatial resolution showed that the flow was highly heterogeneous, with ∼10% of the pores carrying more than 50% of the flow. Structure-flow correlations were found between the local flow velocities through pores and the size and topology (coordination number) associated with these pores. These data show consistent trends with analogous data acquired for flow through a packing of 4-mm-diameter spheres, which may be due to the microstructure of Ketton rock being a consolidation of approximately spherical grains. Using two-dimensional and 3D visualization of coregistered μCT images and velocity maps, complex pore-scale flow patterns were identified. Second, 3D spatially resolved propagators were acquired at 94 μm isotropic spatial resolution. Flow dispersion within the rock was examined by analyzing each of the 331 776 local propagators as a function of observation time. Again, the heterogeneity of flow within the rock was shown. Quantification of the mean and standard deviation of each of the local propagators showed enhanced mixing occurring within the pore space at longer observation times. These spatially resolved measurements also enable investigation of the length scale of a representative elementary volume. It is shown that for a 4-mm-diameter plug this length scale is not reached.

Equivalent Scalar Stress Formulation Taking into Account Non-Resolved Turbulent Scales

PURPOSE

Cardiovascular engineering includes flows with fluid-dynamical stresses as a parameter of interest. Mechanical stresses are high-risk factors for blood damage and can be assessed by computational fluid dynamics. By now, it is not described how to calculate an adequate scalar stress out of turbulent flow regimes when the whole share of turbulence is not resolved by the simulation method and how this impacts the stress calculation.

METHODS

We conducted direct numerical simulations (DNS) of test cases (a turbulent channel flow and the FDA nozzle) in order to access all scales of flow movement. After validation of both DNS with literature und experimental data using magnetic resonance imaging, the mechanical stress is calculated as a baseline. Afterwards, same flows are calculated using state-of-the-art turbulence models. The stresses are computed for every result using our definition of an equivalent scalar stress, which includes the influence from respective turbulence model, by using the parameter dissipation. Afterwards, the results are compared with the baseline data.

RESULTS

The results show a good agreement regarding the computed stress. Even when no turbulence is resolved by the simulation method, the results agree well with DNS data. When the influence of non-resolved motion is neglected in the stress calculation, it is underpredicted in all cases.

CONCLUSION

With the used scalar stress formulation, it is possible to include information about the turbulence of the flow into the mechanical stress calculation even when the used simulation method does not resolve any turbulence.

Opportunities for Characterizing Geological Flows Using Magnetic Resonance Imaging

Geological flows-from mudslides to volcanic eruptions-are often opaque and consist of multiple interacting phases. Scaled laboratory geological experiments using analog materials have often been limited to optical imaging of flow exteriors or ex situ measurements. Geological flows often include internal phase transitions and chemical reactions that are difficult to image externally. Thus, many physical mechanisms underlying geological flows remain unknown, hindering model development. We propose using magnetic resonance imaging (MRI) to enhance geosciences via non-invasive, in situ measurements of 3D flows. MRI is currently used to characterize the interior dynamics of multiphase flows, distinguishing between different chemical species as well as gas, liquid, and solid phases, while quantitatively measuring concentration, velocity, and diffusion fields. This perspective describes the potential of MRI techniques to image dynamics within scaled geological flow experiments and the potential of technique development for geological samples to be transferred to other disciplines utilizing MRI.

Measurement in opaque flows: a review of measurement techniques for dispersed multiphase flows

A review is presented of measurement techniques to characterise dispersed multiphase flows, which are not accessible by means of conventional optical techniques. The main issues that limit the accuracy and effectiveness of optical techniques are briefly discussed: cross-talk, a reduced signal-to-noise ratio, and (biased) data drop-out. Extensions to the standard optical techniques include the use of fluorescent tracers, refractive index matching, ballistic imaging, structured illumination, and optical coherence tomography. As the first non-optical technique, a brief discussion of electrical capacitance tomography is given. While truly non-invasive, it suffers from a low resolving power. Ultrasound-based techniques have rapidly evolved from Doppler-based profiling to recent 2D approaches using feature tracking. The latter is also suitable for time-resolved flow studies. Magnetic resonance velocimetry can provide time-averaged velocity fields in 3D for the continuous phase. Finally, X-ray imaging is demonstrated to be an important tool to quantify local gas fractions. While potentially very powerful, the impact of the techniques will depend on the development of acquisition and measurement protocols for fluid mechanics, rather than for clinical imaging. This requires systematic development, aided by careful validation experiments. As theoretical predictions for multiphase flows are sparse, it is important to formulate standardised 'benchmark' flows to enable this validation.

Elucidating the ionic liquid distribution in monolithic SILP hydroformylation catalysts by magnetic resonance imaging

Monolithic silicon carbide supported ionic liquid-phase (SILP) Rh-catalysts have very recently been introduced for gas-phase hydroformylation as an important step toward industrial upscaling. This study investigates the monolithic catalyst system in combination with different impregnation procedures with non-invasive magnetic resonance imaging (MRI). The findings were supported by X-ray microtomography (micro-CT) data of the monolithic pore structure and a catalytic performance test of the catalyst system for 1-butene gas-phase hydroformylation. MRI confirmed a homogeneous impregnation of the liquid phase throughout the full cross-section of the cylindrical monoliths. Consistent impregnations from one side to the other of the monoliths were achieved with a stabilizer in the system that helped preventing inhomogeneous rim formation. External influences relevant for industrial application, such as long-term storage and temperature exposure, did not affect the homogeneous liquid-phase distribution of the catalyst. The work elucidates important parameters to improve liquid-phase catalyst impregnation to obtain efficient monolithic catalysts for industrial exploitation in gas-phase hydroformylation as well as other important industrial processes.

The action of the liquid catalyst phase in monolithic silicon carbide supported ionic liquid-phase (SILP) Rh-catalysts provide important insight toward industrial upscaling for gas-phase hydroformylation.

Applications of magnetic resonance imaging in chemical engineering

While there are many techniques to study phenomena that occur in chemical engineering applications, magnetic resonance imaging (MRI) receives increasing scientific interest. Its non-invasive nature and wealth of parameters with the ability to generate functional images and contrast favors the use of MRI for many purposes, in particular investigations of dynamic phenomena, since it is very sensitive to motion. Recent progress in flow-MRI has led to shorter acquisition times and enabled studies of transient phenomena. Reactive systems can easily be imaged if NMR parameters such as relaxation change along the reaction coordinate. Moreover, materials and devices can be examined, such as batteries by mapping the magnetic field around them.

Identification of sampling patterns for high-resolution compressed sensing MRI of porous materials: 'learning' from X-ray microcomputed tomography data

There exists a strong motivation to increase the spatial resolution of magnetic resonance imaging (MRI) acquisitions so that MRI can be used as a microscopy technique in the study of porous materials. This work introduces a method for identifying novel data sampling patterns to achieve undersampling schemes for compressed sensing MRI (CS-MRI) acquisitions, enabling 3D spatial resolutions of 17.6 µm to be achieved. A data-driven learning approach is used to derive k-space undersampling schemes for 3D MRI acquisitions from 3D X-ray microcomputed tomography (µCT) datasets acquired at a higher spatial resolution than can be acquired using MRI. The performance of the new sampling approach was compared to other, well-established sampling strategies using simulated MRI data obtained from high-resolution µCT images of rock core plugs. These simulations were performed for a range of different k-space sampling fractions (0.125-0.375) using images of Ketton limestone. The method was then extended to consideration of imaging Estaillades limestone and Fontainebleau sandstone. The results show that the new sampling approach performs as well as or better than conventional variable density sampling and without need for time-consuming parameter optimisation. Further, a bespoke sampling pattern is produced for each rock type. The novel undersampling strategy was employed to acquire 3D magnetic resonance images of a Ketton limestone rock at spatial resolutions of 35 and 17.6 µm. The ability of the k-space sampling scheme produced using the new approach in enabling reconstruction of the pore space characteristics of the rock was then demonstrated by benchmarking against the pore space statistics obtained from high-resolution µCT data. The MRI data acquired at 17.6 µm resolution gave excellent agreement with the pore size distribution obtained from the X-ray microcomputed tomography dataset, while the pore coordination number distribution obtained from the MRI data was slightly skewed to lower coordination numbers. This approach provides a method of producing a k-space undersampling pattern for MRI acquisition at a spatial resolution for which a fully sampled acquisition at that spatial resolution would be impractically long. The approach can be easily extended to other CS-MRI techniques, such as spatially resolved flow and relaxation time mapping. LAY DESCRIPTION: Magnetic resonance imaging (MRI) is widely used to study the microstructure of, and fluid transport phenomena in porous media relevant for engineering applications. A major application is the study of water and hydrocarbon transport in porous sedimentary rocks, which typically have pore sizes smaller than 100 µm. The spatial resolution of routine MRI acquisitions, however, is limited to several hundred µm due to the relatively low sensitivity of the magnetic resonance method. Therefore, there exists a strong motivation to increase the spatial resolution of MRI by one to two orders of magnitude to be able to study these rocks at a pore scale. This work reports the initial step towards achieving this. Three-dimensional images of rock pore structure are acquired at both 35 and 17.6 µm spatial resolution. In ongoing work, these methods are now being incorporated into magnetic resonance velocity imaging methods, thereby enabling imaging of both pore structure and hydrodynamics at these much higher spatial resolutions than were hitherto possible. Although X-ray microcomputed tomography (µCT) produces high spatial resolution images, it is far more limited in being able to spatially map transport processes (i.e. flow) in porous media. This work reports a strategy for accelerating the image acquisition time such that sufficient signal-to-noise ratio (SNR) is achieved to increase the spatial resolution, that is, the voxel size within which there is sufficient SNR within the resulting image. To achieve this, a technique known as compressed sensing is used which exploits undersampling of the acquired data relative to the standard fully sampled image. In MRI, data are acquired in so-called k-space and Fourier transformed to yield the real space image. The challenge, when undersampling, is to optimise the specific points in k-space that are acquired because these will influence the quality of the resulting image. This work reports a straightforward, robust strategy for identifying the optimal sets of k-space points to acquire. The method introduced uses simulated MRI images calculated from high-resolution µCT images of the rocks of interest, from which optimised MRI sampling patterns are obtained. The method does not require any optimisation of parameters for its implementation, which is a significant advantage compared to other strategies. Moreover, we show that the pore space characteristics of the acquired MRI images are in excellent agreement with the same characteristics obtained from a high-resolution µCT image.

Single-shot velocity mapping by rewinding of velocity encoding with Echo-Planar Imaging

The ability of single-shot NMR imaging methods to follow the time evolution of a velocity distribution within an object is strongly limited by the phase errors accumulated as velocity maps are acquired. In the particular case of Carr-Purcell based sequences combined with Echo Planar Imaging acquisition, phase accumulates through subsequent images, hampering the possibility to acquire several velocity maps, which would be useful to determine transient behavior. In this work, we propose the use of a rewinding velocity encoding module applied after the acquisition of each image during the CPMG echo train. In this way, the first velocity module imparts a velocity dependent phase prior to the image acquisition and the second pair cancels this phase out before the next refocusing radiofrequency pulse is applied. The performance and limits of this method are studied by acquiring 100 images of a co-rotating Couette cell over a period of 1.6 s as a function of the rotation speed. The method is applied to determine the kinematic viscosity of a water/alcohol mixture, which is a relevant topic in many physical, chemical and biological processes.

Magnetohydrodynamic flow imaging of ionic solutions using electrical current injection and MR phase measurements

In this study, a method is proposed to image magnetohydrodynamic (MHD) flow of ionic solutions, which is caused by externally injected electrical current to an imaging media, during MRI scans. A multi-physics (MP) model is created by using the electrical current, laminar flow, and MR equations. The conventional spoiled gradient echo MRI pulse sequence with bipolar flow encoding gradients is utilized to encode the MHD flow. Using the MP model and the MRI pulse sequence, relationship between the MHD flow related phase in the acquired MR signal, the injection current, and the MRI pulse sequence parameters is stated. Numerical simulations and physical experiments are performed to validate the proposed method. The simulation and experimental results are in agreement and show that the MHD flow related MR phase depends on the amplitude and duration of the flow encoding gradient and the injected current. This method may be used to evaluate the MHD flow of conductive liquid media during MRI scans with simultaneous electrical current injections. The MHD flow related MR phase is 1.5 radian for an injected current of 1 mA amplitude, 30 ms duration and a flow encoding gradient amplitude of 24 mT/m. This large MR phase range exhibits potential use of this method for clinical applications such as investigation of highly conductive cerebrospinal fluid (CSF) during clinical use of electrical current based neuromodulation in MRI. However, very high and time varying velocities of typical CSF flow compared to the MHD flow velocities should also be considered.

Magnetic resonance imaging of catalytically relevant processes

The main aim of this article is to provide a state-of-the-art review of the magnetic resonance imaging (MRI) utilization in heterogeneous catalysis. MRI is capable to provide very useful information about both living and nonliving objects in a noninvasive way. The studies of an internal heterogeneous reactor structure by MRI help to understand the mass transport and chemical processes inside the working catalytic reactor that can significantly improve its efficiency. However, one of the serious disadvantages of MRI is low sensitivity, and this obstacle dramatically limits possible MRI application. Fortunately, there are hyperpolarization methods that eliminate this problem. Parahydrogen-induced polarization approach, for instance, can increase the nuclear magnetic resonance signal intensity by four to five orders of magnitude; moreover, the obtained polarization can be stored in long-lived spin states and then transferred into an observable signal in MRI. An in-depth account of the studies on both thermal and hyperpolarized MRI for the investigation of heterogeneous catalytic processes is provided in this review as part of the special issue emphasizing the research performed to date in Russia/USSR.

Diffusion weighted magnetic resonance imaging for temperature measurements in catalyst supports with an axial gas flow

In situ thermometry of catalytic gas phase reactions allows to determine temperature profiles in catalyst beds. Diffusion weighted MRI is proposed as an alternative method for temperature measurements using capillaries filled with different liquids.

Magnetic Resonance Imaging in Situ Visualization of an Electrochemical Reaction under Forced Hydrodynamic Conditions

Magnetic resonance imaging (MRI) has proven to be a powerful tool for the characterization and investigation of in situ chemical reactions. This is more relevant when dealing with complex systems, where the spatial distribution of the species, partition equilibrium, flow patterns, among other factors have a determining effect over mass transport and therefore over the reaction rate. The advantage of MRI is that it provides spatial information in a noninvasive way and does not require any molecular sensor or sample extraction. In this work, MRI is used to fully characterize an electrochemical reaction under forced hydrodynamic conditions. Reaction rates, flow patterns, and quantitative concentration of the chemical species involved are spatially monitored in situ in a complex system that involves metallic pieces and a heterogeneous cementation reaction. Experimental data are compared with numerical simulations.

Acquisition of spatially-resolved displacement propagators using compressed sensing APGSTE-RARE MRI

A method is presented for accelerating the acquisition of spatially-resolved displacement propagators via under-sampling of an Alternating Pulsed Gradient Stimulated Echo - Rapid Acquisition with Relaxation Enhancement (APGSTE-RARE) data acquisition with compressed sensing image reconstruction. The method was demonstrated with respect to the acquisition of 2D spatially-resolved displacement propagators of water flowing through a packed bed of hollow cylinders. The q,k-space was under-sampled according to variable-density pseudo-random sampling patterns. The quality of compressed sensing reconstructions of spatially-resolved propagators at a range of sampling fractions was assessed using the peak signal-to-noise ratio (PSNR) as a quality metric. Propagators of good quality (PSNR 33.2 dB) were reconstructed from only 6.25% of all data points in q,k-space, resulting in a reduction in the data acquisition time from 4 h to 14 min. The spatially-resolved propagators were reconstructed using both the total variation and nuclear norm sparsifying transforms; use of total variation resulted in a slightly higher quality of the reconstructed image in most cases. To illustrate the power of this method to characterise heterogeneous flow in porous media, the method is applied to the characterisation of flow in a vuggy carbonate rock.

Motion-sensitized SPRITE measurements of hydrodynamic cavitation in fast pipe flow

The pressure variations experienced by a liquid flowing through a pipe constriction can, in some cases, result in the formation of a bubble cloud (i.e., hydrodynamic cavitation). Due to the nature of the bubble cloud, it is ideally measured through the use of non-optical and non-invasive techniques; therefore, it is well-suited for study by magnetic resonance imaging. This paper demonstrates the use of Conical SPRITE (a 3D, centric-scan, pure phase-encoding pulse sequence) to acquire time-averaged void fraction and velocity information about hydrodynamic cavitation for water flowing through a pipe constriction.

Magnetic resonance imaging of granular materials

Magnetic Resonance Imaging (MRI) has become one of the most important tools to screen humans in medicine; virtually every modern hospital is equipped with a Nuclear Magnetic Resonance (NMR) tomograph. The potential of NMR in 3D imaging tasks is by far greater, but there is only "a handful" of MRI studies of particulate matter. The method is expensive, time-consuming, and requires a deep understanding of pulse sequences, signal acquisition, and processing. We give a short introduction into the physical principles of this imaging technique, describe its advantages and limitations for the screening of granular matter, and present a number of examples of different application purposes, from the exploration of granular packing, via the detection of flow and particle diffusion, to real dynamic measurements. Probably, X-ray computed tomography is preferable in most applications, but fast imaging of single slices with modern MRI techniques is unmatched, and the additional opportunity to retrieve spatially resolved flow and diffusion profiles without particle tracking is a unique feature.

Adapted MR velocimetry of slow liquid flow in porous media

MR velocimetry of liquid flow in opaque porous filters may play an important role in better understanding the mechanisms of deep bed filtration. With this knowledge, the efficiency of separating the suspended solid particles from the vertically flowing liquid can be improved, and thus a wide range of industrial applications such as wastewater treatment and desalination can be optimized. However, MR velocimetry is challenging for such studies due to the low velocities, the severe B₀ inhomogeneity in porous structures, and the demand for high spatial resolution and an appropriate total measurement time during which the particle deposition will change velocities only marginally. In this work, a modified RARE-based MR velocimetry method is proposed to address these issues for velocity mapping on a deep bed filtration cell. A dedicated RF coil with a high filling factor is constructed considering the limited space available for the vertical cell in a horizontal MR magnet. Several means are applied to optimize the phase contrast RARE MRI pulse sequence for accurately measuring the phase contrast in a long echo train, even in the case of a low B₁ homogeneity. Two means are of particular importance. One uses data acquired with zero flow to correct the phase contrast offsets from gradient imperfections, and the other combines the phase contrast from signals of both odd and even echoes. Results obtained on a 7T preclinical MR scanner indicate that the low velocities in the heterogeneous system can be correctly quantified with high spatial resolution and an adequate total measurement time, enabling future studies on flow during the filtration process.

PFG NMR and Bayesian analysis to characterise non-Newtonian fluids

Many industrial flow processes are sensitive to changes in the rheological behaviour of process fluids, and there therefore exists a need for methods that provide online, or inline, rheological characterisation necessary for process control and optimisation over timescales of minutes or less. Nuclear magnetic resonance (NMR) offers a non-invasive technique for this application, without limitation on optical opacity. We present a Bayesian analysis approach using pulsed field gradient (PFG) NMR to enable estimation of the rheological parameters of Herschel-Bulkley fluids in a pipe flow geometry, characterised by a flow behaviour index n, yield stress τ₀, and consistency factor k, by analysis of the signal in q-space. This approach eliminates the need for velocity image acquisition and expensive gradient hardware. We investigate the robustness of the proposed Bayesian NMR approach to noisy data and reduced sampling using simulated NMR data and show that even with a signal-to-noise ratio (SNR) of 100, only 16 points are required to be sampled to provide rheological parameters accurate to within 2% of the ground truth. Experimental validation is provided through an experimental case study on Carbopol 940 solutions (model Herschel-Bulkley fluids) using PFG NMR at a ¹H resonance frequency of 85.2MHz; for SNR>1000, only 8 points are required to be sampled. This corresponds to a total acquisition time of <60s and represents an 88% reduction in acquisition time when compared to MR flow imaging. Comparison of the shear stress-shear rate relationship, quantified using Bayesian NMR, with non-Bayesian NMR methods demonstrates that the Bayesian NMR approach is in agreement with MR flow imaging to within the accuracy of the measurement. Furthermore, as we increase the concentration of Carbopol 940 we observe a change in rheological characteristics, probably due to shear history-dependent behaviour and the different geometries used. This behaviour highlights the need for online, or inline, rheological characterisation in industrial process applications.

In situ nuclear magnetic resonance microimaging of live biofilms in a microchannel

The firstin situnuclear magnetic resonance microimaging of live biofilms in a transferrable microfluidic platform.

Quantitative velocity distributions via nuclear magnetic resonance flow metering

We demonstrate the use of Tikhonov regularisation as a data inversion technique to determine the velocity distributions of flowing liquid streams. Regularisation is applied to the signal produced by a nuclear magnetic resonance (NMR) flow measurement system consisting of a pre-polarising permanent magnet located upstream of an Earth's magnetic field NMR detection coil. A simple free induction decay (FID) NMR signal is measured for the flowing stream in what is effectively a 'time-of-flight' measurement. The FID signal is then modelled as a function of fluid velocity and acquisition time, enabling determination of the velocity probability distributions via regularisation. The mean values of these velocity distributions were successfully validated against in-line rotameters. The ability to quantify multi-modal velocity distributions was also demonstrated using a two-pipe system.

Diagnostic Efficacy of 3-T MRI for Knee Injuries Using Arthroscopy as a Reference Standard: A Meta-Analysis

OBJECTIVE

The objectives of our study were to assess the evidence for the diagnostic efficacy of 3-T MRI for meniscal and anterior cruciate ligament (ACL) injuries in the knee using arthroscopy as the reference standard and to compare these results with the results of a previous meta-analysis assessing 1.5-T MRI.

MATERIALS AND METHODS

The online Cochrane Library, MEDLINE, and PubMed databases were searched using the following terms: MRI AND ((3 OR three) AND (Tesla OR T)) AND knee AND arthroscopy AND (menisc* OR ligament). Patient demographics, patient characteristics, MRI scanning details, and diagnostic results were investigated. The methodologic quality of the included studies was assessed using the revised Quality Assessment of Diagnostic Accuracy Studies (QUADAS-2) tool. A meta-analysis of studies using 3-T MRI was performed, and the results were compared with a previous meta-analysis of studies using 1.5-T MRI.

RESULTS

One hundred one studies were identified by the search strategy, and 13 studies were included in our review. Twelve studies were considered to have level 1b evidence, and one study was considered to have level 2b evidence. All 13 studies had high methodologic integrity and low risk of bias using the QUADAS-2 tool. The studies included 1197 patients with a mean age of 41.9 years. Ten of the 13 studies were eligible for meta-analysis. The mean sensitivity and mean specificity of 3-T MRI for knee injuries by location were as follows: medial meniscus, 0.94 (95% CI, 0.91-0.96) and 0.79 (95% CI, 0.75-0.83), respectively; lateral meniscus, 0.81 (95% CI, 0.75-0.85) and 0.87 (95% CI, 0.84-0.89); and ACL, 0.92 (95% CI, 0.83-0.96) and 0.99 (95% CI, 0.96-1.00). The specificity of 3-T MRI for injuries of the lateral meniscus was significantly lower than that of 1.5-T MRI (p = 0.0013).

CONCLUSION

This study does not provide evidence that 3-T scanners have superior diagnostic efficacy for meniscal damage and ACL integrity when compared with previous studies of 1.5-T machines.

11-interval PFG pulse sequence for improved measurement of fast velocities of fluids with high diffusivity in systems with short T2(∗)

Magnetic resonance (MR) was used to measure SF6 gas velocities in beds filled with particles of 1.1 mm and 0.5 mm in diameter. Four pulse sequences were tested: a traditional spin echo pulse sequence, the 9-interval and 13-interval pulse sequence of Cotts et al. (1989) and a newly developed 11-interval pulse sequence. All pulse sequences measured gas velocity accurately in the region above the particles at the highest velocities that could be achieved (up to 0.1 ms(-1)). The spin echo pulse sequence was unable to measure gas velocity accurately in the bed of particles, due to effects of background gradients, diffusivity and acceleration in flow around particles. The 9- and 13-interval pulse sequence measured gas velocity accurately at low flow rates through the particles (expected velocity <0.06 ms(-1)), but could not measure velocity accurately at higher flow rates. The newly developed 11-interval pulse sequence was more accurate than the 9- and 13-interval pulse sequences at higher flow rates, but for velocities in excess of 0.1 ms(-1) the measured velocity was lower than the expected velocity. The increased accuracy arose from the smaller echo time that the new pulse sequence enabled, reducing selective attenuation of signal from faster moving nuclei.

Measurement with microscopic MRI and simulation of flow in different aneurysm models

PURPOSE

The impact and the development of aneurysms depend to a significant degree on the exchange of liquid between the regular vessel and the pathological extension. A better understanding of this process will lead to improved prediction capabilities. The aim of the current study was to investigate fluid-exchange in aneurysm models of different complexities by combining microscopic magnetic resonance measurements with numerical simulations. In order to evaluate the accuracy and applicability of these methods, the fluid-exchange process between the unaltered vessel lumen and the aneurysm phantoms was analyzed quantitatively using high spatial resolution.

METHODS

Magnetic resonance flow imaging was used to visualize fluid-exchange in two different models produced with a 3D printer. One model of an aneurysm was based on histological findings. The flow distribution in the different models was measured on a microscopic scale using time of flight magnetic resonance imaging. The whole experiment was simulated using fast graphics processing unit-based numerical simulations. The obtained simulation results were compared qualitatively and quantitatively with the magnetic resonance imaging measurements, taking into account flow and spin-lattice relaxation.

RESULTS

The results of both presented methods compared well for the used aneurysm models and the chosen flow distributions. The results from the fluid-exchange analysis showed comparable characteristics concerning measurement and simulation. Similar symmetry behavior was observed. Based on these results, the amount of fluid-exchange was calculated. Depending on the geometry of the models, 7% to 45% of the liquid was exchanged per second.

CONCLUSIONS

The result of the numerical simulations coincides well with the experimentally determined velocity field. The rate of fluid-exchange between vessel and aneurysm was well-predicted. Hence, the results obtained by simulation could be validated by the experiment. The observed deviations can be caused by the noise in the measurement and by the limited resolution of the simulation. The resulting differences are small enough to allow reliable predictions of the flow distribution in vessels with stents and for pulsed blood flow.

Real-time fluid transport characterization through direct acquisition of the averaged propagator

Brownian motion (diffusion) and coherent flow are fundamental for many natural and artificial processes. Therefore, its accurate measurement and description is highly desirable in many areas of science, engineering, and technology. Currently available methods commonly involve multiple experiments and substantial processing of acquired data. This contribution proposes a theoretical and experimental framework that enables one to directly examine the dynamics of fluid matter subject to diffusion and flow through the acquisition of the so-called averaged propagator. This statistical function holds all information on particle mobility due to flow and diffusion averaged over the observed fluid. The proposed method is based on a single instantaneous nuclear magnetic resonance measurement event. It also removes the need for data postprocessing by capturing the averaged propagator directly as the acquired signal, which enables the monitoring of diffusion and flow in real time.

Characterising the rheology of non-Newtonian fluids using PFG-NMR and cumulant analysis

Conventional rheological characterisation using nuclear magnetic resonance (NMR) typically utilises spatially-resolved measurements of velocity. We propose a new approach to rheometry using pulsed field gradient (PFG) NMR which readily extends the application of MR rheometry to single-axis gradient hardware. The quantitative use of flow propagators in this application is challenging because of the introduction of artefacts during Fourier transform, which arise when realistic sampling strategies are limited by experimental and hardware constraints and when particular spatial and temporal resolution are required. The method outlined in this paper involves the cumulant analysis of the acquisition data directly, thereby preventing the introduction of artefacts and reducing data acquisition times. A model-dependent approach is developed to enable the pipe-flow characterisation of fluids demonstrating non-Newtonian power-law rheology, involving the use of an analytical expression describing the flow propagator in terms of the flow behaviour index. The sensitivity of this approach was investigated and found to be robust to the signal-to-noise ratio (SNR) and number of acquired data points, enabling an increase in temporal resolution defined by the SNR. Validation of the simulated results was provided by an experimental case study on shear-thinning aqueous xanthan gum solutions, whose rheology could be accurately characterised using a power-law model across the experimental shear rate range of 1-100 s(-1). The flow behaviour indices calculated using this approach were observed to be within 8% of those obtained using spatially-resolved velocity imaging and within 5% of conventional rheometry. Furthermore, it was shown that the number of points sampled could be reduced by a factor of 32, when compared to the acquisition of a volume-averaged flow propagator with 128 gradient increments, without negatively influencing the accuracy of the characterisation, reducing the acquisition time to only 3% of its original value.

Accelerated three-dimensional cine phase contrast imaging using randomly undersampled echo planar imaging with compressed sensing reconstruction

The aim of this study was to implement and evaluate an accelerated three-dimensional (3D) cine phase contrast MRI sequence by combining a randomly sampled 3D k-space acquisition sequence with an echo planar imaging (EPI) readout. An accelerated 3D cine phase contrast MRI sequence was implemented by combining EPI readout with randomly undersampled 3D k-space data suitable for compressed sensing (CS) reconstruction. The undersampled data were then reconstructed using low-dimensional structural self-learning and thresholding (LOST). 3D phase contrast MRI was acquired in 11 healthy adults using an overall acceleration of 7 (EPI factor of 3 and CS rate of 3). For comparison, a single two-dimensional (2D) cine phase contrast scan was also performed with sensitivity encoding (SENSE) rate 2 and approximately at the level of the pulmonary artery bifurcation. The stroke volume and mean velocity in both the ascending and descending aorta were measured and compared between two sequences using Bland-Altman plots. An average scan time of 3 min and 30 s, corresponding to an acceleration rate of 7, was achieved for 3D cine phase contrast scan with one direction flow encoding, voxel size of 2 × 2 × 3 mm(3) , foot-head coverage of 6 cm and temporal resolution of 30 ms. The mean velocity and stroke volume in both the ascending and descending aorta were statistically equivalent between the proposed 3D sequence and the standard 2D cine phase contrast sequence. The combination of EPI with a randomly undersampled 3D k-space sampling sequence using LOST reconstruction allows a seven-fold reduction in scan time of 3D cine phase contrast MRI without compromising blood flow quantification.