Year: 2022
Author: R. Chauvin, S. Guisset, B. Manach-Perennou, L. Martaud
Communications in Computational Physics, Vol. 31 (2022), Iss. 1 : pp. 293–330
Abstract
A positivity-preserving, conservative and entropic numerical scheme is presented for the three-temperature grey diffusion radiation hydrodynamics model. More precisely, the dissipation matrices of the colocalized semi-Lagrangian scheme are defined in order to enforce the entropy production on each species (electron or ion) proportionally to its mass as prescribed in [34]. A reformulation of the model is then considered to enable the derivation of a robust convex combination based scheme. This yields the positivity-preserving property at each sub-iteration of the algorithm while the total energy conservation is reached at convergence. Numerous pure hydrodynamics and radiation hydrodynamics test cases are carried out to assess the accuracy of the method. The question of the stability of the scheme is also addressed. It is observed that the present numerical method is particularly robust.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.OA-2021-0059
Communications in Computational Physics, Vol. 31 (2022), Iss. 1 : pp. 293–330
Published online: 2022-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 38
Keywords: Colocalized Lagrangian scheme radiation hydrodynamics grey diffusion discrete entropy production plasma physics simulations.