Computation of shock wave/turbulent boundary layer interactions using a two-equation model with compressibility corrections

In this paper computational results for two different types of shock wave / turbulent boundary layer interaction flows are presented. It is shown that upstream effects of the shock induced separation cannot be reproduced by Wilcox's (1991) k--model, whereas downstream of the interaction, predictions of pressure distribution and skin friction are acceptable. The inclusion of the compressible part of the dissipation rate and the pressure dilatation in the model has noticeable, but not dramatic effects on wall pressure and skin friction in the selected flow cases.     

Investigation of flux formulae in transonic shock wave/turbulent boundary layer interaction

The aim of the present study is to examine the accuracy and improvement of various numerical methods in the solution of the transonic shock/turbulent boundary layer interaction problem and to show that a significant source of numerical inaccuracies in turbulent flows is not only the inadequacy of the turbulence model but also the numerical discretization. Comparisons between a Riemann solver and a flux-vector-splitting method as well as between various numerical high-order extrapolation schemes with corresponding experimental results are presented.     

Investigation of a hypersonic crossing shock wave/turbulent boundary layer interaction

A combined theoretical and experimental study is presented for the interaction between crossing shock waves generated by (10°, 10°) sharp fins and a flat plate turbulent boundary layer at Mach 8.3. The theoretical model is the full 3-D mean compressible Reynolds-averaged Navier-Stokes RANS) equations incorporating the algebraic turbulent eddy viscosity model of Baldwin and Lomax. A grid refinement study indicated that adequate resolution of the flowfield has been achieved. Computed results agree well with experiment for surface pressure and surface flow patterns and for pitot pressure and yaw angle profiles in the flowfield. The computations, however, significantly overpredict surface heat transfer. Analysis of the computed flowfield results indicates the formation of complex streamline and wave structures within the interaction region.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

The structures of turbulent boundary layer in the vicinity of separation

Large coherent structures of turbulent boundary layer in the vicinity of separation were observed in a water channel by the hydrogen bubble method. Motion pictures of the de views were taken. The features of the instantaneous velocity profiles, the large transverse and streamwise vortices were discussed.     

FLAT-PLATE BOUNDARY-LAYER FLOWS INDUCED BY DUSTY SHOCK WAVE

ＦＬＡＴ－ＰＬＡＴＥＢＯＵＮＤＡＲＹ－ＬＡＹＥＲＦＬＯＷＳＩＮＤＵＣＥＤＢＹＤＵＳＴＹＳＨＯＣＫＷＡＶＥ（王柏懿）（陶锋）ＦＬＡＴ－ＰＬＡＴＥＢＯＵＮＤＡＲＹ－ＬＡＹＥＲＦＬＯＷＳＩＮＤＵＣＥＤＢＹＤＵＳＴＹＳＨＯＣＫＷＡＶＥ￥ＷａｎｇＢｏｙｉ；Ｔａ...     

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0∼4.7)×10⁷/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30° for a sweepback angle of 67.6°. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10° and secondary separation were detected at deflection of ϑ≥20°. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements atM≥6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers. The project supported by China Academy of Launch Vehicle Technology     

Near-wall damping in model predictions of separated flows

Different near-wall scalings are reviewed by the use of data from direct numerical simulations (DNS) of attached and separated adverse pressure gradient turbulent boundary layers. The turbulent boundary layer equation is analysed in order to extend the validity of existing wall damping functions to turbulent boundary layers under severe adverse pressure gradients. A proposed near-wall scaling is based on local quantities and the wall distance, which makes it applicable for general computational fluid dynamics (CFD) methods. It was found to have a similar behaviour as the pressure-gradient corrected analytical y^* scaling and avoids the inconsistencies present in the y⁺ scaling. The performance of the model is illustrated by model computations using explicit algebraic Reynolds stress models with near-wall damping based on different scalings.     

Incipient separation in shock wave/boundary layer interactions as induced by sharp fin

The incipient separation induced by the shock wave/turbulent boundary layer interaction at the sharp fin is the subject of the present study. Existing theories for the prediction of incipient separation, such as those put forward by McCabe (1966) and Dou and Deng (1992), can thus far only predic the direction of surface streamline and tend to overpredict the incipient separation condition based on the Stanbrook’s criterion. In this paper, the incipient separation is first predicted with Dou and Deng (1992)’s theory and then compared with Lu and Settles’ (1990) experimental data. The physical mechanism of the incipient separation as induced by the shock wave/turbulent boundary layer interactions at sharp fin is explained via surface flow pattern analysis. Furthermore, the reason for the observed discrepancy between the predicted and experimental incipient separation conditions is clarified. It is found that when the wall-limiting streamlines behind the shock wave becomes aligned with one ray from the virtual origin as the strength of the shock wave increases, the incipient separation line is formed at which the wall-limiting streamline becomes perpendicular to the local pressure gradient. The formation of this incipient separation line is the beginning of the separation process. The effects of Reynolds number and Mach number on incipient separation are also discussed. Finally, a correlation for the correction of the incipient separation angle as predicted by the theory is also given.

---

     

Structure of supersonic turbulent flows in the vicinity of inclined backward-facing steps

Supersonic (M_∞ = 2−5) turbulent flows in the vicinity of a two-dimensional backward-facing step with an inclined leeward side are considered by methods of mathematical modeling. The wave structure of the flow with a varied angle of inclination of the leeward side of the step and a varied free-stream Mach number is considered. __________ Translated From Prikladnaya Mekhanika I Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 48–58, November–December, 2006.     

Lagrangian model on the turbulent motion of small solid particle in turbulent boundary layer flows

.Intr0ductionSurfaceerosionofmaterialbysolid-particleimpactisanimportantprobleminmultiphaseflowindustriaIdevicesandthecharacteristicsoftheparticIe'smotioninaturbulentboundarylayerflowisthebaseofthestudyofthematerialsurfaceerosion.Manycalculationmodelshave…     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A calculating method of shock wave oscillating frequency due to turbulent shear layer fluctuations in supersonic flow

One of the more severe fluctuating pressure environments encountered in supersonic orhypersonic flows is the shock wave oscillation driven by interaction of a shock wave withboundary layer.The high intensity oscillating shock wave may induce structure resonanceof a high speed vehicle.The research for the shock oscillation used to adopt empirical orsemiempirical methods because the phenomenon is very complex.In this paper atheoretical solution on shock oscillating frequency due to turbulent shear layer fluctuationshas been obtained from basic conservation equations.Moreover,we have attained theregularity of the frequency of oscillating shock varying with incoming flow Mach numbersM_∞and turning angleθ.The calculating results indicate excellent agreement withmeasurements.This paper has supplied a valuable analytical method to study aeroelasticproblems produced by shock wave oscillation.     

Conical flow breakdown in non-free shock-wave/boundary-layer interaction

The non-free interaction between a shock wave and the boundary layer on a swept plate set at incidence in the undisturbed flow is studied using different experimental methods including special laser techniques for visualizing supersonic conical gas flows. It is shown that under shock-layer conditions the non-free interaction can lead to conical flow breakdown before the incident shock reaches the leading edge of the plate.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 45–58. Original Russian Text Copyright © 2004 by Zubin and Ostapenko.     

On the influence of elevated surface temperatures on hypersonic shock wave/boundary layer interaction at a heated ramp model

Although important flow parameters as Mach number, Reynolds number and total enthalpy can be reproduced in most hypersonic experiments quite well, due to different surface temperature effects in wind tunnel and flight, scaling as well as specific flow properties of shock wave/boundary layer interactions are different. This especially holds for short-duration facilities like, e.g. shock tunnels where due to short running times the models remain more or less at ambient temperature. To overcome this shortcoming, an experimental study has been conducted using a preheatable ramp model with 15° ramp angle. This allowed us to adjust the surfaces to an arbitrary temperature just before the experiment started. Pressure and heat flux measurements clearly showed the effect of varying surface and free stream temperatures. These results are supported by schlieren pictures and infrared measurements. The comparison of the measurements with theoretical and numerical results shows a good agreement. Separation bubble scaling laws proposed by Katzer and Davis have been applied and partially confirmed using the local conditions of the boundary layer at separation.     

Wall temperature effects on shock unsteadiness in a reattaching compressible turbulent shear layer

In this study, the effect of heat transfer on the compressible turbulent shear layer and shockwave interaction in a scramjet has been investigated. To this end, highly resolved Large Eddy Simulations (LES) are performed to explore the effect of wall thermal conditions on the behavior of a reattaching free shear layer interacting with an oblique shock in compressible turbulent flows. Various wall-to-recovery temperature ratios are considered, and results are compared to the adiabatic wall. It is found that the wall temperature affects the reattachment location and the shock behavior in the interaction region. Furthermore, fluctuating heat flux exhibits a strong intermittent behavior with severe heat transfer compared to the mean, characterized by scattered spots. The distribution of the Stanton number shows a strong heat transfer and complex pattern within the interaction, with the maximum thermal (heat transfer rates) and dynamic loads (root-mean-square wall pressure) found for the case of the cold wall. The analysis of LES data reveals that the thermal boundary condition can significantly impact the wall pressure fluctuations level. The primary mechanism for changes in the flow unsteadiness due to the wall thermal condition is linked to the reattaching shear layer, which agrees with the compressible turbulent boundary layer theory.     

On the boundary layer/riblets interaction mechanisms and the prediction of turbulent drag reduction

Turbulent drag reduction experienced by ribletted surfaces is the result of both (1) the interaction between riblet peaks and the coherent structures that characterize turbulent near-wall flows, and (2) the laminar sublayer flow modifications caused by the riblet shape, which can balance, under appropriate conditions, the drag penalty due to the increased wetted surface. The latter “viscous” mechanism is investigated by means of an analytical model of the laminar sublayer, which removes geometrical restrictions and allows us to take into account “real” shapes of riblet contours, affected by manufacturing inaccuracies, and to compute even for such cases a parameter, called protrusion height, related to the longitudinal mean flow. By considering real geometries, riblet effectiveness is clearly shown to be related to the difference between the longitudinal and the transversal protrusion heights. A simple method for the prediction of the performances of ribletted surfaces is then devised. The predicted and measured drag reduction data, for different riblet geometries and flow characteristics, are in close agreement with each other. The soundness of the physical interpretation underlying this prediction method is consequently confirmed.     

Effect of the Turbulent Energy Gradient and the Presence of a Wall on the Turbulent Process of Momentum Transport

A mathematical model of turbulent transport processes is modified to make allowance for the turbulent energy gradient and the presence of walls. The modification consists in making the variance tensor in the Gaussian probability density distribution for the initial mole velocities anisotropic for nonzero turbulent energy gradient and a ratio of the turbulence scale to the distance from the wall of the order of unity. Formulas for the variance tensor components are derived and the empirical coefficients of these formulas are determined. The expression for the dimensionless turbulent friction stress is compared with experimental data for three boundary-layer-type flows, namely, in the wake of a cylinder, in the boundary layer on a flat plate, and in a channel with parallel walls.     

Aerodynamic heating in the region of shock and turbulent boundary layer interaction induced by a cylinder

Detailed distributions of heat flux in the region of shock wave and turbulent boundary layer interaction induced by a cylinder were measured in the shock tunnel. Oil flow patterns and Schlieren photographs were taken. Empirical relations were given for determining separation shock angle, peaks of heat flux and their locations on both cylinder leading edge and flat plate surface, and other characteristic parameters of the interaction region.     

Simulation of supersonic turbulent flow in the vicinity of an inclined backward-facing step

Large eddy simulation (LES) is a viable and powerful tool to analyse unsteady three-dimensional turbulent flows. In this article, the method of LES is used to compute a plane turbulent supersonic boundary layer subjected to different pressure gradients. The pressure gradients are generated by allowing the flow to pass in the vicinity of an expansion–compression ramp (inclined backward-facing step with leeward-face angle of 25°) for an upstream Mach number of 2.9. The inflow boundary condition is the main problem for all turbulent wall-bounded flows. An approach to solve this problem is to extract instantaneous velocities, temperature and density data from an auxiliary simulation (inflow generator). To generate an appropriate realistic inflow condition to the inflow generator itself the rescaling technique for compressible flows is used. In this method, Morkovin's hypothesis, in which the total temperature fluctuations are neglected compared with the static temperature fluctuations, is applied to rescale and generate the temperature profile at inlet. This technique was successfully developed and applied by the present author for an LES of subsonic three-dimensional boundary layer of a smooth curved ramp. The present LES results are compared with the available experimental data as well as numerical data. The positive impact of the rescaling formulation of the temperature is proven by the convincing agreement of the obtained results with the experimental data compared with published numerical work and sheds light on the quality of the developed compressible inflow generator.     

Comparison of kinetic and continuum approaches for simulation of shock wave/boundary layer interaction

The present paper, which is a collaboration between three different research groups, analyzes the efficiency of various numerical approaches to describe the complex problem of shock wave/boundary layer interaction. Computations were carried out based on a kinetic approach (Direct Simulation Monte Carlo method) and on two continuum approaches (Navier-Stokes equations and quasigasdynamic equations), which are validated by comparison with experimental results obtained in the R5Ch blowdown Hypersonic Wind Tunnel in ONERA. The influence of the slip boundary conditions for two continuum approaches are also studied. The results obtained by all models display the good prediction of the main structure of the flow and the levels of the flux coefficients are very close to those measured. The implementation of the slip boundary condition for the continuum approaches improves the agreement with the experimental data. Received 12 July 2001 / Accepted 24 May 2002 /Published online 4 December 2002 Correspondence to: D. Zeitoun (e-mail: David.Zeitoun@polytech.univ-mrs.fr) An abridged version of this paper was presented at the 23rd Int. Symposium on Shock Waves at Fort Worth, Texas, from July 22 to 27, 2001