Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization

Proper initialization of the solar corona magnetic field is important for easing the iterative process of realistic and efficient global magnetohydrodynamics (MHD) simulations. This study introduces a high-order flux reconstruction (FR) framework for solving the Poisson equation, a necessary step fo...

Full description

Saved in:
Bibliographic Details
Main Authors: Rayan Dhib, Firas Ben Ameur, Vatsalya Sharma, Andrea Lani, Stefaan Poedts
Format: Article
Language:English
Published: IOP Publishing 2025-01-01
Series:The Astrophysical Journal
Subjects:
Online Access:https://doi.org/10.3847/1538-4357/adace5
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1823856583992410112
author Rayan Dhib
Firas Ben Ameur
Vatsalya Sharma
Andrea Lani
Stefaan Poedts
author_facet Rayan Dhib
Firas Ben Ameur
Vatsalya Sharma
Andrea Lani
Stefaan Poedts
author_sort Rayan Dhib
collection DOAJ
description Proper initialization of the solar corona magnetic field is important for easing the iterative process of realistic and efficient global magnetohydrodynamics (MHD) simulations. This study introduces a high-order flux reconstruction (FR) framework for solving the Poisson equation, a necessary step for computing a potential-field source-surface (PFSS) to initialize the magnetic field for global solar corona simulations with MHD. By hyperbolizing the elliptic Poisson equation into a set of hyperbolic equations, we develop an efficient and robust high-order PFSS solver. Our contributions include developing a Q 2 (i.e., quadratic) geometrical representation using prismatic elements for the computational domain, which enables high-order mesh generation. Such a hyperbolized Poisson solver effectively relaxes magnetic fields extrapolated from solar magnetograms, producing scalar potential fields that align well with theoretical expectations. Extensive verification was conducted on the high-order FR solver for polynomial orders up to P 3, achieving fourth-order spatial accuracy. The hyperbolized solver demonstrates comparable accuracy to reference solutions (both analytical and numerical) while offering efficient performance, particularly on coarser meshes, making it competitive with state-of-the-art low-order finite volume solvers, which are mostly used for solar MHD simulations. The described developments are a milestone for enabling high-order global solar corona simulations on 3D unstructured grids.
format Article
id doaj-art-61abdf2c8d6f44639462b4fd851653c3
institution Kabale University
issn 1538-4357
language English
publishDate 2025-01-01
publisher IOP Publishing
record_format Article
series The Astrophysical Journal
spelling doaj-art-61abdf2c8d6f44639462b4fd851653c32025-02-12T06:40:36ZengIOP PublishingThe Astrophysical Journal1538-43572025-01-01980216310.3847/1538-4357/adace5Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field InitializationRayan Dhib0https://orcid.org/0000-0003-3670-4678Firas Ben Ameur1https://orcid.org/0000-0001-5097-2885Vatsalya Sharma2https://orcid.org/0000-0003-2781-6563Andrea Lani3https://orcid.org/0000-0003-4017-215XStefaan Poedts4https://orcid.org/0000-0002-1743-0651Centre for Mathematical Plasma-Astrophysics , KU Leuven, Celestijnenlaan 200-B, 3001 Leuven, Belgium ; rayan.dhib@kuleuven.be, vatsalya.sharma@kuleuven.be, andrea.lani@kuleuven.be, stefaan.poedts@kuleuven.beKing Abdullah University of Science and Technology , KSA-23955 Thuwal, Saudi Arabia ; firas.benameur@kaust.edu.saCentre for Mathematical Plasma-Astrophysics , KU Leuven, Celestijnenlaan 200-B, 3001 Leuven, Belgium ; rayan.dhib@kuleuven.be, vatsalya.sharma@kuleuven.be, andrea.lani@kuleuven.be, stefaan.poedts@kuleuven.beCentre for Mathematical Plasma-Astrophysics , KU Leuven, Celestijnenlaan 200-B, 3001 Leuven, Belgium ; rayan.dhib@kuleuven.be, vatsalya.sharma@kuleuven.be, andrea.lani@kuleuven.be, stefaan.poedts@kuleuven.beCentre for Mathematical Plasma-Astrophysics , KU Leuven, Celestijnenlaan 200-B, 3001 Leuven, Belgium ; rayan.dhib@kuleuven.be, vatsalya.sharma@kuleuven.be, andrea.lani@kuleuven.be, stefaan.poedts@kuleuven.be; Institute of Physics, University of Maria Curie-Skłodowska , ul. Radziszewskiego 10, PL-20-031 Lublin, PolandProper initialization of the solar corona magnetic field is important for easing the iterative process of realistic and efficient global magnetohydrodynamics (MHD) simulations. This study introduces a high-order flux reconstruction (FR) framework for solving the Poisson equation, a necessary step for computing a potential-field source-surface (PFSS) to initialize the magnetic field for global solar corona simulations with MHD. By hyperbolizing the elliptic Poisson equation into a set of hyperbolic equations, we develop an efficient and robust high-order PFSS solver. Our contributions include developing a Q 2 (i.e., quadratic) geometrical representation using prismatic elements for the computational domain, which enables high-order mesh generation. Such a hyperbolized Poisson solver effectively relaxes magnetic fields extrapolated from solar magnetograms, producing scalar potential fields that align well with theoretical expectations. Extensive verification was conducted on the high-order FR solver for polynomial orders up to P 3, achieving fourth-order spatial accuracy. The hyperbolized solver demonstrates comparable accuracy to reference solutions (both analytical and numerical) while offering efficient performance, particularly on coarser meshes, making it competitive with state-of-the-art low-order finite volume solvers, which are mostly used for solar MHD simulations. The described developments are a milestone for enabling high-order global solar corona simulations on 3D unstructured grids.https://doi.org/10.3847/1538-4357/adace5Computational methodsSolar coronaSpace weather
spellingShingle Rayan Dhib
Firas Ben Ameur
Vatsalya Sharma
Andrea Lani
Stefaan Poedts
Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization
The Astrophysical Journal
Computational methods
Solar corona
Space weather
title Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization
title_full Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization
title_fullStr Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization
title_full_unstemmed Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization
title_short Toward High-order Solar Corona Simulations: A High-order Hyperbolized Poisson Approach for Magnetic Field Initialization
title_sort toward high order solar corona simulations a high order hyperbolized poisson approach for magnetic field initialization
topic Computational methods
Solar corona
Space weather
url https://doi.org/10.3847/1538-4357/adace5
work_keys_str_mv AT rayandhib towardhighordersolarcoronasimulationsahighorderhyperbolizedpoissonapproachformagneticfieldinitialization
AT firasbenameur towardhighordersolarcoronasimulationsahighorderhyperbolizedpoissonapproachformagneticfieldinitialization
AT vatsalyasharma towardhighordersolarcoronasimulationsahighorderhyperbolizedpoissonapproachformagneticfieldinitialization
AT andrealani towardhighordersolarcoronasimulationsahighorderhyperbolizedpoissonapproachformagneticfieldinitialization
AT stefaanpoedts towardhighordersolarcoronasimulationsahighorderhyperbolizedpoissonapproachformagneticfieldinitialization