Simulation techniques for spatially evolving instabilities in compressible flow over a flat plate

Abstract

Abstract-In this paper we present numerical techniques suitable for a direct numerical simulation in the spatial setting. We demonstrate the application to the simulation of compressible flat plate flow instabilities. We compare second and fourth order accurate spatial discretization schemes in combination with explicit multistage time stepping for the simulation of the 2D Navier-Stokes equations. We consider Mach numbers 0.5 and 4.5. In the vicinity of the outflow boundary, an efficient buffer domain treatment is introduced, which is suitable in conjunction with an explicit time integration scheme. This treatment requires only a short buffer domain to damp wave reflections at the outflow boundary. Results for the instability of Tollmien-Schlichting (T-S) waves are compared with two instability theories, linear stability theory (LST) and linear parabolized stability equations (PSE). The growth rates of T-S waves for parallel base flow at both Mach numbeis compare well with LST results. Moreover, the growth rates of T-S waves for nonparallel base flow compare well with results obtained by solving the PSE at Mach number 0.5. The second order discretization scheme requires, however, considerably higher grid resolution than the fourth order method to achieve accurate results. High amplitude disturbances were also considered to activate nonlinear terms. The nonlinearity strongly affects the form of the T-S waves and the growth rate of the disturbance:. The results obtained here support the use of these numerical techniques in flow simulations with increasing complexity such as flat plate flow simulations up to the turbulent regime and with separation regions in 3D. The results also encourage the use of perturbations derived from the compressible PSE as inlet perturbations for nonparallel flow. 0 1997

Investigation of Three-Dimensional Shock Wave/Turbulent-Boundary-Layer Interaction Initiated by a Single Fin

“…Alternatively, a time sequence of a fully developed turbulent flow at the same flow condition was generated by running another independent LES of a flat plate turbulent boundary layer with a small domain size along the homogeneous spanwise direction, but still large enough to accommodate the near-wall turbulent streaks and the outer-layer wave numbers. The generated flow slices were then introduced to the main LES via the inlet plane by replicating itself along the spanwise direction and using the hybrid supersonic/subsonic inflow condition [53] . A similar inflow turbulence generation method was adopted by Gao et al [54] for an LES of turbulence with a full 3D geometry and good agreement of the LES with measured data was shown.…”

Large Eddy Simulation of stratified flows over structures

“…At the outlet we applied the Neumann boundary conditions with zero derivative and a buffer layer where all fluctuations were damped [11]. At the lateral boundaries and at the top, the free slip boundary conditions were prescribed.…”

Towards a transparent boundary condition for compressible Navier–Stokes equations

“…In order to obtain a solution of the equations, a buffer zone adjacent to the computational domain in used. Inside this buffer, the mesh is stretched in the horizontal direction [23], and the solution progressively brought back to a blasius flow along the dihedral [29], which can be considered as a local solution. Consequently, no BC is needed on the (d) boundary.…”

“…For it, Street and Macaraeg [26] and Liu et al [27] for incompressible flows, and latter Pruett et al [28] for compressible ones use a parabolization procedure. Wasistho et al [29] bring the flow to a reference one, using a rather crude method. A second type of poblems is to make a vortex leaving the computational domain with as little acoustic reflections as possible.…”