In silico evaluation of the thermal stress induced by MRI switched gradient fields in patients with metallic hip implant

Simplified modeling of implanted medical devices with metallic filamentary closed loops exposed to low or medium frequency magnetic fields

BACKGROUND AND OBJECTIVES

Electric currents are induced in implanted medical devices with metallic filamentary closed loops (e.g., fixation grids, stents) when exposed to time varying magnetic fields, as those generated during certain diagnostic and therapeutic biomedical treatments. A simplified methodology to efficiently compute these currents, to estimate the altered electromagnetic field distribution in the biological tissues and to assess the consequent biological effects is proposed for low or medium frequency fields.

METHODS

The proposed methodology is based on decoupling the handling of the filamentary wire and the anatomical body. To do this, a circuital solution is adopted to study the metallic filamentary implant and this solution is inserted in the electromagnetic field solution involving the biological tissues. The Joule losses computed in the implant are then used as a forcing term for the thermal problem defined by the bioheat Pennes' equation. The methodology is validated against a model problem, where a reference solution is available.

RESULTS

The proposed simplified methodology is proved to be in good agreement with solutions provided by alternative approaches. In particular, errors in the amplitude of the currents induced in the wires result to be always lower than 3%. After the validation, the methodology is applied to check the interactions between the magnetic field generated by different biomedical devices and a skull grid, which represents a complex filamentary wire implant.

CONCLUSIONS

The proposed simplified methodology, suitable to be applied to closed loop wires in the low to intermediate frequency range, is found to be sufficiently accurate and easy to apply in realistic exposure scenarios. This modeling tool allows analyzing different types of small implants, from coronary and biliary duct stents to orthopedic grids, under a variety of exposure scenarios.

Computational dosimetry in MRI in presence of hip, knee or shoulder implants: do we need accurate surgery models?

Objective.To quantify the effects of different levels of realism in the description of the anatomy around hip, knee or shoulder implants when simulating, numerically, radiofrequency and gradient-induced heating in magnetic resonance imaging. This quantification is needed to define how precise the digital human model modified with the implant should be to get realistic dosimetric assessments.Approach. The analysis is based on a large number of numerical simulations where four 'levels of realism' have been adopted in modelling human bodies carrying orthopaedic implants.Main results. Results show that the quantification of the heating due to switched gradient fields does not strictly require a detailed local anatomical description when preparing the digital human model carrying an implant. In this case, a simple overlapping of the implant CAD with the body anatomy is sufficient to provide a quite good and conservative estimation of the heating. On the contrary, the evaluation of the electromagnetic field distribution and heating caused by the radiofrequency field requires an accurate description of the tissues around the prosthesis.Significance. The results of this paper provide hints for selecting the 'level of realism' in the definition of the anatomical models with embedded passive implants when performing simulations that should reproduce, as closely as possible, thein vivoscenarios of patients carrying orthopaedic implants.

Classification Scheme of Heating Risk during MRI Scans on Patients with Orthopaedic Prostheses

Due to the large variety of possible clinical scenarios, a reliable heating-risk assessment is not straightforward when patients with arthroplasty undergo MRI scans. This paper proposes a simple procedure to estimate the thermal effects induced in patients with hip, knee, or shoulder arthroplasty during MRI exams. The most representative clinical scenarios were identified by a preliminary frequency analysis, based on clinical service databases, collecting MRI exams of 11,658 implant carrier patients. The thermal effects produced by radiofrequency and switching gradient fields were investigated through 588 numerical simulations performed on an ASTM-like phantom, considering four prostheses, two static field values, seven MR sequences, and seven regions of imaging. The risk assessment was inspired by standards for radiofrequency fields and by scientific studies for gradient fields. Three risk tiers were defined for the radiofrequency, in terms of whole-body and local SAR averages, and for GC fields, in terms of temperature elevation. Only 50 out of 588 scenarios require some caution to be managed. Results showed that the whole-body SAR is not a self-reliant safety parameter for patients with metallic implants. The proposed numerical procedure can be easily extended to any other scenario, including the use of detailed anatomical models.

A contribution to MRI safety testing related to gradient-induced heating of medical devices

PURPOSE

To theoretically investigate the feasibility of a novel procedure for testing the MRI gradient-induced heating of medical devices and translating the results into clinical practice.

METHODS

The concept of index of stress is introduced by decoupling the time waveform characteristics of the gradient field signals from the field spatial distribution within an MRI scanner. This index is also extended to consider the anisotropy of complex bulky metallic implants. Merits and drawbacks of the proposed index of stress are investigated through virtual experiments. In particular, the values of the index of stress evaluated for realistic orthopedic implants placed within an ASTM phantom are compared with accurate heating simulations performed with 2 anatomic body models (a man and a woman) implanted through a virtual surgery procedure.

RESULTS

The manipulation of the proposed index of stress allows to identify regions within the MRI bore where the implant could affect the safety of the examinations. Furthermore, the conducted analysis shows that the power dissipated into the implant by the induced eddy currents is a dosimetric quantity that estimates well the maximum temperature increase in the tissues surrounding the implant.

CONCLUSION

The results support the adoption of an anisotropic index of stress to regulate the gradient-induced heating of geometrically complex implants. They also pave the way for a laboratory characterization of the implants based on electrical measurements, rather than on thermal measurements. The next step will be to set up a standardized experimental procedure to evaluate the index of stress associated with an implant.

Parameters Affecting Worst-Case Gradient-Field Heating of Passive Conductive Implants

BACKGROUND

Testing MRI gradient-induced heating of implanted medical devices is required by regulatory organizations and others. A gradient heating test of the ISO 10974 Technical Specification (TS) for active implants was adopted for this study of passive hip implants. All but one previous study of hip implants used nonuniform gradient exposure fields in clinical scanners and reported heating of less than 5 °C. This present study adapted methods of the TS, addressing the unmet need for identifying worst-case heating via exposures to uniform gradient fields.

PURPOSE

To identify gradient-field parameters affecting maximum heating in vitro for a hip implant and a cylindrical titanium disk.

STUDY TYPE

Computational simulations and experimental validation of induced heating.

PHANTOM

Tissue-simulating gel.

FIELD STRENGTH

42 T/s RMS, sinusoidal, continuous B fields with high spatial uniformity ASSESSMENT: Hip implant heating at 1-10 kHz, via computational modeling, validated by limited point measurements. Experimental measurements of exposures of an implant at 42 T/s for 4, 6, and 9 kHz, analyzed at 50, 100, and 150 seconds.

STATISTICAL TESTS

One sample student's t-test to assess difference between computational and experimental results. Experimental vs. computational results were not significantly different (p < 0.05).

RESULTS

Maximum simulated temperature rise (10-minute exposure) was 10 °C at 1 kHz and 0.66 °C at 10 kHz. The ratio of the rise for 21 T/s vs. 42 T/s RMS was 4, after stabilizing at 50 seconds (dB/dt ratio squared).

DATA CONCLUSIONS

Heating of an implant is proportional to the frequency of the B field and the implant's cross-sectional area and is greater for a thickness on the order of its skin depth. Testing with lower values of dB/dt RMS with lower cost amplifiers enables prediction of heating at higher values for dB/dt squared (per ISO TS) with identical frequency components and waveforms, once thermal equilibrium occurs.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

A fast tool for the parametric analysis of human body exposed to LF electromagnetic fields in biomedical applications

A numerical procedure for analyzing electromagnetic (EM) fields interactions with biological tissues is presented. The proposed approach aims at drastically reducing the computational burden required by the repeated solution of large scale problems involving the interaction of the human body with EM fields, such as in the study of the time evolution of EM fields, uncertainty quantification, and inverse problems. The proposed volume integral equation (VIE), focused on low frequency applications, is a system of integral equations in terms of current density and scalar potential in the biological tissues excited by EM fields and/or electrodes connected to the human body. The proposed formulation requires the voxelization of the human body and takes advantage of the regularity of such discretization by speeding-up the computational procedure. Moreover, it exploits recent advancements in the solution of VIE by means of iterative preconditioned solvers and ad hoc parametric Model Order Reduction techniques. The efficiency of the proposed tool is demonstrated by applying it to a couple of realistic model problems: the assessment of the peripheral nerve stimulation, performed in terms of evaluation of the induced electric field, due to the gradient coils of a magnetic resonance imaging scanner during a clinical examination and the assessment of the exposure to environmental fields at 50 Hz of live-line workers with uncertain properties of the biological tissues. Thanks to the proposed method, uncertainty quantification analyses and time domain simulations are possible even for large scale problems and they can be performed on standard computers and reasonable computation time. Sample implementation of the method is made publicly available at https://github.com/UniPD-DII-ETCOMP/BioMOR.

Gradient coil and radiofrequency induced heating of orthopaedic implants in MRI: influencing factors

Patients with implanted orthopaedic devices represent a growing number of subjects undergoing magnetic resonance imaging (MRI) scans each year. MRI safety labelling is required for all implants under the EU Medical Device Regulations to ensure regulatory compliance, with each device assessed through standardised testing procedures. In this paper, we employ parametric studies to assess a range of clinically relevant factors that cause tissue heating, performing simulations with both radiofrequency (RF) and gradient coil (GC) switching fields, the latter of which is often overlooked in the literature. A series of 'worst-case' scenarios for both types of excitation field is discussed. In the case of GC fields, large volume implants and large plate areas with the field orientated perpendicular to the plane cause the highest heating levels, along with sequences with high rates of field switching. Implant heating from RF fields is driven primarily from the 'antenna effect', with thin, linear implants of resonant length resulting in the highest temperature rises. In this work, we show that simplifications may be made to the field sequence and in some cases the device geometry without significantly compromising the accuracy of the simulation results, enabling the possibility for generic estimates of the implant heating for orthopaedic device manufacturers and opportunities to simplify the safety compliance process.

In silico assessment of collateral eddy current heating in biocompatible implants subjected to magnetic hyperthermia treatments

Purpose: Bearing partially or fully metallic passive implants represents an exclusion criterion for patients undergoing a magnetic hyperthermia procedure, but there are no specific studies backing this restrictive decision. This work assesses how the secondary magnetic field generated at the surface of two common types of prostheses affects the safety and efficiency of magnetic hyperthermia treatments of localized tumors. The paper also proposes the combination of a multi-criteria decision analysis and a graphical representation of calculated data as an initial screening during the preclinical risk assessment for each patient.Materials and methods: Heating of a hip joint and a dental implant during the treatment of prostate, colorectal and head and neck tumors have been assessed considering different external field conditions and exposure times. The Maxwell equations including the secondary field produced by metallic prostheses have been solved numerically in a discretized computable human model. The heat exchange problem has been solved through a modified version of the Pennes' bioheat equation assuming a temperature dependency of blood perfusion and metabolic heat, i.e. thermorregulation. The degree of risk has been assessed using a risk index with parameters coming from custom graphs plotting the specific absorption rate (SAR) vs temperature increase, and coefficients derived from a multi-criteria decision analysis performed following the MACBETH approach.Results: The comparison of two common biomaterials for passive implants - Ti6Al4V and CoCrMo - shows that both specific absorption rate (SAR) and local temperature increase are found to be higher for the hip prosthesis made by Ti6Al4V despite its lower electrical and thermal conductivity. By tracking the time evolution of temperature upon field application, it has been established that there is a 30 s delay between the time point for which the thermal equilibrium is reached at prostheses and tissues. Likewise, damage may appear in those tissues adjacent to the prostheses at initial stages of treatment, since recommended thermal thresholds are soon surpassed for higher field intensities. However, it has also been found that under some operational conditions the typical safety rule of staying below or attain a maximum temperature increase or SAR value is met.Conclusion: The current exclusion criterion for implant-bearing patients in magnetic hyperthermia should be revised, since it may be too restrictive for a range of the typical field conditions used. Systematic in silico treatment planning using the proposed methodology after a well-focused diagnostic procedure can aid the clinical staff to find the appropriate limits for a safe treatment window.

Heating of hip joint implants in MRI: The combined effect of RF and switched-gradient fields

PURPOSE

To investigate how the simultaneous exposure to gradient and RF fields affects the temperature rise in patients with a metallic hip prosthesis during an MRI session.

METHODS

In silico analysis was performed with an anatomically realistic human model with CoCrMo hip implant in 12 imaging positions. The analysis was performed at 1.5 T and 3 T, considering four clinical sequences: turbo spin-echo, EPI, gradient-echo, and true fast imaging sequence with steady precession. The exposure to gradient and RF fields was evaluated separately and superposed, by adopting an ad hoc computational algorithm. Temperature increase within the body, rather than specific absorption rate, was used as a safety metric.

RESULTS

With the exception of gradient-echo, all investigated sequences produced temperature increases higher than 1 K after 360 seconds, at least for one body position. In general, RF-induced heating dominates the turbo spin-echo sequence, whereas gradient-induced heating prevails with EPI; the situation with fast imaging sequence with steady precession is more diversified. The RF effects are enhanced when the implant is within the RF coil, whereas the effects of gradient fields are maximized if the prosthesis is outside the imaging region. Cases for which temperature-increase thresholds were exceeded were identified, together with the corresponding amount of tissue mass involved and the exposure time needed to reach these limits.

CONCLUSION

The analysis confirms that risky situations may occur when a patient carrying a hip implant undergoes an MRI exam and that, in some cases, the gradient field heating may be significant. In general, exclusion criteria only based on whole-body specific absorption rate may not be sufficient to ensure patients' safety.

Heating of Hip Arthroplasty Implants During Metal Artifact Reduction MRI at 1.5- and 3.0-T Field Strengths

OBJECTIVES

The aim of this study was to quantify the spatial temperature rises that occur during 1.5- and 3.0-T magnetic resonance imaging (MRI) of different types of hip arthroplasty implants using different metal artifact reduction techniques.

MATERIALS AND METHODS

Using a prospective in vitro study design, we evaluated the spatial temperature rises of 4 different total hip arthroplasty constructs using clinical metal artifact reduction techniques including high-bandwidth turbo spin echo (HBW-TSE), slice encoding for metal artifact correction (SEMAC), and compressed sensing SEMAC at 1.5 and 3.0 T. Each MRI protocol included 6 pulse sequences, with imaging planes, parameters, and coverage identical to those in patients. Implants were immersed in standard American Society for Testing and Materials phantoms, and fiber optic sensors were used for temperature measurement. Effects of field strength, radiofrequency pulse polarization at 3.0 T, pulse protocol, and gradient coil switching on heating were assessed using nonparametric Friedman and Wilcoxon signed-rank tests.

RESULTS

Across all implant constructs and MRI protocols, the maximum heating at any single point reached 13.1°C at 1.5 T and 1.9°C at 3.0 T. The temperature rises at 3.0 T were similar to that of background in the absence of implants (P = 1). Higher temperature rises occurred at 1.5 T compared with 3.0 T (P < 0.0001), and circular compared with elliptical radiofrequency pulse polarization (P < 0.0001). Compressed sensing SEMAC generated equal or lower degrees of heating compared with HBW-TSE at both field strengths (P < 0.0001).

CONCLUSIONS

Magnetic resonance imaging of commonly used total hip arthroplasty implants is associated with variable degrees of periprosthetic tissue heating. In the absence of any perfusion effects, the maximum temperature rises fall within the physiological range at 3.0 T and within the supraphysiologic range at 1.5 T. However, with the simulation of tissue perfusion effects, the heating at 1.5 T also reduces to the upper physiologic range. Compressed sensing SEMAC metal artifact reduction MRI is not associated with higher degrees of heating than the HBW-TSE technique.