Magnetic shielding of an MRI magnet using the Thin plate boundary condition

Magnetic flux density in shielded room wall and in the vicinity of the MRI magnet

Magnetic flux density in shielded room wall and in the vicinity of the MRI magnet

Many electromagnetic simulations include ferrous geometry with extreme aspect ratios. A good example of this is modelling the magnetic shielding in a calculation of the stray field from a shielded room surrounding a representative superconducting coil set of an MRI magnet. It is an important legislative requirement that public areas are only subjected to DC field levels below a limit that is deemed to be not harmful to humans (usually 5 gauss).

However, even though active shielding is used in modern MRI magnets (superconducting coils that oppose the main field of the magnet), most installations also require a steel shield attached to the structural walls, and possibly to the floor and ceiling. The shield designer needs to ensure that it is effective, but also wishes to minimize the weight (and hence cost). The MRI magnet designer is concerned that the magnetization of the shield does not perturb the homogeneous field at the centre of the magnet to such an extent that it cannot be corrected on-site using shims.

Opera-3d already has useful features, such as mesh layering and hexahedral / prism elements, that enable these structures to be meshed with high aspect ratio volume elements. But even these tools reach a limit when several plates meet – such as around windows and doors. Either the junction area becomes considerably over-meshed, increasing throughput time for simulations without improving accuracy, or the model building requires considerable effort to create “bevelled” edges where the plates join.

In Opera version 18R1, a method to improve this has been introduced where the three dimensional volume structure of the plate is replaced with a two dimensional surface representation, accompanied by a Thin plate boundary condition [1]. Hence, any surfaces can meet without introducing extra elements or requiring special modelling.

The data requirements for the boundary condition are very simple. A material label is used to specify from which magnetic material the plate is constructed, and the thickness of the plate is given. The boundary condition can use either magnetically linear or non-linear materials.

Although the Thin plate boundary condition is applied to a two dimensional surface, the field on the surface is multivalued. For example, if the magnetic flux enters the plate normally and then largely turns to flow in the plate, the normal component will be discontinuous across the plate.

5 gauss iso-surface

5 gauss iso-surface

As can be seen, using 4 cm thick mild steel shielding plates, fields exceeding the 5 gauss limit are found above and below the shielded room, at floor and ceiling level just outside the shield and near the viewing window. Consequently, the shield needs to be made more effective, ceiling and floor plates may also be needed, if the installation is in a multi-story building, and the size of the viewing window reviewed. Increasing the shield effectiveness can be accomplished by increasing thickness or by using a better quality magnetic steel. This is simply achieved with the new Thin Plate boundary condition, since, as previously mentioned, its parameters include both thickness and material characteristics.

Reference

[1] Christopher S. Biddlecombe and Christopher P. Riley, “Improvements to finite element meshing for magnetic signature simulations”, presented at MARELEC 2015 conference, Philadelphia, PA, USA, June 2015