Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5–1 GHz)

A theoretical and experimental investigation on a volume coil with slotted end-rings for rat MRI at 7 T

OBJECTIVE

A volume coil with squared slots-end ring was developed to attain improved sensitivity for imaging of rat's brain at 7 T.

MATERIAL AND METHODS

The principles of the high cavity resonator for the low-pass case and the law of Biot-Savart were used to derive a theoretical expression of [Formula: see text]. The slotted-end ring resonator showed a theoretical 2.22-fold sensitivity improvement over the standard birdcage coil with similar dimensions. Numerical studies were carried out for the electromagnetic fields and specific absorption rates for our coil and a birdcage coil loaded with a saline-filled spherical phantom and a digital brain of a rat.

RESULTS

An improvement of the signal-to-noise ratio (SNR) can be observed for the slotted volume coil over the birdcage regardless of the load used in the electromagnetic simulations. The specific absorption rate simulations show a decrement for the digital brain and quite similar values with the saline solution phantom. Phantom and rat's brain images were acquired at 7 T to prove the viability of the coil design. The experimental noise figure of our coil design was four times less than the standard birdcage with similar dimensions, which showed a 44.5% increase in experimental SNR.

DISCUSSION

There is remarkable agreement among the theoretical, numerical and experimental sensitivity values, which all demonstrate that the coil performance for MR imaging of small rodents can be improved using slotted end-rings.

Multiphysics and Thermal Response Models to Improve Accuracy of Local Temperature Estimation in Rat Cortex under Microwave Exposure

The rapid development of wireless technology has led to widespread concerns regarding adverse human health effects caused by exposure to electromagnetic fields. Temperature elevation in biological bodies is an important factor that can adversely affect health. A thermophysiological model is desired to quantify microwave (MW) induced temperature elevations. In this study, parameters related to thermophysiological responses for MW exposures were estimated using an electromagnetic-thermodynamics simulation technique. To the authors’ knowledge, this is the first study in which parameters related to regional cerebral blood flow in a rat model were extracted at a high degree of accuracy through experimental measurements for localized MW exposure at frequencies exceeding 6 GHz. The findings indicate that the improved modeling parameters yield computed results that match well with the measured quantities during and after exposure in rats. It is expected that the computational model will be helpful in estimating the temperature elevation in the rat brain at multiple observation points (that are difficult to measure simultaneously) and in explaining the physiological changes in the local cortex region.

The specific absorption rate of tissues in rats exposed to electromagnetic plane waves in the frequency range of 0.05-5 GHz and SARwb in free-moving rats

As electromagnetic exposure experiments can only be performed on small animals, usually rats, research on the characteristics of specific absorption rate (SAR) distribution in the rat has received increasing interest. A series of calculations, which simulated the SAR in a male rat anatomical model exposed to electromagnetic plane waves ranging from 0.05 to 5 GHz with different incidence and polarization, were conducted. The whole-body-averaged SAR (SARwb) and the tissue-averaged SAR (SARavg) in 20 major tissues were determined. Results revealed that incidence has great impact on SAR in the rat at higher frequencies owing to the skin effect and the effect on SARavg in tissues is much more apparent than that on SARwb; while polarization plays an important role under lower frequencies. Not only the incidence, but also the polarization in the rat keeps changing when the rat is in free movement. Thus, this article discussed a convenient way to obtain relatively accurate SARwb in a free-moving rat.

Development of computational small animal models and their applications in preclinical imaging and therapy research

The development of multimodality preclinical imaging techniques and the rapid growth of realistic computer simulation tools have promoted the construction and application of computational laboratory animal models in preclinical research. Since the early 1990s, over 120 realistic computational animal models have been reported in the literature and used as surrogates to characterize the anatomy of actual animals for the simulation of preclinical studies involving the use of bioluminescence tomography, fluorescence molecular tomography, positron emission tomography, single-photon emission computed tomography, microcomputed tomography, magnetic resonance imaging, and optical imaging. Other applications include electromagnetic field simulation, ionizing and nonionizing radiation dosimetry, and the development and evaluation of new methodologies for multimodality image coregistration, segmentation, and reconstruction of small animal images. This paper provides a comprehensive review of the history and fundamental technologies used for the development of computational small animal models with a particular focus on their application in preclinical imaging as well as nonionizing and ionizing radiation dosimetry calculations. An overview of the overall process involved in the design of these models, including the fundamental elements used for the construction of different types of computational models, the identification of original anatomical data, the simulation tools used for solving various computational problems, and the applications of computational animal models in preclinical research. The authors also analyze the characteristics of categories of computational models (stylized, voxel-based, and boundary representation) and discuss the technical challenges faced at the present time as well as research needs in the future.

Monitoring local heating around an interventional MRI antenna with RF radiometry

PURPOSE

Radiofrequency (RF) radiometry uses thermal noise detected by an antenna to measure the temperature of objects independent of medical imaging technologies such as magnetic resonance imaging (MRI). Here, an active interventional MRI antenna can be deployed as a RF radiometer to measure local heating, as a possible new method of monitoring device safety and thermal therapy.

METHODS

A 128 MHz radiometer receiver was fabricated to measure the RF noise voltage from an interventional 3 T MRI loopless antenna and calibrated for temperature in a uniformly heated bioanalogous gel phantom. Local heating (ΔT) was induced using the antenna for RF transmission and measured by RF radiometry, fiber-optic thermal sensors, and MRI thermometry. The spatial thermal sensitivity of the antenna radiometer was numerically computed using a method-of-moment electric field analyses. The gel's thermal conductivity was measured by MRI thermometry, and the localized time-dependent ΔT distribution computed from the bioheat transfer equation and compared with radiometry measurements. A "H-factor" relating the 1 g-averaged ΔT to the radiometric temperature was introduced to estimate peak temperature rise in the antenna's sensitive region.

RESULTS

The loopless antenna radiometer linearly tracked temperature inside a thermally equilibrated phantom up to 73 °C to within ±0.3 °C at a 2 Hz sample rate. Computed and MRI thermometric measures of peak ΔT agreed within 13%. The peak 1 g-average temperature was H = 1.36 ± 0.02 times higher than the radiometric temperature for any media with a thermal conductivity of 0.15-0.50 (W/m)/K, indicating that the radiometer can measure peak 1 g-averaged ΔT in physiologically relevant tissue within ±0.4 °C.

CONCLUSIONS

Active internal MRI detectors can serve as RF radiometers at the MRI frequency to provide accurate independent measures of local and peak temperature without the artifacts that can accompany MRI thermometry or the extra space needed to accommodate alternative thermal transducers. A RF radiometer could be integrated in a MRI scanner to permit "self-monitoring" for assuring device safety and/or monitoring delivery of thermal therapy.

Surface coil with reduced specific absorption rate for rat MRI at 7 T

OBJECTIVE

A scaled-down slotted surface radio frequency (RF) coil was built, and the specific absorbance rate (SAR) in 100 mg of tissue (SAR100 mg) produced in a rat brain phantom was computed at 7 T.

MATERIALS AND METHODS

A slotted coil 2-cm in diameter with six circular slots was developed. Its theoretical and experimental performance was computed and compared using the signal-to-noise ratio (SNR) expression and phantom images obtained with a spin echo sequence. Electromagnetic simulations were performed using the finite integral method with saline sphere and rat brain phantoms. SAR100 mg was computed for the circular coil, by varying its radius, and was also computed for the slotted coil.

RESULTS

The slotted coil quality factor gave a twofold increment over the circular coil, and noise was reduced by 17%. The experimental SNR of the slotted coil produced a 30% improvement for points near the coil plane. The theoretical and experimental results showed substantial agreement. Axial map histograms and profiles showed greater SAR100 mg values for the circular coil than for the slotted coil.

CONCLUSIONS

The slotted surface coil offers improved performance and low SAR100 mg for rat brain imaging at 7 T. This approach may be used with new RF coils to investigate SAR in humans.

Consideration of the effects of intense tissue heating on the RF electromagnetic fields during MRI: simulations for MRgFUS in the hip

Due to the strong dependence of tissue electrical properties on temperature, it is important to consider the potential effects of intense tissue heating on the RF electromagnetic fields during MRI, as can occur in MR-guided focused ultrasound surgery. In principle, changes of the RF electromagnetic fields could affect both efficacy of RF pulses, and the MRI-induced RF heating (SAR) pattern. In this study, the equilibrium temperature distribution in a whole-body model with 2 mm resolution before and during intense tissue heating up to 60 °C at the target region was calculated. Temperature-dependent electric properties of tissues were assigned to the model to establish a temperature-dependent electromagnetic whole-body model in a 3T MRI system. The results showed maximum changes in conductivity, permittivity, [absolute value]B(1)(+)[absolute value] and SAR of about 25%, 6%, 2%, and 20%, respectively. Though the B1 field and SAR distributions are both temperature-dependent, the potential harm to patients due to higher SARs is expected to be minimal and the effects on the B1 field distribution should have minimal effect on images from basic MRI sequences.

A virtual rat for simulating environmental and exertional heat stress

Severe cases of environmental or exertional heat stress can lead to varying degrees of organ dysfunction. To understand heat-injury progression and develop efficient management and mitigation strategies, it is critical to determine the thermal response in susceptible organs under different heat-stress conditions. To this end, we used our previously published virtual rat, which is capable of computing the spatiotemporal temperature distribution in the animal, and extended it to simulate various heat-stress scenarios, including 1) different environmental conditions, 2) exertional heat stress, 3) circadian rhythm effect on the thermal response, and 4) whole body cooling. Our predictions were consistent with published in vivo temperature measurements for all cases, validating our simulations. We observed a differential thermal response in the organs, with the liver experiencing the highest temperatures for all environmental and exertional heat-stress cases. For every 3°C rise in the external temperature from 40 to 46°C, core and organ temperatures increased by ∼0.8°C. Core temperatures increased by 2.6 and 4.1°C for increases in exercise intensity from rest to 75 and 100% of maximal O2 consumption, respectively. We also found differences as large as 0.8°C in organ temperatures for the same heat stress induced at different times during the day. Even after whole body cooling at a relatively low external temperature (1°C for 20 min), average organ temperatures were still elevated by 2.3 to 2.5°C compared with normothermia. These results can be used to optimize experimental protocol designs, reduce the amount of animal experimentation, and design and test improved heat-stress prevention and management strategies.

A 3-D mathematical model to identify organ-specific risks in rats during thermal stress

Early prediction of the adverse outcomes associated with heat stress is critical for effective management and mitigation of injury, which may sometimes lead to extreme undesirable clinical conditions, such as multiorgan dysfunction syndrome and death. Here, we developed a computational model to predict the spatiotemporal temperature distribution in a rat exposed to heat stress in an attempt to understand the correlation between heat load and differential organ dysfunction. The model includes a three-dimensional representation of the rat anatomy obtained from medical imaging and incorporates the key mechanisms of heat transfer during thermoregulation. We formulated a novel approach to estimate blood temperature by accounting for blood mixing from the different organs and to estimate the effects of the circadian rhythm in body temperature by considering day-night variations in metabolic heat generation and blood perfusion. We validated the model using in vivo core temperature measurements in control and heat-stressed rats and other published experimental data. The model predictions were within 1 SD of the measured data. The liver demonstrated the greatest susceptibility to heat stress, with the maximum temperature reaching 2°C higher than the measured core temperature and 95% of its volume exceeding the targeted experimental core temperature. Other organs also attained temperatures greater than the core temperature, illustrating the need to monitor multiple organs during heat stress. The model facilitates the identification of organ-specific risks during heat stress and has the potential to aid in the development of improved clinical strategies for thermal-injury prevention and management.

Potassium Titanyl Phosphate Laser Tissue Ablation: Development and Experimental Validation of a New Numerical Model

The photoselective vaporization of prostate (PVP) green light (532 nm) laser is increasingly being used as an alternative to the transurethral resection of prostate (TURP) for treatment of benign prostatic hyperplasia (BPH) in older patients and those who are poor surgical candidates. In order to achieve the goals of increased tissue removal volume (i.e., “ablation” in the engineering sense) and reduced collateral thermal damage during the PVP green light treatment, a two dimensional computational model for laser tissue ablation based on available parameters in the literature has been developed and compared to experiments. The model is based on the control volume finite difference and the enthalpy method with a mechanistically defined energy necessary to ablate (i.e., physically remove) a volume of tissue (i.e., energy of ablation Eab). The model was able to capture the general trends experimentally observed in terms of ablation and coagulation areas, their ratio (therapeutic index (TI)), and the ablation rate (AR) (mm3/s). The model and experiment were in good agreement at a smaller working distance (WD) (distance from the tissue in mm) and a larger scanning speed (SS) (laser scan speed in mm/s). However, the model and experiment deviated somewhat with a larger WD and a smaller SS; this is most likely due to optical shielding and heat diffusion in the laser scanning direction, which are neglected in the model. This model is a useful first step in the mechanistic prediction of PVP based BPH laser tissue ablation. Future modeling efforts should focus on optical shielding, heat diffusion in the laser scanning direction (i.e., including 3D effects), convective heat losses at the tissue boundary, and the dynamic optical, thermal, and coagulation properties of BPH tissue.

Exposure of workers to pulsed gradients in MRI

PURPOSE

To numerically evaluate the electric field/current density magnitudes and spatial distributions in healthcare workers when they are standing close to the gradient coil windings near the magnetic resonance imaging (MRI) scanner ends.

MATERIALS AND METHODS

Anatomically realistic, whole-body male and female voxel phantoms are engaged to model the workers at various positions near the ends of three cylindrical gradient coils (x-, y-, and z-axis gradients). The numerical calculations of induced fields are based on an efficient, quasistatic finite-difference method.

RESULTS

The simulations show that it is possible to induce electric fields/current densities above levels recommended by the International Commission for Non-ionizing Radiation Protection (ICNIRP) and the Institute of Electrical and Electronics Engineers (IEEE) standards when the workers are standing close to the gradient coils and when two or three gradients are switched simultaneously, as is often the case.

CONCLUSION

The longitudinal gradient tends to induce more fields in workers than the transverse coils. The strongest levels of field exposure are observed when all three gradients are operated simultaneously and can be above regulations when the healthcare worker is close to the gradient coils. Other postures such as bending into the magnet shall be investigated in further studies.

Radiofrequency heating at 9.4T: in vivo temperature measurement results in swine

In vivo temperatures were correlated to the whole head average specific absorption rate (SAR(avg)) at 9.4T using 12 anesthetized swine (mean animal weight = 52 kg, standard deviation = 6.7 kg). Correlating the temperatures and SAR(avg) is necessary to ensure safe levels of human heating during ultra-high field MR exams. The temperatures were measured at three depths inside the brain, in the rectum, and at the head-skin of swine. A 400 MHz, continuous wave RF power was deposited to the head using a volume coil. The SAR(avg) values were varied between 2.7-5.8 W/kg. The RF power exposure durations were varied between 1.4-3.7 hr. To differentiate the temperature response caused by the RF from that of the anesthesia, the temperatures were recorded in four unheated swine. To study the effect of the spatial distribution of the RF and tissue properties, the temperature probes were placed at two brain locations (n = 4 swine for each location). Results showed that the in vivo brain temperatures correlated to the SAR(avg) in a geometry-dependent manner. Additionally, 1) the skin temperature change was not the maximum temperature change; 2) the RF heating caused an inhomogeneous brain temperature distribution; and 3) the maximum temperature occurred inside the brain.

Electric field measurements and computational modeling at ultrahigh-field MRI

While magnetic resonance images essentially contain a map of the both circularly polarized components of the RF transverse magnetic fields (B(1) field), the thermal heat and electromagnetic power deposition is generated by the associated electric fields. Measurement of electric field distributions/intensities across a sample yields an indirect indication of possible cause of heating within the sample and potentially enables the detection of "hot spots," which can be present within inhomogeneous radiofrequency (RF) fields, such as the case with magnetic resonance imaging at high field strength. As a result, establishing a valid technique for direct measurements of the electric field and its correlation, obtained using computational electromagnetics, is essential in assessing (1) the safety of the RF coil designs and (2) the validity of the calculations. In this work, a probe was built and used to measure the transverse electric field (E(1) field) distributions within an empty 8 T (tuned to 340 MHz) RF head coil and within a saline water phantom loaded in the same coil. The measured E(1) field distributions were favorably compared to the distributions obtained utilizing a finite difference time domain in-house package.

Transient temperature rise in a mouse due to low-frequency regional hyperthermia

A refined nonlinear heat transfer model of a mouse has been developed to simulate the transient temperature rise in a neoplastic tumour and neighbouring tissue during regional hyperthermia using a 150 kHz inductive coil. In this study, we incorporate various bio-energetic enhancements to the heat transfer equation and numerical validations based on experimental findings for the mouse, in terms of nonlinear metabolic heat production, homeothermy, blood perfusion parameters, thermoregulation, psychological and physiological effects. The discretized bio-heat transfer equation has been validated with the commercial software FEMLAB on a canonical multi-sphere object before applying the scheme to the inhomogeneous mouse voxel phantom. The time-dependent numerical results of regional hyperthermia of mouse thigh have been compared with the available experimental temperature results with only a few small disparities. During the first 20 min of local unfocused heating, the temperature in the tumour and the surrounding tissue increased by around 7.5 degrees C. The objective of this preliminary study was to develop a validated electrothermal numerical scheme for inductive hyperthermia of a small mammal with the intention of expanding the model into a complete numerical solution involving ferromagnetic nanoparticles for targeted heating of tumours at low frequencies. In addition, the numerical scheme herein could assist in optimizing and tailoring of focused electromagnetic fields for hyperthermia.

The welfare and scientific advantages of non-invasive imaging of animals used in biomedical research

At present, animal experimentation remains central to our understanding of human disease-related processes and of the biological effects of many substances. Traditional experiments have relied heavily on invasive techniques to monitor changes in blood biochemistry, tissue structure or function, or to phenotype or genotype genetically modified animals. In some cases, a proportion or all of the animals used during the course of a study may be sacrificed for histopathological assessment. In most cases, this is to track the progression or regression of a disease over time, or to determine the levels of toxicity evident in specific organs or tissues. However, many of these techniques fail to provide details of how a disease develops or how a substance elicits its effects. In recent years there has been a gradual increase in the application of imaging techniques that were originally developed and used in fundamental research or in medicine. These non-invasive techniques allow diseases, and responses to exogenous substances, to be monitored in a temporal and spatial manner, therefore allowing a greater amount of information to be derived from smaller numbers of animals, which in turn, increases the statistical validity of the data by reducing the level of experimental variation. Non-invasive imaging also allows more informative and humane endpoints to be used and, perhaps most importantly, allows functional details to be studied in the context of a living animal. Some of the recent developments within the field of non-invasive imaging and their significance with respect to animal welfare and the understanding of human physiology are discussed.

A theoretical comparison of two optimization methods for radiofrequency drive schemes in high frequency MRI resonators

In this paper, numerical simulations are used in an attempt to find optimal source profiles for high frequency radiofrequency (RF) volume coils. Biologically loaded, shielded/unshielded circular and elliptical birdcage coils operating at 170 MHz, 300 MHz and 470 MHz are modelled using the FDTD method for both 2D and 3D cases. Taking advantage of the fact that some aspects of the electromagnetic system are linear, two approaches have been proposed for the determination of the drives for individual elements in the RF resonator. The first method is an iterative optimization technique with a kernel for the evaluation of RF fields inside an imaging plane of a human head model using pre-characterized sensitivity profiles of the individual rungs of a resonator; the second method is a regularization-based technique. In the second approach, a sensitivity matrix is explicitly constructed and a regularization procedure is employed to solve the ill-posed problem. Test simulations show that both methods can improve the B(1)-field homogeneity in both focused and non-focused scenarios. While the regularization-based method is more efficient, the first optimization method is more flexible as it can take into account other issues such as controlling SAR or reshaping the resonator structures. It is hoped that these schemes and their extensions will be useful for the determination of multi-element RF drives in a variety of applications.