Over the past two decades, there has been much development in discontinuous Galerkin methods for incompressible flows and for compressible flows with a positive Mach number, but almost no attention has been paid to variable-density flows at low speeds. This paper presents a pressure-based discontinuous Galerkin method for flow in the low-Mach number limit. We use a variable-density pressure correction method, which is simplified by solving for the mass flux instead of the velocity. The fluid properties do not depend significantly on the pressure, but may vary strongly in space and time as a function of the temperature. We pay particular attention to the temporal discretization of the enthalpy equation, and show that the specific enthalpy needs to be ‘offset’ with a constant in order for the temporal finite difference method to be stable. We also show how one can solve for the specific enthalpy from the conservative enthalpy transport equation without needing a predictor step for the density. These findings do not depend on the spatial discretization. A series of manufactured solutions with variable fluid properties demonstrate full second-order temporal accuracy, without iterating the transport equations within a time step. We also simulate a Von Kármán vortex street in the wake of a heated circular cylinder, and show good agreement between our numerical results and experimental data.

Accurate methods to solve the Reynolds-Averaged Navier-Stokes (RANS) equations coupled to turbulence models are still of great interest, as this is often the only computationally feasible approach to simulate complex turbulent flows in large engineering applications. In this work, we present a novel discontinuous Galerkin (DG) solver for the RANS equations coupled to the k−ϵ model (in logarithmic form, to ensure positivity of the turbulence quantities). We investigate the possibility of modeling walls with a wall function approach in combination with DG. The solver features an algebraic pressure correction scheme to solve the coupled RANS system, implicit backward differentiation formulae for time discretization, and adopts the Symmetric Interior Penalty method and the Lax-Friedrichs flux to discretize diffusive and convective terms respectively. We pay special attention to the choice of polynomial order for any transported scalar quantity and show it has to be the same as the pressure order to avoid numerical instability. A manufactured solution is used to verify that the solution converges with the expected order of accuracy in space and time. We then simulate a stationary flow over a backward-facing step and a Von Kármán vortex street in the wake of a square cylinder to validate our approach.

We present a new discretization of the mono-energetic Fokker–Planck equation. We build on previous work (Kópházi and Lathouwers, 2015) where we devised an angular discretization for the Boltzmann equation, allowing for both heterogeneous and anisotropic angular refinement. The angular discretization is based on a discontinuous finite element method on the unit sphere. Here we extend the methodology to include the effect of the Fokker–Planck scatter operator describing small angle particle scatter. We describe the construction of an interior penalty method on the sphere surface. Results are provided for a variety of test cases, ranging from purely angular to fully three-dimensional. The results show that the scheme can resolve highly forward-peaked flux distributions with forward-peaked scatter.