A generic bioheat transfer thermal model for a perfused tissue

Towards an integrated radiofrequency safety concept for implant carriers in MRI based on sensor-equipped implants and parallel transmission

To protect implant carriers in MRI from excessive radiofrequency (RF) heating it has previously been suggested to assess that hazard via sensors on the implant. Other work recommended parallel transmission (pTx) to actively mitigate implant-related heating. Here, both ideas are integrated into one comprehensive safety concept where native pTx safety (without implant) is ensured by state-of-the-art field simulations and the implant-specific hazard is quantified in situ using physical sensors. The concept is demonstrated by electromagnetic simulations performed on a human voxel model with a simplified spinal-cord implant in an eight-channel pTx body coil at 3 T . To integrate implant and native safety, the sensor signal must be calibrated in terms of an established safety metric (e.g., specific absorption rate [SAR]). Virtual experiments show that E -field and implant-current sensors are well suited for this purpose, while temperature sensors require some caution, and B 1 probes are inadequate. Based on an implant sensor matrix Q s , constructed in situ from sensor readings, and precomputed native SAR limits, a vector space of safe RF excitations is determined where both global (native) and local (implant-related) safety requirements are satisfied. Within this safe-excitation subspace, the solution with the best image quality in terms of B 1 + magnitude and homogeneity is then found by a straightforward optimization algorithm. In the investigated example, the optimized pTx shim provides a 3-fold higher mean B 1 + magnitude compared with circularly polarized excitation for a maximum implant-related temperature increase ∆ T imp ≤ 1 K . To date, sensor-equipped implants interfaced to a pTx scanner exist as demonstrator items in research labs, but commercial devices are not yet within sight. This paper aims to demonstrate the significant benefits of such an approach and how this could impact implant-related RF safety in MRI. Today, the responsibility for safe implant scanning lies with the implant manufacturer and the MRI operator; within the sensor concept, the MRI manufacturer would assume much of the operator's current responsibility.

Predicting brain temperature in humans using bioheat models: Progress and outlook

Brain temperature, regulated by the balance between blood circulation and metabolic heat generation, is an important parameter related to neural activity, cerebral hemodynamics, and neuroinflammation. A key challenge for integrating brain temperature into clinical practice is the lack of reliable and non-invasive brain thermometry. The recognized importance of brain temperature and thermoregulation in both health and disease, combined with limited availability of experimental methods, has motivated the development of computational thermal models using bioheat equations to predict brain temperature. In this mini-review, we describe progress and the current state-of-the-art in brain thermal modeling in humans and discuss potential avenues for clinical applications.

Computational modeling of microwave ablation with thermal accelerants

PURPOSE

To develop a computational model of microwave ablation (MWA) with a thermal accelerant gel and apply the model toward interpreting experimental observations in ex vivo bovine and in vivo porcine liver.

METHODS

A 3D coupled electromagnetic-heat transfer model was implemented to characterize thermal profiles within ex vivo bovine and in vivo porcine liver tissue during MWA with the HeatSYNC thermal accelerant. Measured temperature dependent dielectric and thermal properties of the HeatSYNC gel were applied within the model. Simulated extents of MWA zones and transient temperature profiles were compared against experimental measurements in ex vivo bovine liver. Model predictions of thermal profiles under in vivo conditions in porcine liver were used to analyze thermal ablations observed in prior experiments in porcine liver in vivo.

RESULTS

Measured electrical conductivity of the HeatSYNC gel was ∼83% higher compared to liver at room temperature, with positive linear temperature dependency, indicating increased microwave absorption within HeatSYNC gel compared to tissue. In ex vivo bovine liver, model predicted ablation zone extents of (31.5 × 36) mm with the HeatSYNC, compared to (32.9 ± 2.6 × 40.2 ± 2.3) mm in experiments (volume differences 4 ± 4.1 cm3). Computational models under in vivo conditions in porcine liver suggest approximating the HeatSYNC gel spreading within liver tissue during ablations as a plausible explanation for larger ablation zones observed in prior in vivo studies.

CONCLUSION

Computational models of MWA with thermal accelerants provide insight into the impact of accelerant on MWA, and with further development, could predict ablations with a variety of gel injection sites.

Comparisons of healthy human brain temperature predicted from biophysical modeling and measured with whole brain MR thermometry

Brain temperature is an understudied parameter relevant to brain injury and ischemia. To advance our understanding of thermal dynamics in the human brain, combined with the challenges of routine experimental measurements, a biophysical modeling framework was developed to facilitate individualized brain temperature predictions. Model-predicted brain temperatures using our fully conserved model were compared with whole brain chemical shift thermometry acquired in 30 healthy human subjects (15 male and 15 female, age range 18-36 years old). Magnetic resonance (MR) thermometry, as well as structural imaging, angiography, and venography, were acquired prospectively on a Siemens Prisma whole body 3 T MR scanner. Bland-Altman plots demonstrate agreement between model-predicted and MR-measured brain temperatures at the voxel-level. Regional variations were similar between predicted and measured temperatures (< 0.55 °C for all 10 cortical and 12 subcortical regions of interest), and subcortical white matter temperatures were higher than cortical regions. We anticipate the advancement of brain temperature as a marker of health and injury will be facilitated by a well-validated computational model which can enable predictions when experiments are not feasible.

Tuning the Pennes Perfusion Rate to Model Large Vessel Cooling Effects in Hepatic Radiofrequency Ablation

Radio-frequency ablation (RFA) has become a popular method for the minimally invasive treatment of liver cancer. However, the success rate of these treatments depends heavily on the amount of experience the clinician possesses. Mathematical modelling can help mitigate this problem by providing an indication of the treatment outcome. Thermal lesions in RFA are affected by the cooling effect of both fine-scale and large-scale blood vessels. The exact model for large-scale blood vessels is advection-diffusion, i.e. a model capable of producing directional effects, which are known to occur in certain cases. In previous research, in situations where directional effects do not occur, the advection term in the blood vessel model has been typically replaced with the Pennes perfusion term, albeit with a higher-than usual perfusion rate. Whether these values of the perfusion rate appearing in literature are optimal for the particular vessel radii in question, has not been investigated so far. The present work aims to address this issue. An attempt has been made to determine, for values of vessel radius between 0.55 mm and 5 mm, best estimates for the perfusion rate which minimize the error in thermal lesion volumes between the perfusion-based model and the advection-based model. The results for the best estimate of the perfusion rate presented may be used in existing methods for fast estimation of RFA outcomes. Furthermore, the possible improvements to the presented methodology have been highlighted.

Measuring radiofrequency field-induced temperature variations in brain MRI exams with motion compensated MR thermometry and field monitoring

PURPOSE

An MR thermometry (MRT) method with motion and field fluctuation compensation is proposed to measure non-invasively sub-degree brain temperature variations occurring through radiofrequency (RF) power deposition during MR exams.

METHODS

MRT at 7T with a multi-slice echo planar imaging (EPI) sequence and concurrent field monitoring was first tested in vitro to assess accuracy in the presence of external field perturbations, an optical probe being used for ground truth. In vivo, this strategy was complemented by a motion compensation scheme based on a dictionary pre-scan, as reported in some previous work, and was adapted to the human brain. Precision reached with this scheme was assessed on eight volunteers with a 5 minute-long low specific absorption rate (SAR) scan. Finally, temperature rise in the brain was measured twice on the same volunteers and with the same strategy, this time by employing a 20-minutes scan at the maximum SAR delivered with a commercial volume head coil.

RESULTS

In vitro, the root mean square (RMS) error between optical probe and MRT measurements was 0.02°C with field sensor correction. In vivo, the low SAR scan returned a precision in temperature change measurement with field monitoring and motion compensation of 0.05°C. The 20-minutes maximum SAR scan returned a temperature rise throughout the inner-brain in the range of 0-0.2°C. Brain periphery remained too sensitive with respect to motion to lead to equally conclusive results.

CONCLUSION

Sub-degree temperature rise in the inner human brain was characterized experimentally throughout RF exposure. Potential applications include improvement of human thermal models and revision of safety margins.

ESHO benchmarks for computational modeling and optimization in hyperthermia therapy

BACKGROUND

The success of cancer hyperthermia (HT) treatments is strongly dependent on the temperatures achieved in the tumor and healthy tissues as it correlates with treatment efficacy and safety, respectively. Hyperthermia treatment planning (HTP) simulations have become pivotal for treatment optimization due to the possibility for pretreatment planning, optimization and decision making, as well as real-time treatment guidance.

MATERIALS AND METHODS

The same computational methods deployed in HTP are also used for in silico studies. These are of great relevance for the development of new HT devices and treatment approaches. To aid this work, 3 D patient models have been recently developed and made available for the HT community. Unfortunately, there is no consensus regarding tissue properties, simulation settings, and benchmark applicators, which significantly influence the clinical relevance of computational outcomes.

RESULTS AND DISCUSSION

Herein, we propose a comprehensive set of applicator benchmarks, efficacy and safety optimization algorithms, simulation settings and clinical parameters, to establish benchmarks for method comparison and code verification, to provide guidance, and in view of the 2021 ESHO Grand Challenge (Details on the ESHO grand challenge on HTP will be provided at https://www.esho.info/).

CONCLUSION

We aim to establish guidelines to promote standardization within the hyperthermia community such that novel approaches can quickly prove their benefit as quickly as possible in clinically relevant simulation scenarios. This paper is primarily focused on radiofrequency and microwave hyperthermia but, since 3 D simulation studies on heating with ultrasound are now a reality, guidance as well as a benchmark for ultrasound-based hyperthermia are also included.

Progress in Imaging the Human Torso at the Ultrahigh Fields of 7 and 10.5 T

Especially after the launch of 7 T, the ultrahigh magnetic field (UHF) imaging community achieved critically important strides in our understanding of the physics of radiofrequency interactions in the human body, which in turn has led to solutions for the challenges posed by such UHFs. As a result, the originally obtained poor image quality has progressed to the high-quality and high-resolution images obtained at 7 T and now at 10.5 T in the human torso. Despite these tremendous advances, work still remains to further improve the image quality and fully capitalize on the potential advantages UHF has to offer.

Simulation study of the cooling effect of blood vessels and blood coagulation in hepatic radio-frequency ablation

PURPOSE

Computer simulations of hepatic radio-frequency ablation (RFA) were performed to: (i) determine the dependence of the vessel wall heat transfer coefficient on geometrical parameters; (ii) study the conditions required for the occurrence of the directional effect of blood; and (iii) classify blood vessels according to their effect on the thermal lesion while considering blood coagulation. The information thus obtained supports the development of a multi-scale bio-heat model tailored for more accurate prediction of hepatic RFA outcomes in the vicinity of blood vessels.

MATERIALS AND METHODS

The simulation geometry consisted of healthy tissue, tumor tissue, a mono-polar RF-needle, and a single cylindrical blood vessel. The geometrical parameters of interest were the RF-needle active length and those describing blood vessel configuration. A simple, novel method to incorporate the effects of blood coagulation into the simulation was developed and tested.

RESULTS

A closed form expression giving the dependence of the vessel wall heat transfer coefficient on geometrical parameters was obtained. Directional effects on the thermal lesion were found to occur for blood vessel radii between 0.4 mm and 0.5 mm. Below 0.4 mm blood coagulation blocked the flow.

CONCLUSIONS

The closed form expression for the heat transfer coefficient can be used in models of RFA to speed up computation. The conditions on vessel radii required for the occurrence of directional effects on the thermal lesion were determined. These conditions allow the classification of blood vessels. Different approximations to the thermal equation can thus be used for these vessel classes.

A Brief on Biological Thermodynamics for Human Physiology

Thermodynamics, the science of energy interactions, governs the direction of processes found in nature. While the subject finds wide applications in science and technology, its connection to biological sciences and in particular to bio-engineering is becoming increasingly important. In this work, after a brief introduction to the fundamental concepts in thermodynamics, we focus on its application in human physiology. A review of application of thermodynamics to the interaction between human body and environment is presented. Research works on biological systems such as the nervous system and the cardiovascular systems are summarized. The thermodynamics of metabolism is reviewed, and finally, the role of the subject in understanding and combating diseases is highlighted.

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.

Validation of a coupled electromagnetic and thermal model for estimating temperatures during magnetic nanoparticle hyperthermia

PURPOSE

Alternating magnetic field (AMF) tissue interaction models are generally not validated. Our aim was to develop and validate a coupled electromagnetic and thermal model for estimating temperatures in large organs during magnetic nanoparticle hyperthermia (MNH).

MATERIALS AND METHODS

Coupled finite element electromagnetic and thermal model validation was performed by comparing the results to experimental data obtained from temperatures measured in homogeneous agar gel phantoms exposed to an AMF at fixed frequency (155 ± 10 kHz). The validated model was applied to a three-dimensional (3D) rabbit liver built from computed tomography (CT) images to investigate the contribution of nanoparticle heating and nonspecific eddy current heating as a function of AMF amplitude.

RESULTS

Computed temperatures from the model were in excellent agreement with temperatures calculated using the analytical method (error < 1%) and temperatures measured in phantoms (maximum absolute error <2% at each probe location). The 3D rabbit liver model for a fixed concentration of 5 mg Fe/cm3 of tumor revealed a maximum temperature ∼44 °C in tumor and ∼40 °C in liver at AMF amplitude of ∼12 kA/m (peak).

CONCLUSION

A validated coupled electromagnetic and thermal model was developed to estimate temperatures due to eddy current heating in homogeneous tissue phantoms. The validated model was successfully used to analyze temperature distribution in complex rabbit liver tumor geometry during MNH. In future, model validation should be extended to heterogeneous tissue phantoms, and include heat sink effects from major blood vessels.

Feasibility and characterization of a safe susceptibility-matched endorectal coil for MR spectroscopy

When using endorectal coils, local radiofrequency (RF) heating may occur in the surrounding tissue. Furthermore, most endorectal coils create a susceptibility artifact detrimental to both anatomical magnetic resonance imaging (MRI) and spectroscopy (MRS) acquisitions. We aimed at assessing the safety and MRS performance of a susceptibility-matched endorectal coil for further rectal wall analysis. Experiments were performed on a General Electric MR750 3 T scanner. A variable number of miniaturized passive RF traps were incorporated in the reception cable. The assessment of RF heating and coil sensitivity was conducted on a 1.5% agar-agar phantom doped with NaCl. Several susceptibility-matched materials such as Ultem, perfluorocarbon and barium sulfate were then compared with an external coil. Finally, Ultem was used as a solid support for an endorectal coil and compared with a reference coil. Phantom experiments exhibited a complete suppression of both the RF heating phenomenon and the coil sensitivity artifact. Ultem was the material that produced the smallest image distortion. The full width at half maximum of MR spectra acquired using the susceptibility-matched endorectal coil showed at least 30% narrowing compared with a reference endorectal coil. A susceptibility-matched endorectal coil with RF traps incorporated was validated on phantoms. This coil appears to be a promising device for future in vivo experiments.

RF heating measurement using MR thermometry and field monitoring: Methodological considerations and first in vivo results

PURPOSE

A MR thermometry (MRT) method with field monitoring is proposed to improve the measurement of small temperature variations induced in brain MRI exams.

METHODS

MR thermometry experiments were performed at 7 Tesla with concurrent field monitoring and RF heating. Images were reconstructed with nominal k-space trajectories and with first-order spherical harmonics correction. Experiments were performed in vitro with deliberate field disturbances and on an anesthetized macaque in 2 different specific absorption rate regimes, that is, at 50% and 100% of the maximal specific absorption rate level allowed in the International Electrotechnical Commission normal mode of operation. Repeatability was assessed by running a second separate session on the same animal.

RESULTS

Inclusion of magnetic field fluctuations in the reconstruction improved temperature measurement accuracy in vitro down to 0.02°C. Measurement precision in vivo was on the order of 0.15°C in areas little affected by motion. In the same region, temperature increase reached 0.5 to 0.8°C after 20 min of heating at 100% specific absorption rates and followed a rough factor of 2 with the 50% specific absorption rate scans. A horizontal temperature plateau, as predicted by Pennes bioheat model with thermal constants from the literature and constant blood temperature assumption, was not observed.

CONCLUSION

Inclusion of field fluctuations in image reconstruction was beneficial for the measurement of small temperature rises encountered in standard brain exams. More work is needed to correct for motion-induced field disturbances to extract reliable temperature maps.

Wire-based sternal closure: MRI-related heating at 1.5 T/64 MHz and 3 T/128 MHz based on simulation and experimental phantom study

PURPOSE

The paper investigates factors that affect the RF-induced heating for commonly used wire-based sternal closure under 1.5 T and 3 T MRI systems and clarifies the heating mechanisms.

METHODS

Numerical simulations based on the finite-difference time-domain method and experimental measurements in ASTM (American Society for Testing and Materials) phantom were used in the study. Various configurations of the wire-based sternal closure in the phantom were studied based on parameter sweeps to understand key factors related to the RF-induced heating. In vivo simulations were further performed to explore the RF-induced heating in computational human phantoms for clinically relevant scenarios.

RESULTS

The wire-based sternal closure can lead to peak 1-g averaged spatial absorption ratio of 106.3 W/kg and 75.2 W/kg in phantom and peak 1-g averaged specific absorption rate of 32.1 W/kg and 62.1 W/kg in computational human models near the device at 1.5 T and 3 T, respectively. In phantom, the simulated maximum temperature rises for 15-minute RF exposure are 9.4°C at 1.5 T and 5.8°C at 3 T. Generally, the RF-induced heating will be higher when the electrical length of the device is close to the resonant length or when multiple components are spaced closely along the longitudinal direction.

CONCLUSION

The RF-induced heating related to wire-based sternal closure can be significant due to the antenna effect and capacitive mutual coupling effect related to the specific geometries of devices.

Evolution of UHF Body Imaging in the Human Torso at 7T: Technology, Applications, and Future Directions

The potential value of ultrahigh field (UHF) magnetic resonance imaging (MRI) and spectroscopy to biomedical research and in clinical applications drives the development of technologies to overcome its many challenges. The increased difficulties of imaging the human torso compared with the head include its overall size, the dimensions and location of its anatomic targets, the increased prevalence and magnitude of physiologic effects, the limited availability of tailored RF coils, and the necessary transmit chain hardware. Tackling these issues involves addressing notoriously inhomogeneous transmit B1 (B1) fields, limitations in peak B1, larger spatial variations of the static magnetic field B0, and patient safety issues related to implants and local RF power deposition. However, as research institutions and vendors continue to innovate, the potential gains are beginning to be realized. Solutions overcoming the unique challenges associated with imaging the human torso are reviewed as are current studies capitalizing on the benefits of UHF in several anatomies and applications. As the field progresses, strategies associated with the RF system architecture, calibration methods, RF pulse optimization, and power monitoring need to be further integrated into the MRI systems making what are currently complex processes more streamlined. Meanwhile, the UHF MRI community must seize the opportunity to build upon what have been so far proof of principle and feasibility studies and begin to further explore the true impact in both research and the clinic.

Journal of Biomechanical Engineering Legacy Paper 2018

The Journal of Biomechanical Engineering has contributed to biomechanical engineering field since 1977. To honor papers published at least 30 years that have had a long-lasting impact on the field, the Editors now recognize "Legacy Papers." The journal is pleased to present the following paper as this year's Legacy Paper: "A New Simplified Bioheat Equation for the Effect of Blood Flow on Local Average Tissue Temperature" by S. Weinbaum and L. Jiji, ASME Journal of Biomechanical Engineering 107: 131-139, 1985.

Mathematical and numerical analysis of thermal distribution in cancerous tissues under the local heat therapy

The main purpose of this study is to investigate the thermal behavior of living tissues in the presence of spatial external heat source. An effort has been made to formulate the mathematical model to study the temperature distribution in in vivo tissues of the human body. The mathematical formulation is governed by bio-heat equation together with appropriate initial, boundary and interface conditions. The solution of the model was carried out using variational finite element method and computational simulations. The model describes the exchange of heat between the internal biological tissues and other surrounding media. The effect of external heat source under different conditions of atmospheric temperature and as a local hyperthermic method provides an important information to the temperature regulation in biological tissues under normal and malignant conditions. Thermal fluctuations at the targeted regions were obtained with respect to various time-dependent heating sources and scattering coefficients. The results obtained may be helpful for clinical purposes especially in the treatment of cancerous tumors through radiotherapy and other local hyperthermic approaches.

New Frontiers for Applications of Thermal Infrared Imaging Devices: Computational Psychopshysiology in the Neurosciences

Thermal infrared imaging has been proposed, and is now used, as a tool for the non-contact and non-invasive computational assessment of human autonomic nervous activity and psychophysiological states. Thanks to a new generation of high sensitivity infrared thermal detectors and the development of computational models of the autonomic control of the facial cutaneous temperature, several autonomic variables can be computed through thermal infrared imaging, including localized blood perfusion rate, cardiac pulse rate, breath rate, sudomotor and stress responses. In fact, all of these parameters impact on the control of the cutaneous temperature. The physiological information obtained through this approach, could then be used to infer about a variety of psychophysiological or emotional states, as proved by the increasing number of psychophysiology or neurosciences studies that use thermal infrared imaging. This paper presents a review of the principal achievements of thermal infrared imaging in computational psychophysiology, focusing on the capability of the technique for providing ubiquitous and unwired monitoring of psychophysiological activity and affective states. It also presents a summary on the modern, up-to-date infrared sensors technology.

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.

SAR Simulations & Safety

At ultra-high fields, the assessment of radiofrequency (RF) safety presents several new challenges compared to low-field systems. Multi-channel RF transmit coils in combination with parallel transmit techniques produce time-dependent and spatially varying power loss densities in the tissue. Further, in ultra-high-field systems, localized field effects can be more pronounced due to a transition from the quasi stationary to the electromagnetic field regime. Consequently, local information on the RF field is required for reliable RF safety assessment as well as for monitoring of RF exposure during MR examinations. Numerical RF and thermal simulations for realistic exposure scenarios with anatomical body models are currently the only practical way to obtain the requisite local information on magnetic and electric field distributions as well as tissue temperature. In this article, safety regulations and the fundamental characteristics of RF field distributions in ultra-high-field systems are reviewed. Numerical methods for computation of RF fields as well as typical requirements for the analysis of realistic multi-channel RF exposure scenarios including anatomical body models are highlighted. In recent years, computation of the local tissue temperature has become of increasing interest, since a more accurate safety assessment is expected because temperature is directly related to tissue damage. Regarding thermal simulation, bio-heat transfer models and approaches for taking into account the physiological response of the human body to RF exposure are discussed. In addition, suitable methods are presented to validate calculated RF and thermal results with measurements. Finally, the concept of generalized simulation-based specific absorption rate (SAR) matrix models is discussed. These models can be incorporated into local SAR monitoring in multi-channel MR systems and allow the design of RF pulses under constraints for local SAR.

Validating subject-specific RF and thermal simulations in the calf muscle using MR-based temperature measurements

PURPOSE

Ongoing discussions occur to translate the safety restrictions on MR scanners from specific absorption rate (SAR) to thermal dose. Therefore, this research focuses on the accuracy of thermal simulations in human subjects during an MR exam, which is fundamental information in that debate.

METHODS

Radiofrequency (RF) heating experiments were performed on the calves of 13 healthy subjects using a dedicated transmit-receive coil while monitoring the temperature with proton resonance frequency shift (PRFS) thermometry. Subject-specific models and one generic model were used for electromagnetic and thermal simulations using Pennes' bioheat equation, with the blood equilibration constant equaling zero. The simulations were subsequently compared with the experimental results.

RESULTS

The mean B1+ equaled 15 µT in the center slice of all volunteers, and 95% of the voxels had errors smaller than 2.8 µT between the simulation and measurement. The intersubject variation in RF power to achieve the required B1+ was 11%. The resulting intersubject variation in median temperature rise was 14%. Thermal simulations underestimated the median temperature increase on average, with 34% in subject-specific models and 28% in the generic model.

CONCLUSIONS

Although thermal measures are directly coupled to tissue damage and therefore suitable for RF safety assessment, insecurities in the applied thermal modeling limit their estimation accuracy. Magn Reson Med 77:1691-1700, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Modeling of electrodes and implantable pulse generator cases for the analysis of implant tip heating under MR imaging

PURPOSE

The authors' purpose is to model the case of an implantable pulse generator (IPG) and the electrode of an active implantable medical device using lumped circuit elements in order to analyze their effect on radio frequency induced tissue heating problem during a magnetic resonance imaging (MRI) examination.

METHODS

In this study, IPG case and electrode are modeled with a voltage source and impedance. Values of these parameters are found using the modified transmission line method (MoTLiM) and the method of moments (MoM) simulations. Once the parameter values of an electrode/IPG case model are determined, they can be connected to any lead, and tip heating can be analyzed. To validate these models, both MoM simulations and MR experiments were used. The induced currents on the leads with the IPG case or electrode connections were solved using the proposed models and the MoTLiM. These results were compared with the MoM simulations. In addition, an electrode was connected to a lead via an inductor. The dissipated power on the electrode was calculated using the MoTLiM by changing the inductance and the results were compared with the specific absorption rate results that were obtained using MoM. Then, MRI experiments were conducted to test the IPG case and the electrode models. To test the IPG case, a bare lead was connected to the case and placed inside a uniform phantom. During a MRI scan, the temperature rise at the lead was measured by changing the lead length. The power at the lead tip for the same scenario was also calculated using the IPG case model and MoTLiM. Then, an electrode was connected to a lead via an inductor and placed inside a uniform phantom. During a MRI scan, the temperature rise at the electrode was measured by changing the inductance and compared with the dissipated power on the electrode resistance.

RESULTS

The induced currents on leads with the IPG case or electrode connection were solved for using the combination of the MoTLiM and the proposed lumped circuit models. These results were compared with those from the MoM simulations. The mean square error was less than 9%. During the MRI experiments, when the IPG case was introduced, the resonance lengths were calculated to have an error less than 13%. Also the change in tip temperature rise at resonance lengths was predicted with less than 4% error. For the electrode experiments, the value of the matching impedance was predicted with an error less than 1%.

CONCLUSIONS

Electrical models for the IPG case and electrode are suggested, and the method is proposed to determine the parameter values. The concept of matching of the electrode to the lead is clarified using the defined electrode impedance and the lead Thevenin impedance. The effect of the IPG case and electrode on tip heating can be predicted using the proposed theory. With these models, understanding the tissue heating due to the implants becomes easier. Also, these models are beneficial for implant safety testers and designers. Using these models, worst case conditions can be determined and the corresponding implant test experiments can be planned.

Efficient Lattice Boltzmann Solver for Patient-Specific Radiofrequency Ablation of Hepatic Tumors

Radiofrequency ablation (RFA) is an established treatment for liver cancer when resection is not possible. Yet, its optimal delivery is challenged by the presence of large blood vessels and the time-varying thermal conductivity of biological tissue. Incomplete treatment and an increased risk of recurrence are therefore common. A tool that would enable the accurate planning of RFA is hence necessary. This manuscript describes a new method to compute the extent of ablation required based on the Lattice Boltzmann Method (LBM) and patient-specific, pre-operative images. A detailed anatomical model of the liver is obtained from volumetric images. Then a computational model of heat diffusion, cellular necrosis, and blood flow through the vessels and liver is employed to compute the extent of ablated tissue given the probe location, ablation duration and biological parameters. The model was verified against an analytical solution, showing good fidelity. We also evaluated the predictive power of the proposed framework on ten patients who underwent RFA, for whom pre- and post-operative images were available. Comparisons between the computed ablation extent and ground truth, as observed in postoperative images, were promising (DICE index: 42%, sensitivity: 67%, positive predictive value: 38%). The importance of considering liver perfusion while simulating electrical-heating ablation was also highlighted. Implemented on graphics processing units (GPU), our method simulates 1 minute of ablation in 1.14 minutes, allowing near real-time computation.

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

PURPOSE

Present a novel method for rapid prediction of temperature in vivo for a series of pulse sequences with differing levels and distributions of specific energy absorption rate (SAR).

THEORY AND METHODS

After the temperature response to a brief period of heating is characterized, a rapid estimate of temperature during a series of periods at different heating levels is made using a linear heat equation and impulse-response (IR) concepts. Here the initial characterization and long-term prediction for a complete spine exam are made with the Pennes' bioheat equation where, at first, core body temperature is allowed to increase and local perfusion is not. Then corrections through time allowing variation in local perfusion are introduced.

RESULTS

The fast IR-based method predicted maximum temperature increase within 1% of that with a full finite difference simulation, but required less than 3.5% of the computation time. Even higher accelerations are possible depending on the time step size chosen, with loss in temporal resolution. Correction for temperature-dependent perfusion requires negligible additional time and can be adjusted to be more or less conservative than the corresponding finite difference simulation.

CONCLUSION

With appropriate methods, it is possible to rapidly predict temperature increase throughout the body for actual MR examinations.

Feasibility of measuring thermoregulation during RF heating of the human calf muscle using MR based methods

PURPOSE

One of the main safety concerns in MR is heating of the subject due to radiofrequency (RF) exposure. Recently was shown that local peak temperatures can reach dangerous values and the most prominent parameter for accurate temperature estimations is thermoregulation. Therefore, the goal of this research is testing the feasibility of measuring thermoregulation in vivo using MR methods.

THEORY AND METHODS

The calves of 13 volunteers were scanned at 3 tesla. A Proton Resonance Frequency Shift method was used for temperature measurement. Arterial Spin Labeling and phase contrast scans were used for perfusion and flow measurements respectively. The calves were monitored during extreme RF exposure (20 W/kg, 16 min) and after physical exercise.

RESULTS

Temperature increases due to RF absorption (range of the 90th percentile of all volunteers: 1.1-2.5°C) matched with the reference skin temperature changes. Increases in perfusion and flow were defined on the whole leg and normalized to baseline. Perfusion showed a significant increase due to RF heating (ratio compared with baseline: 1.28 ± 0.37; P < 0.05), the influence of exercise was much greater, however (2.97 ± 2.45, P < 0.01).

CONCLUSION

This study represents a first exploration of measuring thermoregulation, which will become essential when new safety guidelines are based on thermal dose.

Modeling the effect of adverse environmental conditions and clothing on temperature rise in a human body exposed to radio frequency electromagnetic fields

This study considers the computationally determined thermal profile of a fully clothed, finely discretized, heterogeneous human body model, subject to the maximum allowable reference level for a 1-GHz radio frequency electromagnetic field for a worker, and also subject to adverse environmental conditions, including high humidity and high ambient temperature. An initial observation is that while electromagnetic fields at the occupational safety limit will contribute an additional thermal load to the tissues, and subsequently, cause an elevated temperature, the magnitude of this effect is far outweighed by that due to the conditions including the ambient temperature, relative humidity, and the type of clothing worn. It is envisaged that the computational modeling approach outlined in this paper will be suitably modified in future studies to evaluate the thermal response of a body at elevated metabolic rates, and for different body shapes and sizes including children and pregnant women.

In vivo radiofrequency heating in swine in a 3T (123.2-MHz) birdcage whole body coil

PURPOSE

To study in vivo radiofrequency (RF) heating produced due to power deposition from a 3T (Larmour frequency = 123.2 MHz), birdcage, whole body coil.

METHODS

The RF heating was simulated in a digital swine by solving the mechanistic generic bioheat transfer model (GBHTM) and the conventional, empirical Pennes bioheat transfer equation for two cases: 1) when the swine head was in the isocenter and 2) when the swine trunk was in the isocenter. The simulation results were validated by making direct fluoroptic temperature measurements in the skin, brain, simulated hot regions, and rectum of 10 swine (case 1: n = 5, mean animal weight = 84.03 ± 6.85 kg, whole body average SAR = 2.65 ± 0.22 W/kg; case 2: n = 5, mean animal weight = 81.59 ± 6.23 kg, whole body average SAR = 2.77 ± 0.26 W/kg) during 1 h of exposure to a turbo spin echo sequence.

RESULTS

The GBHTM simulated the RF heating more accurately compared with the Pennes equation. In vivo temperatures exceeded safe temperature thresholds with allowable SAR exposures. Hot regions may be produced deep inside the body, away from the skin.

CONCLUSION

SAR exposures that produce safe temperature thresholds need reinvestigation.

Measurement of SAR-induced temperature increase in a phantom and in vivo with comparison to numerical simulation

PURPOSE

To compare numerically simulated and experimentally measured temperature increase due to specific energy absorption rate from radiofrequency fields.

METHODS

Temperature increase induced in both a phantom and in the human forearm when driving an adjacent circular surface coil was mapped using the proton resonance frequency shift technique of magnetic resonance thermography. The phantom and forearm were also modeled from magnetic resonance image data, and both specific energy absorption rate and temperature change as induced by the same coil were simulated numerically.

RESULTS

The simulated and measured temperature increase distributions were generally in good agreement for the phantom. The relative distributions for the human forearm were very similar, with the simulations giving maximum temperature increase about 25% higher than measured.

CONCLUSION

Although a number of parameters and uncertainties are involved, it should be possible to use numerical simulations to produce reasonably accurate and conservative estimates of temperature distribution to ensure safety in magnetic resonance imaging. R01 EB006563

NEW COMPUTATIONAL APPROACHES FOR BIOPHYSICAL HEAT TRANSFER IN TISSUE UNDER ULTRASONIC WAVES: THE VARIATIONAL ITERATION AND CHEBYSCHEV SPECTRAL SIMULATIONS

A mathematical and numerical study is presented for simulating temperature distribution in a two-dimensional tissue medium using Pennes bioheat transfer equation, when the tissue is subjected to ultrasonic waves. Following nondimensionalization of the governing partial differential equation, a novel variational iteration method (VIM) solution is developed. This excellent technique introduced by He [Variational iteration method — a kind of non-linear analytical technique: Some examples, Int J Non-Linear Mech.34:699–708, 1999] employs Lagrange multipliers which can be identified optimally via variational theory. The space and time distributions of temperature are studied and solutions visualized via Mathematica. The influence of thermal conductivity and relaxation time are also examined. Excellent stability and convergence characteristics of VIM are demonstrated. Validation is achieved with a Chebyschev spectral collocation method (CSCM). The present work demonstrates the excellent potential of this powerful semi-numerical method in nonlinear biological heat transfer and furthermore provides an alternative strategy to conventional finite element and finite difference computational simulations. The model finds applications in minimally-invasive spinal laser treatments, glaucoma therapy in ophthalmology and thermoradiotherapy for malignant tumors.

MRI at 7 Tesla and above: demonstrated and potential capabilities

With more than 40 installed MR systems worldwide operating at 7 Tesla or higher, ultra-high-field (UHF) imaging has been established as a platform for clinically oriented research in recent years. Along with technical developments that, in part, have also been successfully transferred to lower field strengths, MR imaging and spectroscopy at UHF have demonstrated capabilities and potentials for clinical diagnostics in a variety of studies. In terms of applications, this overview article focuses on already achieved advantages for in vivo imaging, i.e., in imaging the brain and joints of the musculoskeletal system, but also considers developments in body imaging, which is particularly challenging. Furthermore, new applications for clinical diagnostics such as X-nuclei imaging and spectroscopy, which only really become feasible at ultra-high magnetic fields, will be presented.

MR safety assessment of potential RF heating from cranial fixation plates at 7 T

PURPOSE

The increasing number of clinically oriented MRI studies at 7 T motivates the safety assessment of implants, since many 7 T research sites conservatively exclude all subjects with metallic implants, regardless of type or location. The purpose of this study was to investigate potential RF-induced heating during a 7 T MRI scan using a self-built transmit/receive RF coil in patients with implants used for refixation of the bone flap after craniotomy. Going beyond standard ASTM safety tests, a comprehensive test procedure for safety assessments at 7 T is presented which takes into account the more complex coupling of the electromagnetic field with the human body and the implant as well as polarization effects.

METHODS

The safety assessment consisted of three main investigations using (1) numerical simulations in simplified models, (2) electric and magnetic field measurements and validation procedures in homogeneous phantoms, and (3) analysis of exposure scenarios in a heterogeneous human body model including thermal simulations. Finally, 7 T in vivo images show the degree of image artifact around the implants.

RESULTS

The simulations showed that the field distortions remain localized within the direct vicinity of the implants. A parallel E-field polarization was found to be the most relevant component in creating local SAR deviations, resulting in a 10% increase in 10-g-averaged SAR and 53% in 1-g-averaged SAR. Using a heterogeneous human head model, the implants caused field distortions and SAR elevations in the numerical simulations which were distinctly lower than the maximum local SAR value caused by the RF coil alone. Also, the position of the maximum 10-g-averaged SAR remained unchanged by the presence of the implants. Similarly, the maximum absolute local temperature remained below 39 °C in the thermal simulations. Only minor artifacts from the implants were observed in the in vivo images that would not likely affect the diagnostic image quality in patients.

CONCLUSIONS

The findings suggested no evidence for noteworthy RF-related heating in humans after craniotomy using the described implants and for the particular RF coil that was used in this study. Here, identical transmit power restrictions apply with or without the implants. For other RF coils, the maximum permissible input power should be reduced by 10% until further simulations may indicate otherwise.

Heating induced near deep brain stimulation lead electrodes during magnetic resonance imaging with a 3 T transceive volume head coil

Heating induced near deep brain stimulation (DBS) lead electrodes during magnetic resonance imaging with a 3 T transceive head coil was measured, modeled, and imaged in three cadaveric porcine heads (mean body weight = 85.47 ± 3.19 kg, mean head weight = 5.78 ± 0.32 kg). The effect of the placement of the extra-cranial portion of the DBS lead on the heating was investigated by looping the extra-cranial lead on the top, side, and back of the head, and placing it parallel to the coil's longitudinal axial direction. The heating was induced using a 641 s long turbo spin echo sequence with the mean whole head average specific absorption rate of 3.16 W kg(-1). Temperatures were measured using fluoroptic probes at the scalp, first and second electrodes from the distal lead tip, and 6 mm distal from electrode 1 (T(6 mm)). The heating was modeled using the maximum T(6 mm) and imaged using a proton resonance frequency shift-based MR thermometry method. Results showed that the heating was significantly reduced when the extra-cranial lead was placed in the longitudinal direction compared to the other placements (peak temperature change = 1.5-3.2 °C versus 5.1-24.7 °C). Thermal modeling and MR thermometry may be used together to determine the heating and improve patient safety online.

Radiofrequency heating in porcine models with a "large" 32 cm internal diameter, 7 T (296 MHz) head coil

Temperatures were measured in vivo in four pigs (mean animal weight = 110.75 kg and standard deviation = 6.13 kg) due to a continuous wave radiofrequency (RF) power irradiation with a 31.75 cm internal diameter and a 15.24 cm long, 7 T (296 MHz), eight channel, transverse electromagnetic head coil. The temperatures were measured in the subcutaneous layer of the scalp, 5, 10, 15, and 20 mm deep in the brain, and rectum using fluoroptic temperature probes. The RF power was delivered to the pig's head for ∼3 h (mean deposition time = 3.14 h and standard deviation = 0.06 h) at the whole head average specific absorption rate of ∼3 W kg(-1) (mean average specific absorption rate = 3.08 W kg(-1) and standard deviation = 0.09 W kg(-1)). Next, simple bioheat transfer models were used to simulate the RF power induced temperature changes. Results show that the RF power produced uniform temperature changes in the pigs' heads (mean temperature change = 1.68°C and standard deviation = 0.13°C) with no plateau achieved during the heating. No thermoregulatory alterations were detected due to the heating because the temperature responses of the pre-RF and post-RF epochs were not statistically significantly different. Simple, validated bioheat models may provide accurate temperature changes.

Effect of the extracranial deep brain stimulation lead on radiofrequency heating at 9.4 Tesla (400.2 MHz)

PURPOSE

To study the effect of the extracranial portion of a deep brain stimulation (DBS) lead on radiofrequency (RF) heating with a transmit and receive 9.4 Tesla head coil.

MATERIALS AND METHODS

The RF heating was studied in four excised porcine heads (mean animal head weight = 5.46 +/- 0.14 kg) for each of the following two extracranial DBS lead orientations: one, parallel to the coil axial direction; two, perpendicular to the coil axial direction (i.e., azimuthal). Temperatures were measured using fluoroptic probes at four locations: one, scalp; two, near the second DBS lead electrode-brain contact; three, near the distal tip of the DBS lead; and four, air surrounding the head. A continuous wave RF power was delivered to each head for 15 min using the coil. Net, delivered RF power was measured at the coil (mean whole head average specific absorption rate = 2.94 +/- 0.08 W/kg).

RESULTS

RF heating was significantly reduced when the extracranial DBS lead was placed in the axial direction (temperature change = 0-5 degrees C) compared with the azimuthal direction (temperature change = 1-27 degrees C).

CONCLUSION

Development of protocols seems feasible to keep RF heating near DBS electrodes clinically safe during ultra-high field head imaging.

Model-based planning and real-time predictive control for laser-induced thermal therapy

In this article, the major idea and mathematical aspects of model-based planning and real-time predictive control for laser-induced thermal therapy (LITT) are presented. In particular, a computational framework and its major components developed by authors in recent years are reviewed. The framework provides the backbone for not only treatment planning but also real-time surgical monitoring and control with a focus on MR thermometry enabled predictive control and applications to image-guided LITT, or MRgLITT. Although this computational framework is designed for LITT in treating prostate cancer, it is further applicable to other thermal therapies in focal lesions induced by radio-frequency (RF), microwave and high-intensity-focused ultrasound (HIFU). Moreover, the model-based dynamic closed-loop predictive control algorithms in the framework, facilitated by the coupling of mathematical modelling and computer simulation with real-time imaging feedback, has great potential to enable a novel methodology in thermal medicine. Such technology could dramatically increase treatment efficacy and reduce morbidity.