A new high-performance computation method of the Eddington tensor in radiation hydrodynamics simulation
Abstract:
Radiative hydrodynamics models the coupling between the dynamics of a hypersonic hot plasma and the radiation it produces or external radiation. Almost every numerical code uses simplified models, which are often either limited or inaccurate. To accurately model photon transport, the HADES 2D code was specifically developed. Such a code is indispensable for studying astrophysical objects, in which optically intermediate regions are still poorly modeled but commonly encountered within such phenomena. This code couples hydrodynamics with the M1-multigroup model for radiation transfer to accurately represent the spectral behavior of light, involving the partitioning of the electromagnetic spectrum into groups. However, simulating radiative hydrodynamics flows remains highly time-consuming, constraining our capacity to conduct comprehensive numerical studies within this field.
The most computationally expensive part of the M1-multigroup simulations is the calculation of the closure relation, which relates the radiative pressure to the radiative energy and the radiative flux via the Eddington factor. This is due to the lack of an analytical solution. Consequently, two methods exist: One method is accurate but costly, relying on expensive search algorithms implemented in HADES. Another method is quicker but incorrect, utilizing the analytical grey case closure relation for each group, implemented in HERACLES.
To mitigate these challenges, we've pioneered an inventive approach intertwining neural networks with simplified models. This innovative method dramatically reduces computation time while maintaining an acceptable precision, revolutionizing the efficiency of these calculations within M1-multigroup simulations.
To affirm the efficiency of our approach, we conducted validation simulations, beginning with the renowned benchmark simulation of a 1D radiative shock, wherein we used up to five groups. Additionally, we undertook a radial test to assess the efficiency of our method in a 2D situation.