Simulating the interaction of a non-magnetized planet with the stellar wind produced by a sun-like star using the FLASH Code
Pith reviewed 2026-05-24 22:18 UTC · model grok-4.3
The pith
The FLASH code can be adapted to place rigid bodies in MHD simulations of a non-magnetized planet interacting with stellar wind from a sun-like star.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors adapt the existing rigid body placement resources in FLASH to enable their use in MHD scenarios. With this adaptation they simulate the interaction of a non-magnetized planet that has no significant atmosphere with the stellar wind produced by a sun-like star. They examine the patterns of density, magnetic field and velocity around the planet and the influence of viscosity on those patterns, providing an improved methodological approach for other users.
What carries the argument
Adaptation of the rigid body placement feature from hydrodynamic to MHD cases within the FLASH code.
If this is right
- Density, magnetic field, and velocity patterns around the planet become quantifiable under the adapted scheme.
- Viscosity is shown to modify those patterns in a measurable way.
- The adapted method supports simulation of other solid objects in MHD flows.
- The improved approach is made available for reuse by other researchers.
Where Pith is reading between the lines
- The same adaptation could be tested on other non-magnetized bodies such as moons or asteroids.
- Adding a thin atmosphere layer in follow-on runs would isolate its separate contribution to the interaction.
- Direct comparison of the simulated magnetic draping with spacecraft data from similar encounters would test the model's realism.
- Varying the viscosity parameter across a wider range could map how dissipation affects bow-shock standoff distance.
Load-bearing premise
The rigid body placement feature, tested only for pure hydrodynamics, can be extended to MHD without introducing major unaccounted errors or instabilities.
What would settle it
Perform the adapted MHD simulation of the planet in stellar wind and check whether the magnetic field and flow patterns remain stable and free of numerical artifacts that would not appear in the pure hydrodynamic case.
Figures
read the original abstract
The study of the interaction between solid objects and magnetohydrodynamic (MHD) fluids is of great importance in physics as consequence of the significant phenomena generated, such as planets interacting with stellar wind produced by their host stars. There are several computational tools created to simulate hydrodynamic and MHD fluids, such as the FLASH code. In this code there is a feature which permits the placement of rigid bodies in the domain to be simulated. However, it is available and tested for pure hydrodynamic cases only. Our aim here is to adapt the existing resources of FLASH to enable the placement of a rigid body in MHD scenarios and, with such a scheme, to produce the simulation of a non-magnetized planet interacting with the stellar wind produced by a sun-like star. Besides, we consider that the planet has no significant atmosphere. We focus our analysis on the patterns of the density, magnetic field and velocity around the planet, as well as the influence of the viscosity on such patterns. At last, an improved methodological approach is available to other interested users.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper describes an adaptation of the rigid-body placement feature in the FLASH code (previously tested only for pure hydrodynamics) to MHD simulations. This is applied to model the interaction of a non-magnetized planet with no significant atmosphere and the stellar wind from a sun-like star, with analysis focused on resulting patterns in density, magnetic field, and velocity, plus the role of viscosity. The authors state that an improved methodological approach is made available to other users.
Significance. If the MHD adaptation of the rigid-body module is correctly implemented and validated, the work would provide a practical extension of an existing open code for space-physics applications involving solid-body boundaries in magnetized flows. The topic is relevant to planetary magnetospheres and stellar-wind interactions, but the absence of any reported tests leaves the reliability of the produced patterns undetermined.
major comments (2)
- [Abstract] Abstract: the central claim that the rigid-body feature can be adapted to MHD scenarios without major errors rests on an unverified extension; the text supplies no description of the modified boundary conditions for the magnetic field, no adjustments to divergence cleaning at the interface, and no benchmark or convergence tests against known MHD problems.
- [Abstract] The reported density, magnetic-field and velocity patterns around the planet cannot be assessed for physical fidelity versus numerical artifacts because no validation results, error analysis, or comparison with analytic or other-code solutions are presented.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. The points raised highlight the need for greater detail on the MHD adaptation and for explicit validation of the results. We agree with these assessments and will revise the manuscript to address them.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that the rigid-body feature can be adapted to MHD scenarios without major errors rests on an unverified extension; the text supplies no description of the modified boundary conditions for the magnetic field, no adjustments to divergence cleaning at the interface, and no benchmark or convergence tests against known MHD problems.
Authors: We agree that the manuscript does not describe the MHD-specific modifications to the rigid-body module. The extension was implemented by applying the existing rigid-body boundary routines to the additional MHD variables (magnetic field components and related quantities), with the magnetic field inside the planet set to zero to represent a non-magnetized body. However, no explicit discussion of boundary conditions for B or modifications to divergence cleaning at the interface is provided. We will add a new methods subsection detailing these adaptations. We also acknowledge the absence of MHD-specific benchmarks or convergence tests; the original hydro implementation was tested, but this work did not include equivalent MHD tests. We will incorporate benchmark cases (e.g., MHD flow past a sphere) and convergence studies in the revised version. revision: yes
-
Referee: [Abstract] The reported density, magnetic-field and velocity patterns around the planet cannot be assessed for physical fidelity versus numerical artifacts because no validation results, error analysis, or comparison with analytic or other-code solutions are presented.
Authors: The patterns presented are the direct output of applying the adapted code to the planet-stellar wind interaction scenario. We recognize that without dedicated validation, it is not possible to fully separate physical features from numerical effects. In the revision we will add a validation section that includes comparisons with analytic expectations for magnetic draping around a conducting obstacle and, where feasible, with published results from other MHD codes for analogous stellar-wind/planet interactions. This will allow readers to evaluate the fidelity of the reported density, magnetic-field, and velocity structures. revision: yes
Circularity Check
No significant circularity; direct code adaptation and simulation description.
full rationale
The paper describes adapting the existing rigid-body module in FLASH (noted as tested only for hydrodynamics) to MHD and running a simulation of a non-magnetized planet in stellar wind. No derivation chain, equations, or predictions are presented that reduce to inputs by construction. No fitted parameters are relabeled as predictions, no self-citation chains support central claims, and no ansatz or uniqueness theorems are invoked. The work is a methods-and-results description of numerical setup, making it self-contained against external benchmarks with no circular steps.
Axiom & Free-Parameter Ledger
free parameters (1)
- viscosity
axioms (2)
- domain assumption MHD equations govern the stellar wind plasma and its interaction with the planet.
- domain assumption The planet can be modeled as a rigid body with no significant magnetic field or atmosphere.
Reference graph
Works this paper leans on
-
[1]
T. Takahashi, H. Ueki, A. Kunimatsu, H. Fujii, in ACM SIGGRAPH 2002 Conference Abstracts and Applications (ACM, New York-USA, 2002), p.266.https://doi.org/10. 1145/1242073.1242279
-
[2]
N. Takashi, T.J.R. Hughes, Comput. Method. Appl. M. 95, 115–138 (1992). https://doi.org/10.1016/ 0045-7825(92)90085-X
work page 1992
-
[3]
K.G. Powell, P.L. Roe, T.J. Linde, T.I. Gombosi, D.L. De Zeeuw, J. Comp. Phys. 154, 284–309 (1999). https: //doi.org/10.1006/jcph.1999.6299
-
[4]
D.G.E. Grigoriadis, I.E. Sarris, S.C. Kassinos, Comput. Fluids 39, 345–258 (2010). https://doi.org/10.1016/j. compfluid.2009.09.012
work page doi:10.1016/j 2010
-
[5]
C.P. Johnstone, et al., Astron. Astrophys. 577, A122 (2015). https://doi.org/10.1051/0004-6361/201425134
-
[6]
Y. Vernisse, et al., Planet. Space Sci. 84, 37–47 (2013). https://doi.org/10.1016/j.pss.2013.04.004
-
[7]
O. Cohen, et al., Astrophys. J. 806(1), 41 (2015). https: //doi.org/10.1088/0004-637X/806/1/41
-
[8]
J.D. Nichols, S.E. Milan, Mon. Not. R. Astron. Soc. 461, 2353–2366 (2016). https://doi.org/10.1093/mnras/ stw1430
-
[9]
V. Bourrier, et al., Astron. Astrophys. 591, A121 (2016). https://doi.org/10.1051/0004-6361/201628362
-
[10]
J.R. Spreiter, A.L. Summers, A.W. Rizzi, Planet. Space Sci. 18, 1281–1299 (1970). https://doi.org/10.1016/ 0032-0633(70)90139-X
work page 1970
-
[11]
M. Dryer, A.W. Rizzi, Wen-Wu Shen, Astrophys. Space Sci. 22, 329–351 (1973). https://doi.org/10.1007/ BF00647431
work page 1973
-
[12]
J.D. Alvarado-G´ omez, et al., Astron. Astrophys. 594, A95 (2016). https://doi.org/10.1051/0004-6361/ 201628988
-
[13]
Bateman, MHD Instabilities (The MIT Press, Massachusetts-USA, 1978)
G. Bateman, MHD Instabilities (The MIT Press, Massachusetts-USA, 1978)
work page 1978
-
[14]
Biskamp, Magnetohydrodynamic Turbulence (Cam- bridge University Press, Cambridge-UK, 2003)
D. Biskamp, Magnetohydrodynamic Turbulence (Cam- bridge University Press, Cambridge-UK, 2003)
work page 2003
-
[15]
J.P. Goedbloed, S. Poedts, Principles of Magnetohy- drodynamics (Cambridge University Press, Cambridge- UK,2004)
work page 2004
-
[16]
A.E. Lifschitz, Developments in Electromagnetic The- ory and Applications: Magnetohydrodynamics and Spec- tral Theory (Kluwer Academic Publishers, Dordrecht-The Netherlands, 1989)
work page 1989
-
[17]
P. MacNeice, et al., Comput. Phys. Commun. 126(3), 330–354 (2000). https://doi.org/10.1016/ S0010-4655(99)00501-9
work page 2000
-
[18]
R. L¨ ohner, Comp. Meth. App. Mech. Eng.61(3), 323–338 (1987). https://doi.org/10.1016/0045-7825(87)90098-3
-
[19]
B. Fryxell, et al., Astrophys. J. Suppl. S. 131(1), 273–334 (2000). https://doi.org/10.1086/317361
-
[20]
P.L. Roe, J. Comp. Phys. 43(2), 357–372 (1981). https: //doi.org/10.1016/0021-9991(81)90128-5
-
[21]
S. Li, J. Comput. Phys. 203(1), 344–357 (2005). https: //doi.org/10.1016/j.jcp.2004.08.020
-
[22]
B. Einfeldt, C.D. Munz, P.L. Roe, B. Sj¨ ogreen, J. Com- put. Phys. 92(2), 273–295 (1991). https://doi.org/10. 1016/0021-9991(91)90211-3
work page 1991
-
[23]
B. van Leer, SIAM J. Sci. Stat. Comp. 5(1), 1–20 (1984). https://doi.org/10.1137/0905001
-
[24]
D. Lee, A.E. Deane, J. Comput. Phys. 228(4), 952–975 (2009). https://doi.org/10.1016/j.jcp.2008.08.026
-
[25]
R. Courant, K. Friedrichs, H. Lewy, IBM J. Res. Dev. 11(2), 215–234 (1967). https://doi.org/10.1147/rd.112. 0215 16
-
[26]
E.N. Parker, Astrophys. J. 128, 664–676 (1958). https: //doi.org/10.1086/146579
-
[27]
M.G. Kivelson, C.T. Russell (ed.), Introduction to Space Physics (Cambridge University Press, Cambridge- UK,1995)
work page 1995
-
[28]
R.C. Tautz, A. Shalchi, A. Dosch, J. Geophys. Res.: Space Phys. 116(A2), A02102 (2011). https://doi.org/ 10.1029/2010JA015936
-
[29]
H. Korth, et al., Planet. Space Sci. 59, 2075–2085 (2011). https://doi.org/10.1016/j.pss.2010.10.014
-
[30]
P. Subramanian, A. Lara, A. Borgazzi, Geophys. Res. Lett. 39(19), L19107 (2012). https://doi.org/10.1029/ 2012GL053625
work page 2012
-
[31]
H. P´ erez-de-Tejada, Astrophys. J. Lett. 618(2), L145– L148 (2005). https://doi.org/10.1086/425864
-
[32]
E. Marsch, Living Rev. Solar Phys. 3(1), 1 (2006). https: //doi.org/10.12942/lrsp-2006-1
-
[33]
S. Livi, E. Marsch, H. Rosenbauer, J. Geophys. Res.: Space Phys. 91(A7), 8045–8050 (1986). https://doi.org/ 10.1029/JA091iA07p08045
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.