Entropy method for calculating separated-flow parameters

A variational method of calculating two-dimensional supersonic turbulent separated flows, based on the Prigogine principle (minimum entropy production), the relations of irreversible-process thermodynamics, and the integral energy equation of a viscous gas, is suggested. An application of this method is illustrated by an example of the calculations of the parameters of separated flows with a fixed separation point. The outer flow in the inviscid region is calculated by known methods. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 122–132, January–February, 1997.     

离心叶轮正冲角流场计算的边界层方法

本文使用stanitz快速分析法计算主流,用积分法计算带旋转和曲率的二维不可压湍流边界层.通过流量方程将两者紧密联系起来、同时求解,能够连续计算至分离区.可方便快速地用于估算离心叶轮在正冲角工况下的流场.     

Fast numerical evaluation of flow fields with vortex cells

A vortex cell (in this paper) is an aerodynamically shaped cavity in the surface of a body, for example a wing, designed specially to trap the separated vortex within it, thus preventing large-scale unsteady vortex shedding from the wing. Vortex stabilisation can be achieved either by the special geometry, as has already been done experimentally, or by a system of active control. In realistic conditions the boundary and mixing layers in the vortex cell are always turbulent. In the present study a model for calculating the flow in a vortex cell was obtained by replacing the laminar viscosity with the turbulent viscosity in the known high-Reynolds-number asymptotic theory of steady laminar flows in vortex cells. The model was implemented numerically and was shown to be faster than solving the Reynolds-averaged Navier–Stokes equations. An experimental facility with a vortex cell was built and experiments performed. Comparisons of the experimental results with the predictions of the model are reasonably satisfactory. The results also indicate that at least for flows in near-circular vortex cells it is sufficient to have accurate turbulence models only in thin viscous layers, while outside the viscosity should only be small enough to make the flow effectively inviscid.     

Time resolved flow-field measurements of a turbulent mixing layer over a rectangular cavity

High Reynolds number, low Mach number, turbulent shear flow past a rectangular, shallow cavity has been experimentally investigated with the use of dual-camera cinematographic particle image velocimetry (CPIV). The CPIV had a 3 kHz sampling rate, which was sufficient to monitor the time evolution of large-scale vortices as they formed, evolved downstream and impinged on the downstream cavity wall. The time-averaged flow properties (velocity and vorticity fields, streamwise velocity profiles and momentum and vorticity thickness) were in agreement with previous cavity flow studies under similar operating conditions. The time-resolved results show that the separated shear layer quickly rolled-up and formed eddies immediately downstream of the separation point. The vortices convect downstream at approximately half the free-stream speed. Vorticity strength intermittency as the structures approach the downstream edge suggests an increase in the three-dimensionality of the flow. Time-resolved correlations reveal that the in-plane coherence of the vortices decays within 2–3 structure diameters, and quasi-periodic flow features are present with a vortex passage frequency of ~1 kHz. The power spectra of the vertical velocity fluctuations within the shear layer revealed a peak at a non-dimensional frequency corresponding to that predicted using linear, inviscid instability theory.     

Dynamics of large-scale eddies in turbulent flows

A model equation based on the equipartition of the turbulent dissipation is proposed for describing the dynamics of large-scale eddies in turbulent flows. The equation is reducible to the equation of motion of an inviscid fluid, so that the motion of the large-scale eddies can be described in terms of inviscid fluid dynamics. It is found that the large-scale eddies are always weakened by the background turbulence and their evolution is slowed down compared with the corresponding inviscid motion. In the case of turbulent mixing layer, its linear growth in downstream direction is accounted for by the exponential growth in time of the perturbation in an inviscid plane vortex sheet.     

Turbulent flow calculations using unstructured and adaptive meshes

A method of efficiently computing turbulent compressible flow over complex two-dimensional configurations is presented. The method makes use of fully unstructured meshes throughout the entire flow field, thus enabling the treatment of arbitrarily complex geometries and the use of adaptive meshing techniques throughout both viscous and inviscid regions of the flow field. Mesh generation is based on a locally mapped Delaunay technique in order to generate unstructured meshes with highly stretched elements in the viscous regions. The flow equations are discretized using a finite element Navier-Stokes solver, and rapid convergence to steady state is achieved using an unstructured multigrid algorithm. Turbulence modelling is performed using an inexpensive algebraic model, implemented for use on unstructured and adaptive meshes. Compressible turbulent flow solutions about multiple-element aerofoil geometries are computed and compared with experimental data.     

Thin weakly curved wings with maximum aerodynamic quality

Lifting wings that only slightly disturb the supersonic gas flow are considered. The plan shape and thickness distribution of the wing and the free-stream parameters are given. The flow problem is solved within the framework of the Prandtl model. The outer potential flow is determined in accordance with the linear theory. The turbulent boundary layer is found by the method of plane sections with allowance for the three-dimensional inviscid flow pattern. A numerical model of the flow is constructed in the class of piecewise-constant functions on characteristic calculation grids [1]. The variational problem of finding the weakly curved middle surface of the wing giving maximum aerodynamic quality is reduced, by analogy with [2], to a problem of nonlinear programming and is solved by the gradient projection method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 165–168, July–August, 1991.     

Isolating curvature effects in computing wall-bounded turbulent flows

An adjoint optimization method is utilized to design an inviscid outer wall shape required for a turbulent flow field solution of the So–Mellor convex curved wall experiment using the Navier–Stokes equations. The associated cost function is the desired pressure distribution on the inner wall. Using this optimized wall shape with a Navier–Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient are evaluated. The one-equation Spalart–Allmaras (SA) turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model (EASM) possesses no heuristic functions or additional constants. It slightly lowers the computed skin friction coefficient and the turbulent stress levels for this case, in better agreement with experiment. The effect on computed velocity profiles is minimal.     

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.     

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

Shallow cavity flow tone experiments: onset of locked-on states

Fully turbulent inflow past a shallow cavity is investigated for the configuration of an axisymmetric cavity mounted in a pipe. Emphasis is on conditions giving rise to coherent oscillations, which can lead to locked-on states of flow tones in the pipe–cavity system. Unsteady surface pressure measurements are interpreted using three-dimensional representations of amplitude–frequency, and velocity; these representations are constructed for a range of cavity depth. Assessment of these data involves a variety of approaches. Evaluation of pressure gradients on plan views of the three-dimensional representations allows extraction of the frequencies of the instability (Strouhal) modes of the cavity oscillation. These frequency components are correlated with traditional models originally formulated for cavities in a free-stream. In addition, they are normalized using two length scales: inflow boundary-layer thickness and pipe diameter. These scales are consistent with those employed for the hydrodynamic instability of the separated shear layer, and are linked to the large-scale mode of the shear layer oscillation, which occurs at relatively long cavity length. In fact, a simple scaling based on pipe diameter can correlate the frequencies of the dominant peaks over a range of cavity depth.The foregoing considerations provide evidence that pronounced flow tones can be generated from a fully turbulent inflow at very low Mach number, including the limiting case of fully developed turbulent flow in a pipe. These tones can arise even for the extreme case of a cavity having a length over an order of magnitude longer than its depth. Suppression of tones is generally achieved if the cavity is sufficiently shallow.     

Chaos and direct numerical simulation in turbulence

We show that direct numerical simulation will yield turbulent flowfields which are strongly dependent upon computer hardware and software. A computed flow trajectory is apparently uncorrelated to the true solution of a flowfield if it is allowed to evolve over a long time, and hence is called a pseudo-orbit. This is due to the trajectory instability of chaotic turbulent flows. All is not lost, however; a long-time average of flow quantities can now be computed using a pseudo-orbit by invoking the shadowing lemma. For the inviscid flow, this time average tends to approach asymptotically the phase average as predicted by the classical ergodic theorem. Although the inviscid two-dimensional flow has no real physical importance, the existence of canonical (equilibrium) distribution permits us to examine the accuracy of time averaging based on the pseudo-orbit and its inherent limitations.This work was supported by AFOSR task 2304N1.     

Experiments on turbulent natural convection in an enclosed tall cavity

Experiments have been undertaken to investigate the natural convection of air in a tall differentially heated rectangular cavity (2.18 m high by 0.076 m wide by 0.52 m in depth). They were performed with temperature differentials between the vertical plates of 19.6°C and 39.9°C, giving Rayleigh numbers based on the width of 0.86×10⁶ and 1.43×10⁶. Under these conditions the flow in the core of the cavity is fully turbulent and property variations with temperature are comparatively small. A previously used experimental rig has been modified, by fitting partially conducting top and bottom walls and outer guard channels, to provide boundary conditions which avoid the inadequately defined sharp changes in temperature gradient and other problems associated with insufficient insulation on nominally adiabatic walls. Mean and turbulent temperature and velocity variations within the cavity have been measured, together with heat fluxes and turbulent shear stresses. The temperature and flow fields were found to be closely two-dimensional, except close to the front and back walls, and anti-symmetric across the diagonal of the cavity. The partially conducting roof and floor provide locally unstable thermal stratification in the wall jet flows there, which enhances the turbulence as the flow moves towards the temperature controlled plates. The results provide a greatly improved benchmark for the testing of turbulence models in this low turbulence Reynolds number flow.     

On the Computation of Compressible Turbulent Flows on Unstructured Grids

An accurate, fast, matrix-free implicit method has been developed to solve compressible turbulent How problems using the Spalart and Allmaras one equation turbulence model on unstructured meshes. The mean-flow and turbulence-model equations are decoupled in the time integration in order to facilitate the incorporation of different turbulence models and reduce memory requirements. Both mean flow and turbulent equations are integrated in time using a linearized implicit scheme. A recently developed, fast, matrix-free implicit method, GMRES+LU-SGS, is then applied to solve the resultant system of linear equations. The spatial discretization is carried out using a hybrid finite volume and finite element method, where the finite volume approximation based on a containment dual control volume rather than the more popular median-dual control volume is used to discretize the inviscid fluxes, and the finite element approximation is used to evaluate the viscous flux terms. The developed method is used to compute a variety of turbulent flow problems in both 2D and 3D. The results obtained are in good agreement with theoretical and experimental data and indicate that the present method provides an accurate, fast, and robust algorithm for computing compressible turbulent flows on unstructured meshes.     

Flow and heat transfer in a rotating cavity with a radial inflow of fluid Part 1: The flow structure

Flow visualization has been conducted in a rotating cavity, comprising two steel discs and a peripheral polycarbonate shroud, for dimensionless flow rates of air up to |C_w|8000 and rotational Reynolds number up to Re_φ10⁶. For all the experiments, the ratio of the inner to outer radii of the discs was 0.1 and the ratio of the axial clearance between the discs to their outer radii was 0.133; five different shroud geometries were tested. The flow visualization has confirmed that the flow structure comprises a source region near the shroud, laminar or turbulent Ekman layers on the discs, a sink layer near the centre of the cavity, and an interior core of rotating fluid. Above a certain flow rate, this structure was found to be unstable; heating one disc tended to stabilize the flow. For isothermal flow, measurements of the size of the source region were in good agreement with values predicted from a simple theoretical model.     

Velocity field in a cylindrical hydrocyclone

The flow in a cylindrical hydrocyclone is investigated at moderate Reynolds numbers. Such a flow regime is realized in a hydrocyclone when, for example, certain solid particles are separated. The flow is calculated in the frame-work of the theory of an inviscid incompressible fluid, since the influence of turbulent pulsations on the flow structure in the investigated flow regime (Reynolds numbers of the order of a few thousand) is slight. Comparison of the results with experimental data indicates good qualitative agreement.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 12–20, November–December, 1980.I thank Yu. P. Gupalo and Yu. S. Ryazantsev for discussing the work and for valuable comments.     

Compressibility effects due to turbulent fluctuations

This paper deals with intrinsic effects of compressibility, i.e. with dilatation fluctuations in response to pressure fluctuations. Three different types of turbulent flows are considered in more detail: homogeneous turbulent shear flow, wall-bounded turbulent shear flow and shock/turbulence interaction. A survey of the present knowledge in this field, mainly based on DNS data, is given. Using the linear inviscid perturbation equations a direct link between fluctuations of dilatation and of velocity in the direction of mean shear is presented for homogeneous shear flow. This relation might form the basis for a more universal pressure-dilatation model. It is conjectured that the insignificance of intrinsic compressibility effects in wall-bounded supersonic shear flow is mainly due to the impermeability constraint of the wall. To this end, a linear stability analysis of supersonic channel flow along cooled, but permeable walls has been performed based on Coleman et al.'s [5] mean flow data. It shows an increase in the moduli of eigenfunctions related to compressibility, like pressure, and in moduli of quantities derived from eigenfunctions such as ‘pressure dilatation’ and squared dilatation. Although these results do not prove our hypothesis they provide hints in this direction. Shock/turbulence interaction is viewed as a source of compressibility. Former DNS data of Hannappel and Friedrich [10] for shock/isotropic turbulence interaction showing the effect of compressibility on the amplification of fluctuations are interpreted based on linear perturbation equations.     

Problem of Deceleration of a Supersonic Conducting Channel Flow by a Magnetic Field

Problems of the deceleration of a supersonic conducting flow by a magnetic field are investigated. A conducting gas flow in a circular tube is considered in the presence of an axisymmetric magnetic field induced by a unit current loop or solenoid of finite length. The analysis is carried out on the basis of both the Euler equations (inviscid gas) and the complete system of Navier-Stokes equations for laminar viscous gas flow and turbulent flow using a one-parameter turbulence model. The numerical simulation is based on an implicit relaxation finite-difference scheme which is a modification of the Godunov method. The total pressure losses are determined for various values of the magnetohydrodynamic (MHD) interaction, the initial Mach number, and different magnetic field geometries and it is shown that the irreversible losses are significant in MHD supersonic flow deceleration.     

Extension du modèle d'échelles mixtes à la diffusivité de sous-maille

In the large-eddy simulation frame for non-isothermal turbulent flow, the Mixed Scale Model is extended to the subgrid diffusivity, in order to dissociate the computation of subgrid viscosity and diffusivity. The identification of the subgrid thermal dissipation term in the subgrid flux transport equation leads to an algebraic expression of the subgrid diffusivity. This diffusive model, as the Smagorinsky one, is weighted by a model based on scale similarity. This expression leads to satisfactory results when applied to a buoyant turbulent flow in a differentially heated cavity.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.