Arbitrary order spline representation of cohomology generators for isogeometric analysis of eddy current problems

Bernard Kapidani,Melina Merkel,Sebastian Schöps,Rafael Vázquez

DOI: https://doi.org/10.48550/arXiv.2310.14313

2023-10-22

Abstract:The eddy current problem has many relevant practical applications in science, ranging from non-destructive testing to magnetic confinement of plasma in fusion reactors. It arises when electrical conductors are immersed in an external time-varying magnetic field operating at frequencies for which electromagnetic wave propagation effects can be neglected. Popular formulations of the eddy current problem either use the magnetic vector potential or the magnetic scalar potential. They have individual advantages and disadvantages. One challenge is related to differential geometry: Scalar potential based formulations run into trouble when conductors are present in non-trivial topology, as approximation spaces must be then augmented with generators of the first cohomology group of the non-conducting domain. For all existing algorithms based on lowest order methods it is assumed that the extension of the graph-based algorithms to high-order approximations requires hierarchical bases for the curl-conforming discrete spaces. However, building on insight on de Rham complexes approximation with splines, we will show in the present submission that the hierarchical basis condition is not necessary. Algorithms based on spanning tree techniques can instead be adapted to work on an underlying hexahedral mesh arising from isomorphisms between spline spaces of differential forms and de Rham complexes on an auxiliary control mesh.

Numerical Analysis,Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

What this paper attempts to solve is how to deal with the challenges brought by non - trivial topological structures when solving eddy - current problems within the framework of Isogeometric Analysis (IGA). Specifically, when a conductor has a complex topological structure, traditional discretization methods need to expand the generators of the first homology group in non - conductive regions to ensure the accuracy of numerical solutions.

### Main problems
1. **Background of eddy - current problems**
- Eddy - current problems occur when a conductor is immersed in an external time - varying magnetic field and the working frequency is low, so that the electromagnetic wave propagation effect can be ignored.
- Traditional methods usually use the magnetic vector potential \( \mathbf{A} \) or the magnetic scalar potential \( \phi \) to describe eddy - current problems.
2. **Limitations of existing methods**
- When a conductor has a non - trivial topological structure (such as the existence of holes or complex connections), methods based on the magnetic scalar potential need to expand the discrete space by adding the generators of the first homology group.
- Most existing algorithms are based on the lowest - order method and assume that high - order approximations require hierarchical basis functions, which increases the computational complexity and implementation difficulty.
3. **Research objectives**
- This paper proposes a new method that uses the spline representation method to approximate the generators of the homology group, thus avoiding the need for hierarchical basis functions.
- By applying the spanning - tree technique on the auxiliary hexahedral mesh, high - order spline spaces can be efficiently computed while maintaining topological properties.

### Key points of the solution
- **Spline space and de Rham complex**: The paper shows how to approximate the de Rham complex of differential forms through the construction of spline spaces, thereby simplifying the handling of non - trivial topological structures.
- **Spanning - tree technique**: Through the spanning - tree technique, a low - order finite - element space can be defined on the hexahedral control mesh and applied to the high - order spline space through the exchange isomorphism mapping.
- **Maintaining topological properties**: Through the exchange isomorphism mapping, the topological properties between the control mesh and the high - order spline space are ensured to be consistent, thus guaranteeing the accuracy of the numerical solution.

### Conclusion
The main contribution of this paper is to provide a new method that enables efficient high - order spline approximation when dealing with eddy - current problems with complex topological structures without relying on hierarchical basis functions. This provides new ideas and tools for the application of Isogeometric Analysis in electromagnetic field simulation.