Open access peer-reviewed chapter

MRI RF-Induced Heating in Heterogeneous Human Body with Implantable Medical Device

By Ran Guo, Jianfeng Zheng and Ji Chen

Submitted: June 7th 2017Reviewed: October 2nd 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71384

Downloaded: 352

Abstract

Magnetic resonance imaging (MRI) radio frequency (RF)-induced heating is one of the most important concerns of MRI safety for patients, especially with orthopaedic healthcare products. In this chapter, numerical modelling and simulations were conducted to study the RF-induced heating within a 1.5T and 3T magnetic resonance (MR) environment. Numerical simulations were firstly conducted to study the difference between the cases of implantable medical devices inside the phantom and the human body. Then, numerical modelling were conducted to describe the difference of electromagnetic behaviours between the homogeneous phantom and heterogeneous human tissues. The MRI RF-induced heating due to an implantable medical device behaves significantly different in homogeneous media and in heterogeneous human body. For typical orthopaedic medical devices, such as healthcare products applied to shoulder, humerus, hip, femur and tibia, the properties of the RF-induced heating are different in general phantom and in human body. The hot spot location, as well as worst case configuration were evaluated and it was found that they were determined by the incident field and electromagnetic properties of medium. With further scaling, the RF-induced heating in human body for the orthopedic devices can be assessed by phantom studies.

Keywords

  • MRI
  • RF-induced heating
  • orthopedic implant
  • phantom

1. Introduction

Many of the MR-related injuries and the few fatalities that have occurred were the apparent result of failure to follow safety guidelines or of the use of inappropriate information related to the safety aspects of biomedical implants and devices [1, 2, 3, 4, 5, 6, 7]. The preservation of a safe MR environment requires constant attention to the care of patients and individuals with metallic implants and devices, because the variety and complexity of these objects constantly changes [5, 6, 7]. Therefore, to guard against accidents in the MR environment, it is important to understand the risk associated with implantable medical devices which may cause potential problems.

The radiofrequency coils could send energy, in the form of electromagnetic radiation, into the human body. Since the energy is in the radio frequency range, the radiation is not ionizing. But it still can influence biological tissue. During MR procedures, the majority of the RF power transmitted for imaging or spectroscopy (especially for carbon decoupling) is transformed into heat within the patient’s tissue as a result of resistive losses, through convection, conduction, radiation or evaporation [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. Thus, a potential concern in MRI is the heating of the body during image acquisition.

To evaluate the RF-induced heating, the specific absorption rate (SAR) is applied to determine how much electromagnetic energy is absorbed by the body. SAR is typically expressed in unites of watts per kilogram, or W/kg. So the SAR could be defined as:

SARr=σ2ρE2rE1

where E is the total electric field and σ and ρ are the conductivity and density of biological tissue, respectively. The temperature rise in human body or phantom could be calculated by the total SAR according to the bio heat equation. SAR depends on the pulse sequence and the size, geometry, and conductivity of the absorbing object. To ensure participant safety, SAR in MRI studies is limited to minimize temperature increases.

The first study of human thermal responses to RF radiation-induced heating during an MR procedure was conducted by Schaefer et al. [19]. Temperature changes and other physiological parameters were assessed in volunteer subjects exposed to relatively high, whole-body averaged SARs (approximately 4.0 W/kg). The data indicated that there were no excessive temperature elevations or other deleterious physiological consequences related to the exposure to RF radiation [19].

However, for patients with medical implants, MRI-related RF induced heating is potentially problematic. The evaluation of heating for an implant or device is particularly challenging because of the many factors that affect temperature increase in these items. Variables that impact heating include the following: the specific type of implant or device; the electrical characteristic of the implant or device; the RF wavelength of the MR system; the type of transmit RF coil that is used (i.e., transmit head versus transmit body RF coil); the amount of RF energy delivered (i.e., the specific absorption rate, SAR); the landmark position or body part undergoing MRI relative to the transmit RF coil; and the orientation or configuration of the implant or device relative to the source of transmit RF coil.

In this chapter, it shows the importance of evaluation the MRI-related RF induced heating issues for patient with implantable medical devices. Generally, the estimation and measurement is based on in-vitro numerical simulation and experiment. And assessment methods could be separated into active and passive medical implants, respectively, due to the configuration difference of these devices. With the help of the in-vitro evaluation methods, it provides a highly possible way to estimate the temperature increase for patient with implants or devices during MRI examination.

2. In-vivo and in-vitro

MRI may be contraindicated for a given patient primarily due to its potential risks associated with a metallic implant or device. Although many investigations have been performed using laboratory animals to determine thermoregulatory reactions to tissue heating associated with exposure to RF radiation, these experiments do not directly apply to the conditions that occur during MR procedures, nor can they be extrapolated to provide useful information for various reasons [20, 21]. For example, the pattern of RF absorption or the coupling of radiation to biological tissues is primarily dependent on the organism’s size, anatomical features, duration of exposure, sensitivity of the involved tissues (e.g., some tissues are more “thermal sensitive” than others), and a myriad of other variables [14, 21, 22]. Furthermore, there is no laboratory animal that sufficiently mimics or simulates the thermoregulatory responses of an organism with the dimensions and specific responses to that of a human subject. Therefore, experimental results obtained in laboratory animals cannot be simply “scaled” or extrapolated to predict thermoregulatory or other physiological changes in human subjects exposed to RF radiation-induced heating during MR procedures [14, 15, 22], and. In consideration of the above, in-vitro testing is performed to assess the various MRI issues for implants and devices in order to properly characterize the possible risks.

One of in-vitro methods is to use standard American Society for Testing and Materials (ASTM) phantom. ASTM F2182-11A depicts the guideline to measure the RF heating induced by implanted medical devices in a standard phantom filled with gelled-saline which mimic the muscles [23]. Studies have been conducted to evaluate the RF heating induced by orthopedic implants. Commonly a phantom or homogenous media is used to mimic the environments as the implants locate in human body in experiments and/or numerical simulations [24, 25, 26, 27, 28, 29, 30, 31, 32].

Although the RF-induced heating evaluating method using the phantom filled with gelled-saline is widely used, it is obvious that the RF environment of a human body and a phantom filled with gelled-saline are quite different. The power deposition due to an implant for a given incident RF field is a function of the physical properties of the implant and electrical properties of the surrounding medium. Compared with homogeneous gelled-saline in phantom, human body is an inhomogeneous circumstance which includes different tissues with various permittivity and conductivity in a wide range. Hence, it is necessary to study a feasible guide with in-vitro phantom to assess the RF-induced heating in heterogeneous human body.

2.1. Human body: heterogeneous medium

With the development of computational electromagnetics, anatomical computer models of the human body have been used for nearly four decades for dosimetric applications in electromagnetics (EM) [33] and in medical radiation physics [34]. The most prominent numerical methods used in computational dosimetry of electromagnetic fields are based on finite-difference formulations of the underlying differential equations. For the simulation of both RF fields and induced tissue heating, the finite-difference time-domain (FDTD) method in its formulations by Yee [35] and Patankar [36] is applied to rectilinear grids to optimally handle large voxel models. The reconstructed human model used in this Chapter is from the Virtual Family [37]. It is based on high resolution magnetic resonance images of healthy volunteers. Seventy seven different tissue types were distinguished during the segmentation. Currently, the models are being widely applied in several studies on electromagnetic exposure, device optimization and medical applications. Table 1 shows the characteristics of the anatomical model. Duke model is an anatomical model of adult male which is shown in Figures 1 and 2 . And Table 2 shows the segmented tissues and organs of the model, as well as the electromagnetic properties.

NameAge (years)GenderHeight (m)Mass (kg)BMI (kg/m2)
Duke34Male1.747023.1

Table 1.

The characteristics of the anatomical Duke model.

Figure 1.

The segmented tissues and organs of anatomic body.

Figure 2.

The segmented tissues and organs of anatomic brain.

Tissue or organElectric conductivity (S/m)Relative permittivityDensity (kg/m3)
1.5 T/64 MHz3 T/128 MHz1.5 T/64 MHz3 T/128 MHz
Adrenal gland0.7783050.80416673.947266.78391027.5
Air internal00111.2
Artery1.206671.2486386.444173.1591049.75
Bladder0.2873520.29801424.594321.86071035
Blood vessel1.206671.2486386.444173.1591049.75
Bone0.05952550.067352416.681214.71711908
Brain gray matter0.5108680.5867397.429473.52041044.5
Brain white matter0.2915040.34215167.835852.53381041
Bronchi0.5284150.55934658.889650.57141101.5
Bronchi lumen00111.2
Cartilage0.4521030.48837562.914552.92421099.5
Cerebellum0.7190030.829397116.3579.73771045
Cerebrospinal fluid2.065972.1430197.312484.04061007
Commissure anterior0.2915040.34215167.835852.53381041
Commissure posterior0.2915040.34215167.835852.53381041
Connective tissue0.4743310.49872759.489251.85681525
Cornea1.000581.0587487.377971.45661050.5
Diaphragm0.6882130.71923572.234763.49481090.4
Ear cartilage0.4521030.48837562.914552.92421099.5
Ear skin0.435750.52270492.167965.4371109
Epididymis0.8848710.92640484.527272.12791082
Esophagus0.8778420.91280785.820474.8951040
Esophagus lumen00111.2
Eye lens0.285880.31268450.339242.79111075.5
Eye sclera0.8826730.91766575.299864.99911032
Eye vitreous humor1.503151.5053669.126469.06191004.5
Fat0.06615580.069729913.643612.3711911
Gall bladder1.481791.5764105.44388.8995928
Heart lumen1.206671.2486386.444173.1591049.75
Heart muscle0.6784230.766108106.51484.25731080.8
Hippocampus0.5108680.5867397.429473.52041044.5
Hypophysis0.7783050.80416673.947266.78391053
Hypothalamus0.7783050.80416673.947266.78391053
Intervertebral disc0.4521030.48837562.914552.92421099.5
Kidney cortex0.7413160.852313118.55689.61681049
Kidney medulla0.7413160.852313118.55689.61681044
Large intestine0.6381520.70521494.663976.57221088
Large intestine lumen0.6882130.71923572.234763.49481045.2
Larynx0.4521030.48837562.914552.92421099.5
Liver0.4479840.51089780.559564.25071078.75
Lung0.2889770.31561637.102229.4677394
Mandible0.05952550.067352416.681214.71711908
Marrow0.1543350.16202116.435513.53771028.5
Medulla oblongata0.7190030.829397116.3579.73771045.5
Meniscus0.4521030.48837562.914552.92421099.5
Midbrain0.7190030.829397116.3579.73771045.5
Mucosa0.4880390.54420276.723361.58521102
Muscle0.6882130.71923572.234763.49481090.4
Nerve0.3121740.35380255.062144.06531075
Pancreas0.7783050.80416673.947266.78391086.5
Patella0.05952550.067352416.681214.71711908
Penis0.4293110.47893468.636855.98881101.5
Pharynx00111.2
Pineal body0.7783050.80416673.947266.78391053
Pons0.7190030.829397116.3579.73771045.5
Prostate0.8848710.92640484.527272.12791045
SAT0.06615580.069729913.643612.3711911
Skin0.435750.52270492.167965.4371109
Skull0.05952550.067352416.681214.71711908
Small intestine1.591451.69285118.36387.97251030
Small intestine lumen0.6882130.71923572.234763.49481045.2
Spinal cord0.3121740.35380255.062144.06531075
Spleen0.7439140.835186110.55982.89171089
Stomach0.8778420.91280785.820474.8951088
Stomach lumen0.6882130.71923572.234763.49481045.2
Teeth0.05952550.067352416.681214.71712180
Tendon ligament0.4743310.49872759.489251.85681142
Testis0.8848710.92640484.527272.12791082
Thalamus0.5108680.5867397.429473.52041044.5
Thymus0.7783050.80416673.947266.78391023
Thyroid gland0.7783050.80416673.947266.78391050
Tongue0.6521450.68713775.299864.99911090.4
Trachea0.5284150.55934658.889650.57141080
Trachea lumen00111.2
Ureter Urethra0.4293110.47893468.636855.98881101.5
Vein1.206671.2486386.444173.1591049.75
Vertebrae0.05952550.067352416.681214.71711908

Table 2.

The electromagnetic properties of the segmented tissues and organs.

2.2. ASTM phantom: in-vitro measurement

The standard F2182 describe a test method for measurement of RF induced heating on or near passive implants and its surrounding during MRI procedure. A design of phantom container is introduced in the standard with its dimension shown in Figure 3 . The material of phantom container are electrical insulators and non-magnetic and non-metallic. The phantom container is filled with a gelled-saline which has a relative permittivity εr = 80.4 and conductivity of σ = 0.47 S/m. In order to have a great conductivity and viscosity, a suitable gelled saline should be made with 1.32 g/L NaCl and 10 g/L polyacrylic acid (PAA) in water. Numerical simulations indicate that the maximum electric field inside the ASTM phantom is at mid-axial plane about 2 cm away from the vertical phantom side wall. To maximize the heating, and thereby maximizing the signal-to-noise ratio, we placed the implants at this location.

Figure 3.

The structure and dimension of standard ASTM phantom.

A generic RF transmit body coil is developed and shown in Figure 4 . The upper plots represent a 1.5 T RF coil and the lower two plots represent a 3 T RF coil. A physical coil is usually difficult to model and it takes much longer simulation time to reach the steady state of the simulation. It has been shown that using a non-physical coil could significantly reduce the simulation time while providing the same result as that from a physical coil. Thus, rather than modeling the exact physical coil, the non-physical coils were modeled in this study. The two coils have the same dimensions, and both have 8 rungs. The diameter of the RF coil is 63 cm, and the height of the RF coil is 65 cm. The eight parallel lines or the rungs are one dimensional line current excitation. The end rings on top and bottom of the RF coils are tuning capacitors which are also modelled as one dimensional line segments.

Figure 4.

The generic coil model of 1.5-T RF coil (top) and 3-T RF coil (bottom) in SEMCAD X.

The capacitance value is determined from several broadband simulations so that the second highest resonant frequency was adjusted to 64 MHz for 1.5-T and 128 MHz for 3-T systems. The detailed steps are: set an initial capacitance value for all capacitors on end rings and add a broadband pulse signal on one single rung. The other seven rungs are modeled as zero ohm resistors. After the simulation is finished, the power spectrum is extracted. If the second highest resonant frequency is not at appropriate resonant frequencies, the capacitance needs to be adjusted. From this study, the capacitance for the end ring tuning capacitor values is 7.2 pF for 1.5-T RF coil and 1.3 pF for 3-T RF coil.

3. Passive implantable medical device

Any device intended to be totally or partially introduced into the human body through surgical intervention and intended to remain in place after the procedure for at a long-term duration is considered as an implantable device. Passive devices in terms of their form of operation can be classified as device used for transportation and storage of pharmaceutical liquid, device for alteration of blood, body fluids, medical dressing, surgical instruments; reusable surgical instruments, disposable aseptic device, implantable device, device for contraception and birth control, device for sterilization and cleaning, patient care device, in vitro diagnostic reagent, as well as other passive contacting device or passive supplementary device.

In this chapter, three typical categories of orthopedic implantable devices, bone plate system, hip prostheses and tibia intramedullary nails, are selected for MRI related RF induced heating study which are shown in Figures 5 7 . The configuration of each implantable device is shown in figure. For bone fragment compression plate, it is designed to offer multiple compression and reconstruction plating options for the treatment of bone fractures. The application of hip prostheses is related to hip revision and arthroplasty. As for intramedullary nails, they are characterized by the anatomic shape, which is intended to replicate the natural anatomic shape of the bones. They have been designed to help restore the shape of the bone and treat the fractured bones.

Figure 5.

The bone plant system of AxSOS system from Stryker®.

Figure 6.

The hip prostheses of Excia® T from Aesculap®.

Figure 7.

The tibia intramedullary nails of PROTect™ from Depuy Synthes.

4. Numerically evaluate RF-induced heating

4.1. FDTD method

In this numerical investigation, we use the finite difference time domain (FDTD) based SEMCAD X 14.8 (SPEAG) simulation platform. Graphics processing unit (GPU) hardware acceleration was achieved using the SPEAG CUDA library with Tesla C2075 graphic card which is can handle up to 90 million cells. To assure convergence of the numerical simulations, the simulation time was set to 20 periods for each simulation. Additionally, the convergence was checked after the simulations were finished. The material of orthopedic devices is set to perfect electric conductor (PEC), and all the numerical results are normalized to a whole body average SAR of 2 W/kg. The SAR distribution is studied for each case.

4.2. Bone plate system

To ensure a comprehensive comparison, the 1g local average peak SAR value at device is extracted for each configuration of femur and humerus system. Tables 3 and 4 show the value for femur and humerus system. For each numerical result, whole-body average SAR is normalized to 2 W/kg. Since the interaction between RF induced field and implant is dependent on the physical structure of device, the heating effect variations related to the length of plate and screw are studied separately. For femur system, the plate length varies from 100 to 300 mm, and the screw length changes from 10 to 32 mm. For humerus system, the screw dimension is the same as femur system. But the plate length varies only from 100 to 250 mm due to the limit of bone structure. The plate length is studied at first for minimum and maximum screw length. Then the screw length is investigated under the worst case of plate length study which has the highest 1g average peak SAR value for in-vivo simulation. Figures 8 13 show the results which are corresponding to femur and humerus plate system. The solid and dash curve and indicate the in-vivo and in-vitro results, respectively.

Plate length (mm)Screw length (mm)1.5 T/64 MHz3 T/128 MHz
In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)
1001064.20125.9779.7580.90
1501094.82178.6274.5264.01
17510107.00190.4468.5050.20
20010116.65188.8763.3744.69
22510117.00185.0463.1037.72
25010123.00169.7561.2337.27
27510117.00149.8153.8035.33
30010105.00134.9142.3738.53
1003285.02100.2288.2671.90
15032108.17135.9055.4847.74
20032104.91140.2751.5641.59
25032111.73123.1748.9437.09
3003279.94128.9327.7740.70
250(1.5 T)100(3 T)15120.00150.9179.3072.19
250(1.5 T)100(3 T)20121.00135.7176.2068.00
250(1.5 T)100(3 T)25123.00132.1564.266.52

Table 3.

Peak 1g averaged SAR of femur system for in-vivo and in-vitro cases.

Plate length (mm)Screw length (mm)1.5 T/64 MHz3 T/128 MHz
In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)
1001038.02135.6165.4790.51
1501063.74192.1494.6859.90
2001069.47204.5781.8445.54
25010104.57193.6186.7341.50
1003230.7092.1763.4664.6
1503254.10124.7685.3953.70
2003255.90161.9775.1943.39
25032109.00156.38108.9242.51
2501595.53161.6477.8235.33
2502091.02163.5572.9138.94
2502588.27164.5468.4041.32

Table 4.

The peak 1g average SAR value of humerus system for in-vivo and in-vitro cases.

Figure 8.

The femur bone plate length study of 10 mm screw for 1.5 T (left) and 3 T (right).

Figure 9.

The femur bone plate length study of 32 mm screw for 1.5 T (left) and 3 T (right).

Figure 10.

The femur screw length study for 1.5 T (left) and 3 T (right).

Figure 11.

The humerus bone plate length study of 10 mm screw for 1.5 T (left) and 3 T (right).

Figure 12.

The humerus bone plate length study of 32 mm screw for 1.5 T (left) and 3 T (right).

Figure 13.

The humerus screw length study for 1.5 T (left) and 3 T (right).

4.3. Hip prostheses

For hip prostheses, the 1g average local peak SAR value at device is also extracted for each configuration. Table 5 shows the value for hip system of various dimensions. The height of hip prostheses stem ranges from 100 to 170 mm. For in-vivo simulation, the stem is inserted into the bone marrow. And for the trochanter region, the hip prostheses is touching with soft tissue and muscle. Figure 14 represents the results of hip prostheses. The solid and dash curve and indicate the in-vivo and in-vitro results, respectively.

Stem height (mm)1.5 T/64 MHz3 T/128 MHz
In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)
10083.2764270.84343.474655.4143
11081.5782260.50235.824756.3111
12077.8568248.11624.462757.0498
13074.9259247.63521.909360.352
14064.5142237.5116.305961.7026
15062.8772205.32813.148263.2055
16059.6273221.78111.559963.9114
17055.0469213.60211.187763.8915

Table 5.

Peak 1g average SAR of hip prostheses system for in-vivo and in-vitro cases.

Figure 14.

The hip prostheses stem length study for 1.5 T (left) and 3 T (right).

4.4. Tibia intramedullary nails

The 1g average local peak SAR value at device is also extracted for each configuration of tibia intramedullary nails. The length of stem ranges from 255 to 360 mm. The entire nail is penetrated into the bone marrow. The four screws are inserted perpendicularly through the nail and bone. Table 6 shows the value for nail system of various dimensions. Figure 15 represents the results of tibia intramedullary nails. The solid and dash curve and indicate the in-vivo and in-vitro results, respectively.

Nail length (mm)1.5 T/64 MHz3 T/128 MHz
In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)In-vivo SAR1g (W/kg)In-vitro SAR1g (W/kg)
25577.6968136.33193.553255.4068
27077.6346129.24993.234153.6038
28578.9892122.2889.176749.8946
30082.6751115.74591.564846.8385
31581.7577109.2892.734644.0529
33074.7444103.31988.998941.6544
34566.609597.323290.779539.3275
36067.346991.816890.180937.1875

Table 6.

Peak 1g average SAR of tibia nails system for in-vivo and in-vitro cases.

Figure 15.

The nail length study for 1.5 T (left) and 3 T (right).

5. Summary

From the comparison between in-vitro and in-vivo simulations, the RF-induced heating are different because of the variance of incident electric field and surrounding medium. For incident field study, the antenna resonance effect would mainly lead to a heating issue for both in-vitro and in-vivo situation. Although the wavelength of human muscle and gelled-saline nearly equals to each other, due to the variance of incident RF field, the device dimension causing the resonance would be different. Hence, the trend of peak 1g average SAR value along with plate length is unlike from in-vitro to in-vivo circumstance. Additionally, when the screw is inserted across the human bone into the muscle, a huge amount of power would dissipated to the human tissue through the screw tip so that induce a large peak SAR value.

Based on the comparison result, conservatively, the in-vitro method, such as ASTM phantom, could be used to assess RF-induced heating. However, to accurately assess the RF-induced heating in heterogeneous human body with implantable medical device, due to the limit of homogeneous ASTM phantom, it still needs some improvement to handle several particular cases, especially, when the implantable device is penetrating through various human tissues and organs.

Disclaimer

The mention of commercial products, their sources, or their use in connection with material reported herein is not to be construed as either an actual or implied endorsement of such products by the Department of Health and Human Services.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Ran Guo, Jianfeng Zheng and Ji Chen (December 20th 2017). MRI RF-Induced Heating in Heterogeneous Human Body with Implantable Medical Device, High-Resolution Neuroimaging - Basic Physical Principles and Clinical Applications, Ahmet Mesrur Halefoğlu, IntechOpen, DOI: 10.5772/intechopen.71384. Available from:

chapter statistics

352total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Advanced Brain Tumour Segmentation from MRI Images

By Kavitha Angamuthu Rajasekaran and Chellamuthu Chinna Gounder

Related Book

First chapter

Advances in Cardiac Computed Tomography

By Karthik Ananthasubramaniam, Nishtha Sareen and Gjeka Rudin

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us