Computational modelling of shock wave/boundary layer interaction with a cell-vertex scheme and transport models of turbulence

Abstract A calculation procedure for modelling the interaction between shock waves and attached or separated turbulent boundary layers is introduced. The numerical framework, applicable to general curved grids, combines cell-vertex storage, a Lax-Wendroff time-marching scheme and multigrid convergence acceleration. The main numerical ingredients of the procedure are documented in some detail, with particular emphasis placed on the inclusion of viscous and turbulence transport within the cell-vertex framework, which was originally formulated for inviscid flow. An investigation of the predictive performance of alternative transport models of turbulence has been the primary objective of the present work. Particular attention has been focused on a comparison between variants of low Reynolds number k-ε models and an algebraic variant of a Reynolds-stress transport closure in strong interaction situations, including shock-induced separation. The turbulence models are introduced, and important numerical issues affecting their stable implementation are discussed. The calculation procedure is then applied to two confined transonic flows over bumps — one incipiently and the other extensively separated (Delery Cases A and C) — and to the transonic flow around the RAE 2822 aerofoil at two angles of incidence (Cases 9 and 10). The investigation demonstrates that the eddy-viscosity models tend to seriously underestimate the strength of interaction, particularly when separation is extensive. The performance of the Reynolds-stress model is not entirely consistent across the range of conditions examined. In the case of bump flows, the model displays strong sensitivity to the shock, predicting excessive boundary layer displacement in Case A, a broadly correct separation pattern in Case C and insufficient rate of post shock recovery in both cases. The aerofoil flows are either attached or incipiently separated, and the benefits arising from Reynolds-stress modelling are modest. Neither the k-ε model nor the Reynolds-stress closure is able to return a satisfactory representation of the most challenging RAE 2822 Case 10; at least not with the recommended windtunnel corrections for freestream Mach number and angle of incidence.

[1]  T. J. Coakley,et al.  An implicit Navier-Stokes code for turbulent flow modeling , 1992 .

[2]  M. Wolfshtein.  The velocity and temperature distribution in one-dimensional flow with turbulence augmentation and pressure gradient , 1969 .

[3]  K. W. Morton,et al.  Finite volume solutions of convection-diffusion test problems , 1993 .

[4]  Brian Launder,et al.  Second-moment closure: present… and future? , 1989 .

[5]  K. Chien,et al.  Predictions of Channel and Boundary-Layer Flows with a Low-Reynolds-Number Turbulence Model , 1982 .

[6]  Antony Jameson,et al.  Numerical solution of the Euler equation for compressible inviscid fluids , 1985 .

[7]  W. Jones,et al.  The prediction of laminarization with a two-equation model of turbulence , 1972 .

[8]  N. Ron-Ho,et al.  A Multiple-Grid Scheme for Solving the Euler Equations , 1982 .

[9]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[10]  Brian Launder,et al.  Numerical methods in laminar and turbulent flow , 1983 .

[11]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[12]  Jubaraj Sahu,et al.  Navier-Stokes computations of transonic flows with a two-equation turbulence model , 1986 .

[13]  K. W. Morton,et al.  Numerical methods for fluid dynamics , 1987 .

[14]  W. Rodi.  A new algebraic relation for calculating the Reynolds stresses , 1976 .

[15]  Implicit Navier-Stokes Calculations of Transonic Shock/Turbulent Boundary-Layer Interactions , 1986 .

[16]  Thomas J. R. Hughes,et al.  Numerical Methods in Laminar and Turbulent Flows , 1979 .

[17]  B. Launder,et al.  THE NUMERICAL COMPUTATION OF TURBULENT FLOW , 1974 .

[18]  Jean Delery,et al.  A Study of Turbulence Modelling in Transonic Shock-Wave Boundary-Layer Interactions , 1989 .

[19]  Michael A. Leschziner,et al.  Modelling engineering flows with Reynolds stress turbulence closure , 1990 .

[20]  Budugur Lakshminarayana,et al.  Turbulence modeling for complex shear flows , 1985 .

[21]  M. Leschziner,et al.  Multilevel convergence acceleration for viscous and turbulent transonic flows computed with a cell-vertex method , 1989 .

[22]  D. Wilcox.  Reassessment of the scale-determining equation for advanced turbulence models , 1988 .

[23]  B. Launder,et al.  Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc , 1974 .

[24]  F. Lien,et al.  Modelling shock/turbulent-boundary-layer interaction with second-moment closure within a pressure-velocity strategy , 1993 .

[25]  Jean-Antoine Désidéri,et al.  Numerical methods for the Euler equations of fluid dynamics , 1985 .

[26]  A. J. Peace,et al.  Turbulent flow predictions for afterbody/nozzle geometries includingbase effects , 1989 .

[27]  M. Nallasamy,et al.  Turbulence models and their applications to the prediction of internal flows: a review , 1987 .

[28]  B. Launder,et al.  Ground effects on pressure fluctuations in the atmospheric boundary layer , 1978, Journal of Fluid Mechanics.

[29]  C. C. Horstman,et al.  On the use of wall functions as boundary conditions for two-dimensional separated compressible flows , 1985 .

[30]  A. Sugavanam.  Evaluation of low Reynolds number turbulence models for attached and separated flows , 1985 .

[31]  K. W. Morton,et al.  A finite volume scheme with shock fitting for the steady euler equations , 1989 .

[32]  J. Délery.  Turbulent Shear-Layer/Shock-Wave Interactions , 1986 .