pith. sign in

arxiv: 1206.6159 · v2 · pith:VMF6NWN3new · submitted 2012-06-27 · 🌌 astro-ph.IM · astro-ph.GA· astro-ph.SR· physics.comp-ph

Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics

classification 🌌 astro-ph.IM astro-ph.GAastro-ph.SRphysics.comp-ph
keywords nablacleaningconstrainedconstraintdivergenceschemespmhdenergy
0
0 comments X
read the original abstract

We present a constrained formulation of Dedner et al's hyperbolic/parabolic divergence cleaning scheme for enforcing the \nabla\dot 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 \nabla\dot B and \nabla{\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 & 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 |\nabla\dot B| / |B| 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.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. An Adjoint Projection Formulation for Enforcing the divergence-free Constraint in Smoothed Particle Magnetohydrodynamics

    astro-ph.SR 2026-06 unverdicted novelty 7.0

    Introduces an adjoint projection method using volume-weighted metric to control numerical ∇·B errors in SPMHD, achieving roundoff-level accuracy with 1-10% overhead and consistent physical structures versus divergence...