A structure preserving H-curl algebraic multigrid method for the eddy current equations
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866912422342688768 |
|---|---|
| author | Tuminaro, Raymond Glusa, Christian |
| author_facet | Tuminaro, Raymond Glusa, Christian |
| contents | A new algebraic multigrid method (AMG) is presented for solving the linear systems associated with the eddy current approximation to the Maxwell equations. This AMG method extends an idea proposed by Reitzinger and Schoberl. The main feature of the Reitzinger and Schoberl algorithm (RSAMG) is that it maintains null-space properties of the Curl-Curl operator throughout all levels of the AMG hierarchy. It does this by enforcing a commuting relationship involving grid transfers and the discrete gradient operator. This null-space preservation property is critical to the algorithm's success, however enforcing this commuting relationship is non-trivial except in the special case where one leverages a piece-wise constant nodal interpolation operator. For this reason, mesh independent convergence rates are generally not observed for RSAMG due to its reliance on sub-optimal piece-wise constant interpolation. We present a new AMG algorithm that enforces the same commuting relationship. The main advance is that the new structure preserving H-curl algorithm (SpHcurlAMG) does not rely on piece-wise constant interpolation and can leverage fairly general and more sophisticated nodal interpolation operators. The key idea is to employ energy minimization AMG (EAMG) to construct edge interpolation grid transfers and to enforce the commuting relationship by embedding it as constraints within an EAMG procedure. While it might appear that solving such a constrained energy minimization is costly, we illustrate how this is not the case in our context. Numerical results are then given demonstrating mesh independent convergence over a range of test problems. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2506_08284 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A structure preserving H-curl algebraic multigrid method for the eddy current equations Tuminaro, Raymond Glusa, Christian Numerical Analysis A new algebraic multigrid method (AMG) is presented for solving the linear systems associated with the eddy current approximation to the Maxwell equations. This AMG method extends an idea proposed by Reitzinger and Schoberl. The main feature of the Reitzinger and Schoberl algorithm (RSAMG) is that it maintains null-space properties of the Curl-Curl operator throughout all levels of the AMG hierarchy. It does this by enforcing a commuting relationship involving grid transfers and the discrete gradient operator. This null-space preservation property is critical to the algorithm's success, however enforcing this commuting relationship is non-trivial except in the special case where one leverages a piece-wise constant nodal interpolation operator. For this reason, mesh independent convergence rates are generally not observed for RSAMG due to its reliance on sub-optimal piece-wise constant interpolation. We present a new AMG algorithm that enforces the same commuting relationship. The main advance is that the new structure preserving H-curl algorithm (SpHcurlAMG) does not rely on piece-wise constant interpolation and can leverage fairly general and more sophisticated nodal interpolation operators. The key idea is to employ energy minimization AMG (EAMG) to construct edge interpolation grid transfers and to enforce the commuting relationship by embedding it as constraints within an EAMG procedure. While it might appear that solving such a constrained energy minimization is costly, we illustrate how this is not the case in our context. Numerical results are then given demonstrating mesh independent convergence over a range of test problems. |
| title | A structure preserving H-curl algebraic multigrid method for the eddy current equations |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2506.08284 |