Abstract:The numerical solution of problems in nonlinear magnetostatics is typically based on a variational formulation in terms of magnetic potentials, the discretization by finite elements, and iterative solvers like the Newton method. The vector potential approach aims at minimizing a certain energy functional and, in three dimensions, requires the use of edge elements and appropriate gauging conditions. The scalar potential approach, on the other hand, seeks to maximize the negative coenergy and can be realized by standard Lagrange finite elements, thus reducing the number of degrees of freedom and simplifying the implementation. The number of Newton iterations required to solve the governing nonlinear system, however, has been observed to be usually higher than for the vector potential formulation. In this paper, we propose a method that combines the advantages of both approaches, i.e., it requires as few Newton iterations as the vector potential formulation while involving the magnetic scalar potential as the primary unknown. We discuss the variational background of the method, its well-posedness, and its efficient implementation. Numerical examples are presented for illustration of the accuracy and the gain in efficiency compared to other approaches.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in the numerical solution of nonlinear magnetostatic problems, how to combine the advantages of the scalar potential and vector potential methods to propose an improved simplified scalar potential method, thereby reducing the number of Newton iterations and improving computational efficiency. Specifically, the two existing main methods - the vector potential method and the scalar potential method - each have their own advantages and disadvantages: - **Vector potential method**: Solved by minimizing the energy functional. When applied to three - dimensional problems, edge elements need to be used and appropriate gauge conditions need to be imposed. Although this method usually requires a smaller number of Newton iterations, it is complex to implement and has more degrees of freedom. - **Scalar potential method**: Solved by maximizing the negative co - energy. Standard Lagrange finite elements can be used, with fewer degrees of freedom and simple implementation. However, this method usually requires more Newton iterations. For this reason, the author proposes a new method, aiming at: 1. **Reducing the number of Newton iterations**: Achieving the same number of iterations as the vector potential method. 2. **Maintaining a low number of degrees of freedom**: Similar to the scalar potential method, using standard Lagrange finite elements. 3. **Simplifying the implementation**: Avoiding complex gauge conditions. By combining the advantages of these two methods, the author hopes to significantly improve the efficiency of solving nonlinear magnetostatic problems while ensuring computational accuracy. ### Summary of mathematical formulas The main formulas involved in the paper are as follows: - The relationship between the magnetic field intensity \( \mathbf{h} \) and the magnetic induction intensity \( \mathbf{b} \): \[ \mathbf{h} = \frac{\partial w(\mathbf{b})}{\partial \mathbf{b}} \] where \( w(\mathbf{b}) \) is the magnetic energy density. - The relationship between the scalar potential \( \psi \) and the source field \( \mathbf{h}_s \): \[ \mathbf{h} = \mathbf{h}_s - \nabla \psi \] - The nonlinear scalar potential formula: \[ \nabla \cdot \left( \frac{\partial w^*(\mathbf{h}_s - \nabla \psi)}{\partial \mathbf{h}} \right) = 0 \quad \text{in } \Omega \] \[ \mathbf{b} \cdot \mathbf{n} = 0 \quad \text{on } \partial \Omega \] where \( w^*(\mathbf{h}) \) is the co - energy density. - The discretization problem in weak form: \[ \langle \frac{\partial w(\mathbf{b}_h), \mathbf{b}'_h \rangle_h + \langle \nabla \psi'_h, \mathbf{b}'_h \rangle_h = \langle \mathbf{h}_s, \mathbf{b}'_h \rangle_h \quad \forall \mathbf{b}'_h \in Q_h \] \[ \langle \mathbf{b}_h, \nabla \psi_h \rangle_h = 0 \quad \forall \psi'_h \in V_h \] Through these formulas, the author proposes an improved scalar potential method, which can reduce the number of Newton iterations and maintain a low number of degrees of freedom and a simple implementation method when solving nonlinear magnetostatic problems.

A reduced scalar potential approach for magnetostatics avoiding the coenergy

Stability and convergence of the penalty formulation for nonlinear magnetostatics

An efficient mixed finite element method for nonlinear magnetostatics and quasistatics

On the use of mixed potential formulation for finite-element analysis of large-scale magnetization problems with large memory demand

Magnetic Field Analysis of Electric Current in Coplanar Coils Using Equivalent Permanent Magnetization and Magnetic Scalar Potential

Magnetic Scalar and Vector Potential Approaches in Dynamic Response Studies of Magneto-Electro-Elastic Plates

On the convergence of higher order finite element methods for nonlinear magnetostatics

Comparison of magnetic vector and total scalar potential formulations for finite-element modeling of dipole magnet with COMSOL Multiphysics

Efficient Reduced Magnetic Vector Potential Formulation for the Magnetic Field Simulation of Accelerator Magnets

A combined vector potential-scalar potential method for FE computation of 3D magnetic fields in electrical devices with iron cores

Global convergence of iterative solvers for problems of nonlinear magnetostatics

Variational methods for solving numerically magnetostatic systems

On some energy‐based variational principles in non‐dissipative magneto‐mechanics using a vector potential approach

Physics aware machine learning for micromagnetic energy minimization: recent algorithmic developments

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

On nonlinear magnetic field solvers using local Quasi-Newton updates

A Material Point Method for Nonlinearly Magnetized Materials

Decoupled, linear, unconditionally energy stable and charge-conservative finite element method for an inductionless magnetohydrodynamic phase-field model

Magnetostatics and micromagnetics with physics informed neural networks

An Adaptive Edge Element Method and Its Convergence for a Saddle‐Point Problem from Magnetostatics

Error analysis of the linearized Crank-Nicolson FEM for the incompressible vector potential magnetohydrodynamic system