Beat phenomena of oscillating readouts

PURPOSE

To demonstrate slowly varying, erroneous magnetic field gradients for oscillating readouts due to the mechanically resonant behavior of gradient systems.

METHODS

Projections of a static phantom were acquired using a one-dimensional (1D) EPI sequence with varying EPI frequencies ranging from 1121 to 1580 Hz on clinical 3T systems (30 mT/m, 200 T/m/s). Phase due to static B0 inhomogeneities was eliminated by a complex division of two separate scans with different polarities of the EPI readout. The temporal evolution of phase was evaluated and related to the mechanical resonances of the gradient systems derived from the gradient modulation transfer function. Additionally, the impact of temporally varying mechanical resonance effects on EPI was evaluated using an echo-planar spectroscopic imaging sequence.

RESULTS

A beat phenomenon resulting in a slowly varying phase was observed. Its temporal frequency was given by the difference between the EPI frequency and the mechanical resonance frequency of the activated gradient axis. The maximum erroneous, oscillating phase during phase encoding was ±0.5 rad for an EPI frequency of 1281 Hz. Echo-planar spectroscopic imaging images showed the resulting time-dependent stretching/compression of the FOV.

CONCLUSION

Oscillating readouts such as those used in EPI can result in low-frequency, erroneous phase contributions, which are explained by the beat phenomenon. Therefore, EPI phase-correction approaches may need to include beat effects for accurate image reconstruction.

Ramping down a clinical 3 T scanner: a journey into MRI and MRS at 0.75 T

OBJECT

Lower-field MR is reemerging as a viable, potentially cost-effective alternative to high-field MR, thanks to advances in hardware, sequence design, and reconstruction over the past decades. Evaluation of lower field strengths, however, is limited by the availability of lower-field systems on the market and their considerable procurement costs. In this work, we demonstrate a low-cost, temporary alternative to purchasing a dedicated lower-field MR system.

MATERIALS AND METHODS

By ramping down an existing clinical 3 T MRI system to 0.75 T, proton signals can be acquired using repurposed ¹³C transmit/receive hardware and the multi-nuclei spectrometer interface. We describe the ramp-down procedure and necessary software and hardware changes to the system.

RESULTS

Apart from presenting system characterization results, we show in vivo examples of cardiac cine imaging, abdominal two- and three-point Dixon-type water/fat separation, water/fat-separated MR Fingerprinting, and point-resolved spectroscopy. In addition, the ramp-down approach allows unique comparisons of, e.g., gradient fidelity of the same MR system operated at different field strengths using the same receive chain, gradient coils, and amplifiers.

DISCUSSION

Ramping down an existing MR system may be seen as a viable alternative for lower-field MR research in groups that already own multi-nuclei hardware and can also serve as a testing platform for custom-made multi-nuclei transmit/receive coils.

Fourier transform temporal diffusion spectroscopy

Temporal diffusion spectroscopy (TDS) currently uses the oscillating gradient spin echo (OGSE) experiment to measure the spectral density of translational velocity autocorrelation at single frequencies. Due to timing restrictions imposed by the transverse relaxation, the frequency selectivity and the sampling density of OGSE are limited, especially at low frequencies. We propose to overcome this problem by adopting the principles of Fourier transform spectroscopy. The new method of Fourier transform TDS (FTDS) uses two broadband gradient waveforms with different relative delays to make the spin echo attenuation sensitive to a broad range of diffusion frequencies with different harmonic modulations and calculates the spectrum by discrete Fourier transform. The method was validated by a measurement of diffusion spectra in highly restrictive tissues of a celery stalk and provided results consistent with OGSE, however, on a denser frequency grid.

Direct comparison of gradient Fidelity and acoustic noise of the same MRI system at 3 T and 0.75 T

PURPOSE

To analyze the difference between gradient fidelity and acoustic noise of the same MRI scanner operated at product field strength (3 T) and lower field strength (0.75 T).

METHODS

Gradient modulation transfer functions (GMTFs) were measured using a four-slice 2D phase-encoded chirp-based sequence on the same scanner operated at 3 T and, following ramp-down, at 0.75 T with identical gradient specifications (40 mT/m, 200 T/m/s). Calibrated audio measurements were performed at both field strengths to correlate audio spectra with GMTFs.

RESULTS

While eddy currents were independent of field strength, mechanical resonances were substantially decreased at lower field, resulting in a reduction of GMTF distortions by up to 95% (88% on average) at the mechanical resonances of the gradient system. Audio spectra amplitudes were reduced by up to 87% when comparing 0.75 T versus 3 T.

CONCLUSION

Lower static fields lead to reduced Lorentz forces on the gradient coil and, in turn, to reduced mechanical resonances, thereby improving gradient fidelity. Simultaneously, the reduction of acoustic noise may help to improve patient comfort.

Fundamentals of turbulent flow spectrum imaging

PURPOSE

To introduce a mathematical framework and in-silico validation of turbulent flow spectrum imaging (TFSI) of stenotic flow using phase-contrast MRI, evaluate systematic errors in quantitative turbulence parameter estimation, and propose a novel method for probing the Lagrangian velocity spectra of turbulent flows.

THEORY AND METHODS

The spectral response of velocity-encoding gradients is derived theoretically and linked to turbulence parameter estimation including the velocity autocorrelation function spectrum. Using a phase-contrast MRI simulation framework, the encoding properties of bipolar gradient waveforms with identical first gradient moments but different duration are investigated on turbulent flow data of defined characteristics as derived from computational fluid dynamics. Based on theoretical insights, an approach using velocity-compensated gradient waveforms is proposed to specifically probe desired ranges of the velocity autocorrelation function spectrum with increased accuracy.

RESULTS

Practical velocity-encoding gradients exhibit limited encoding power of typical turbulent flow spectra, resulting in up to 50% systematic underestimation of intravoxel SD values. Depending on the turbulence level in fluids, the error due to a single encoding gradient spectral response can vary by 20%. When using tailored velocity-compensated gradients, improved quantification of the Lagrangian velocity spectrum on a voxel-by-voxel basis is achieved and used for quantitative correction of intravoxel SD values estimated with velocity-encoding gradients.

CONCLUSION

To address systematic underestimation of turbulence parameters using bipolar velocity-encoding gradients in phase-contrast MRI of stenotic flows with short correlation times, tailored velocity-compensated gradients are proposed to improve quantitative mapping of turbulent blood flow characteristics.

Toward an accurate estimation of wall shear stress from 4D flow magnetic resonance downstream of a severe stenosis

PURPOSE

First, to investigate the agreement between velocity, velocity gradient, and Reynolds stress obtained from four-dimensional flow magnetic resonance (4D flow MRI) measurements and direct numerical simulation (DNS). Second, to propose and optimize based on DNS, 2 alternative methods for the accurate estimation of wall shear stress (WSS) when the resolution of the flow measurements is limited. Thirdly, to validate the 2 methods based on 4D flow MRI data.

METHODS

In vitro 4D MRI has been conducted in a realistic rigid stenosed aorta model under a constant flow rate of 12 L/min. A DNS of transitional stenotic flow has been performed using the same geometry and boundary conditions.

RESULTS

Time-averaged velocity and Reynolds stresses are in good agreement between in vitro 4D MRI data and DNS (errors between 2% and 8% of the reference downsampled data). WSS estimation based on the 2 proposed methods applied to MRI data provide good agreement with DNS for slice-averaged values (maximum error is less than 15% of the mean reference WSS for the first method and 25% for the second method). The performance of both models is not strongly sensitive to spatial resolution up to 1.5 mm voxel size. While the performance of model 1 deteriorates appreciably at low signal-to-noise ratios, model 2 remains robust.

CONCLUSIONS

The 2 methods for WSS magnitude give an overall better agreement than the standard approach used in the literature based on direct calculation of the velocity gradient close to the wall (relative error of 84%).

On the limitations of echo planar 4D flow MRI

PURPOSE

To compare EPI and GRE readout in high-flow velocity regimes and evaluate their impact on measurement accuracy in silico and in vitro.

THEORY AND METHODS

Phase-contrast sequences for EPI and GRE were simulated using CFD velocity data to assess displacement artifacts as well as effective spatial resolution. In silico findings were validated experimentally using a steady flow phantom.

RESULTS

For EPI factor 5 and simulated stenotic flow with peak velocity of 2.2 m s - 1 , displacement artifacts resulted in misregistration of 7.3 mm at echo time and the effective resolution was locally reduced by factors 5 and 8 compared to GRE for flow along phase and frequency encoding directions, respectively. In vitro, a maximum velocity difference between EPI factor 5 and GRE of 0.97 m s - 1 was found.

CONCLUSIONS

Four-dimensional flow MRI using EPI readout results not only in considerable velocity misregistration but also in spatially varying degradation of resolution. The proposed work indicates that EPI is inferior to standard GRE for 4D flow MRI.

Multipoint 5D flow cardiovascular magnetic resonance - accelerated cardiac- and respiratory-motion resolved mapping of mean and turbulent velocities

BACKGROUND

Volumetric quantification of mean and fluctuating velocity components of transient and turbulent flows promises a comprehensive characterization of valvular and aortic flow characteristics. Data acquisition using standard navigator-gated 4D Flow cardiovascular magnetic resonance (CMR) is time-consuming and actual scan times depend on the breathing pattern of the subject, limiting the applicability of the method in a clinical setting. We sought to develop a 5D Flow CMR framework which combines undersampled data acquisition including multipoint velocity encoding with low-rank image reconstruction to provide cardiac- and respiratory-motion resolved assessment of velocity maps and turbulent kinetic energy in fixed scan times.

METHODS

Data acquisition and data-driven motion state detection was performed using an undersampled Cartesian tiny Golden angle approach. Locally low-rank (LLR) reconstruction was implemented to exploit correlations among heart phases and respiratory motion states. To ensure accurate quantification of mean and turbulent velocities, a multipoint encoding scheme with two velocity encodings per direction was incorporated. Velocity-vector fields and turbulent kinetic energy (TKE) were obtained using a Bayesian approach maximizing the posterior probability given the measured data. The scan time of 5D Flow CMR was set to 4 min. 5D Flow CMR with acceleration factors of 19 .0 ± 0.21 (mean ± std) and velocity encodings (VENC) of 0.5 m/s and 1.5 m/s per axis was compared to navigator-gated 2x SENSE accelerated 4D Flow CMR with VENC = 1.5 m/s in 9 subjects. Peak velocities and peak flow were compared and magnitude images, velocity and TKE maps were assessed.

RESULTS

While net scan time of 5D Flow CMR was 4 min independent of individual breathing patterns, the scan times of the standard 4D Flow CMR protocol varied depending on the actual navigator gating efficiency and were 17.8 ± 3.9 min on average. Velocity vector fields derived from 5D Flow CMR in the end-expiratory state agreed well with data obtained from the navigated 4D protocol (normalized root-mean-square error 8.9 ± 2.1%). On average, peak velocities assessed with 5D Flow CMR were higher than for the 4D protocol (3.1 ± 4.4%).

CONCLUSIONS

Respiratory-motion resolved multipoint 5D Flow CMR allows mapping of mean and turbulent velocities in the aorta in 4 min.