Synopsis
We evaluated locations of maximum temperature rise (max(ΔT)) and the dependence of max(ΔT)
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
max(Δ
T) 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 V
peak-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
max(Δ
T).
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
max(Δ
T),
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
max(Δ
T); 3)
essentially smaller Δ
Ttip compared to
max(Δ
T) (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,
DeNeCoRReferences
[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.