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

Engineering notes

Engineering notes will be added by the Power for Solar editorial team.

Chinese explanation / 中文解读

中文解读待补充：本站会优先为光伏效率、钙钛矿太阳能电池、储能技术、太阳能热利用、BIPV、并网技术等高价值论文补充中文说明。

Original abstract

We present a projection method for controlling numerical \(\nabla\cdot\B\) errors in smoothed particle magnetohydrodynamics (SPMHD). The method corrects the magnetic field after an MHD update by solving an elliptic projection problem constructed from the same discrete divergence operator used to measure the error. A key ingredient is to use the adjoint gradient associated with a volume-weighted metric. With this choice, the projection gives an energy-minimizing correction, does not increase the discrete magnetic energy, and leads to a symmetric positive semidefinite linear system that can be solved by the conjugate-gradient method without explicitly assembling the matrix. We test the method using two-dimensional Dedner-type divergence tests and three-dimensional magnetized collapse calculations. With sufficiently many iterations, the projection reduces the divergence error to the floating-point roundoff level in both test problems. In realistic collapse runs, practical stopping criteria designed to reduce the divergence error generated by the underlying SPMHD update suppress the normalized divergence error well below that obtained in the divergence-cleaning run, with a projection cost of only about \(1\)--\(10\%\) of the SPMHD update cost. The density and plasma-\(β\) structures remain consistent when the projection interval is varied, whereas the divergence-cleaning run shows quantitative differences. These results indicate that the projection method is a robust and attractive alternative to divergence cleaning for controlling \(\nabla\cdot\B\) errors in SPMHD and related particle or meshless MHD schemes.

Links and sources

• Official / arXiv page
• PDF from original source

Comments