Utilizza questo identificativo per citare o creare un link a questo documento:
http://hdl.handle.net/20.500.12386/23591
Titolo: | Astrophysical fluid simulations of thermally ideal gases with non-constant adiabatic index: numerical implementation | Autori: | Vaidya, B. Mignone, A. BODO, Gianluigi Massaglia, S. |
Data pubblicazione: | 2015 | Rivista: | ASTRONOMY & ASTROPHYSICS | Numero: | 580 | Da pagina:: | A110 | Abstract: | Context. An equation of state (EoS) is a relation between thermodynamic state variables and it is essential for closing the set of equations describing a fluid system. Although an ideal EoS with a constant adiabatic index Γ is the preferred choice owing to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Aims: In the present work we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation, and temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation-equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes. Methods: We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS. The first employs root-finder methods and it is best suited for EoS in analytical form. The second is based on lookup tables and interpolation and results in a more computationally efficient approach, although care must be taken to ensure thermodynamic consistency. Results: A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is 500-10<SUP>4</SUP> K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although the employment of lookup table methods (when possible) can significantly reduce the overhead by a factor of ~ 3-4. | Acknowledgments: | We would like to thank the referee, Prof. A. Raga, for his constructive suggestions on the paper. B.V. is grateful to the funding support from the University of Torino under the contract: “Progetto di Ateneo-Compagnia di SanPaolo”. | URI: | http://hdl.handle.net/20.500.12386/23591 | URL: | https://www.aanda.org/articles/aa/abs/2015/08/aa26247-15/aa26247-15.html | ISSN: | 0004-6361 | DOI: | 10.1051/0004-6361/201526247 | Bibcode ADS: | 2015A&A...580A.110V | Fulltext: | open |
È visualizzato nelle collezioni: | 1.01 Articoli in rivista |
File in questo documento:
File | Descrizione | Dimensioni | Formato | |
---|---|---|---|---|
SM_paper_AandA_format.pdf | postprint | 1.37 MB | Adobe PDF | Visualizza/apri |
aa26247-15.pdf | pdf editoriale | 1.56 MB | Adobe PDF | Visualizza/apri |
Page view(s)
56
controllato il 29-mar-2024
Download(s)
23
controllato il 29-mar-2024
Google ScholarTM
Check
Altmetric
Altmetric
Tutti i documenti in DSpace sono pubblicati ad Accesso Aperto, salvo diversa indicazione per alcuni documenti specifici.