Numerical generation of vector potentials from specified magnetic fields

Abstract

Many codes have been developed to study highly relativistic, magnetized flows around and inside compact objects. Depending on the adopted formalism, some of these codes evolve the vector potential $\mathbf{A}$, and others evolve the magnetic field $\mathbf{B}=\nabla \times \mathbf{A}$ directly. Given that these codes possess unique strengths, sometimes it is desirable to start a simulation using a code that evolves $\mathbf{B}$ and complete it using a code that evolves $\mathbf{A}$. Thus transferring the data from one code to another would require an inverse curl algorithm. This paper describes two new inverse curl techniques in the context of Cartesian numerical grids: a cell-by-cell method, which scales approximately linearly with the numerical grid, and a global linear algebra approach, which has worse scaling properties but is generally more robust (e.g., in the context of a magnetic field possessing some nonzero divergence). We demonstrate these algorithms successfully generate smooth vector potential configurations in challenging special and general relativistic contexts.

20 pages, 9 figures

References in corpus (7)

• HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
• GRHydro: A new open source general-relativistic magnetohydrodynamics code for the Einstein Toolkit
• IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes
• Magnetic Helicity and Large Scale Magnetic Fields: A Primer
• Variational Integrators for Inertial Magnetohydrodynamics
• Inverse Vector Operators
• Method of Relative Magnetic Helicity Computation II: Boundary Conditions for the Vector Potentials

Cited by in corpus (6)

• H-AMR FORGE'd in FIRE I: Magnetic state transitions, jet launching and radiative emission in super-Eddington, highly magnetized quasar disks formed from cosmological initial conditions
• Curl-Flow: Boundary-Respecting Pointwise Incompressible Velocity Interpolation for Grid-Based Fluids
• Discrete Vector Calculus and Helmholtz Hodge Decomposition for Classical Finite Difference Summation by Parts Operators
• Inverting the discrete curl operator: a novel graph algorithm to find a vector potential of a given vector field
• Turbulence in Magnetised Neutron Stars
• BHaHAHA: A Fast, Robust Apparent Horizon Finder Library for Numerical Relativity