The effect of wall cooling on compressible Görtler vortices

It is known that in adiabatic boundary layer flow over a curved surface the detailed structure of the spanwise periodic Görtler vortex instability varies markedly over the range of spanwise wavelength. At short wavelengths the modes tend to be concentrated in a well-defined thin zone located within the boundary layer. As the vortex wavenumber diminishes so the region of vortex activity is first driven to the bounding wall but subsequently expands to cover the entire boundary layer at which stage the modes take on a principally inviscid form. At yet longer wavelengths the vortices are given by the solution of an interactive multi-deck structure which has some similarities with that for Tollmien–Schlichting waves.In this work we investigate how the application of wall cooling affects the above scenario. It is shown how cooling both restricts the range of mode types and gives rise to two new structures. The first, for moderate cooling and which relates to longer wavelengths, is interactive in nature. Here the viscous–inviscid interaction between an essentially inviscid Görtler problem, albeit for an effective basic flow which in its general form has a non-standard near-wall structure, and a viscous sublayer is provided by novel boundary conditions. Shorter wavelength vortices are largely unaffected by wall cooling unless this is quite severe. However when this degree of cooling is applied, the vortices take on a fully viscous form and are confined to a thin region next to the bounding wall wherein the basic flow assumes an analytic form. Numerical solutions are obtained and we provide evidence as to how the two new structures are related both to each other and to the previously known uncooled results.     

ANALYSIS/DESIGN CALCULATIONS OF TRANSONIC FLOWS

SUMMARY

Analysis/design calculations of transonic flow are discussed and several improvements are made. The nonisentropic potential method is used to calculate the inviscid transonic flow analysis problem instead of the traditional potential method. An inverse integral 3D boundary layer method is used to calculate the boundary layer in the viscous transonic flow analysis problem. The viscous/inviscid interaction calculations are carried out by a semi-inverse coupling scheme. In design problem calculations, an improved residual-correction method is used. Three individual methods are combined in a global algorithm and computing code. The improvements speed up the convergence, increase applicability and computational efficiency. Some numerical results are given to illustrate that the present method provides an effective engineering tool of high accuracy and efficiency in three dimensional transonic analysis and design situations.     

A HYBRID SPECTRAL/FINITE VOLUME METHOD FOR VISCOUS FLOWS

Abstract

This paper presents a hybrid spectral/finite volume method for steady-state compressible viscous flows. The method is evaluated for accuracy via test cases for various Mach numbers. The domain is divided into a viscous region and an inviscid region. The viscous region uses the full Navier-Stokes equations, while the inviscid region employs the Euler equations. A high order Chebyshev collocation spectral method is developed for the viscous region to resolve boundary layers. This method avoids the dense grids needed by finite-volume methods to resolve the viscous areas. A low order finite-volume method based on a Lax-Wendroff type scheme is employed for the inviscid region. A special interface formulation is developed for coupling the spectral with the finite-volume method. Comparisons with analytic results as well as convergence histories are presented.     

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.     

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.     

Analytical study of heat transfer from circular cylinder in liquid metals

In this study the influence of a thin hydrodynamic boundary layer on the heat transfer from a single circular cylinder in liquid metals having low Prandtl number (0.004–0.03) is investigated under isothermal and isoflux boundary conditions. Two separate analytical heat transfer models, viscous and inviscid, are developed to clarify the discrepancy between previous results. For both models, integral approach of the boundary layer analysis is employed to derive closed form expressions for the calculation of the average heat transfer coefficients. For an inviscid model, the energy equation is solved using potential flow velocity only whereas for a viscous model, a fourth-order velocity profile is used in the hydrodynamic boundary layer and potential flow velocity is used outside the boundary layer. The third-order temperature profile is used inside the thermal boundary layer for both models. It is shown that the inviscid model gives higher heat transfer coefficients whereas viscous flow model gives heat transfer results in a fairly good agreement with the previous experimental/numerical results.     

翼身组合体气动力计算研究

采用求解Euler方程结合附面层修正的方法在结构网格上对翼身组合体跨音速流场进行了数值模拟.附面层方程的求解应用Whitfield提出的动量积分方程和平均流动能积分方程,为了保持Euler方程求解过程中计算网格的固定性,用加在物面上的溢出速度来模拟附面层效应.针对传统的近场方法计算阻力,计算精度较低、误差较大并且不能给出各阻力分量值的缺点,将基于动量定理的远场方法用于飞机的阻力估算,采用远场法将阻力分解为:粘性阻力,激波阻力,诱导阻力,并对各个分量分别进行了求解,将计算结果与近场法以及风洞实验值做了比较.以DLR-F4翼身组合体为考核算例,对所述方法进行了验证,结果显示远场法的计算结果与风洞实验值吻合的很好.     

Unsteady flow calculation in a radial flow centrifugal pump with spiral casing

The prediction of the two-dimensional unsteady flow established in a radial flow centrifugal pump is considered. Assuming the fluid incompressible and inviscid, the velocity field is represented by means of source and vorticity surface distributions as well as a set of point vortices. Using this representation, a grid-free (Lagrangian) numerical method is derived based on the coupling of the boundary element and vortex particle methods. In this context the source and vorticity surface distributions are determined through the non-entry boundary condition together with the unsteady Kutta condition. In order to satisfy Kelvin's theorem, vorticity is shed at the trailing edges of the impeller blades. Then the vortex particle method is used to approximate the convection of the free vorticity distribution. Results are given for a pump configuration experimentally tested by Centre Technique des Industries Mécaniques (CETIM). Comparisons between predictions and experimental data show the capability of the proposed method to reproduce the main features of the flow considered.     

Treatment of vortex sheets for the transonic full-potential equation

This paper summarizes a combined analytical-computational technique which models vortex sheets in transonic potential-flow methods. In this approach, the inviscid nature of discontinuities across vortex sheets is preserved by employing the step function to remove singularities at these surfaces. The location and strength of the vortex sheets are determined by satisfying the flow-tangency boundary condition and the vorticity transport equation. The theory is formulated for the general three-dimensional case, but its application is confined to the problem of computing slipstreams behind propellers with free-vortex blading in axisymmetric flows.     

Application of a vortex method to free surface flows

Vortex methods have found wide applications in various practical problems. The use of vortex methods in free surface flow problems, however, is still very limited. This paper demonstrates a vortex method for practical computation of non-linear free surface flows produced by moving bodies. The method is a potential flow formulation which uses the exact non-linear free surface boundary condition at the exact location of the instantaneous free surface. The position of the free surface, on which vortices are distributed, is updated using a Lagrangian scheme following the fluid particles on the free surface. The vortex densities are updated by the non-linear dynamic boundary condition, derived from the Euler equations, with an iterative Lagrangian numerical scheme. The formulation is tested numerically for a submerged circular cylinder in unsteady translation. The iteration is shown to converge for all cases. The results of the unsteady simulations agree well with classical linearized solutions. The stability of the method is also discussed.     

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.     

Viscous and inviscid interactions of an oblique shock with a flexible panel

The complex self-sustained oscillations arising from the interaction of an oblique shock with a flexible panel in both the inviscid and viscous regimes have been investigated numerically. The aeroelastic interactions are simulated using either the Euler or the full compressible Navier–Stokes equations coupled to the nonlinear von Karman plate equations. Results demonstrate that for a sufficiently strong shock limit-cycle oscillations emerge from either subcritical or supercritical bifurcations even in the absence of viscous separated flow effects. The critical dynamic pressure diminishes with increasing shock strength and can be much lower than that corresponding to standard panel flutter. Significant changes in panel dynamics were also found as a function of the shock impingement point and cavity pressure. For viscous laminar flow above the panel without a shock, high-frequency periodic oscillations appear due to the coupling of boundary-layer instabilities with high-mode flexural deflections. For a separated shock laminar boundary layer interaction, non-periodic self-excited oscillations arise which can result in a significant reduction in the extent of the time-averaged separation region. This finding suggests the potential use of an aeroelastically tailored flexible panel as a means of passive flow control. Forced panel oscillations, induced by a specified variable cavity pressure underneath the panel, were also found to be effective in reducing separation. For both inviscid and viscous interactions, the significant unsteadiness generated by the fluttering panel propagates along the complex reflected expansion/recompression wave system.     

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.     

Asymptotic Analysis of Viscous Fluctuations in Turbulent Boundary Layers

A two-dimensional boundary layer of an incompressible viscous fluid is investigated in the presence of velocity and pressure fluctuations. The characteristic Reynolds number is high and, as a consequence, the unsteady (turbulent) boundary layer is thin. An asymptotic approach is used to analyze the complete unsteady Navier–Stokes equations, which makes it possible to separate out the characteristic viscous and inviscid flow zones in the boundary layer and to solve the corresponding problems. The analytical expressions for the viscous fluctuations governed by the Hamel equation with a large value of the parameter are derived.     

Asymptotic model of the development of separations inside a boundary layer under the action of a traveling pressure wave

A model of the nonlinear interaction between a pressure perturbation traveling at a constant velocity and an incompressible boundary layer is constructed when its near-wall part is described by the “inviscid boundary layer” equations. A steady-state solution is managed to obtain in the finite form under the assumption that it exists in a moving coordinate system. It is shown that the boundary layer can easily overcome pressure perturbations whose amplitude is not higher than the dynamic pressure calculated from the velocity of the pressure perturbation. At the higher pressure perturbation amplitudes a vortex sheet sheds from the body surface to the boundary layer. The vortex sheet represents an unstable surface of the tangential discontinuity which separates the regions of the direct and reverse separation flows. In the case of an arbitrary shape of the pressure perturbation the surface of the tangential discontinuity sheds from the body surface at a finite angle with the formation of a stagnation point. An example of the pressure perturbation in which the vortex sheet sheds from the body surface along the tangent is constructed.     

Linear stability of swirling flows computed as solutions to the quasi-cylindrical equations of motion

The linear stability of numerical solutions to the quasi-cylindrical equations of motion for swirling flows is investigated. Initial conditions are derived from Batchelor's similarity solution for a trailing line vortex. The stability calculations are performed using a second-order-accurate finite-difference scheme on a staggered grid, with the accuracy of the computed eigenvalues enhanced through Richardson extrapolation. The streamwise development of both viscous and inviscid instability modes is presented. The possible relationship to vortex breakdown is discussed.     

Unsteady Viscous Flows about Bodies: Vorticity Release and Forces

The force (drag and lift) exerted on a body moving in a viscous fluid is expressed via the free and bound vorticity moments, and the role of vortex shedding is discussed. The formulation encompasses classical, inviscid flows, and leads to efficient computational methods. Numerical results for a few prototype flows are presented.     

Hypersonic flow over a six-finned configuration

Experimental data of the flow field over a six-finned configuration are compared with computations from an inviscid Euler code and a parabolized Navier-Stokes (PNS) code to assess their accuracy. Surface oil flow patterns, surface pressures, and flow field pilot surveys are presented. The PNS computations are in good agreement with the measurements in the inviscid regions. Differences are noted in the viscous regions possibly due to nose tip bluntness effects and the associated entropy layer.     

Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

A combined Lagrangian method for simulation of axisymmetric gas-particle vortex flows

A combined fully Lagrangian approach for meshless modeling of unsteady axisymmetric vortex flows of a gas-particle mixture with an incompressible carrier phase is developed. The approach proposed is based on the combination of a meshless vortex method for calculating axisymmetric flows of the carrier phase described by the Navier–Stokes (or Euler) equations and the full Lagrangian method for calculating the parameters of the dispersed phase. The combination of these methods reduces the problem of modeling the two phase flows to the solution of a high-order system of ordinary differential equations for the coordinates of toroidal vortex elements in the carrier phase and the particle trajectories, the velocity components, and the components of the Jacobian of transformation from the Eulerian to the Lagrangian variables in the dispersed phase. The application of the method is illustrated by modeling the behavior of an admixture of inertial Stokes particles with a small mass concentration in unsteady flows like solitary vortex rings in a viscous carrier phase and groups of vortex rings in an effectively inviscid carrier phase.