1
|
Modestov M, Khomenko E, Vitas N, de Vicente A, Navarro A, González-Morales PA, Collados M, Felipe T, Martínez-Gómez D, Hunana P, Luna M, Koll Pistarini M, Popescu Braileanu B, Perdomo García A, Liakh V, Santamaria I, Gomez Miguez MM. Mancha3D Code: Multipurpose Advanced Nonideal MHD Code for High-Resolution Simulations in Astrophysics. Sol Phys 2024; 299:23. [PMID: 38390515 PMCID: PMC10879365 DOI: 10.1007/s11207-024-02267-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 01/31/2024] [Indexed: 02/24/2024]
Abstract
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic (MHD) processes in solar/stellar atmospheres. The code includes nonideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high-quality realistic simulations of magnetoconvection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D (Multifluid (-purpose -physics -dimensional) Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics 3D) code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and a memory-saving Runge-Kutta (RK) scheme. It handles nonideal effects through super-time-stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a nonreflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up to several thousands of processors (CPUs).
Collapse
Affiliation(s)
- M. Modestov
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - E. Khomenko
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - N. Vitas
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - A. de Vicente
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - A. Navarro
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - P. A. González-Morales
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - M. Collados
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - T. Felipe
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - D. Martínez-Gómez
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - P. Hunana
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - M. Luna
- Departament de Física, Universitat de les Illes Balears, E-07122 Palma, Spain
- Institute of Applied Computing and Community Code (IAC3), UIB, Palma, Spain
| | - M. Koll Pistarini
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | | | - A. Perdomo García
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - V. Liakh
- Centre for Mathematical Plasma Astrophysics, KU Leuven, 3001 Leuven, Belgium
| | - I. Santamaria
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| | - M. M. Gomez Miguez
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife Spain
- Dpto. de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife Spain
| |
Collapse
|