Receptivity of the inclined-cylinder attachment-line boundary layer to vortex perturbations

The receptivity of the boundary layer in the neighborhood of the attachment line of a cylinder inclined to the flow with respect to periodic vortex perturbations frozen into the stream is investigated. The problem considered simulates the interaction between external turbulence and the leading-edge swept wing boundary layer. It is shown that if the direction of the external perturbation vector is almost parallel to the leading edge, then the external perturbations are considerably strengthened at the outer boundary layer edge. This effect can cause laminar-turbulent transition on the attachment line at subcritical Reynolds numbers.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 72–85. Original Russian Text Copyright © 2004 by Ustinov.     

Supersonic boundary layer transition induced by roughness on the attachment line of a yawed cylinder

The effect of a single two-dimensional irregularity and sandy roughness on boundary layer transition in supersonic flow over a yawed cylinder (M = 6)-has been experimentally investigated. The characteristic flow regimes beyond the roughness are identified, and their limits are determined as a function of the Reynolds number and the ratio of the height of the roughness to the characteristic thickness of the boundary layer. A qualitative comparison is made with the flow regimes induced by roughness on the attachment line in incompressible flow over a cylinder [1–3]. The thermal indicator coating method is used to measure the heat fluxes along the attachment line and a comparison is made with calculations carried out in accordance with the methods of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 28–35, November–December, 1991.The authors are grateful to A. F. Kiselev for helping to calculate the heat flux in the turbulent boundary layer.     

Receptivity of a boundary layer to external sonic waves

The effect of external acoustic perturbations on a two-dimensional laminar boundary layer is studied within the framework of the asymptotic theory. Essentially nonparallel regimes of the basic flow when flows with a separation zone develop are investigated.     

采用不同黏性处理方法的宽无叶扩压器不稳定流动研究

径向无叶扩压器的全局稳定性可能受到核心主流失稳,出口回流与壁面边界层分离等因素影响,对于宽无叶扩压器,无黏核心主流与壁面边界层流动对不稳定扰动诱发的作用机理是当前研究的重点.本文首先通过数值计算获得了大宽度比孤立无叶扩压器平均流动,然后基于小扰动理论和周向均质假设,分别对欧拉方程与 Navier-Stokes 方程进行线性化,建立了基于无黏核心流动的稳定性分析方法,以及基于涡黏性与分子黏性的混合稳定性分析方法;通过与实验结果的对比,验证了混合稳定性分析方法预测所得流动失稳频率和全局直接模态的准确性;最后基于伴随方法获得了特征值的结构敏感性,揭示了不同黏性处理条件下宽无叶扩压器内全局不失稳扰动的源发区域.在只考虑核心主流的无黏条件下,宽无叶扩压器内流动不稳定扰动来源于流场中部,为二维的离心失稳;在同时考虑核心主流与边界层的作用时,宽无叶扩压器不稳定扰动不仅来源于扩压器流场中部的核心主流,壁面回流对于不稳定扰动的产生了重要影响.      

一种鲁棒的HLLC格式及其稳定性分析

精确捕捉接触波和剪切波的Godunov型数值方法,如流行的HLLC格式,在模拟高超声速流动问题时会出现激波异常现象。对HLLC格式进行稳定性分析发现,流体主流方向的扰动都能有效衰减,但是横向的密度与剪切速度的扰动不会衰减。具有特殊对称性的二维Sedov爆轰波问题证明了横向通量和不稳定现象之间的密切联系。利用压力比和马赫数来探测数值激波层亚声速区的横向网格界面,并且在该界面的数值通量上增加熵波粘性和剪切波粘性来构造一种激波稳定的HLLC格式。分析表明,在熵波粘性和剪切波粘性的作用下,横向的所有扰动都会衰减。一系列数值测试证明了新格式不仅可以成功地抑制各类激波异常现象,还保留了原HLLC格式低耗散性的优点。     

Effect of Heat Supply on the Stability of the Supersonic Swept Atiachment-Line Boundary Layer

The stability of a boundary layer with volume heat supply on the attachment line of a swept wing is investigated within the framework of the linear theory at supersonic inviscid-free-stream Mach numbers. The results of numerical calculations of the flow stability and neutral curves are presented for the flow on the leading edge of a swept wing with a swept angle χ=60° at various free-stream Mach numbers. The effect of volume heat supply on the characteristics of boundary layer stability on the attachment line is studied at a surface temperature equal to the temperature of the external inviscid flow. It is shown that in the case of a supersonic external inviscid flow volume heat supply may result in an increase in the critical Reynolds number and stabilization of disturbances corresponding to large wave numbers. For certain energy supply parameters the situation is reversed, the unstable disturbances corresponding to the main flow-instability zone are stabilized but another zone of flow-instability with small wave numbers and a significantly lower critical Reynolds number appears.     

A finite volume Navier-Stokes algorithm for adaptive grids

A compact, finite volume, time-marching scheme for the two-dimensional Navier-Stokes equations of viscous fluid flow is presented. The scheme is designed for unstructured (locally refined) quadrilateral meshes. An earlier inviscid equation (Euler) scheme is employed for the convective terms and the emphasis is on treatment of the viscous terms. An essential feature of the algorithm is that all necessary operations are restricted to within each cell, which is very important when dealing with unstructured grids. Numerical issues which have to be addressed when developing a Navier-Stokes scheme are investigated. These issues are not limited to the particular Navier-Stokes scheme developed in the present work but are general problems. Specifically, the extent of the numerical molecule, which is related to the compactness of the scheme and to its suitability for unstructured grids, is examined. An approach which considers suppression of odd-even mode decoupling of the solution when designing a scheme is presented. In addition, accuracy issues related to grid stretching as well as boundary layer solution contamination due to artificial dissipation are addressed. Although the above issues are investigated with respect to the specific scheme presented, the conclusions are valid for an entire class of finite volume algorithms. The Navier-Stokes solver is validated through test cases which involve comparisons with analytical, numerical and experimental results. The solver is coupled to an adaptive algorithm for high-Reynolds-number aerofoil flow computations.     

Interferometry and reconstruction of strongly refracting fields in two-dimensional boundary layer flow

 Optical methods like interferometry as non-intrusive experimental techniques are used for fine analysis of flowfields. The accuracy of the measuring method is very important for the applicability of its results to CFD validation. A common evaluation method to reconstruct interferograms is based on the assumption, that the object ray propagates along a straight line. But the strong bending of rays that occurs, e.g. in supersonic boundary layer flow, cannot be neglected without losses in reconstruction quality. Since the reconstruction of a two-dimensional boundary layer flow can be considered as an one-dimensional problem, the phase difference of the object and reference ray at the interferogram can be related analytically to the refractive-index distribution using a Taylor series expansion. The resulting interferometric equation is an ordinary non-linear second-order differential equation, which can be integrated by numerical methods. By application of this interferometric equation on the one hand, the error in the “classical” interferometry resulting from the ray bending neglection can be estimated. On the other hand, the accuracy in evaluation of interferograms of two-dimensional boundary layer flow can be improved by solving this equation. Received: 23 May 1996 / Accepted: 21 September 1996     

Stability of a nonisothermal boundary layer in the presence of periodic perturbations of the exterior flow velocity

The stability of a laminar boundary layer in the presence of high-frequency time-periodic perturbations of the exterior flow velocity, in particular, acoustic vibrations, is investigated in a series of papers which are reviewed in detail in monograph [1]. The mechanisms by which such perturbations influence the stability and transition to the turbulent flow regime can vary. For example, they can lead to the deformation of the averaged field of the basic flow. However, there was good reason not to discuss the effect of this type of perturbation earlier, as it was considered that the change in the basic flow was very small even for perturbations of great amplitude. The aim of the present paper is to demonstrate how perturbations or pulsations in the exterior flow velocity can, by changing the basic flow, have a strong influence on the stability of the laminar boundary layer of a gas under appreciably nonisothermal conditions. Examples of calculations that support this assertion are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 183–186, September–October, 1984.The authors wish to thank A. A. Maslov for his help with the calculations.     

On asymptotic theory of separated flows past bodies

Symmetric two-dimensional steady flow past a body in a homogeneous incompressible fluid stream at high Reynolds numbers is considered. A slow motion in the reverse flow zone is investigated and the solution for the flow in the external region is obtained in the second approximation. Additional considerations of the fact that the flow in the closure region of the separation zone and in the wake behind this zone is turbulent are presented. The laminar-turbulent transition in the mixing layer is analyzed and an analogy between this process and the propagation of perturbations upstream of the boundary layer interaction regions is revealed.     

A numerical simulation of boundary layer effects in a shock tube

A numerical scheme is used to investigate boundary layer effects in a shock tube. The method consists of a mixture of Roe's approximate Riemann solver and central differences for the convective fluxes and central differences for the viscous fluxes and is implicit in one space dimension. Comparisons are made with experimental data and with solutions obtained via boundary layer equations. Examination of the calculated flow field explains the observed behaviour and highlights the approximate nature of boundary layer solutions.     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

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.     

Selection of convective rolls in a thin evaporating-fluid layer in an air flow

The process of selection of longitudinal convective rolls in a thin layer of evaporating fluid immersed in an air turbulent boundary layer flow is studied numerically. The dependence of the two-dimensional flow patterns on the Rayleigh number and boundary conditions is analyzed. Calculations with account for the thermocapillary effect are carried out. The numerical results are compared with experimental data.     

Free interaction of a three-dimensional boundary layer with an exterior irrotational flow

Asymptotic equations describing the unsteady free interaction of a three-dimensional boundary layer with an exterior flow are derived. The orders of the independent variables and perturbations of the flow parameters are chosen in such a way that the pressure gradient that occurs in the equation of the wall layer is due to the displacement of streamlines situated near the surface of the body. The Fourier method is used to construct a solution to the linearized problem. A class of perturbations satisfying homogeneous boundary conditions on the surface of the body is found.     

Finite analytic boundary condition for a higher‐order finite analytic scheme

A higher‐order finite analytic scheme based on one‐dimensional finite analytic solutions is used to discretize three‐dimensional equations governing turbulent incompressible free surface flow. In order to preserve the accuracy of the numerical scheme, a new, finite analytic boundary condition is proposed for an accurate numerical solution of the partial differential equation. This condition has higher‐order accuracy. Thus, the same order of accuracy is used for the boundary. Boundary conditions were formulated and derived for fluid inflow, outflow, impermeable surfaces and symmetry planes. The derived boundary conditions are treated implicitly and updated with the solution of the problem. The basic idea for the derivation of boundary conditions was to use the discretized form of the governing equations for the fluid flow simplified on the boundaries and flow information. To illustrate the influence of the higher‐order effects at the boundaries, another, lower‐order finite analytic boundary condition, is suggested. The simulations are performed to demonstrate the validity of the present scheme and boundary conditions for a Wigley hull advancing in calm water. Copyright © 2005 John Wiley & Sons, Ltd.     

A spectral collocation method for compressible,non-similar boundary layers

An efficient and highly accurate algorithm based on a spectral collocation method is developed for numerical solution of the compressible, two-dimensional and axisymmetric boundary layer equations. The numerical method incorporates a fifth-order, fully implicit marching scheme in the streamwise (timelike) dimension and a spectral collocation method based on Chebyshev polynomial expansions in the wall-normal (spacelike) dimension. The discrete governing equations are cast in residual form and the residuals are minimized at each marching step by a preconditioned Richardson iteration scheme which fully couples energy, momentum and continuity equations. Preconditioning on the basis of the finite difference analogues of the governing equations results in a computationally efficient iteration with acceptable convergence properties. A practical application of the algorithm arises in the area of compressible linear stability theory, in the investigation of the effects of transverse curvature on the stability of flows over axisymmetric bodies. The spectral collocation algorithm is used to derive the non-similar mean velocity and temperature profiles in the boundary layer of a ‘fuselage’ (cylinder) in a high-speed (Mach 5) flow parallel to its axis. The stability of the flow is shown to be sensitive to the gradual streamwise evolution of the mean flow and it is concluded that the effects of transverse curvature on stability should not be ignored routinely.     

Stability of a compressible boundary layer

The problem of stability in a compressible boundary layer, as opposed to an incompressible layer, involves many parameters and requires consideration of three-dimensional perturbations. The transverse component of the velocity, the thermal regime at the wall, etc., take on great significance. Investigation of all aspects of this problem requires systematic calculations performed by electronic computers. There do exist a few calculations of stability of a compressible boundary layer with respect to three-dimensional disturbances for particular cases. It follows from those studies (see, for example, [1]) that consideration of three-dimensional perturbations and of the transverse component of the basic flow velocity is important. Many aspects of this problem remain uninvestigated. Aside from the sheer cumbersomeness of the problem, there exist purely mathematical difficulties connected with the presence of a small parameter with higher derivatives in the differential equations for the perturbations, which causes losses in accuracy of calculation. In this present study an algorithm will be developed for solution of the problem of stability of a compressible boundary layer relative to three-dimensional disturbances with consideration of the transverse component of the basic velocity. Calculations are performed for a boundary layer on a plane thermally insulating plate, and the effects of the transverse velocity component and the three-dimensionality of the perturbations on stability at various Mach numbers are demonstrated.     

碟式离心机层流边界层方程的对称性和不变解

为运用边界层理论进行碟式离心机的碟片型线设计,首先对于碟式离心机内部多个碟片浅层固液两相界面流动模型,从流体力学基本的N-S方程组出发,得出碟片间隙斜窄流二维层流边界层方程;其次,采用李群演绎算法,推导出系统偏微分方程组允许的无穷小对称和不变解,更重要的是给出了其详细的解析表达式和有关性质。结果发现,碟式离心机层流边界层厚度与碟片结构(如碟片尺寸和间距、转速)以及流体物性(如液相粘性系数)密切相关。本文研究为考虑边界和初始条件的流场边界层精确解析解提供了有效途径,也为边界层分离控制模型的各种数值算法提供了验证。     

THE NEW SOLUTION FOR THE AXIALLY SYMMETRICAL LAMINAR BOUNDARY LAYER EQUATIONS BETWEEN TWO PARALLEL SPHERICAL SURFACES

The transformations, which are similar to Mangler's transformations, are given in this paper. They change the entrance region flow of axially symmetrical laminar boundary layer between two parallel spherical surfaces into the flow of two-dimensional boundary layer, and simplify the problems. The simplified equattons can be solved by the two-dimensional boundary layer theory and numerical methods. Therefore, a new way is opened up to solve the diffusive laminar fiow in the entrance region between two parallel spherical surfaces.