Numerical simulation of axisymmetric turbulent jets

The flow in axisymmetric turbulent jets is numerically simulated with the use of a semi-empirical second-order turbulence model including differential transport equations for the normal Reynolds stresses. Calculated results are demonstrated to agree with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 55–60, September–October, 2008.     

Verification of parabolized stability equations for its application to compressible boundary layers

Parabolized stability equations (PSE) were used to study the evolution of disturbances in compressible boundary layers.The results were compared with those ob- tained by direct numerical simulations (DNS),to check if the results from PSE method were reliable or not.The results of comparison showed that no matter for subsonic or supersonic boundary layers,results from both the PSE and DNS method agreed with each other reasonably well,and the agreement between temperatures was better than those between velocities.In addition,linear PSE was used to calculate the neutral curve for small amplitude disturbances in a supersonic boundary layer.Compared with those obtained by linear stability theory (LST),the situation was similar to those for incom- pressible boundary layer.     

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

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

Similarity solution of thermal boundary layers for laminar narrow axisymmetric jets

In this work, the similarity equation describing the thermal boundary layers of laminar narrow axisymmetric jets is derived based on boundary layer assumptions. The equation is solved exactly. Some properties of the thermal jet are discussed. By introducing new-defined non-dimensional coordinates, the similarity solution results in a “universal” format. The results can also be applied in the boundary layer problem of species diffusion.     

PSE as applied to problems of secondary instability in supersonic boundary layers

Parabolized stability equations (PSE) approach is used to investigate prob-lems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether the 2-D fundamental disturbance is of the first mode or second mode T-S wave. The variation of the growth rates of the 3-D sub-harmonic wave against its span-wise wave number and the amplitude of the 2-D fundamental wave is found to be similar to those found in incompressible boundary layers. But even as the amplitude of the 2-D wave is as large as the order 2%, the maximum growth rate of the 3-D sub-harmonic is still much smaller than the growth rate of the most unstable second mode 2-D T-S wave. Consequently, secondary instability is unlikely the main cause leading to transition in supersonic boundary layers.     

Stability of an axisymmetric swirled flow in a supersonic cocurrent stream for volume energy supply in the viscous vortex core

The results of numerical calculations of the stability of axisymmetric swirled flows in a viscous vortex embedded in a supersonic cocurrent stream with a constant circulation of the azimuthal velocity component are presented. The stability characteristics of the swirled three-dimensional viscous flow in the streamwise vortex are determined on the basis of the linearized system of Navier-Stokes equations for a viscous heat-conducting gas under the assumption that the basic undisturbed flow is locally plane-parallel. The disturbed flow stability is studied in the temporal formulation with respect to both symmetric and asymmetric three-dimensional waves traveling along the vortex axis and corresponding to both positive and negative values of the azimuthal wavenumber. It is shown that at external inviscid flow Mach numbers M = 2 and 3 thermal energy supply in a small region near the vortex axis leads to considerable restructuring of the basic undisturbed flow in the vicinity of the vortex core, the growth of the adverse pressure gradient along the vortex axis, and a significant change in the small perturbation stability and behavior.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 71–80. Original Russian Text Copyright © 2005 by Kazakov.     

Pseudospectral analysis of the interfacial stability of free jets with different velocity profiles

A double fluid model for a liquid jet surrounded by a coaxial gas stream was constructed. The interfacial stability of the model was studied by Chebyshev pseudospectral method for different basic velocity profiles. The physical variables were mapped into computational space using a nonlinear coordinates transformation. The general eigenvalues of the dispersion relation obtained are solved by QZ method, and the basic characteristics and their dependence on the flow parameters are analyzed.     

Numerical calculations of supersonic underexpanded jets in a cocurrent flow with the use of parabolized Navier-Stokes equations

Results of a numerical study of three-dimensional supersonic jets propagating in a cocurrent flow are described. Averaged parabolized Navier-Stokes equations are solved numerically on the basis of a developed scheme, which allows calculations in supersonic and subsonic flow regions to be performed in a single manner. A jet flow with a cocurrent flow Mach number 0.05 ⩽ M_∞ ⩽ 7.00 is studied, and its effect on the structure of the mixing layer is demonstrated. The calculated results are compared with available experimental and numerical data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 54–63, May–June, 2008.     

Vortex structure simulation for supersonic mixing layers using nonlinear PSE method

The method of nonlinear parabolized stability equations (PSE) is applied in the simulation of vortex structures in compressible mixing layer. The spatially-evolving unstable waves, which dominate the vortex structure, are investigated through spatial marching method. The instantaneous flow field is obtained by adding the harmonic waves to basic flow. The results show that T-S waves do not keep growing exponentially as the linear evolution, the energy transfer to high order harmonic modes, and that finally all harmonic modes get saturated due to nonlinear interaction. The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones, makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up, vortex pairing, and Λ vortex.     

The stability of viscous liquid jets in a swirling gas

Based on the linear analysis of stability, a dispersion equation is deduced which delineates the evolution of a general 3-dimensional disturbance on the free surface of an incompressible viscous liquid jet injected into a gas with swirl. Here, the dimensionless parameterJ _e is again introduced, in the meantime, another dimensionless parameterE called as circulation is also introduced to represent the relative swirling intensity. With respect to the spatial growing disturbance mode, the numerical results obtained from solving the dispersion equation reveal the following facts. First, at the same value ofE, in pace with the changing ofJ _e , the variation of disturbance and the critical disturbance mode still keep the same characters. Second, the present results are the same as that of S.P. Lin whenJ _e >1; but in the range ofJ _e <1, it's no more the case, the swirl decreases the axisymmetric disturbance, yet increases the asymmetric disturbance, furthermore the swirl may make the character of the most unstable disturbance mode changed (axisymmetric or asymmetric); the above action of the swirl becomes much stronger whenJ _e ≪1. The project supported by the National Natural Science Foundation of China     

Stress function method in shakedown analysis of axisymmetric shells

A self-equilibrated stress obtained from the stress functions of thin shells is used for the static shakedown theorem as a residual stress. In combination with the finite element method, a linear programming formulation of the shakedown analysis of axisymmetric shells is derived. The physical meaning of the stress function method is clear and its computing amount is small. Some examples of the plates and shells show that the method is reasonable and efficient.The project was supported by the National Natural Science Foundation of China     

Predictions of axisymmetric and two-dimensional impinging turbulent jets

The k − turbulence model and a version of a second-moment closure, modified to include the effect of pressure reflections from a solid surface, have been used as the basis of predictions of the flow that results from the orthogonal impingement of circular and two-dimensional (2-D) jets on a flat surface. Comparison of model predictions has been made with velocity measurements obtained in the stagnation and wall jet regions of the impinging flows. Results, in general, confirm the superiority of the Reynolds stress transport equation model for predicting mean and fluctuating velocities within the latter regions of such flows. In particular, modifications to the second-moment closure to account for the influence of the surface in distorting the fluctuating pressure field away from the wall successfully predict the damping of normal-to-wall velocity fluctuations throughout the impinging flows. In contrast, results derived from the eddy-viscosity-based approach do not, in general, accurately reproduce experimental observations.     

Pitot pressures of correctly-expanded and underexpanded free jets from axisymmetric supersonic nozzles

The structures of the axisymmetric free jets from supersonic nozzles with the exit Mach numbers of 1.5 and 2.0 are studied with special attention to the decay of the Pitot pressures downstream of the Mach disk. The Pitot pressure probe and schlieren method are used in the experiments to diagnose the flowfield. A TVD numerical method is also applied to the Euler equations, and the computed jet structures are compared with experiments. In the underexpanded jet, the experimentally obtained Pitot pressure near the jet centerline is found to substantially recover downstream of the Mach disk. By comparing the numerical computation, this phenomenon is thought to be caused by the turbulent momentum transfer to the central region from the region outside the slip line where the stagnation pressure loss is small.     

Linear and non‐linear stability analysis for finite difference discretizations of high‐order Boussinesq equations

This paper considers a method of lines stability analysis for finite difference discretizations of a recently published Boussinesq method for the study of highly non‐linear and extremely dispersive water waves. The analysis demonstrates the near‐equivalence of classical linear Fourier (von Neumann) techniques with matrix‐based methods for formulations in both one and two horizontal dimensions. The matrix‐based method is also extended to show the local de‐stabilizing effects of the non‐linear terms, as well as the stabilizing effects of numerical dissipation. A comparison of the relative stability of rotational and irrotational formulations in two horizontal dimensions provides evidence that the irrotational formulation has significantly better stability properties when the deep‐water non‐linearity is high, particularly on refined grids. Computation of matrix pseudospectra shows that the system is only moderately non‐normal, suggesting that the eigenvalues are likely suitable for analysis purposes. Numerical experiments demonstrate excellent agreement with the linear analysis, and good qualitative agreement with the local non‐linear analysis. The various methods of analysis combine to provide significant insight into the numerical behaviour of this rather complicated system of non‐linear PDEs. Copyright © 2004 John Wiley & Sons, Ltd.     

低速轴对称层流射流流动形态和失稳机制的实验研究

本文应用染色液和悬浮粒子显示方法,进一步实验研究轴对称层流射流的流动形态及其失稳机制。首次成功地在从一定口径的喷嘴流出的低速轴对称层流射流中观察到环形回流流动。给出了射流随速度演化及实验容器边界对其流动形态影响的显示照片,发现实际射流的轴对称波动及失稳过程正好对应射流在容器底部产生的环形旋涡的生长和破碎过程。本文认为由于实验空间有限尺度对流动的限制改变了原来射流的流动形态和流场空间的拓扑性质,射流与实际边界的相互作用对实际射流的失稳和转捩有重要的影响。     

Behavior of the vorticity vector in supersonic axisymmetric swirl flows behind a detonation wave

The behavior of the vorticity vector on a discontinuity surface arising in a supersonic nonuniform combustible gas flow with the formation of a shock or detonation wave is studied. In the general case, it is a vortex flow with prescribed distributions of parameters. It is demonstrated that the ratio of the tangential component of vorticity to density remains continuous in passing through the discontinuity surface, while the quantities proper become discontinuous. Results calculated for flow vorticity behind a steady-state detonation wave in an axisymmetric supersonic flow of a combustible mixture of gases are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 15–21, November–December, 2007     

Stability analysis method considering non-parallelism: EPSE method and its application

The e-N method is widely used in transition prediction. The amplitude growth rate used in the e-N method is usually provided by the linear stability theory (LST) based on the local parallel hypothesis. Considering the non-parallelism effect, the parabolized stability equation (PSE) method lacks local characteristic of stability analysis. In this paper, a local stability analysis method considering non-parallelism is proposed, termed as EPSE since it may be considered as an expansion of the PSE method. The EPSE considers variation of the shape function in the streamwise direction. Its local characteristic is convenient for stability analysis. This paper uses the EPSE in a strong non-parallel flow and mode exchange problem. The results agree well with the PSE and the direct numerical simulation (DNS). In addition, it is found that the growth rate is related to the normalized method in the non-parallel flow. Different results can be obtained using different normalized methods. Therefore, the normalized method must be consistent.     

Computation of axisymmetric confined jets in a diffuser

The paper reports on a numerical study of turbulent confined jets in a conical duct with a 5° divergence. The flow has a large ratio of jet to ambient velocities at the entrance so that it gives rise to strong recirculation. The calculations are carried out with a general finite volume method designed for calculating incompressible elliptic flows with complex boundaries. Turbulence is simulated by the standard κ–? model. The sensitivity of the solution to numerical discretization errors is examined using three convection schemes, i.e. hybrid central/upwind differencing, QUICK and SOUCUP, on two grids consisting of 68 × 50 and 102 × 82 points respectively. An examination is also made of the influence of inlet boundary conditions on the predicted flow field. The computed results are compared with experimental data for mean axial velocity, turbulent shear stress and turbulent kinetic energy profiles. It is shown that the calculations reproduce the essential features of the flow observed in the experiments.