Computation of the two-dimensional viscous incompressible fluid flow with separation at the nlet edge around a cascade

A complex flow consisting of an outer inviscid stream, a dead-water separation domain, and a boundary layer, which interact strongly, is formed in viscous fluid flows with separation at the streamlined profile with high Re numbers. Different jet and vortex models of separation flow are known for an inviscid fluid; numerical, asymptotic, and integral methods [1–3] are used for a viscous fluid. The plane, stationary, turbulent flow through a turbine cascade by a constant-density fluid without and with separation from the inlet edge of the profile and subsequent attachment of the stream to the profile (a short, slender separation domain) is considered in this paper.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 34–44, May–June, 1978.     

Method of calculation of interaction of wall flows

A flow of viscous compressible fluid in the neighborhood of the line of interaction of wall flows is considered. A method of calculating the line of interaction and the direction of the self-induced secondary flow is developed. Papers [1–3] are devoted to the simulation of a separation flow with singularities in the neighborhood of singular lines and points, where boundary-layer equations are invalid. However, the theories of local separation used at present have mainly been developed only for two-dimensional problems, while the models of viscous-inviscid interaction have restrictions in application for turbulent flows with developed separation. The interaction of three-dimensional wall turbulent flows is considered below. It is assumed that the thickness of the boundary layers and the scales of the interaction zones are small in comparison with the characteristic dimension of the system, while the line of discontinuity of the solutions of the three-dimensional boundary layer equations is the same as the line of interaction of the wall flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 53–59, March–April, 1987.The author is grateful to G. Yu. Stepanov and V. N. Ershov for their interest in my work and their valuable remarks.     

Boundary-layer separation on conical bodies

Some characteristics of the variation in the linear dimensions of the flow separation zones on conical bodies with expanding conical skirts and of variation of the pressure within these zones as a function of variation of the Mach number, Reynolds number, and intensity of the disturbance that causes the boundary layer separation are examined. Experiments were conducted in laminar, transitional, and turbulent flows in flow separation regions. The interaction of viscous and nearly inviscid flows is quite common. This phenomenon occurs in flow past a concave corner, when a compression shock impinges on a boundary layer, and in many other cases. The characteristics of this phenomenon in flow about two-dimensional bodies have been investigated experimentally in [1, 2] and other studies. Attempts have been made to analyze the interaction of compression shocks with the boundary layer theoretically. In “free” separated flows, when the points of separation and reattachment of the boundary layer are not fixed (for example, on a flat plate with a long wedge attached to it), theoretical studies are usually made within the framework of the boundary layer theory with use of the approximate integral methods [3, 4]. In this article we examine some results from studies of free separated flows on conical bodies with conical skirts in laminar, transitional, and turbulent flows (Fig. 1).     

Transonic viscous-inviscid interaction by a finite element method

A method is outlined for solving two-dimensional transonic viscous flow problems, in which the velocity vector is split into the gradient of a potential and a rotational component. The approach takes advantage of the fact that for high-Reynolds-number flows the viscous terms of the Navier-Stokes equations are important only in a thin shear layer and therefore solution of the full equations may not be needed everywhere. Most of the flow can be considered inviscid and, neglecting the entropy and vorticity effects, a potential model is a good approximation in the flow core. The rotational part of the flow can then be calculated by solution of the potential, streamfunction and vorticity transport equations. Implementation of the no-slip and no-penetration boundary conditions at the walls provides a simple mechanism for the interaction between the viscous and inviscid solutions and no extra coupling procedures are needed. Results are presented for turbulent transonic internal choked flows.     

Influence of the mixing of two flows with different total parameters in a laval nozzle on its integral characteristics

A calculation is made of the turbulent zone of mixing of two flows of viscous and heat conducting gas in a Laval nozzle. For such a nozzle of given geometry, a comparison is made of calculations of the integrated characteristics of flows that are nonuniform with respect to the total parameters in the framework of various models: laminar hydraulics, viscous laminar hydraulics, and total mixing without hydraulic losses. The calculations are made for a stationary, nonswlrling flow of a viscous heat conducting gas with nearly discontinuous step distribution of the total parameters at the entrance to an axisymmetric Laval nozzle of given geometry. In this situation, the gas flows with different total parameters at the entrance to the nozzle are separated by a surface near which the profiles of the flow parameters are specified on the basis of boundary-layer theory. In the blocked regime investigated here, the flow in the part where the nozzle becomes narrower and at least at the beginning of the expanding part does not depend on the pressure of the surrounding medium. The integrated characteristics of the nozzle (gas flow rate G, impulse I, specific impulse i = I/G, etc.) depend on the parameter distributions at the entrance to the nozzle, and also on the turbulent mixing of the flows in the mixing zone. To analyze the dependence of the integrated characteristics on the turbulent mixing, the values of these characteristics calculated in the framework of the three models are compared. The model of mixing without hydraulic losses presupposes complete equalization of the parameters of the original inhomogeneous flow in the constant-area chamber in front of the nozzle with conservation of the mass, energy, and momentum fluxes. The model of laminar hydraulics is described in detail in [1, 2]. The model of viscous laminar hydraulics will be described in Sec. 1.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 114–119, July–August, 1979.I thank A. N. Kraiko for supervising the work, A. N. Sekundov for helpful discussions, and I. P. Smirnova and A. B. Lebedev for making available the computer program.     

Experimental study of boundary layer separation in truncated ideal contour nozzles

DLR Lampoldshausen carried out a cold flow test series to study the boundary layer separation and the related flow field in a truncated ideal contour nozzle. A special focus was set on low nozzle pressure ratios to identify the origin of a locally re-attached flow condition that was detected in previous test campaigns. A convex shaped Mach disc was found for nozzle pressure ratios less than 10 and a slight concave one for nozzle pressure ratios more than 20. Due to boundary layer transition at low nozzle pressure ratios the convex Mach disc is temporary tilted and redirects the flow towards the nozzle wall. A simple separation criterion for turbulent nozzle flows is presented that fits well for both cold and hot flows. It is shown that the oblique separation shock recompresses the flow to 90% of the ambience. The separation zone of the presented film cooled nozzle is compared with a conventional one around 40% longer. Furthermore a relation between shear layer shape and forced side loads is described.      

Cavity detachment on a hydrofoil with the inclusion of surface tension effects

The non-linear problem of cavity flow past a hydrofoil is considered with taking into account fluid viscosity in the cavity closure region and surface tension, which affect the cavity detachment. The theoretical model is based on the concept of viscous–inviscid interaction between the outer inviscid cavity flow and the inner turbulent separated flow downstream of the cavity. The outer inviscid flow is solved by constructing the complex flow potential, and the wake model is based on the method of integral relationships for separated turbulent flows. The obtained numerical results are compared with experimental data.     

Dynamic stall modelling on airfoils based on strong viscous–inviscid interaction coupling

The prediction of the aerodynamic performance of pitching airfoils in stall conditions is considered in the context of strong viscous–inviscid interaction modelling. The aim of the work is to demonstrate the capabilities of a low‐cost dynamic stall model well suited for engineering applications. The model is formulated on the basis of a standard panel method combined with a vortex blob approximation of the wake. The development of the boundary layer over the airfoil and the evolution of the shear layer in the wake are taken into account by means of strong viscous–inviscid interaction coupling. To this end a transpiration layer is added to the inviscid formulation which represents the displacement effect viscosity results in the flow while the non‐linear coupled equations are solved simultaneously. Separation is modelled by introducing a second wake originating from the separation point (‘double‐wake’ concept) which is provided as part of the boundary layer solution. The theoretical presentation of the model is supported with favourable comparisons to four sets of wind tunnel measurements. Copyright © 2007 John Wiley & Sons, Ltd.     

Study of laminar separation in the vicinity of the point of piecewise-constant pressure distribution over the wall

A family of two-dimensional divergent channels with piecewise-constant velocity and pressure distributions over the wall is considered. The method of matched asymptotic expansions is applied to study the two-dimensional viscous incompressible flow at high but subcritical Reynolds numbers in the vicinity of a pressure jump point on the channel wall. It is shown that if the pressure difference is of the order O(Re^?1/4), then in the vicinity of this point a classical region of interaction between the viscous boundary layer on the wall and the outer inviscid flow occurs. The problem formulated for the interaction region is solved numerically. The asymptotic values of the pressure difference corresponding to separationless flow are determined and the separation flow patterns are constructed.     

粘流-无粘干扰流动(IF)理论

对不可压缩层流二维干扰流动,本文提出一个干扰流动(IF)理论。IF理论要点为:1)干扰流动沿主流的法向被分为三层即粘性层、干扰层和无粘层,引进了法向动量交换为主导过程的干扰层概念。2)利用力学守恒律、三层匹配关系及文中引进的干扰模型,把三层的空间尺度及惯性-粘性诸力的数置级表示为单参数m的函数,m<1/2·3)导出描述各层流动的控制方程、导出描述全城流动的控制方程为简化Navie-Stokes(SNS)方程。IF理论适用于不存在分离的附着干扰流动以及存在分离的大范围干扰流动,经典边界层(CBL)理论和流动分离局部区域Triple-Deck(TD)理论分别是本文理论在参数m=O和1/4时的两个特例,本文理论容易推广到可压缩、三维及湍流流动。     

Calculation of hysteresis and flow-rate oscillations of base pressure in supersonic annular nozzles

The interaction of the turbulent axisymmetric near wake behind the face of the central body of an annular nozzle with the supersonic annular jet discharging from this nozzle is analyzed. The flow in the monoparametric near wake is calculated by the integral method [1] while the flow in the nonviscous jet is calculated by the method of through calculation using a monotonic explicit difference system of the first order of accuracy [2]. The interaction between the nonviscous and turbulent streams is determined by the displacement thickness of the wake. The initial conditions of the wake are determined from the integral conditions of attachment with the mixing flow in the isobaric base region. The interaction flow is described by the particular solution of the equations which passes through the singular saddle point — the throat of the wake. The near wake and base pressure in different modes of discharge from an annular nozzle at the exit cross section of which the ratio of outer and inner radii is y₂/y₁ = 1.3 and the Mach number is M = 2.54 are calculated as an example. The region of hysteresis of the base pressure, connected with the ambiguity of the interaction flow owing to the formation of the throat of the wake within the first or second barrel of the jet, and the parameters of the low-frequency flow-rate oscillations of base pressure in this region are determined. The results of the calculations are in satisfactory agreement with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 125–130, January–February, 1977.     

The effect of streamwise vortices on the turbulence structure of a separating boundary layer

A high Reynolds number flat plate turbulent boundary layer is investigated in a wind-tunnel experiment. The flow is subjected to an adverse pressure gradient which is strong enough to generate a weak separation bubble. This experimental study attempts to shed some new light on separation control by means of streamwise vortices with emphasize on the change in the boundary layer turbulence structure. In the present case, counter-rotating and initially non-equidistant streamwise vortices become and remain equidistant and confined within the boundary layer, contradictory to the prediction by inviscid theory. The viscous diffusion cause the vortices to grow, the swirling velocity component to decrease and the boundary layer to develop towards a two-dimensional state. At the position of the eliminated separation bubble the following changes in the turbulence structure were observed. The anisotropy state in the near-wall region is unchanged, which indicates that it is determined by the presence of the wall rather than the large scale vortices. However, the turbulence in the outer part of the boundary layer becomes overall more isotropic due to an increased wall-normal mixing and a significantly decreased production of streamwise fluctuations. The turbulent kinetic energy is decreased as a consequence of the latter. Despite the complete change in mean flow, the spatial turbulence structure and the anisotropy state, the process of transfer of turbulent kinetic energy to the spanwise fluctuating component seems to be unchanged. Local regions of anisotropy are strongly connected to maxima in the turbulent production. For example, at spanwise positions in between those of symmetry, the spanwise gradient of the streamwise velocity cause significant production of turbulent fluctuations. Transport of turbulence in the spanwise direction occurs in the same direction as the rotation of the vortices.     

The structure of the supersonic turbulent boundary layer in its interaction with a shock wave

An experimental study is made of the turbulent boundary layer in its interaction with a shock wave, the purpose being to clarify questions connected with the increase in the fullness of the velocity profiles. New systematic data are obtained on the development of the boundary layer, and its structure and asymptotic behavior beyond the interaction region. These results are for axisymmetric flow in the range of Mach numbers M=2–4 and angles of rotation of the flow 10–25°. Conditions of developed separation are included. Extensive information about the general properties of flows with separation has been obtained in a number of studies. A survey of these may be found, for example, in [1, 2]. Certain questions about the separation and reattachment of the boundary layer are clarified. The dimensions of the separation region are determined and its structure studied in detail for various shapes of the surface around which the flow takes place. Nevertheless it has not yet proved possible to reach a complete understanding of this complex phenomenon. Usually plane models have been used for the investigations, but in this case it is evidently impossible to exclude completely the influence of end effects on the flow in the interaction zone. Therefore it is preferable to study such flows in axisymmetric models; this considerably eases the task of analyzing and interpreting the results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 75–82, September–October, 1985.     

Advances in the study of separated flows

This paper consists of three parts. The first deals with the separation conditions for three dimensional steady viscous separated flows, in which the behaviour of separated flow described by Navier-Stokes equations and boundary layer equations is studied. The second part involves an application of differential topology to qualitative analysis of flow fields. Here the distribution rule of singular points on the separation line is studied. The last part discusses the numerical method solving Navier-Stokes equations for separated flows. The obtained computational results are analysed by the above mentioned theories and methods.     

Unsteady Separation Processes at High Reynolds Number and Their Control

Physical situations where a viscous boundary layer breaks down and interacts strongly with an effectively inviscid external flow are common place. For large Reynolds numbers, viscous effects are normally confined to thin boundary layers on all solid surfaces for the majority of any observation time. In most practical situations, exposure of such layers to an adverse pressure gradient is inevitable and in this circumstance, a sequence of events commences near the wall that culminates in an eruption and a strong viscous-inviscid interaction with the external flow. The events leading up to eruption are known as the Van Dommelen–Shen process and the eruption itself is referred to as boundary-layer separation; here the term ‘separation’ denotes the first process of interaction between a hitherto thin boundary layer and the external flow. The event is sufficiently complicated that extraordinary measures are needed to compute its evolution. In most situations, the onset of separation is subtle and hard to detect and thus development of rational control procedures is a challenging task. Here recent calculations of unsteady separation events are discussed for two- and three-dimensional flows. The phenomena involved are generic but leading-edge separation on airfoils and rotorcraft blades is emphasized. Recent studies on various control mechanisms are described, which are found to have the effect of slowing down and/or weakening the separation process. For some control processes, it has proved possible to eliminate separation entirely.     

Quasihomogeneous model of cavitation flows in diffuser channels

Starting with the experiments carried out by Reynolds in 1894, the flow in Venturi tubes has traditionally been used to study and demonstrate various forms of cavitation. Numerous authors have carried out experimental research on the various flow regimes in diffuser channels [1–7] or have investigated theoretical models of such flows [6, 8]. The occurrence and development of cavitation is closely associated with the phenomenon of turbulent separation complicated by the presence of two-phase flow in the dissipation zone. For a long time these effects were considered separately, until Gogish and Stepanov [9] proposed a single model of cavitation and separation based on the theory of intense interaction of an incompressible potential flow and a turbulent cavitation layer of variable density and embracing the various stages of cavitation. The object of this study is to demonstrate the possibilities of this model with reference to the simple example of flows accompanied by cavitation and separation in plane and axisymmetric diffuser channels of the Venturi tube type with straight and curved walls. The dissipative flow near the walls is described by a quasihomogeneous model of turbulent two-phase flow, in which the presence of two phases is taken into account only by varying the mean density. The potential core of the flow is considered in the one-dimensional formulation. The displacement thickness serves as the flow interaction parameter. The conditions of ocurrence and development of circulatory flows are determined. Examples of symmetrical and nonsymmetrical flows are presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 47–54, September–October, 1986.     

The Fascination of Vortex Rings

We present inviscid and viscous models for the formation and propagation of single, and co-axial pairs of, vortex rings. Inviscid flows are based on both thin rings, and thick rings treated by a contour dynamics approach, whilst viscous flows are determined from numerical solutions of the Navier–Stokes equations. A kaleidoscope of different flow behaviours for these axisymmetric flows is presented.     

High-Velocity Laminar and Turbulent Flow in Porous Media

We model high-velocity flow in porous media with the multiple scale homogenization technique and basic fluid mechanics. Momentum and mechanical energy theorems are derived. In idealized porous media inviscid irrotational flow in the pores and wall boundary layers give a pressure loss with a power of 3/2 in average velocity. This model has support from flow in simple model media. In complex media the flow separates from the solid surface. Pressure loss effects of flow separation, wall and free shear layers, pressure drag, flow tube velocity and developing flow are discussed by using phenomenological arguments. We propose that the square pressure loss in the laminar Forchheimer equation is caused by development of strong localized dissipation zones around flow separation, that is, in the viscous boundary layer in triple decks. For turbulent flow, the resulting pressure loss due to average dissipation is a power 2 term in velocity.     

The turbulent boundary layer on the moving surface of a cylindrical body

A study is made of the problem of a two-dimensional turbulent boundary layer on the moving surface of a cylindrical body (a Rankine oval with a relative elongation of four) moving at constant velocity in an incompressible fluid. For the numerical simulation of the turbulent flow of the fluid, the boundary layer is divided into exterior and interior regions in accordance with a two-layer model, using different expressions for the coefficients of turbulent transfer for each region. A study was nade of the development of the boundary layer on the body at different speeds of the body surface and different Reynolds numbers. The following integral characteristics were found by numerical calculation: the work of friction as the body is displaced; the work expended on the movement of its surface; and, for a flow regime with separation, the work of the pressure force. In this case the following model of separation flow is assumed: beyond the singular point in the solution of the boundary layer equations that indicates the appearance of a region of reverse flow, the pressure and friction stress on the wall are constant and are determined by their values at the singular point.Translated from Izvestiya Akademii Nauk SSSH, Mekhanika Zhidkosti i Gaza, No. 5, pp. 61–67, September–October, 1984.Finally, the author would like to thank G. G. Chernyi and Yu. D. Shevelev for useful discussions and for their interest in this work.     

Effect of the turbulent viscosity relaxation time on the modeling of nozzle and jet flows

Certain modifications of three-equation turbulence models are proposed. They are intended for increasing the accuracy of the calculations of turbulent flows in nozzles with boundary layer separation and in supersonic jets with complicated shock wave structures. Basing on the idea of the inclusion of flow prehistory in terms of an additional relaxation equation for nonequilibrium turbulent viscosity we propose three modifications of the k-ω-µ _t model based on the k-ω model and a version of the k- ?-µ _t turbulence model. In these modifications we introduce an additional dependence of the nonequilibrium turbulent viscosity relaxation time on different physical parameters which can be important near the point of boundary layer separation from the nozzle wall, such as viscous effects and effects of large gradients of the mean velocity and the kinetic energy of turbulence (turbulent pressure). The comparison of the results of the calculations with the experimental data shows that all the proposed versions of the three-equation models make it possible to improve the accuracy of the calculations of turbulent flows in nozzles and jets.