In this article, we discuss numerical aspects related to the accuracy and the computational efficiency of numerical dosimetric simulations, performed in the context of human exposure to static inductive charging systems of electric vehicles. Two alternative numerical methods based on electric vector potential and electric scalar potential formulations, respectively, are here considered for the electric field computation in highly detailed anatomical human models. The results obtained by the numerical implementation of both approaches are discussed in terms of compliance assessment with ICNIRP guidelines limits for human exposure to electromagnetic fields. In particular, different strategies for smoothing localized unphysical outliers are compared, including novel techniques based on statistical considerations. The outlier removal is particularly relevant when comparison with basic restrictions is required to define the safety of electromagnetic fields exposure. The analysis demonstrates that it is not possible to derive general conclusions about the most robust method for dosimetric solutions. Nevertheless, the combined use of both formulations, together with the use of an algorithm for outliers removal based on a statistical approach, allows to determine final results to be compared with reference limits with a significant level of reliability.

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.

Magnetic resonance imaging scanners can provide, by B1-mapping techniques, measurements of the transmit sensitivity magnitude, whereas the transmit sensitivity phase can only be approximated under symmetry assumptions of the coil and the sample. Nonetheless, up to now many techniques have been proposed to retrieve electric properties maps at radiofrequency by elaborating the complex transmit sensitivity. These techniques are usually referred to as magnetic resonance-based electric properties tomography. In order to perform imaging of the elecric properties also in cases which do not fulfill the assumptions for transmit sensitivity phase estimation, a variation of the phaseless contrast source inversion technique is proposed in this paper to retrieve the electric properties by elaborating only the magnitude of the transmit sensitivity. The novel technique is described in a generalised framework based on a functional viewpoint and it is specialised to deal with two-dimensional boundary-free problems. Moreover, a phaseless back propagation solution is proposed as initial guess. The proposed technique and initial guess have been tested on realistic two-dimensional model problems, which show that the numerical minimisation converges towards the global minimum when multiple linearly independent measurements are provided and that it is quite robust with respect to noisy inputs.

Electric properties tomography (EPT) is a novel magnetic resonance imaging–based method to estimate non-invasively the distribution of the electric properties in the human body. In this paper, EPTlib, an open-source extensible C++ library collecting ready-to-use algorithms for electric properties tomography, is presented. Currently, EPTlib implements three techniques, named Helmholtz-EPT, convection-reaction-EPT and gradient-EPT, whose derivation and implementation is deeply discussed. Moreover, the configuration files needed by the terminal application included in EPTlib to apply the implemented techniques are outlined. The three techniques are applied to a couple of model problems in order to highlight their main features and the effects of the tunable parameters.

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.

Contrast source inversion electric properties tomography (CSI-EPT) is a recently developed technique for the electric properties tomography that recovers the electric properties distribution starting from measurements performed by magnetic resonance imaging scanners. This method is an optimal control approach based on the contrast source inversion technique, which distinguishes itself from other electric properties tomography techniques for its capability to recover also the local specific absorption rate distribution, essential for online dosimetry. Up to now, CSI-EPT has only been described in terms of integral equations, limiting its applicability to homogeneous unbounded background. In order to extend the method to the presence of a shield in the domain-as in the recurring case of shielded radio frequency coils-a more general formulation of CSI-EPT, based on a functional viewpoint, is introduced here. Two different implementations of CSI-EPT are proposed for a 2-D transverse magnetic model problem, one dealing with an unbounded domain and one considering the presence of a perfectly conductive shield. The two implementations are applied on the same virtual measurements obtained by numerically simulating a shielded radio frequency coil. The results are compared in terms of both electric properties recovery and local specific absorption rate estimate, in order to investigate the requirement of an accurate modeling of the underlying physical problem.