Modelling solar coronal magnetic fields with physics-informed neural networks

Hubert Baty,Vincent Vigon
2023-10-27
Abstract:We present a novel numerical approach aiming at computing equilibria and dynamics structures of magnetized plasmas in coronal environments. A technique based on the use of neural networks that integrates the partial differential equations of the model, and called Physics-Informed Neural Networks (PINNs), is introduced. The functionality of PINNs is explored via calculation of different magnetohydrodynamic (MHD) equilibrium configurations, and also obtention of exact two-dimensional steady-state magnetic reconnection solutions (Craig & Henton 1995). Advantages and drawbacks of PINNs compared to traditional numerical codes are discussed in order to propose future improvements. Interestingly, PINNs is a meshfree method in which the obtained solution and associated different order derivatives are quasi-instantaneously generated at any point of the spatial domain. We believe that our results can help to pave the way for future developments of time dependent MHD codes based on PINNs
Solar and Stellar Astrophysics,Plasma Physics,Space Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to use Physics - Informed Neural Networks (PINNs) to model the equilibrium state and dynamic structure of magnetized plasma in the solar corona environment. Specifically, the author explores the functions of PINNs in calculating different Magnetohydrodynamic (MHD) equilibrium configurations and obtaining accurate two - dimensional steady - state magnetic field reconnection solutions. The paper discusses the advantages and disadvantages of PINNs compared with traditional numerical methods and proposes directions for future improvement. As a mesh - free method, PINNs can quickly generate solutions and their derivatives of all orders at any point in the spatial domain, which lays the foundation for the future development of PINNs - based time - dependent MHD codes. ### Key issues: 1. **Calculating the equilibrium state of magnetized plasma in the solar corona**: - The paper first explores the use of PINNs to calculate two - dimensional force - free magnetic equilibria, which represent arch and loop structures in the corona. By solving partial differential equations similar to the Grad - Shafranov equation, the author demonstrates the application of PINNs in this area. 2. **Calculating accurate two - dimensional steady - state magnetic field reconnection solutions**: - The author further extends the application of PINNs to solve a more complex set of partial differential equations, namely the incompressible resistive MHD equations, to calculate two - dimensional magnetic field reconnection solutions. Particular attention is paid to reconnective annihilation solutions, which are specific accurate steady - state solutions proposed by Craig and Henton in 1995. ### Specific problem descriptions: - **Force - free magnetic equilibrium**: - By solving the linear force - free equilibrium equation \(\Delta\psi + c^{2}\psi = 0\), the author calculates different types of magnetic field structures, such as dipole, quadrupole, and hybrid dipole/quadrupole configurations. These structures are of great observational significance in the solar corona, especially in the study of solar flares and coronal mass ejections (CME). - **Magnetic field reconnection**: - By solving the incompressible resistive MHD equations, the author calculates pure annihilation solutions (\(\beta = 0\)) and reconnective diffusion solutions (\(\beta = 0.5\)). These solutions describe the process in which two anti - parallel magnetic field regions annihilate each other under the action of an incompressible plasma flow, as well as the accompanying reconnective diffusion process. ### Methods and techniques: - **Basic principles of PINNs**: - PINNs solve partial differential equations by minimizing the values of equation residuals on pre - defined data points (called collocation points). This method combines data loss functions and physical loss functions and can effectively solve complex nonlinear equations. - **Network architecture**: - The author selects a feed - forward neural network structure with multiple hidden layers and neurons. For different problems, the complexity of the network (the number of layers and the number of neurons in each layer) varies to meet the needs of different problems. ### Conclusions: - The author believes that PINNs perform well in calculating the equilibrium state and dynamic structure of magnetized plasma in the solar corona, especially in terms of its mesh - free characteristics and the ability to quickly generate solutions. However, compared with traditional numerical methods, PINNs have slightly lower accuracy in some cases, but the quality of the results can be further improved by optimizing network parameters and loss function weights. Future work will focus on developing PINNs - based time - dependent MHD codes to better simulate dynamic processes in the solar corona.