Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 231, Issue 21, Pages 7214-7236Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2012.06.039
Keywords
Particle methods; Magnetohydrodynamics (MHD); Smoothed particle magnetohydrodynamics (SPMHD); Astrophysics; Magnetic fields; Divergence constraint; Divergence cleaning
Funding
- Endeavour IPRS
- APA
- Australian Research Council [DP1094585]
- Australian Research Council [DP1094585] Funding Source: Australian Research Council
Ask authors/readers for more resources
We present a constrained formulation of Dedner et al.'s hyperbolic/parabolic divergence cleaning scheme for enforcing the del . B = 0 constraint in smoothed particle magnetohydrodynamics (SPMHD) simulations. The constraint we impose is that energy removed must either be conserved or dissipated, such that the scheme is guaranteed to decrease the overall magnetic energy. This is shown to require use of conjugate numerical operators for evaluating del . B and del psi in the SPMHD cleaning equations. The resulting scheme is shown to be stable at density jumps and free boundaries, in contrast to an earlier implementation by Price and Monaghan (2005). Optimal values of the damping parameter are found to be sigma = 0.2-0.3 in 2D and sigma = 0.8-1.2 in 3D. With these parameters, our constrained Hamiltonian formulation is found to provide an effective means of enforcing the divergence constraint in SPMHD, typically maintaining average values of h vertical bar del . B vertical bar / vertical bar B vertical bar to 0.1-1%, up to an order of magnitude better than artificial resistivity without the associated dissipation in the physical field. Furthermore, when applied to realistic, 3D simulations we find an improvement of up to two orders of magnitude in momentum conservation with a corresponding improvement in numerical stability at essentially zero additional computational expense. (C) 2012 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available