On the Development of an Implicit Discontinuous Galerkin Solver for Turbulent Real Gas Flows

• Published in: Fluids (Basel) Vol. 8; no. 4; p. 117
• Main Authors: Mantecca, Edoardo, Colombo, Alessandro, Ghidoni, Antonio, Noventa, Gianmaria, Pasquale, David, Rebay, Stefano
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01.04.2023
• Subjects: Accuracy; Boundary conditions; discontinuous Galerkin; equation of state; Finite element method; Finite volume method; Flow velocity; Fluid dynamics; Galerkin method; Gas flow; Jacobi matrix method; Jacobian matrix; Linearization; Mathematical analysis; NICFD; Numerical analysis; Polynomials; real gas; Real gases; Riemann solver; Temperature dependence; Thermodynamic models; Thermodynamics; Time integration; Transport properties; Turbulence models; Italy
• ISSN: 2311-5521; 2311-5521
• DOI: 10.3390/fluids8040117

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.