Assessment of Radio Frequency Induced Heating On or Near Implants during MRI – some open issues
Mikhail Kozlov1,2 and Gregor Schaefers1,3

1MR:comp GmbH, Gelsenkirchen, Germany, 2Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany, 3Magnetic Resonance Institute for Safety, Technology and Research GmbH, Gelsenkirchen, Germany

Synopsis

We evaluated locations of maximum temperature rise (max(ΔT)) and the dependence of maxT) on RF-induced power deposition (Ptotal) for some generic implants. ΔT spatial and temporal variations were investigated. To fulfill ASTM F2182-11a setup requirements, the temperature probe should be placed with submillimetre precision at location that cannot be predicted by a full wave electromagnetic simulation alone. It is a challenge to validate with small uncertainty Ptotal calculated using EM simulation by only measuring SAR or VLD value at some points in space, if the field probe sensor size is larger than one tenth of the wire diameter.

Introduction

A standard test method for measurement of RF-induced heating (ΔT) on or near passive implants during MRI is described in ASTM F2182-11a [1]. It requires demonstrating the worst-case maximal temperature increase (max(ΔT)) without prescribing a method of worst-case analysis. In some reports worst-case analysis is based on searching cases with the highest value of SAR integrated over 0.1 gram of tissue (SAR0.1g) and estimating maxT) by temperature measurement at the location of the maximum (max(SAR0.1g)). Our numerical simulation goals in this study were: 1) to compare locations of max(ΔT) and max(SAR0.1g); 2) to evaluate the dependence of max(ΔT) on RF-induced power deposition at the hot spot (Ptotal) defined in ISO/TS 10974 [2]; 3) to analyze how ΔT changes along the path normal to the implant end face for two heating time points.

Method

The implants investigated were insulated stainless steel wires of 300 mm in length, of 0.73 and 1.6 mm in diameter, respectively. At one end of the wire, the insulation with thickness of 0.5 mm was removed, at the other side the wire was capped by the same insulation. The wires were placed in the ASTM phantom, which was located inside a 64 MHz high pass 16 rung birdcage coil with a diameter of 704 mm and a length of 650 mm. The ASTM phantom was filled with tissue-simulating medium (TSM) (εr = 78 and σ = 0.47 S/m). All wires were oriented along the Z axis, and shifted by 150 mm in X direction from the coil iso-center. The abbreviation "d×i×l" denotes: wire diameter “d”, outer insulation diameter “i”, length of insulation removed “l”. To evaluate the influence of wire end face shape on max(ΔT), wires with a chamfered wire end were also simulated (0.73×1.46×3_ch and 0.73×1.46×10_ch). ΔT(x,y,z), ΔTtip, and ΔT0.5mm were calculated as the temperature at location with coordinates x, y, z, the maximum wire temperature, and the temperature at location 0.5 mm from the wire tip (150, 0, 150.5), minus the ambient temperature of 22°C, respectively. Ptotal integration volume enveloped an area where volume loss density (VLD) decayed more than 20 dB from the peak value. A path from the wire tip (150, 0, 150) to a point (150, 0, 152) was defined to evaluate ΔT spatial variation for two heating time points (t = 10 and t = 900 s). All results were scaled to incident tangential electrical field Etan = 1 Vpeak-to-peak/m at the wire center.

Results and discussion

The results in Table 1 shows that correlations existed neither between Ptotal and SAR0.1g nor between Ptotal or SAR0.1g and maxT). The area with highest VLD was located in close proximity to the wire tip (Fig. 2). Despite similar VLD distribution in this area, for a wire with l = 10 mm, a larger integration volume resulted in a higher Ptotal. All integration volumes were significantly smaller than the volume of 0.1 gram of TSM. The location of max(SAR0.1g) was far away from the location of max(VLD) and maxT), especially at t = 10 s. Thermal properties of two major thermal sinks - implant wire and TSM – varied considerably. The latter and different contact surface areas between the implant wire and TSM resulted in: 1) significant difference in thermal time constants, for example, the wires 0.73×1.46×3 and 0.73×1.46×3_ch heated up faster than other wires investigated (Fig. 3); 2) varying locations of maxT); 3) essentially smaller ΔTtip compared to maxT) (Fig. 4); 4) temporal change in shape of the thermal hot spot (Fig. 5); 5) less heating of implants with a larger wire diameter.

Conclusion

To fulfill the ASTM F2182-11a setup requirements, the temperature probe should be placed with submillimeter precision at a location that cannot be predicted by an EM simulation alone. Further, in some cases it also is impossible to predict the location of max(ΔT) at a given point in time by either SAR or the temperature measurement at another point in time. It is a challenge to validate with small uncertainty Ptotal calculated using EM simulation by only measuring SAR or VLD value at some points in space, if the field probe sensor size is larger than one tenth of the wire diameter. We doubt that an acceptable uncertainty of Ptotal can be achieved using thermal, SAR or VLD measurement only.

Acknowledgements

This work was supported by the German Federal Ministry of Education and Research (BMBF) and within the European Joint Undertakings ENIAC JU, grant # 16ES0028, DeNeCoR

References

[1] ASTM F2182–11a, “Standard Test Method for Measurement of Radio Frequency Induced Heating On or Near Passive Implants During Magnetic Resonance Imaging”.

[2] ISO/TS 10974 “Assessment of the safety of magnetic resonance imaging for patients with an active implantable medical device” 1st edition 2012.

Figures

Table 1. Results summary.

Figure 1. VLD profiles. a) wire 0.73×1.46×3. b) wire 0.73×1.46×10

Figure 2. max(ΔT) normalized to incident Etan at location (150,0,0)

Figure 3. ΔT along the path. a) t = 10 s. b) t = 900 s.

Figure 4. Temperature profiles scaled at local maximum for 0.73×1.46×10. Dark red is highest temperature, dark blue is 22 °C.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2215