Context. Modeling the dynamics of most astrophysical structures requires an adequate description of the interaction of radiation and matter. Several numerical (magneto-) hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those that used either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to local thermodynamic equilibrium (LTE). This assumption may not be valid in some astrophysical cases. Aims. We present an upgraded version of the LTE radiation-hydrodynamics (RHD) module implemented in the PLUTO code, which we have extended to handle non-LTE regimes. Methods. Starting from the general frequency-integrated comoving-frame equations of RHD, we have justified all the assumptions that were made to obtain the non-LTE equations that are implemented in the module under the FLD approximation. An operator-split method with two substeps was employed: The hydrodynamics part was solved with an explicit method by the solvers that are currently available in PLUTO, and the non-LTE radiation diffusion and energy exchange part was solved with an implicit method. The module was implemented in the PLUTO environment. It uses databases of radiative quantities that can be provided independently by the user: The radiative power loss, and the Planck and Rosseland mean opacities. In our case, these quantities were determined from a collisional-radiative steady-state model, and they are tabulated as functions of temperature and density. Results. Our implementation has been validated through different tests, in particular, radiative shock tests. The agreement with the semi-analytical solutions (when available) is good, with a maximum error of 7%. Moreover, we have proved that a non-LTE approach is of paramount importance to properly model accretion shock structures. Conclusion. Our radiation FLD module represents a step toward a general non-LTE RHD modeling.

Colombo S., Colombo S., Orlando S., Rodriguez R., Espinosa G., Gonzalez M., et al. (2019). Non-LTE radiation hydrodynamics in PLUTO. ASTRONOMY & ASTROPHYSICS, 631, A41 [10.1051/0004-6361/201935991].

Non-LTE radiation hydrodynamics in PLUTO

Colombo S.
;
Colombo S.
;
Peres G.
Membro del Collaboration Group
2019-01-01

Abstract

Context. Modeling the dynamics of most astrophysical structures requires an adequate description of the interaction of radiation and matter. Several numerical (magneto-) hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those that used either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to local thermodynamic equilibrium (LTE). This assumption may not be valid in some astrophysical cases. Aims. We present an upgraded version of the LTE radiation-hydrodynamics (RHD) module implemented in the PLUTO code, which we have extended to handle non-LTE regimes. Methods. Starting from the general frequency-integrated comoving-frame equations of RHD, we have justified all the assumptions that were made to obtain the non-LTE equations that are implemented in the module under the FLD approximation. An operator-split method with two substeps was employed: The hydrodynamics part was solved with an explicit method by the solvers that are currently available in PLUTO, and the non-LTE radiation diffusion and energy exchange part was solved with an implicit method. The module was implemented in the PLUTO environment. It uses databases of radiative quantities that can be provided independently by the user: The radiative power loss, and the Planck and Rosseland mean opacities. In our case, these quantities were determined from a collisional-radiative steady-state model, and they are tabulated as functions of temperature and density. Results. Our implementation has been validated through different tests, in particular, radiative shock tests. The agreement with the semi-analytical solutions (when available) is good, with a maximum error of 7%. Moreover, we have proved that a non-LTE approach is of paramount importance to properly model accretion shock structures. Conclusion. Our radiation FLD module represents a step toward a general non-LTE RHD modeling.
2019
Colombo S., Colombo S., Orlando S., Rodriguez R., Espinosa G., Gonzalez M., et al. (2019). Non-LTE radiation hydrodynamics in PLUTO. ASTRONOMY & ASTROPHYSICS, 631, A41 [10.1051/0004-6361/201935991].
File in questo prodotto:
File Dimensione Formato  
aa35991-19.pdf

accesso aperto

Descrizione: articolo principale
Tipologia: Versione Editoriale
Dimensione 568.34 kB
Formato Adobe PDF
568.34 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10447/389732
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 4
  • ???jsp.display-item.citation.isi??? 4
social impact