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
|
2
|
Khomenko E, Collados M, Vitas N, González-Morales PA. Influence of ambipolar and Hall effects on vorticity in three-dimensional simulations of magneto-convection. Philos Trans A Math Phys Eng Sci 2021; 379:20200176. [PMID: 33342386 PMCID: PMC7780135 DOI: 10.1098/rsta.2020.0176] [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] [Subscribe] [Scholar Register] [Accepted: 09/18/2020] [Indexed: 06/12/2023]
Abstract
This paper presents the results of the analysis of three-dimensional simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation runs are compared: one including both ambipolar diffusion and the Hall effect; one including only ambipolar diffusion and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 h of simulated solar time after the local dynamo reaches the saturation regime. We analyse the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.
Collapse
Affiliation(s)
- E. Khomenko
- Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
- Departamento 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
- Departamento 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
- Departamento 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
- Departamento de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife, Spain
| |
Collapse
|