Benedict Newling

Flow imaging of fluids in porous media by magnetization prepared centric-scan SPRITE

MRI has considerable potential as a non-destructive probe of porous media, permitting rapid quantification of local fluid content and the possibility of local flow visualization and quantification. In this work we explore a general approach to flow velocity measurement in porous media by combining Cotts pulsed field gradient flow encoding with SPRITE MRI. This technique permits facile and accurate flow and dispersion coefficient mapping of fluids in porous media. This new approach has proven to be robust in characterizing fluid behavior. This method is illustrated through measurements of flow in pipes, flow in sand packs and flow in porous reservoir rocks. Spatially resolved flow maps and local fluid velocity distribution were acquired.

SPRITE MRI of bubbly flow in a horizontal pipe.

Bubble flow is characterised by numerous phase interfaces and turbulence, leading to fast magnetic resonance signal decay and artefacts in spin-warp imaging. In this paper, the SPRITE MRI pulse sequence, with its potential for very short encoding times, is demonstrated as an ideal technique for studying such dynamic systems. It has been used to acquire liquid velocity and relative intensity maps of two-phase gas-liquid dispersed bubble flow in a horizontal pipe at a liquid Reynolds number of 14,500. The fluids were air and water and a turbulence grid was used to generate a dispersed bubble flow pattern. The SPRITE technique shows promise for future research in gas-liquid flow.

Mapping B1-induced eddy current effects near metallic structures in MR images: A comparison of simulation and experiment

Magnetic resonance imaging (MRI) in the presence of metallic structures is very common in medical and non-medical fields. Metallic structures cause MRI image distortions by three mechanisms: (1) static field distortion through magnetic susceptibility mismatch, (2) eddy currents induced by switched magnetic field gradients and (3) radio frequency (RF) induced eddy currents. Single point ramped imaging with T1 enhancement (SPRITE) MRI measurements are largely immune to susceptibility and gradient induced eddy current artifacts. As a result, one can isolate the effects of metal objects on the RF field. The RF field affects both the excitation and detection of the magnetic resonance (MR) signal. This is challenging with conventional MRI methods, which cannot readily separate the three effects. RF induced MRI artifacts were investigated experimentally at 2.4 T by analyzing image distortions surrounding two geometrically identical metallic strips of aluminum and lead. The strips were immersed in agar gel doped with contrast agent and imaged employing the conical SPRITE sequence. B1 mapping with pure phase encode SPRITE was employed to measure the B1 field around the strips of metal. The strip geometry was chosen to mimic metal electrodes employed in electrochemistry studies. Simulations are employed to investigate the RF field induced eddy currents in the two metallic strips. The RF simulation results are in good agreement with experimental results. Experimental and simulation results show that the metal has a pronounced effect on the B1 distribution and B1 amplitude in the surrounding space. The electrical conductivity of the metal has a minimal effect.

Measurement of pulsatile flow using MRI and a Bayesian technique of probability analysis

This work shows that complete spatial information of periodic pulsatile fluid flows can be rapidly obtained by Bayesian probability analysis of flow encoded magnetic resonance imaging data. These data were acquired as a set of two-dimensional images (complete two-dimensional sampling of k-space or reciprocal position space) but with a sparse (six point) and nonuniform sampling of q-space or reciprocal displacement space. This approach enables more precise calculation of fluid velocity to be achieved than by conventional two q-sample phase encoding of velocities, without the significant time disadvantage associated with the complete flow measurement required for Fourier velocity imaging. For experimental comparison with the Bayesian analysis applied to nonuniformly sampled q-space data, a Fourier velocity imaging technique was used with one-dimensional spatial encoding within a selected slice and a uniform sampling of q-space using 64 values of the pulsed gradients to encode fluid flow. Because the pulsatile flows were axially symmetric within the resolution of the experiment, the radial variation of fluid velocity, in the direction of the pulsed gradients, was reconstructed from one-dimensional spatial projections of the velocity by exploiting the central slice theorem. Data were analysed for internal consistency using linearised flow theories. The results show that nonuniform q-space sampling followed by Bayesian probability analysis is at least as accurate as the combined uniform q-space sampling with Fourier velocity imaging and projection reconstruction method. Both techniques give smaller errors than a two-point sampling of q-space (the conventional flow encoding experiment).

Local T2 measurement employing longitudinal Hadamard encoding and adiabatic inversion pulses in porous media.

Band selective adiabatic inversion radio frequency pulses were employed for multi-slice T2 distribution measurements in porous media samples. Multi-slice T2 measurement employing longitudinal Hadamard encoding has an inherent sensitivity advantage over slice-by-slice local T2 measurements. The slice selection process is rendered largely immune to B1 variation by employing hyperbolic secant adiabatic inversion pulses, which simultaneously invert spins in several well-defined slices. While Hadamard encoding is well established for local spectroscopy, the current work is the first use of Hadamard encoding for local T2 measurement.

Short, shaped pulses in a large magnetic field gradient

A number of materials MR developments require that measurements be made in a large magnetic field gradient, including unilateral (single-sided) magnet designs for portability and open access. In such cases, all radiofrequency (RF) pulses are slice selective. Typically, little effort is made to tailor the shape of the selected slice, because shaped RF excitations are viewed as too lengthy in duration to be useful in materials MRI, where signal lifetimes are mostly less than 1 ms. We compare measured magnetization responses to various standard shaped pulses under extreme conditions of application (approximately 30 micros duration, offset frequencies up to 0.3 MHz, and in the presence of a 13 T/m permanent magnetic field gradient). We discuss the feasibility of their implementation for materials MRI in a large gradient, including the difficulty of choosing optimized pulse area, and propose viable solutions.

Local diffusion and diffusion-T2 distribution measurements in porous media

Slice-selective pulsed field gradient (PFG) and PFG-T2 measurements are developed to measure spatially-resolved molecular diffusion and diffusion-T2 distributions. A spatially selective adiabatic inversion pulse was employed for slice-selection. The slice-selective pulse is able to select a coarse slice, on the order of 1cm, at an arbitrary position in the sample. The new method can be employed to characterize oil-water mixtures in porous media. The new technique has an inherent sensitivity advantage over phase encoding imaging based methods due to signal being localized from a thick slice. The method will be advantageous for magnetic resonance of porous media at low field where sensitivity is problematic. Experimental CPMG data, following PFG diffusion measurement, were compromised by a transient ΔB0(t) field offset. The off resonance effects of ΔB0(t) were examined by simulation. The ΔB0 offset artifact in D-T2 distribution measurements may be avoided by employing real data, instead of magnitude data.

Shear effects on the diffusive movement of oil in triacylglycerol networks

Oil migration is the foremost contributor to the quality loss of many high fat food products; thus, developing an understanding on how this phenomenon takes place in food systems is crucial for the food industry. Diffusion of triacylglycerols (TAG) through the fat network has often been modeled using simplified solutions to Ficks second law. However, a drawback of the use of diffusion models is the lack of sensitivity toward some microstructural characteristics of the matrix. This work reports the study of molecular and macroscopic oil diffusion coefficients of fat crystal networks using nuclear magnetic resonance measurements and from gravimetric determinations of oil loss, respectively. Blends of fully hydrogenated soybean oil (FHSO) in soybean oil (SO) were crystallized statically, under laminar shear rates of 30 and 240 s−1 at different wall temperatures (−10, 0, 20 °C). Another batch of similar blends prepared with and without emulsifier was crystallized in a scraped surface heat exchanger. The results showed a good correlation between diffusion coefficients obtained using both methodologies only when samples exhibited low oil loss, particularly in blends crystallized statically or under mild-shearing conditions. This work suggests that although the Ziegleder model has great advantages and provides valuable insight into the oil migration in fat matrices, there remains a large need for further evaluating its suitability, principally in cases where crystalline structure is severely affected by processing conditions and where oil migration mechanisms other than diffusion may be involved.

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.

Local T 1 -T 2 distribution measurements in porous media

A novel slice-selective T1-T2 measurement is proposed to measure spatially resolved T1-T2 distributions. An adiabatic inversion pulse is employed for slice-selection. The slice-selective pulse is able to select a quasi-rectangular slice, on the order of 1 mm, at an arbitrary position within the sample.The method does not employ conventional selective excitation in which selective excitation is often accomplished by rotation of the longitudinal magnetization in the slice of interest into the transverse plane, but rather a subtraction based on CPMG data acquired with and without adiabatic inversion slice selection. T1 weighting is introduced during recovery from the inversion associated with slice selection. The local T1-T2 distributions measured are of similar quality to bulk T1-T2 measurements. The new method can be employed to characterize oil-water mixtures and other fluids in porous media. The method is beneficial when a coarse spatial distribution of the components is of interest.