The aim of this work is to describe an efficient implementation of cubic and multiparameter real gas models in an existing discontinuous Galerkin solver to extend its capabilities to the simulation of turbulent real gas flows. The adopted thermodynamic models are van der Waals, Peng–Robinson, and Span–Wagner, which differ from each other in terms of accuracy and computational cost. Convective numerical fluxes across elements interfaces are calculated with a thermodynamic consistent linearized Riemann solver, whereas for boundary conditions, a linearized expression of the generalized Riemann invariants is employed. Transport properties are treated as temperature- and density-dependent quantities through multiparameter correlations. An implicit time integration is adopted; Jacobian matrix and thermodynamic derivatives are obtained with the automatic differentiation tool Tapenade. The solver accuracy is assessed by computing both steady and unsteady real gas test cases available in the literature, and the effect of the mesh size and polynomial degree of approximation on the solution accuracy is investigated. A good agreement with experimental and numerical reference data is observed and specific non-classical phenomena are well reproduced by the solver.
(2023). On the Development of an Implicit Discontinuous Galerkin Solver for Turbulent Real Gas Flows [journal article - articolo]. In FLUIDS. Retrieved from https://hdl.handle.net/10446/243630
On the Development of an Implicit Discontinuous Galerkin Solver for Turbulent Real Gas Flows
Mantecca, Edoardo;Colombo, Alessandro;Ghidoni, Antonio;Rebay, Stefano
2023-01-01
Abstract
The aim of this work is to describe an efficient implementation of cubic and multiparameter real gas models in an existing discontinuous Galerkin solver to extend its capabilities to the simulation of turbulent real gas flows. The adopted thermodynamic models are van der Waals, Peng–Robinson, and Span–Wagner, which differ from each other in terms of accuracy and computational cost. Convective numerical fluxes across elements interfaces are calculated with a thermodynamic consistent linearized Riemann solver, whereas for boundary conditions, a linearized expression of the generalized Riemann invariants is employed. Transport properties are treated as temperature- and density-dependent quantities through multiparameter correlations. An implicit time integration is adopted; Jacobian matrix and thermodynamic derivatives are obtained with the automatic differentiation tool Tapenade. The solver accuracy is assessed by computing both steady and unsteady real gas test cases available in the literature, and the effect of the mesh size and polynomial degree of approximation on the solution accuracy is investigated. A good agreement with experimental and numerical reference data is observed and specific non-classical phenomena are well reproduced by the solver.File | Dimensione del file | Formato | |
---|---|---|---|
fluids-08-00117.pdf
accesso aperto
Versione:
publisher's version - versione editoriale
Licenza:
Creative commons
Dimensione del file
6.1 MB
Formato
Adobe PDF
|
6.1 MB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
Aisberg ©2008 Servizi bibliotecari, Università degli studi di Bergamo | Terms of use/Condizioni di utilizzo