Global stability analysis of axisymmetric boundary layer over a circular cylinder

This paper presents a linear global stability analysis of the incompressible axisymmetric boundary layer on a circular cylinder. The base flow is parallel to the axis of the cylinder at inflow boundary. The pressure gradient is zero in the streamwise direction. The base flow velocity profile is fully non-parallel and non-similar in nature. The boundary layer grows continuously in the spatial directions. Linearized Navier–Stokes (LNS) equations are derived for the disturbance flow quantities in the cylindrical polar coordinates. The LNS equations along with homogeneous boundary conditions forms a generalized eigenvalues problem. Since the base flow is axisymmetric, the disturbances are periodic in azimuthal direction. Chebyshev spectral collocation method and Arnoldi’s iterative algorithm is used for the solution of the general eigenvalues problem. The global temporal modes are computed for the range of Reynolds numbers and different azimuthal wave numbers. The largest imaginary part of the computed eigenmodes is negative, and hence, the flow is temporally stable. The spatial structure of the eigenmodes shows that the disturbance amplitudes grow in size and magnitude while they are moving towards downstream. The global modes of axisymmetric boundary layer are more stable than that of 2D flat-plate boundary layer at low Reynolds number. However, at higher Reynolds number they approach 2D flat-plate boundary layer. Thus, the damping effect of transverse curvature is significant at low Reynolds number. The wave-like nature of the disturbance amplitudes is found in the streamwise direction for the least stable eigenmodes.     

基于扰动方程的超音速轴对称射流马赫波辐射研究

超音速不稳定波是导致剪切流失稳和转捩的主要不稳定模态,这种模态以马赫波的形式辐射到远场,从而产生强烈的声场。采用线性稳定性理论和非线性扰动方程(NLDE)分析,计算超音速轴对称射流不稳定波的扰动演化(Ma=2.1),对马赫波辐射进行研究,包括马赫波辐射方向、辐射源位置,以及随斯特劳哈尔数的变化情况。研究结果表明,在超音速轴对称射流中,马赫波沿固定方向辐射向远方,不稳定波相位沿另一方向传播,这两个方向相互正交;马赫波辐射源位置位于不稳定波压力幅值最大处;斯特劳哈尔数St越大,马赫波辐射的能力越强,辐射区域越集中。     

On the Stability of a Laminar Wall Jet with Heat Transfer

The hydrodynamic stability of a low speed, plane, non-isothermal laminar wall jet at a constant temperature boundary condition was investigated theoretically and experimentally. The mean velocity and temperature profiles used in the stability analysis were obtained by implementing the Illingworth–Stewartson transformation that allows one to extend the classical Glauert solution to a thermally non-uniform flow. The stability calculations showed that the two unstable eigenmodes coexisting at moderate Reynolds numbers are significantly affected by the heat transfer. Heating is destabilizing the flow while cooling is stabilizing it. However, the large-scale instabilities associated with the inflection point of the velocity profile still amplify in spite of the high level of the stabilizing temperature difference. The calculated stability characteristics of the wall jet with heat transfer were compared with experimental data. The comparison showed excellent agreement for small amplitudes of the imposed perturbations. The agreement is less good for the phase velocities of the sub-harmonic wave and this is attributed to experimental difficulties and to nonlinear effects.     

Asymptotic theory of neutral stability of the Couette flow of a vibrationally excited gas

An asymptotic theory of the neutral stability curve for a supersonic plane Couette flow of a vibrationally excited gas is developed. The initial mathematical model consists of equations of two-temperature viscous gas dynamics, which are used to derive a spectral problem for a linear system of eighth-order ordinary differential equations within the framework of the classical linear stability theory. Unified transformations of the system for all shear flows are performed in accordance with the classical Lin scheme. The problem is reduced to an algebraic secular equation with separation into the “inviscid” and “viscous” parts, which is solved numerically. It is shown that the thus-calculated neutral stability curves agree well with the previously obtained results of the direct numerical solution of the original spectral problem. In particular, the critical Reynolds number increases with excitation enhancement, and the neutral stability curve is shifted toward the domain of higher wave numbers. This is also confirmed by means of solving an asymptotic equation for the critical Reynolds number at the Mach number M ≤ 4.     

Flow instability of nanofuilds in jet

The flow instability of nanofluids in a jet is studied numerically under various shape factors of the velocity profile, Reynolds numbers, nanoparticle mass loadings,Knudsen numbers, and Stokes numbers. The numerical results are compared with the available theoretical results for validation. The results show that the presence of nanoparticles enhances the flow stability, and there exists a critical particle mass loading beyond which the flow is stable. As the shape factor of the velocity profile and the Reynolds number increase, the flow becomes more unstable. However, the flow becomes more stable with the increase of the particle mass loading. The wavenumber corresponding to the maximum of wave amplification becomes large with the increase of the shape factor of the velocity profile, and with the decrease of the particle mass loading and the Reynolds number. The variations of wave amplification with the Stokes number and the Knudsen number are not monotonic increasing or decreasing, and there exists a critical Stokes number and a Knudsen number with which the flow is relatively stable and most unstable,respectively, when other parameters remain unchanged. The perturbation with the first azimuthal mode makes the flow unstable more easily than that with the axisymmetric azimuthal mode. The wavenumbers corresponding to the maximum of wave amplification are more concentrated for the perturbation with the axisymmetric azimuthal mode.     

Large-scale instability of a fine-grained turbulent jet

The initial growth of a large scale perturbation on a fine-grained turbulent jet is studied via linear stability analysis. The turbulent jet is assumed to be homogeneous and isotropic with zero mean shear, and the inviscid stream outside the jet has a uniform velocity profile. The incremental Reynolds stress caused by the large scale perturbation is modeled by a viscoelastic constitutive equation, following the analysis of Crow (1968). It is found that the jet is always unstable to both sinuous and varicose types of perturbation, with the sinuous mode having a larger growth rate. In particular, short waves are always amplified, in contrast to the short wave stabilization at low speed found by Townsend (1966) for a purely elastic jet. The growth rates of these short waves are finite, and are smaller than those for the classical Kelvin-Helmholtz instability of an inviscid jet, but larger than those for the Hooper-Boyd (1983) instability of a viscous jet with continuous velocity profile.     

Subcritical spatial transition of swept Hiemenz flow

It is known from experimental investigations that the leading-edge boundary layer of a swept wing exhibits transition to turbulence at subcritical Reynolds numbers, i.e. at Reynolds numbers which lie below the critical Reynolds number predicted by linear stability theory. In the present work, we investigate this subcritical transition process by direct numerical simulations of a swept Hiemenz flow in a spatial setting. The laminar base flow is perturbed upstream by a pair of stationary counter-rotating vortex-like disturbances. This perturbation generates high- and low-speed streaks by a non-modal growth mechanism. Further downstream, these streaky structures exhibit a strong instability to secondary perturbations which leads to a breakdown to turbulence.The observed transition mechanism has strong similarities to by-pass transition mechanisms found for two-dimensional boundary layers. It can be shown that transition strongly depends on the amplitude of the primary perturbation as well as on the frequency of the secondary perturbation.     

Spatial structure of two-dimensional wave disturbances in a boundary layer

In [1] on the basis of a numerical integration of the Navier-Stokes equations the authors investigated the nonlinear evolution of two-dimensional disturbances of the traveling wave type in the boundary layer on a flat plate. The process of interaction of two waves with different wave numbers and initial amplitudes was examined. In this article the study of these interactions is continued. Special attention is paid to the spatial structure of the disturbances with respect to the cross-flow coordinate (with respect to the longitudinal coordinate the disturbances are assumed to be periodic) at various moments of time. It is shown that if the initial amplitude of one of the waves is sufficiently large, i.e., exceeds a certain threshold value, an undamped quasisteady regime is established during the interaction process. At lower amplitudes the process degenerates and the waves develop independently. In these two cases the evolution of the spatial distribution of the perturbation amplitudes is qualitatively different. In the first case the shape of the amplitude distribution varies only slightly with time, while in the second it depends importantly on the parameters of the wave numbers and the Reynolds number. When the parameters are such that one of the finite-amplitude waves is damped, its amplitude distribution rapidly evolves into the form characteristic of disturbances of the continuous spectrum of the linear stability problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 19–24, September–October, 1990.     

Stability of wave forms of motion of a viscous fluid layer under tangential stress

We study the stability of wave flow of a viscous incompressible fluid layer subjected to tangential stress and an inclined gravity force with respect to long-wave disturbances.An asymptotic solution is constructed for the equations of the disturbed motion and the problem is reduced to the study of a second-order ordinary differential equation. It is shown that after loss of stability by a Poiseuille flow the laminar nature of the flow is not destroyed, but the form of the free surface acquires a wave-like profile. The Poiseuille regime is stable for low Reynolds numbers. The critical Reynolds number for wave flow is found, and the stability and instability regions are determined.     

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.     

Large eddy simulation of boundary layer flow under cnoidal waves

Water waves in coastal areas are generally nonlinear, exhibiting asymmetric velocity profiles with different amplitudes of crest and trough. The behaviors of the boundary layer under asymmetric waves are of great significance for sediment transport in natural circumstances. While previous studies have mainly focused on linear or symmetric waves, asymmetric wave-induced flows remain unclear, particularly in the flow regime with high Reynolds numbers.Taking cnoidal wave as a typical example of asymmetric waves, we propose to use an infinite immersed plate oscillating cnoidally in its own plane in quiescent water to simulate asymmetric wave boundary layer. A large eddy simulation approach with Smagorinsky subgrid model is adopted to investigate the flow characteristics of the boundary layer. It is verified that the model well reproduces experimental and theoretical results. Then a series of numerical experiments are carried out to study the boundary layer beneath cnoidal waves from laminar to fully developed turbulent regimes at high Reynolds numbers, larger than ever studied before.Results of velocity profile, wall shear stress, friction coefficient, phase lead between velocity and wall shear stress, and the boundary layer thickness are obtained. The dependencies of these boundary layer properties on the asymmetric degree and Reynolds number are discussed in detail.     

An interpretation of the stability of a turbulent shear flow near a rough wall

A conventional, small perturbation, stability analysis has been applied to a fully developed turbulent shear flow in a parallel duct with rough walls. This is an attempt to detect the inherent state of flow stability to quasi-regular disturbances emanating from the surface roughness elements. The surface roughness is represented by the usual roughness Reynolds number; it is fed into the analysis through an assumed mean velocity profile valid between the viscous sublayer and the inner (wall) region. An eddy viscosity model is used to secure the equation closure and the final equation for the perturbation amplitude has been solved numerically using the techniques developed for the Orr-Sommerfeld equation.Within the domain of realistic flow conditions, and for a range of surface roughness amplitudes, a local minimum of stability in terms of the longitudinal wave number has been found. However, it is not implied that it is a minimum minomorum, as only a limited range of surface roughnesses has been tried.     

具有强SORET效应的充分发展的混合流体行进波对流

混合流体Rayleigh-Benard对流是研究对流稳定性,时空结构和非线性特性的典型模型之一。本文利用流体力学扰动方程组的数值模拟,讨论了偏离传导状态具有强SORET效应的混合流体行进波对流的温度场和浓度场的成长过程,分析了充分发展对流情况下的对流振幅,Nusselt数及混合参数与相对瑞利数的关系。并给出了行进波相速度对相对瑞利数的依赖关系。结果说明混合参数的曲线与行进波相速度的分布曲线是类似的。文末,给出了垂直速度,温度和浓度场的分布并讨论了相对瑞利数对场的分布及不同场之间的相位差的影响。     

Stability of the Hypersonic Shock Layer on a Flat Plate

The stability of hypersonic viscous gas flow in a shock layer in the neighborhood of a flat plate is considered. The stability of the velocity, temperature, density, and pressure profiles calculated on the basis of the complete viscous shock layer equations is investigated within the framework of the linear stability theory with allowance for the shock wave relations. The calculated perturbation growth rates and phase velocities are compared with the experimental data obtained by means of electron-beam fluorescence.     

Direct and adjoint global stability analysis of turbulent transonic flows over a NACA0012 profile

In this work, various turbulent solutions of the two‐dimensional (2D) and three‐dimensional compressible Reynolds averaged Navier–Stokes equations are analyzed using global stability theory. This analysis is motivated by the onset of flow unsteadiness (Hopf bifurcation) for transonic buffet conditions where moderately high Reynolds numbers and compressible effects must be considered. The buffet phenomenon involves a complex interaction between the separated flow and a shock wave. The efficient numerical methodology presented in this paper predicts the critical parameters, namely, the angle of attack and Mach and Reynolds numbers beyond which the onset of flow unsteadiness appears. The geometry, a NACA0012 profile, and flow parameters selected reproduce situations of practical interest for aeronautical applications. The numerical computation is performed in three steps. First, a steady baseflow solution is obtained; second, the Jacobian matrix for the RANS equations based on a finite volume discretization is computed; and finally, the generalized eigenvalue problem is derived when the baseflow is linearly perturbed. The methodology is validated predicting the 2D Hopf bifurcation for a circular cylinder under laminar flow condition. This benchmark shows good agreement with the previous published computations and experimental data. In the transonic buffet case, the baseflow is computed using the Spalart–Allmaras turbulence model and represents a mean flow where the high frequency content and length scales of the order of the shear‐layer thickness have been averaged. The lower frequency content is assumed to be decoupled from the high frequencies, thus allowing a stability analysis to be performed on the low frequency range. In addition, results of the corresponding adjoint problem and the sensitivity map are provided for the first time for the buffet problem. Finally, an extruded three‐dimensional geometry of the NACA0012 airfoil, where all velocity components are considered, was also analyzed as a Triglobal stability case, and the outcoming results were compared to the previous 2D limited model, confirming that the buffet onset is well detected. Copyright © 2014 John Wiley & Sons, Ltd.     

涂布流动的非线性阈值问题

本文研究了沿斜面流动薄层液体的非线性稳定性,即涂布流动的非线性稳定性问题。我们将周恒对平面Poiseuille流提出的弱非线性理论应用于涂布流动。文中对自由表面的世界条件提出了一个合理的简化方法,对亚临界时不同Reynolds数及扰动频率,求出了有限扰动的阈值。     

The effect of non-linear interaction between gas and particle velocities on the hydrodynamic stability in the Blasius boundary layer

A weakly non-linear stability analysis of two phase flow in the Blasius boundary layer has been carried out. Two mathematical models have been established based on the perturbation shape preserved assumption and linear stability model of two phase flow proposed by Stuart [On the non-linear mechanics of hydrodynamic stability, J. Fluid Mech. 4 (1958) 1-21] and Saffman [On the stability of laminar flow of dusty gas, J. Fluid Mech. 13 (1962) 120-128], respectively. The perturbation model and the perturbation energy balance equation are solved numerically with Chebyshev spectral method and artificial boundary condition. The numerical program adopted in the present study is verified by comparison with former works. The results show that the non-linear interaction between mean flow and perturbation reduces the growth rate of perturbation, while the non-linear interaction between particle phase and gas phase increases the growth rate of perturbation amplitude. The distortion of the mean flow caused by the Reynolds stress modifies the rate of transfer of energy from the mean flow to disturbance. The existence of particle alleviates the distortedness. The result also indicates that the weakly non-linear stability theory is consistent to linear stability theory, and the addition of fine and coarse particles reduces and increases the critical Reynolds number.     

Linear stability analysis in compressible, flat-plate boundary-layers

The stability problem of two-dimensional compressible flat-plate boundary layers is handled using the linear stability theory. The stability equations obtained from three-dimensional compressible Navier–Stokes equations are solved simultaneously with two-dimensional mean flow equations, using an efficient shoot-search technique for adiabatic wall condition. In the analysis, a wide range of Mach numbers extending well into the hypersonic range are considered for the mean flow, whereas both two- and three-dimensional disturbances are taken into account for the perturbation flow. All fluid properties, including the Prandtl number, are taken as temperature-dependent. The results of the analysis ascertain the presence of the second mode of instability (Mack mode), in addition to the first mode related to the Tollmien–Schlichting mode present in incompressible flows. The effect of reference temperature on stability characteristics is also studied. The results of the analysis reveal that the stability characteristics remain almost unchanged for the most unstable wave direction for Mach numbers above 4.0. The obtained results are compared with existing numerical and experimental data in the literature, yielding encouraging agreement both qualitatively and quantitatively.      

Surface instabilities of thin liquid film flow on a rotating disk

 Steady flow of a liquid jet from a nozzle onto the centre of a rotating disk is studied with a streak line method to determine the superficial velocity of the spreading liquid film. Good agreement is found with an asymptotic analysis of the unperturbed flow field. Experimentally, the liquid surface is always perturbed by surface waves which appear as regular spirals, steady in the laboratory system in the low Reynolds number range. It could be shown that wave formation is very sensitive to entrance conditions. Therefore, it is assumed that wave generation is an entrance effect which acts as periodic forcing on the forming liquid film. Wave velocities outside the entrance region are measured and proved to be in good agreement with the prediction of a linear stability theory, as long as the flow rate and entrance perturbations are small. At higher flow rates or stronger disturbances, the radial development of the wave velocities takes on the characteristics predicted by nonlinear stability theories and is in qualitative agreement with experiments performed on an inclined plane. Received: 15 January 1998/Accepted: 8 June 1998     

粘性金属中的正激波稳定性问题

G.H.Miller等把高压金属中的粘性激波作为强间断面处理,解析推论出:在大粘性系数条件下小扰动激波是不稳定的,物质粘性是导致失稳的因素。本文中针对平面正激波,认为高压金属中的粘性激波的物理量是连续变化的,利用线性稳定性理论,用数值解推论出:在有粘性条件下小扰动激波都是稳定的,物质粘性是致稳的因素。指出G.H.Miller等获得错误结论的原因在于:从无粘流动解推出的小扰动边界条件导致粘性激波小扰动增长。给出实验确定的小扰动速度梯度的边界条件,这样既可以把粘性正激波作为强间断面处理,也能够保证粘性正激波的稳定性。