A Finite Volume scheme for the solution of discontinuous magnetic field distributions on non-orthogonal meshes

Abstract

We present a Finite Volume formulation for determining discontinuous distributions of magnetic fields within non-orthogonal and non-uniform meshes. The numerical approach is based on the discretization of the vector potential variant of the equations governing static magnetic field distribution in magnetized, permeable and current carrying media. After outlining the derivation of the magnetostatic balance equations and its associated boundary conditions, we propose a cell–centered Finite Volume framework for spatial discretization and a Block Gauss–Seidel multi-region scheme for solution. We discuss the structure of the solver, emphasizing its effectiveness and addressing stabilization and correction techniques to enhance computational robustness. We validate the accuracy and efficacy of the approach through numerical experiments and comparisons with the Finite Element method.