Stable, entropy-pressure compatible subsonic Riemann boundary condition for embedded DG compressible flow simulations

One approach to reducing the computational cost of simulating transitional compressible boundary layer flow is to adopt a near body reduced domain with boundary conditions enforced to be compatible with a computationally cheaper three-dimensional RANS simulation. In such an approach it is desirable to enforce a consistent pressure distribution which is not typically the case when using the standard Riemann inflow boundary conditions. We therefore revisit the Riemann problem adopted in many DG based high fidelity formulations. Through analysis of the one-dimensional linearised Euler equations it is demonstrated that maintaining entropy compatibility with the RANS simulation is important for a stable solution. It is also necessary to maintain the invariant for the Riemann outflow boundary condition in a subsonic flow leav-ing one condition that can be imposed at the inflow boundary. Therefore the entropy-pressure enforcement is the only stable boundary condition to enforce a known pressure distribution. We further demonstrate that all the entropy compatible inflow Riemann boundary conditions are stable providing the invariant compatible Riemann outflow boundary condition is also respected. Although the entropy-pressure compatible Riemann inflow boundary condition is stable from the one-dimensional analysis, two-dimensional tests highlight divergence in the inviscid problem and neutrally stable wiggles in the velocity fields in viscous simulations around the stagnation point. A two-dimensional analysis about a non-uniform baseflow assumption provides insight into this stability issue (ill-posedness) and motivate the use of a mix of inflow boundary conditions in this region of the flow. As a validation we apply the proposed boundary conditions to a reduced domain of a wing section normal to the leading-edge of the CRM-NLF model taken out of a full three-dimensional RANS simulation at Mach 0.86 and a Reynolds number of 8.5 million. The results show that the entropy-pressure compatible Riemann inflow boundary condition leads to a good agreement in pressure distribution.

[1]  Spencer J. Sherwin,et al.  Industry-Relevant Implicit Large-Eddy Simulation of a High-Performance Road Car via Spectral/hp Element Methods , 2020, SIAM Rev..

[2]  Emma E. Cooke Modelling the effect of step and roughness features on swept wing boundary layer instabilities , 2020 .

[3]  Robert Michael Kirby,et al.  Nektar++: enhancing the capability and application of high-fidelity spectral/hp element methods , 2019, Comput. Phys. Commun..

[4]  Florian J. Hindenlang,et al.  Stability of Wall Boundary Condition Procedures for Discontinuous Galerkin Spectral Element Approximations of the Compressible Euler Equations , 2019, Lecture Notes in Computational Science and Engineering.

[5]  U. Rist,et al.  Self-similar compressible turbulent boundary layers with pressure gradients. Part 1. Direct numerical simulation and assessment of Morkovin’s hypothesis , 2019, Journal of Fluid Mechanics.

[6]  M. S. Mughal,et al.  Destabilisation of Stationary and Travelling Crossflow Disturbances Due to Steps over a Swept Wing , 2019, AIAA Aviation 2019 Forum.

[7]  Dan S. Henningson,et al.  Control of Instabilities in an Unswept Wing Boundary Layer , 2018 .

[8]  Jan Nordström,et al.  A Roadmap to Well Posed and Stable Problems in Computational Physics , 2016, Journal of Scientific Computing.

[9]  Dan S. Henningson,et al.  Direct numerical simulation of the flow around a wing section at moderate Reynolds number , 2016 .

[10]  Robert Michael Kirby,et al.  Nektar++: An open-source spectral/hp element framework , 2015, Comput. Phys. Commun..

[11]  Spencer J. Sherwin,et al.  A Guide to the Implementation of Boundary Conditions in Compact High-Order Methods for Compressible Aerodynamics , 2014 .

[12]  D. Henningson,et al.  Swept wing boundary-layer receptivity to localized surface roughness , 2012, Journal of Fluid Mechanics.

[13]  Jan Nordström,et al.  Weak and strong wall boundary procedures and convergence to steady-state of the Navier-Stokes equations , 2012, J. Comput. Phys..

[14]  Jan Nordström,et al.  Weak Versus Strong No-Slip Boundary Conditions for the Navier-Stokes Equations , 2010 .

[15]  Jan Nordström,et al.  The Influence of Weak and Strong Solid Wall Boundary Conditions on the Convergence to Steady-State of the Navier-Stokes Equations , 2009 .

[16]  N. SIAMJ.,et al.  WELL-POSED BOUNDARY CONDITIONS FOR THE NAVIER – STOKES EQUATIONS , 2005 .

[17]  W. Saric,et al.  STABILITY AND TRANSITION OF THREE-DIMENSIONAL BOUNDARY LAYERS , 2003 .

[18]  David L. Darmofal,et al.  Eigenmode Analysis of Boundary Conditions for the One-Dimensional Preconditioned Euler Equations , 1999 .

[19]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[20]  D. Arnal,et al.  Linear Stability Theory Applied to Boundary Layers , 1996 .

[21]  M. Giles Eigenmode analysis of unsteady one-dimensional Euler equations , 1983 .