Predicting long-term temperature increase for time-dependent SAR levels with a single short-term temperature response

Steady state temperature rise in multilayered tissue due to arbitrary periodic SAR using finite difference FFT and transfer function method

Steady state (SS) and transient temperature-rise in tissue from radiofrequency exposure forms the underlying basis for limits in international exposure guidelines. Periodically pulsed or intermittent exposures form a special case of having both peak and average levels, producing temperature-rise oscillations in the SS. Presented here is a method for determining tissue temperature-rise for periodic specific absorption rate (SAR) modulation having arbitrary waveform. It involves the finite difference solution of a form of the Pennes Bioheat Transfer equation (BHTE) and uses the concept of the transfer function and the Fast Fourier Transform (FFT). The time-dependent BHTE is converted to a SS harmonic version by assuming that the time-dependent SAR waveform and tissue temperature can both be represented by Fourier series. The transfer function is obtained from solutions of the harmonic BHTE for an assumed SAR waveform consisting of periodic impulses. The temperature versus time response for an arbitrary periodic SAR waveform is obtained from the inverse FFT of the product of the transfer function and the FFT of the actual SAR waveform. This method takes advantage of existing FFT algorithms on most computational platforms and the ability to store the transfer function for later re-use. The transfer function varies slowly with harmonic number, allowing interpolation and extrapolation to reduce the computational effort. The method is highly efficient for the case where repeated temperature-rise calculations for parameter variations in the SAR waveform are sought. Examples are given for a narrow, circularly symmetric beam incident on a planar skin/fat/muscle model with rectangular, triangular and cosine-pulsed SAR modulation waveforms. Calculations of temperature-rise crest factor as a function of rectangular pulse duty factor and pulse repetition frequency for the same exposure/tissue model are also presented as an example of the versatility of the method.

Transceive phase corrected 2D contrast source inversion-electrical properties tomography

PURPOSE

To remove the necessity of the tranceive phase assumption for CSI-EPT and show electrical properties maps reconstructed from measured data obtained using a standard 3T birdcage body coil setup.

METHODS

The existing CSI-EPT algorithm is reformulated to use the transceive phase rather than relying on the transceive phase assumption. Furthermore, the radio frequency (RF)-shield is numerically implemented to accurately model the RF fields inside the MRI scanner. We verify that the reformulated two-dimensional (2D) CSI-EPT algorithm can reconstruct electrical properties maps given 2D electromagnetic simulations. Afterward, the algorithm is tested with three-dimensional (3D) FDTD simulations to investigate if the 2D CSI-EPT can retrieve the electrical properties for 3D RF fields. Finally, an MR experiment at 3T with a phantom is performed.

RESULTS

From the results of the 2D simulations, it is seen that CSI-EPT can reconstruct the electrical properties using MRI accessible quantities. For 3D simulations, it is observed that the electrical properties are underestimated, nonetheless, CSI-EPT has a lower standard deviation than the standard Helmholtz based methods. Finally, the first CSI-EPT reconstructions based on measured data are presented showing comparable accuracy and precision to reconstructions based on simulated data, and demonstrating the feasibility of CSI-EPT.

CONCLUSIONS

The CSI-EPT algorithm was rewritten to use MRI accessible quantities. This allows for CSI-EPT to fully exploit the benefits of the higher static magnetic field strengths with a standard quadrature birdcage coil setup.

Temperature-based MRI safety simulations with a limited number of tissues

PURPOSE

Demonstrate ability to produce reasonable simulations of temperature using numerical models of the human body with a limited number of tissues.

METHODS

For both a male and female human body model, numerical simulations were used to calculate temperature distributions in three different models of the same human body: the original model with 35 tissues for the male model and 76 tissues for the female model, a simplified model having only three tissues (muscle, fat, and lung), and a simplified model having six tissues (muscle, fat, lung, bone, brain, and skin).

RESULTS

Although a three-tissue model gave reasonable specific absorption rate estimates in comparison to an original with many more tissues, because of tissue-specific thermal and physiological properties that do not affect specific absorption rate, such as rate of perfusion by blood, the three-tissue model did not provide temperature distributions similar to those of the original model. Inclusion of a few additional tissues, as in the six-tissue model, produced results in much better agreement with those from the original model.

CONCLUSION

Reasonable estimates of temperature can be simulated with a limited number of tissues, although this number is higher than the number of tissues required to produce reasonable simulations of specific absorption rate. For exposures primarily in the head and thorax, six tissues may be adequate for reasonable estimates of temperature.

7T MR Safety

Magnetic resonance imaging and spectroscopy (MRI/MRS) at 7T represents an exciting advance in MR technology, with intriguing possibilities to enhance image spatial, spectral, and contrast resolution. To ensure the safe use of this technology while still harnessing its potential, clinical staff and researchers need to be cognizant of some safety concerns arising from the increased magnetic field strength and higher Larmor frequency. The higher static magnetic fields give rise to enhanced transient bioeffects and an increased risk of adverse incidents related to electrically conductive implants. Many technical challenges remain and the continuing rapid pace of development of 7T MRI/MRS is likely to present further challenges to ensuring safety of this technology in the years ahead. The recent regulatory clearance for clinical diagnostic imaging at 7T will likely increase the installed base of 7T systems, particularly in hospital environments with little prior ultrahigh-field MR experience. Informed risk/benefit analyses will be required, particularly where implant manufacturer-published 7T safety guidelines for implants are unavailable. On behalf of the International Society for Magnetic Resonance in Medicine, the aim of this article is to provide a reference document to assist institutions developing local institutional policies and procedures that are specific to the safe operation of 7T MRI/MRS. Details of current 7T technology and the physics underpinning its functionality are reviewed, with the aim of supporting efforts to expand the use of 7T MRI/MRS in both research and clinical environments. Current gaps in knowledge are also identified, where additional research and development are required. Level of Evidence 5 Technical Efficacy 2 J. MAGN. RESON. IMAGING 2021;53:333-346.

Optimization of the order and spacing of sequences in an MRI exam to reduce the maximum temperature and thermal dose

PURPOSE

Evaluate the possibility to reduce specific energy absorption rate (SAR)-induced maximum temperature and thermal dose by rearranging the order and spacing of sequences without increasing duration of the MRI examination.

METHODS

Using numerical simulations based on an actual SAR-intensive MRI examination, optimizations to reduce either maximum temperature or thermal dose were performed. For each permutation of groups of sequences having the same patient table position, temperature and thermal dose were computed very rapidly using recently published methods. Disposition of sequences was further adjusted by optimizing the spacing between each sequence without exceeding the original exam duration.

RESULTS

The maximum simulated temperature in the original exam was 42.38°C, and the maximum thermal dose was 3.23 cumulative effective minutes at 43°C (CEM43). After optimization to reduce maximum temperature, it was 41.77°C, and after optimization to minimize the thermal dose, it was 1.42 CEM43.

CONCLUSION

It is possible to reduce maximum temperature and thermal dose in the exam by changing the arrangement and spacing of the sequences without increasing the duration of the exam (by increasing TR or adding delays) or compromising image quality (by reducing flip angles).

Magnetic resonance thermometry of flowing blood

Blood temperature is a key determinant of tissue temperature and can be altered under normal physiological states, such as exercise, in diseases such as stroke or iatrogenically in therapies which modulate tissue temperature, such as therapeutic hypothermia. Currently available methods for the measurement of arterial and venous temperatures are invasive and, for small animal models, are impractical. Here, we present a methodology for the measurement of intravascular and tissue temperature by magnetic resonance imaging (MRI) using the lanthanide agent TmDOTMA^- (DOTMA, tetramethyl-1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid; Tm, thulium). The approach makes use of phase-sensitive imaging measurements, combined with spectrally selective excitation, to monitor the temperature-dependent shift in the resonance of proton nuclei associated with water and with methyl groups of TmDOTMA^- . Measurements were first made in a flow phantom modelling diastolic blood flow in the mouse aorta or inferior vena cava (IVC) and imaged using 7-T preclinical MRI with a custom-built surface coil. Flowing and static fluid temperatures agreed to within 0.12°C for these experiments. Proof-of-concept experiments were also performed on three healthy adult mice, demonstrating temperature measurements in the aorta, IVC and kidney following a bolus injection of contrast agent. A small (0.7-1°C), but statistically significant, higher kidney temperature compared with the aorta (p = 0.002-0.007) and IVC (p = 0.003-0.03) was shown in all animals. These findings demonstrate the feasibility of the technique for in vivo applications and illustrate how the technique could be used to explore the relationship between blood and tissue temperature for a wide range of applications.

Parallel transmission RF pulse design with strict temperature constraints

RF safety in parallel transmission (pTx) is generally ensured by imposing specific absorption rate (SAR) limits during pTx RF pulse design. There is increasing interest in using temperature to ensure safety in MRI. In this work, we present a local temperature correlation matrix formalism and apply it to impose strict constraints on maximum absolute temperature in pTx RF pulse design for head and hip regions. Electromagnetic field simulations were performed on the head and hip of virtual body models. Temperature correlation matrices were calculated for four different exposure durations ranging between 6 and 24 min using simulated fields and body-specific constants. Parallel transmission RF pulses were designed using either SAR or temperature constraints, and compared with each other and unconstrained RF pulse design in terms of excitation fidelity and safety. The use of temperature correlation matrices resulted in better excitation fidelity compared with the use of SAR in parallel transmission RF pulse design (for the 6 min exposure period, 8.8% versus 21.0% for the head and 28.0% versus 32.2% for the hip region). As RF exposure duration increases (from 6 min to 24 min), the benefit of using temperature correlation matrices on RF pulse design diminishes. However, the safety of the subject is always guaranteed (the maximum temperature was equal to 39°C). This trend was observed in both head and hip regions, where the perfusion rates are very different.