pith. sign in

arxiv: 2606.18745 · v1 · pith:O3TXA7FVnew · submitted 2026-06-17 · ⚛️ physics.comp-ph · physics.plasm-ph

Extension of a multi-region free-surface MHD solver beyond the inductionless approximation

Min Ki Jung , Brian Wynne , Francisco Saenz , Yufan Xu , Jabir Al-Salami , Yong-Su Na , Egemen Kolemen This is my paper

Pith reviewed 2026-06-26 19:10 UTC · model grok-4.3

classification ⚛️ physics.comp-ph physics.plasm-ph
keywords free-surface MHDinduced magnetic fieldvector-potential formulationliquid metal flowsmulti-region solverHartmann numbermagnetic Reynolds number
0
0 comments X

The pith

A vector-potential formulation extends a free-surface MHD solver to resolve the induced magnetic field self-consistently.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends an existing transient three-dimensional free-surface MHD solver to include the effects of the induced magnetic field. Prior solvers relied on the inductionless approximation that neglects this field. The extension uses a vector-potential formulation that preserves the original multi-region two-phase structure and enforces the divergence-free condition on the magnetic field by construction. The updated solver is checked against analytical duct-flow solutions and tested against experimental free-surface measurements. This change removes the approximation entirely rather than relaxing it, creating a basis for treating conditions where the induced field becomes significant.

Core claim

The solver is extended beyond the inductionless approximation to resolve the evolution of the induced magnetic field using a vector-potential formulation that enforces the divergence-free condition by construction while preserving the multi-region two-phase framework.

What carries the argument

The vector-potential formulation for the magnetic field, which enforces the divergence-free condition by construction.

If this is right

  • The solver can treat finite magnetic Reynolds number regimes without approximation.
  • It supplies the foundation for modeling transient events in large-scale liquid-metal systems.
  • Verification against analytical duct solutions holds across a range of Hartmann numbers.
  • Validation reproduces experimental free-surface height data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same formulation could be applied to other multi-region two-phase MHD problems that currently use the inductionless limit.
  • It enables direct numerical studies of how induced fields alter flow stability at higher magnetic Reynolds numbers.
  • Coupling this magnetic treatment with time-dependent external fields would test its behavior under realistic transient conditions.

Load-bearing premise

The vector-potential formulation correctly enforces the divergence-free condition on the magnetic field while preserving the original multi-region two-phase framework.

What would settle it

A mismatch between the solver's predicted free-surface heights and measured values in an experiment conducted at magnetic Reynolds number greater than one, where the inductionless approximation would produce a detectable deviation.

Figures

Figures reproduced from arXiv: 2606.18745 by Brian Wynne, Egemen Kolemen, Francisco Saenz, Jabir Al-Salami, Min Ki Jung, Yong-Su Na, Yufan Xu.

Figure 1
Figure 1. Figure 1: At the start of each step the time-step size is updated to satisfy the flow and [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1. Flowchart of the per-time-step solution procedure. Steps highlighted with an orange [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The two cases differ only in the electrical boundary condition imposed on the Hart [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Schematic of the square duct geometry used for the Shercliff and Hunt benchmark cases. [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Normalized velocity profiles for the Shercliff flow at [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Normalized induced magnetic field [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Normalized velocity profiles for the Hunt flow at [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Normalized induced magnetic field [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Schematic of the LMX-U experimental device. Cross-sectional side view ( [PITH_FULL_IMAGE:figures/full_fig_p021_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Steady-state free-surface height profiles along the LMX-U channel at applied transverse [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Steady-state field, flow, and current structure in the LMX-U simulation at [PITH_FULL_IMAGE:figures/full_fig_p023_9.png] view at source ↗
read the original abstract

Free-surface liquid metal flows are a leading candidate for the plasma-facing components of future fusion reactors. Existing transient, three-dimensional, free-surface MHD solvers rely on the inductionless approximation in which the induced magnetic field is neglected. This paper extends the open-source solver FreeMHD [B. Wynne et al., Phys. Plasmas 32, 013907 (2025)] beyond the inductionless approximation to resolve the induced magnetic field self-consistently using a vector-potential formulation that enforces $\nabla\cdot\boldsymbol{B}=0$ by construction while preserving the original multi-region, two-phase framework. The solver is verified against analytical Shercliff and Hunt duct-flow solutions across a range of Hartmann numbers and validated against free-surface height measurements from the LMX-U experiment. To the best of our knowledge, FreeMHD2 is the first open-source, experimentally validated free-surface liquid metal solver to resolve the evolution of the induced magnetic field without invoking the inductionless approximation. By removing this approximation rather than relaxing it, the formulation provides the basis for future modeling of the finite magnetic Reynolds number conditions expected in large-scale, transient fusion events.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript extends the open-source FreeMHD solver (Wynne et al., Phys. Plasmas 2025) to FreeMHD2 by replacing the inductionless approximation with a vector-potential formulation (B = curl A) that enforces ∇·B = 0 by construction. The multi-region two-phase framework is preserved. The new solver is verified against the Shercliff and Hunt analytical duct-flow solutions over a range of Hartmann numbers and validated against free-surface height measurements from the LMX-U experiment. The authors claim that FreeMHD2 is the first open-source, experimentally validated free-surface liquid-metal solver that resolves the induced magnetic field without the inductionless approximation and that the formulation supplies a basis for finite-Rm transient fusion modeling.

Significance. If the implementation details and reported agreements hold, the work is significant because it removes rather than relaxes the inductionless approximation while retaining an existing multi-region free-surface infrastructure. The combination of analytical verification across Hartmann numbers and experimental validation on height data, together with the open-source release, supplies a concrete, reusable tool for modeling liquid-metal plasma-facing components under conditions where induced-field effects become non-negligible.

minor comments (3)
  1. [Verification section] The abstract states that the solver is 'verified against analytical Shercliff and Hunt duct-flow solutions across a range of Hartmann numbers,' but the main text should supply the specific Hartmann-number values, mesh resolutions, and quantitative error measures (e.g., L2 or L∞ norms of velocity or B) so that readers can judge the accuracy of the induced-field resolution.
  2. [Validation section] The validation is performed on free-surface height; the manuscript should clarify whether any induced-field quantities (e.g., surface current or local B perturbation) were compared with experiment or whether the height agreement is taken as indirect evidence that the induced field is correctly captured.
  3. [Introduction or Discussion] A brief statement of the magnetic Reynolds number range accessed by the LMX-U data and by the intended fusion applications would help readers assess how far the new capability extends beyond the inductionless regime.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the work, the recognition of its significance in removing the inductionless approximation while retaining the multi-region free-surface framework, and the recommendation for minor revision. No specific major comments appear in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The derivation extends the prior FreeMHD solver by adopting the standard vector-potential formulation B = curl(A) to resolve the induced field while preserving the multi-region framework; this is a direct, non-circular removal of the inductionless approximation using an established technique that enforces div(B)=0 by construction. Verification relies on independent analytical duct-flow solutions (Shercliff/Hunt) and external LMX-U experimental height data rather than any fitted parameters, self-definitions, or load-bearing self-citations that reduce the central result to its inputs. The novelty claim is external and does not affect the mathematical chain.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the work relies on standard MHD equations and numerical methods from prior literature.

pith-pipeline@v0.9.1-grok · 5753 in / 1102 out tokens · 31543 ms · 2026-06-26T19:10:48.307372+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

35 extracted references · 34 canonical work pages

  1. [1]

    M. K. Jung. Input files for ``Extension of a multi-region free-surface MHD solver beyond the inductionless approximation''. Zenodo. 2026. doi:10.5281/zenodo.20727182

    work page doi:10.5281/zenodo.20727182 2026

  2. [2]

    Alegre and E

    D. Alegre and E. Oyarzabal and D. Tafalla and M. Liniers and A. Soleto and F. Tabar \'e s. Design and testing of advanced liquid metal targets for DEMO divertor: The OLMAT Project. Journal of Fusion Energy. 2020. doi:10.1007/s10894-020-00254-5

    work page doi:10.1007/s10894-020-00254-5 2020

  3. [3]

    Al-Salami

    J. Al-Salami. Numerical simulation of MHD free surface liquid metal flows for nuclear fusion applications. 2022

    2022

  4. [4]

    Brackbill and D

    J. Brackbill and D. Kothe and C. Zemach. A continuum method for modeling surface tension. Journal of Computational Physics. 1992. doi:10.1016/0021-9991(92)90240-Y

    work page doi:10.1016/0021-9991(92)90240-y 1992

  5. [5]

    Caretto and A

    L. Caretto and A. Gosman and S. Patankar and D. Spalding. Two calculation procedures for steady, three-dimensional flows with recirculation. Proceedings of the Third International Conference on Numerical Methods in Fluid Mechanics. 1973. doi:10.1007/BFb0112677

    work page doi:10.1007/bfb0112677 1973

  6. [6]

    Cortazar and N

    J. Cortazar and N. Lazarus. Creating a softer RoboSquid: Liquid-metal-based compliant pumps for pulsed jet propulsion. Advanced Materials Technologies. 2025. doi:10.1002/admt.202401372

    work page doi:10.1002/admt.202401372 2025

  7. [7]

    Deshpande and L

    S. Deshpande and L. Anumolu and M. Trujillo. Evaluating the performance of the two-phase flow solver interFoam. Computational Science & Discovery. 2012. doi:10.1088/1749-4699/5/1/014016

    work page doi:10.1088/1749-4699/5/1/014016 2012

  8. [8]

    Endeve and D

    E. Endeve and D. Stefanski and M.-O. G. Delchini and S. Slattery and C. D. Hauck and B. Turcksin and S. Smolentsev. A full-induction magnetohydrodynamics solver for liquid metal fusion blankets in Vertex-CFD. Fusion Engineering and Design. 2026. doi:10.1016/j.fusengdes.2026.115792

    work page doi:10.1016/j.fusengdes.2026.115792 2026

  9. [9]

    A. E. Fisher and Z. Sun and E. Kolemen. Liquid metal ``divertorlets'' concept for fusion reactors. Nuclear Materials and Energy. 2020. doi:10.1016/j.nme.2020.100855

    work page doi:10.1016/j.nme.2020.100855 2020

  10. [10]

    Haiqi and W

    Y. Haiqi and W. Baofeng and L. Huiqin and L. Jianchao. Influence of electromagnetic brake on flow field of liquid steel in the slab continuous casting mold. Journal of Materials Processing Technology. 2008. doi:10.1016/j.jmatprotec.2007.08.054

    work page doi:10.1016/j.jmatprotec.2007.08.054 2008

  11. [11]

    R. Issa. Solution of the implicitly discretised fluid flow equations by operator-splitting. Journal of Computational Physics. 1986. doi:10.1016/0021-9991(86)90099-9

    work page doi:10.1016/0021-9991(86)90099-9 1986

  12. [12]

    Kawczynski and S

    C. Kawczynski and S. Smolentsev and M. Abdou. An induction-based magnetohydrodynamic 3D code for finite magnetic Reynolds number liquid-metal flows in fusion blankets. Fusion Engineering and Design. 2016. doi:10.1016/j.fusengdes.2016.02.088

    work page doi:10.1016/j.fusengdes.2016.02.088 2016

  13. [13]

    Kawczynski and S

    C. Kawczynski and S. Smolentsev and M. Abdou. Characterization of the lid-driven cavity magnetohydrodynamic flow at finite magnetic Reynolds numbers using far-field magnetic boundary conditions. Physics of Fluids. 2018. doi:10.1063/1.5036775

    work page doi:10.1063/1.5036775 2018

  14. [14]

    Levitt and E

    B. Levitt and E. T. Meier and R. Umstattd and J. R. Barhydt and I. A. M. Datta and C. Liekhus-Schmaltz and D. A. Sutherland and B. A. Nelson. The Zap Energy approach to commerical fusion. Physics of Plasmas. 2023. doi:10.1063/5.0163361

    work page doi:10.1063/5.0163361 2023

  15. [15]

    Miyazawa and T

    J. Miyazawa and T. Goto and T. Murase and T. Ohgo and N. Yanagi and H. Tanaka and H. Tamura and T. Tanaka and S. Masuzaki and R. Sakamoto and J. Yagi and A. Sagara and the FFHR Design Group. Conceptual design of a liquid metal limiter/divertor system for the FFHR-d1. Fusion Engineering and Design. 2017. doi:10.1016/j.fusengdes.2017.07.003

    work page doi:10.1016/j.fusengdes.2017.07.003 2017

  16. [16]

    N. B. Morley and S. Smolentsev and R. Munipalli and M.-J. Ni and D. Gao and M. Abdou. Progress on the modeling of liquid metal, free surface, MHD flows for fusion liquid walls. Fusion Engineering and Design. 2004. doi:10.1016/j.fusengdes.2004.07.013

    work page doi:10.1016/j.fusengdes.2004.07.013 2004

  17. [17]

    Morley and J

    N. Morley and J. Burris and L. Cadwallader and M. Nornberg. GaInSn usage in the research laboratory. Review of Scientific Instruments. 2008. doi:10.1063/1.2930813

    work page doi:10.1063/1.2930813 2008

  18. [18]

    u ller and L. B \

    U. M \"u ller and L. B \"u hler. Magnetofluiddynamics in Channels and Containers. 2001. doi:10.1007/978-3-662-04405-6

    work page doi:10.1007/978-3-662-04405-6 2001

  19. [19]

    Pintsuk and G

    G. Pintsuk and G. Aiello and S. Dudarev and M. Gorley and J. Henry and M. Richou and M. Rieth and D. Terentyev and R. Vila. Materials for in-vessel components. Fusion Engineering and Design. 2022. doi:10.1016/j.fusengdes.2021.112994

    work page doi:10.1016/j.fusengdes.2021.112994 2022

  20. [20]

    Rindt and J

    P. Rindt and J. van den Eijnden and T. Morgan and N. Lopes Cardozo. Conceptual design of a liquid-metal divertor for the European DEMO. Fusion Engineering and Design. 2021. doi:10.1016/j.fusengdes.2021.112812

    work page doi:10.1016/j.fusengdes.2021.112812 2021

  21. [21]

    Saenz and Z

    F. Saenz and Z. Sun and A. E. Fisher and B. Wynne and E. Kolemen. Divertorlets concept for low-recycling fusion reactor divertor: experimental, analytical and numerical verification. Nuclear Fusion. 2022. doi:10.1088/1741-4326/ac6682

    work page doi:10.1088/1741-4326/ac6682 2022

  22. [22]

    Samulyak and J

    R. Samulyak and J. Du and J. Glimm and Z. Xu. A numerical algorithm for MHD of free surface flows at low magnetic Reynolds numbers. Journal of Computational Physics. 2007. doi:10.1016/j.jcp.2007.06.005

    work page doi:10.1016/j.jcp.2007.06.005 2007

  23. [23]

    Sergeev and B

    V. Sergeev and B. Kuteev and A. Bykov and A. Gervash and D. Glazunov and P. Goncharov and A. Dnestrovskij and R. Khayrutdinov and A. Klishchenko and V. Lukash and I. Mazul and P. Molchanov and V. Petrov and V. Rozhansky and Yu. Shpanskiy and A. Sivak and V. Skokov and A. Spitsyn. Conceptual design of divertor and first wall for DEMO-FNS. Nuclear Fusion. 2...

    work page doi:10.1088/0029-5515/55/12/123013 2015

  24. [24]

    Siriano and L

    S. Siriano and L. Melchiorri and S. Pignatiello and A. Tassone. A multi-region and a multiphase MHD OpenFOAM solver for fusion reactor analysis. Fusion Engineering and Design. 2024. doi:10.1016/j.fusengdes.2024.114216

    work page doi:10.1016/j.fusengdes.2024.114216 2024

  25. [25]

    Smolentsev and R

    S. Smolentsev and R. Moreau and L. B \"u hler and C. Mistrangelo. MHD thermofluid issues of liquid-metal blankets: Phenomena and advances. Fusion Engineering and Design. 2010. doi:10.1016/j.fusengdes.2010.02.038

    work page doi:10.1016/j.fusengdes.2010.02.038 2010

  26. [26]

    Smolyanov and O

    I. Smolyanov and O. Zikanov. Exploratory study of liquid-metal response to rapid variation of applied magnetic field. Physical Review E. 2025. doi:10.1103/nm74-f77r

    work page doi:10.1103/nm74-f77r 2025

  27. [27]

    Suarez and E

    D. Suarez and E. Iraola and C. Lamp \'o n and E. de les Valls and L. Batet. Liquid metal MHD flow influence on heat transfer phenomena in fusion reactor blankets. Fusion Engineering and Design. 2021. doi:10.1016/j.fusengdes.2021.112503

    work page doi:10.1016/j.fusengdes.2021.112503 2021

  28. [28]

    Suarez and A

    D. Suarez and A. Khodak and E. de les Valls. A formal verification and validation of a low magnetic Reynolds number MHD code for fusion applications. IEEE Transactions on Plasma Science. 2022. doi:10.1109/TPS.2022.3203801

    work page doi:10.1109/tps.2022.3203801 2022

  29. [29]

    Sun and J

    Z. Sun and J. Al-Salami and A. Khodak and F. Saenz and B. Wynne and R. Maingi and K. Hanada and C. Hu and E. Kolemen. Magnetohydrodynamics in free surface liquid metal flow relevant to plasma-facing components. Nuclear Fusion. 2023. doi:10.1088/1741-4326/acd864

    work page doi:10.1088/1741-4326/acd864 2023

  30. [30]

    H. G. Weller and G. Tabor and H. Jasak and C. Fureby. A tensorial approach to computational continuum mechanics using object-oriented techniques. Computers in Physics. 1998. doi:10.1063/1.168744

    work page doi:10.1063/1.168744 1998

  31. [31]

    Wynne and F

    B. Wynne and F. Saenz and J. Al-Salami and Y. Xu and Z. Sun and C. Hu and K. Hanada and E. Kolemen. FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows. Physics of Plasmas. 2025. doi:10.1063/5.0230242

    work page doi:10.1063/5.0230242 2025

  32. [32]

    Xiang and S

    L. Xiang and S. Yang and Q. Wang and J. Wu. Design and testing of a direct current electromagnetic pump for liquid metal. Thermal Science and Engineering Progress. 2024. doi:10.1016/j.tsep.2024.102560

    work page doi:10.1016/j.tsep.2024.102560 2024

  33. [33]

    Yang and N

    M. Yang and N. Ma and D. Bliss and G. Bryant. Melt motion during liquid-encapsulated Czochralski crystal growth in steady and rotating magnetic fields. International Journal of Heat and Fluid Flow. 2007. doi:10.1016/j.ijheatfluidflow.2006.08.001

    work page doi:10.1016/j.ijheatfluidflow.2006.08.001 2007

  34. [34]

    Zhang and J

    J. Zhang and J. Li and X. Liu. Investigate the impact of different neutron source models on the neutronics performance for CFETR. AIP Advances. 2025. doi:10.1063/5.0257060

    work page doi:10.1063/5.0257060 2025

  35. [35]

    Zhang and Y

    S. Zhang and Y. Liu and Q. Fan and C. Zhang and T. Zhou and K. Kalantar-Zadeh and Z. Guo. Liquid metal batteries for future energy storage. Energy & Environmental Science. 2021. doi:10.1039/D1EE00531F

    work page doi:10.1039/d1ee00531f 2021