Abstract:We extend the co-scaling formalism of Habegger & Heitsch (2021) implemented in Athena++ to magneto-hydrodynamics. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The formalism is fully consistent with the upwind constrained transport method implemented in Athena++ and is accurate to 2nd order in space. Applying our implementation to standard magneto-hydrodynamic test cases leads to improved results and higher efficiency, compared to the fixed-grid solutions.
What problem does this paper attempt to address?
This paper attempts to solve the problem of how to effectively handle large - scale spatial ranges and complex flow symmetries in astrophysical simulations, especially when dealing with magnetohydrodynamics (MHD) problems. Specifically, the authors extended the co - scaling grid method previously implemented in the Athena++ code to adapt to magnetohydrodynamics simulations. The following are the main problems in the paper and their solutions:
### Research Background and Problems
1. **Numerical Simulation of Large - scale Spatial Ranges**
- When simulating astrophysical phenomena such as supernova remnants and kilonova remnants, a large range of spatial scales needs to be covered.
- Traditional fixed - grid methods are less efficient in handling these problems and may lead to loss of accuracy.
2. **Symmetries and Complex Flows in Magnetohydrodynamics**
- Symmetries such as expansion, contraction, and centroid motion are often involved in astrophysical problems.
- The presence of the magnetic field breaks the radial symmetry, causing materials to "sweep" towards the equatorial region and affecting the evolution of the blast wave.
- The change of the magnetic field is crucial for understanding molecular clouds and star formation.
### Solutions
1. **Co - scaling Grid Method**
- The authors extended the co - scaling grid method proposed by Habegger & Heitsch (2021) to make it applicable to magnetohydrodynamics.
- This method can handle large - scale spatial ranges and complex flow symmetries more efficiently by introducing time - varying grids.
2. **Improved Magnetic Flux Update**
- To ensure the accurate update of the magnetic flux, the authors proposed two corrections:
1. **Change in Magnetic Flux Caused by Boundary Length Change**: Equation (8) describes the change in magnetic flux due to the change in the area boundary length.
2. **Change in Magnetic Flux Caused by Cell Area Change**: Equation (9) describes the change in magnetic flux due to the change in the cell area.
3. **High Precision and High Efficiency**
- This method shows higher precision and efficiency in standard magnetohydrodynamics test cases and is superior to the fixed - grid method.
- The effectiveness of the method was verified through two - dimensional and three - dimensional tests, including the 1D Brio - Wu test, the 1D switch shock test, the 2D blast wave test, and the 3D blast wave test.
### Conclusion
The authors successfully applied the co - scaling grid method to magnetohydrodynamics simulations, solving the problems of large - scale spatial ranges and complex flow symmetries. This method not only improves the precision and efficiency of the simulation but also shows its broad application potential in handling astrophysical fluid dynamics problems.
### Formula Summary
- **Time Evolution Equation of Magnetic Flux**
\[
\frac{\partial}{\partial t} \Phi=\int_{\partial \Sigma}(v \times B) \cdot ds-\int_{\partial \Sigma}(w \times B) \cdot ds
\]
where $\Phi$ is the magnetic flux, $v$ is the volume velocity, $B$ is the magnetic field, and $w$ is the grid speed.
- **Cell Area Change Correction**
\[
B_3^{i, j, k}(t_{n + 1})=\frac{A_3^{i, j, k}(t_n)}{A_3^{i, j, k}(t_{n + 1})}\left[B_3^{i, j, k}(t_n)+\int_{t_n}^{t_{n + 1}} dt\left(\frac{d}{dt} B_3^{i, j, k}\right)\right]
\]
These formulas ensure the accurate update of the magnetic flux in the co - scaling grid, thereby improving the stability and precision of the simulation.