Ronald M. Caplan,Craig D. Johnston,Lars K. S. Daldoff,Jon A. Linker
Abstract:Unconditionally stable time stepping schemes are useful and often practically necessary for advancing parabolic operators in multi-scale systems. However, serious accuracy problems may emerge when taking time steps that far exceed the explicit stability limits. In our previous work, we compared the accuracy and performance of advancing parabolic operators in a thermodynamic MHD model using an implicit method and an explicit super time-stepping (STS) method. We found that while the STS method outperformed the implicit one with overall good results, it was not able to damp oscillatory behavior in the solution efficiently, hindering its practical use. In this follow-up work, we evaluate an easy-to-implement method for selecting a practical time step limit (PTL) for unconditionally stable schemes. This time step is used to `cycle' the operator-split thermal conduction and viscosity parabolic operators. We test the new time step with both an implicit and STS scheme for accuracy, performance, and scaling. We find that, for our test cases here, the PTL dramatically improves the STS solution, matching or improving the solution of the original implicit scheme, while retaining most of its performance and scaling advantages. The PTL shows promise to allow more accurate use of unconditionally stable schemes for parabolic operators and reliable use of STS methods.
Computational Engineering, Finance, and Science,Astrophysics of Galaxies,Solar and Stellar Astrophysics,Numerical Analysis,Computational Physics
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the precision problems encountered when using unconditionally stable time - stepping schemes in the thermodynamical magnetohydrodynamics (MHD) model. Specifically:
1. **Time - step selection problem**: When using very large time steps, although unconditionally stable time - stepping methods can be numerically stable, they will lead to the inability to correctly capture fast processes, thus introducing large errors. In particular, in the explicit super - time - stepping (STS) method, an overly large time step will lead to oscillatory behavior of the solution, affecting its practical application.
2. **Improving the practicality of the STS method**: In previous comparative studies, the author found that although the STS method is superior to implicit methods in computational performance, in some cases it will produce solution artifacts, limiting its wide application. Therefore, this paper proposes a simple and easy - to - implement method to dynamically calculate the practical time - step limit (PTL) to improve the accuracy and reliability of the STS method.
3. **Ensuring stability and accuracy in multi - scale systems**: In multi - scale systems, the time scales of different physical processes vary greatly, leading to difficulties in numerical integration. Unconditionally stable time - stepping algorithms can eliminate the limitation of numerical time steps, but a reasonable choice of time steps is required to ensure the accuracy and stability of the solution.
### Solutions
To solve the above problems, this paper proposes the following solutions:
- **Practical time - step limit (PTL)**: By dynamically calculating a suitable time step, such that within each large time step, multiple small - time - step iterations are performed on the split parabolic operator. This can not only reduce solution artifacts but also maintain high computational performance.
- **Testing and verification**: PTL was implemented in the MAS model, and its effect was tested through two practical cases (Test 1 and Test 2). The results show that PTL significantly improves the solution quality of the STS method, making it close to or better than the implicit BE + PCG method while retaining most of the performance advantages.
### Conclusions
By introducing PTL, this paper has successfully improved the accuracy and reliability of the STS method in the thermodynamical MHD model, making it a more competitive option. PTL not only solves the oscillation problem in the STS method but also maintains its computational efficiency to a certain extent, making the STS method more practical and robust in the simulation of multi - scale systems.