An Integral-based Approach for the Vector Potential in Smoothed Particle Magnetohydrodynamics

Terrence S. Tricco,Daniel J. Price

2023-06-27

Abstract:A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge choice. This new discretisation is tested using the Orszag-Tang MHD vortex in a 3D configuration. The corresponding conservative equations of motion are derived, but are not found to solve the MHD equations in the continuum limit. Tests are performed using a hybrid approach instead, whereby the equations of motion based on the magnetic field instead of vector potential are used. Test results experience the same numerical instability as with the Price (2010) formulation. We conclude that this new formulation is non-viable.

Instrumentation and Methods for Astrophysics,Plasma Physics

What problem does this paper attempt to address?

This paper attempts to solve the numerical instability problem in the time evolution using vector potential in smoothed particle magnetohydrodynamics (SPMHD). Specifically, the paper proposes a new volume - integral - based method to discretize the time - evolution equation of the vector potential, hoping to overcome the numerical instability problems existing in previous methods (such as the formula proposed by Price in 2010). However, the test results show that although the new method is theoretically innovative, it still encounters the same numerical instability problems as Price's 2010 method in practical applications, especially more obvious at high resolutions.

### Main contributions of the paper
1. **New discretization method**: By considering the integral form of the induction equation, the paper proposes a new discretization method. In particular, the "inverse convection" term is improved, using the interpolation of the velocities of neighboring particles instead of directly using the particle velocities.
2. **Galilean invariance**: Through an appropriate gauge choice, the new discretization method can achieve Galilean invariance.
3. **Testing and verification**: The paper tests the new method using 3D Orszag - Tang vortices and compares it with Price's 2010 method and the method of directly evolving the magnetic field.

### Main findings
- **Numerical instability**: Although the new method is theoretically improved, it still shows the same numerical instability as Price's 2010 method in actual tests. This instability is especially obvious at high resolutions and usually begins to appear after particles break away from the lattice structure.
- **Energy growth**: The magnetic field energy in the new method starts to grow exponentially at certain times, indicating that there are serious numerical problems in the method when dealing with the evolution of the magnetic field.
- **Conservative equations**: Although an attempt was made to use conservative equations to solve the problem, these equations still have problems such as negative pressure and the direct second - order derivative noise of the kernel function. These problems are not caused by the vector potential evolution equation itself, but originate from the discretization of the equation of motion.

### Conclusion
The main conclusion drawn by the paper is that although the new volume - integral - based discretization method is theoretically innovative, it still cannot overcome the numerical instability problem in practical applications. Therefore, this method is not feasible in its current form. Future research needs to continue to search for an SPMHD method that can both ensure the divergence - free of the magnetic field and be numerically stable.