Stability and convergence of the penalty formulation for nonlinear magnetostatics

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Egger, Herbert, Engertsberger, Felix, Roppert, Klaus
Natura: Preprint
Pubblicazione: 2024
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866913286159597568
author Egger, Herbert
Engertsberger, Felix
Roppert, Klaus
author_facet Egger, Herbert
Engertsberger, Felix
Roppert, Klaus
contents The magnetostatic field distribution in a nonlinear medium amounts to the unique minimizer of the magnetic coenergy over all fields that can be generated by the same current. This is a nonlinear saddlepoint problem whose numerical solution can in principle be achieved by mixed finite element methods and appropriate nonlinear solvers. The saddlepoint structure, however, makes the solution cumbersome. A remedy is to split the magnetic field into a known source field and the gradient of a scalar potential which is governed by a convex minimization problem. The penalty approach avoids the use of artificial potentials and Lagrange multipliers and leads to an unconstrained convex minimization problem involving a large parameter. We provide a rigorous justification of the penalty approach by deriving error estimates for the approximation due to penalization. We further highlight the close connections to the Lagrange-multiplier and scalar potential approach. The theoretical results are illustrated by numerical tests for a typical benchmark problem
format Preprint
id arxiv_https___arxiv_org_abs_2403_18285
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Stability and convergence of the penalty formulation for nonlinear magnetostatics
Egger, Herbert
Engertsberger, Felix
Roppert, Klaus
Numerical Analysis
Optimization and Control
The magnetostatic field distribution in a nonlinear medium amounts to the unique minimizer of the magnetic coenergy over all fields that can be generated by the same current. This is a nonlinear saddlepoint problem whose numerical solution can in principle be achieved by mixed finite element methods and appropriate nonlinear solvers. The saddlepoint structure, however, makes the solution cumbersome. A remedy is to split the magnetic field into a known source field and the gradient of a scalar potential which is governed by a convex minimization problem. The penalty approach avoids the use of artificial potentials and Lagrange multipliers and leads to an unconstrained convex minimization problem involving a large parameter. We provide a rigorous justification of the penalty approach by deriving error estimates for the approximation due to penalization. We further highlight the close connections to the Lagrange-multiplier and scalar potential approach. The theoretical results are illustrated by numerical tests for a typical benchmark problem
title Stability and convergence of the penalty formulation for nonlinear magnetostatics
topic Numerical Analysis
Optimization and Control
url https://arxiv.org/abs/2403.18285