A componentwise spherical Dirichlet-to-Neumann closure for finite-volume MHD at finite magnetic Reynolds number
Abstract
At finite magnetic Reynolds number, the field induced by a moving conducting fluid or solid extends into open space, so a finite mesh needs an outer boundary condition that represents the decaying exterior field. Magnetohydrodynamic (MHD) codes widely use pseudo-vacuum or truncation at a finite radius. Neither condition is physically correct for a vacuum exterior, but the exact alternatives are costly or limited to simple geometries. We present a Dirichlet-to-Neumann (DtN) closure on a fixed enclosing sphere and release it as part of the open-source multi-region finite-volume solver FreeMHD. The closure matches each Cartesian component of the vector potential to its decaying exterior solution, which retains all three decaying vector-Laplace families. A meshed vacuum buffer connects non-spherical conductors to the sphere. We verify the closure with manufactured modes, sphere free decay, sphere and oblate dynamos, and a freely braking rotor. Under refinement, oblate dynamo thresholds and sphere growth rates converge toward independent spectral references, and the rotor torque matches an exact transient reference. On the same sphere, pseudo-vacuum and truncation give the larger errors predicted by their exterior models.
Related
- Truncated Polyakov bootstrap for BCFTs: Neumann-Dirichlet flow and Ising special transition
- Encounter Propagator for Multiple Targets: A Dirichlet-to-Neumann Spectral Formalism
- Residual-Based Adaptive Refinement in Hybrid Variational-Collocation PINNs for Singular Elliptic PDEs with Mixed Dirichlet-Neumann Boundary Conditions
- Critical magnetic Reynolds number of the weakly-collisional turbulent dynamo in supersonic plasmas
- Can componentwise D7-tadpole cancellation yield a chiral Standard Model? A no-go for smooth O7− orientifolds
BibTeX
@misc{jung2026componentwise,
title = {{A componentwise spherical Dirichlet-to-Neumann closure for finite-volume MHD at finite magnetic Reynolds number}},
author = {Min Ki Jung and Brian Wynne and Francisco Saenz and Yufan Xu and Yong-Su Na and Egemen Kolemen},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.07304},
eprint = {2611.07304},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.07304},
}
RIS
TY - UNPB
TI - A componentwise spherical Dirichlet-to-Neumann closure for finite-volume MHD at finite magnetic Reynolds number
AU - Min Ki Jung
AU - Brian Wynne
AU - Francisco Saenz
AU - Yufan Xu
AU - Yong-Su Na
AU - Egemen Kolemen
PY - 2026
DA - 2026/11/05
PB - arXiv
DO - 10.48550/arxiv.2611.07304
UR - https://arxiv.org/abs/2611.07304
AB - At finite magnetic Reynolds number, the field induced by a moving conducting fluid or solid extends into open space, so a finite mesh needs an outer boundary condition that represents the decaying exterior field. Magnetohydrodynamic (MHD) codes widely use pseudo-vacuum or truncation at a finite radius. Neither condition is physically correct for a vacuum exterior, but the exact alternatives are costly or limited to simple geometries. We present a Dirichlet-to-Neumann (DtN) closure on a fixed enclosing sphere and release it as part of the open-source multi-region finite-volume solver FreeMHD. The closure matches each Cartesian component of the vector potential to its decaying exterior solution, which retains all three decaying vector-Laplace families. A meshed vacuum buffer connects non-spherical conductors to the sphere. We verify the closure with manufactured modes, sphere free decay, sphere and oblate dynamos, and a freely braking rotor. Under refinement, oblate dynamo thresholds and sphere growth rates converge toward independent spectral references, and the rotor torque matches an exact transient reference. On the same sphere, pseudo-vacuum and truncation give the larger errors predicted by their exterior models.
ER -