Improvement for expansion of parabolized stability equation method in boundary layer stability analysis

An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber α and its streamwise gradient dα/dx are unknown variables. This eigenvalue problem is solved for the eigenvalue dα/dx with an initial α, and the correction of α is performed with the conservation relation used in the PSE. The i EPSE is validated in several compressible and incompressible boundary layers. The computational results show that the prediction accuracy of the i EPSE is significantly higher than that of the ESPE, and it is in excellent agreement with the PSE which is regarded as the baseline for comparison. In addition, the unphysical multiple eigenmode problem in the EPSE is solved by using the i EPSE. As a local non-parallel stability analysis tool, the i EPSE has great potential application in the eNtransition prediction in general three-dimensional boundary layers.     

Applications of EPSE method for predicting crossflow instability in swept-wing boundary layers

The nth-order expansion of the parabolized stability equation(EPSEn) is obtained from the Taylor expansion of the linear parabolized stability equation(LPSE) in the streamwise direction. The EPSE together with the homogeneous boundary conditions forms a local eigenvalue problem, in which the streamwise variations of the mean flow and the disturbance shape function are considered. The first-order EPSE(EPSE1) and the second-order EPSE(EPSE2) are used to study the crossflow instability in the swept NLF(2)-0415 wing boundary layer. The non-parallelism degree of the boundary layer is strong. Compared with the growth rates predicted by the linear stability theory(LST),the results given by the EPSE1 and EPSE2 agree well with those given by the LPSE.In particular, the results given by the EPSE2 are almost the same as those given by the LPSE. The prediction of the EPSE1 is more accurate than the prediction of the LST, and is more efficient than the predictions of the EPSE2 and LPSE. Therefore, the EPSE1 is an efficient e~N prediction tool for the crossflow instability in swept-wing boundary-layer flows.     

Self-consistent parabolized stability equation (PSE) method for compressible boundary layer

The parabolized stability equation (PSE) method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulation, for nonlinear problems, the complex wave number of each Fourier mode is determined by the so-called phase-locked rule, which results in non-self-consistency in the wave numbers. In this paper, a modification is proposed to make it self-consistent. The main idea is that, instead of allowing wave numbers to be complex, all wave numbers are kept real, and the growth or decay of each mode is simply manifested in the growth or decay of the modulus of its shape function. The validity of the new formulation is illustrated by comparing the results with those from the corresponding direct numerical simulation (DNS) as applied to a problem of compressible boundary layer with Mach number 6.     

Nonparallel stability of the flow in an axially rotating pipe

The linear stability of the developing flow in an axially rotating pipe is analyzed using parabolized stability equations (PSE). The results are compared with those obtained from a near-parallel stability approximation that only takes into account the axial variation of the basic flow. Though the PSE results obviously coincide with the near-parallel ones far downstream, when the flow has reached a Hagen-Poiseuille axial velocity profile with superimposed solid-body rotation, they differ significantly in the developing region. Therefore, the onset of instability strongly depends on the axial evolution of the perturbations. The PSE results are also compared with experimental data from Imao et al. [Exp. Fluids 12 (1992) 277], showing a good agreement in the frequencies and wavelengths of the unstable disturbances, that take the form of spiral waves. Finally, a simple method for detecting one of the conditions to characterize the onset of absolute instability using PSE is given.     

Applications of parabolized stability equation for predicting transition position in boundary layers

The phenomenon of laminar-turbulent transition exists universally in nature and various engineering practice.The prediction of transition position is one of crucial theories and practical problems in fluid mechanics due to the different characteristics of laminar flow and turbulent flow.Two types of disturbances are imposed at the entrance,i.e.,identical amplitude and wavepacket disturbances,along the spanwise direction in the incompressible boundary layers.The disturbances of identical amplitude are consisted of one two-dimensional(2D) wave and two three-dimensional(3D) waves.The parabolized stability equation(PSE) is used to research the evolution of disturbances and to predict the transition position.The results are compared with those obtained by the numerical simulation.The results show that the PSE method can investigate the evolution of disturbances and predict the transition position.At the same time,the calculation speed is much faster than that of the numerical simulation.     

SHPB试验中岩石试件的端面不平行修正

为研究短圆柱体岩石试件端面不平行对岩石动力学特性测试结果的影响,采用有限元分析软件LS-DYNA对9种端面不平行度和5种杨氏模量的岩石试件开展SHPB(split Hopkinson pressure bar)试验数值模拟,对岩石选用HJC(Holmquist-Johnson-Cook)本构模型。数值模拟结果表明,当端面不平行度在0.40%以内时,端面不平行对动态应力测试结果的影响可忽略不计;但对动态应变测试结果的影响较大。当杨氏模量一定时,平均应变率测试误差和峰值应变测试误差随端面不平行度增大呈线性增大;当端面不平行度一定时,平均应变率测试误差和峰值应变测试误差随杨氏模量增大也呈线性增大。对数值模拟得到的平均应变率测试误差和峰值应变测试误差实施二元线性回归分析,提出了SHPB试验中端面不平行岩石试件平均应变率和峰值应变的修正公式。     

Higher order stability effects in a natural convection boundary layer over a vertical heated wall

 The stability of a laminar boundary layer flow under natural convection on a vertical isothermally heated wall is studied analytically. The analysis is performed by using two different two-dimensional linear models: (1) The non-parallel flow model in which the steady mean flow as well as the disturbance amplitude functions can change in the streamwise direction; (2) The parallel flow model in which the effects of the mean flow and disturbance changes in the streamwise direction are neglected. The linear non-parallel stability analysis is based on the so-called parabolised stability equations (PSEs) which have been successfully applied to the stability analysis of forced convection boundary layers. In this study the PSE equations are applied to natural convection boundary layers in order to show the difference between parallel and non-parallel stability analysis. A second part of this study deals with the effects of variable properties, which are always present in natural convection flows. They are analysed by an extended version of the Orr–Sommerfeld equation (EOSE). Received on 31 May 2000     

A minimal composite theory for stability of non-parallel compressible boundary-layer flow

A new minimal composite theory that extends the approach of Govindarajan and Narasimha [1] is proposed here for 2D non-parallel compressible boundary-layer stability subject to 3D disturbances. The mean profiles are obtained from Hortons analysis, which provides a good approximation for a large range of Prandtl numbers at non-zero pressure gradients. In the lowest order, all effects of order lower than O(R^-2/3) anywhere in the boundary-layer are included, R being the local boundary-layer Reynolds number; the resulting non-parallel formulation yields a set of four ordinary differential equations, as compared to the five coupled equations of classical parallel flow theory of Mack [2]. The largest effect on stability of flow non-parallelism is found to be due to the wall-normal advection of velocity and temperature disturbance quantities by the mean flow. The present theory shows that the bulk viscosity, invariably included in compressible stability theories, is irrelevant at the lowest order. In comparison with the full [O(R^-1)] non-parallel theory, the present theory is marginally better than the parallel flow theory. PACS 03.50.De, 04.20-q, 42.65-k     

Sensitivity Analysis Using Adjoint Parabolized Stability Equations for Compressible Flows

An input/output framework is used to analyze the sensitivity of two- and three-dimensional disturbances in a compressible boundary layer for changes in wall and momentum forcing. The sensitivity is defined as the gradient of the kinetic disturbance energy at a given downstream position with respect to the forcing. The gradients are derived using the parabolized stability equations (PSE) and their adjoint (APSE). The adjoint equations are derived in a consistent way for a quasi-two-dimensional compressible flow in an orthogonal curvilinear coordinate system. The input/output framework provides a basis for optimal control studies. Analysis of two-dimensional boundary layers for Mach numbers between 0 and 1.2 show that wall and momentum forcing close to branch I of the neutral stability curve give the maximum magnitude of the gradient. Forcing at the wall gives the largest magnitude using the wall normal velocity component. In case of incompressible flow, the two-dimensional disturbances are the most sensitive ones to wall inhomogeneity. For compressible flow, the three-dimensional disturbances are the most sensitive ones. Further, it is shown that momentum forcing is most effectively done in the vicinity of the critical layer. This revised version was published online in July 2006 with corrections to the Cover Date.     

二维抛物化稳定性方程的特征和次特征

特征分析表明 :对原始扰动量的抛物化稳定性方程组 (PSE) ,它在亚、超音速区分别具有椭圆和抛物特性 ,给出PSE特征对马赫数的依赖关系 ,阐明PSE仅把信息对流 扩散传播特性抛物化 ,而保留了信息对流 扰动传播特性 ,因此PSE应称为扩散抛物化稳定性方程 (DPSE)     

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations ( PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub- characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic., respectively . The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories , the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time , the methods of removing the remained ellipticity are further obtained from the nonlinear PSE .     

Spatial stability of the incompressible corner flow

The linear spatial stability of the incompressible corner flow under pressure gradient has been studied. A self-similar form has been used for the mean flow, which reduces the related problem to the solution of a two-dimensional problem. The stability problem was formulated using the parabolised stability equations (PSE) and results were obtained for the viscous modes at medium and high frequencies. The related N-factors indicate that the flow is stable at these frequencies, but probably unstable for small frequencies. Furthermore the inviscid mode for each mean flow was obtained and the results indicate that its importance increases considerably with an increase in the adverse pressure gradient. Finally the dependence of the stability characteristics on the extent of the domain is also considered.     

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.     

RBF collocation method and stability analysis for phononic crystals

The mesh-free radial basis function (RBF) collocation method is explored to calculate band structures of periodic composite structures. The inverse multi-quadric (MQ), Gaussian, and MQ RBFs are used to test the stability of the RBF collocation method in periodic structures. Much useful information is obtained. Due to the merits of the RBF collocation method, the general form in this paper can easily be applied in the high dimensional problems analysis. The stability is fully discussed with different RBFs. The choice of the shape parameter and the effects of the knot number are presented.     

Studies on nonlinear stability of three-dimensional H-type disturbance

The three-dimensional H-type nonlinear evolution process for the problem of boundary layer stability is studied by using a newly developed method called parabolic stability equations (PSE). The key initial conditions for sub-harmonic disturbances are obtained by means of the secondaryinstability theory. The initial solutions of two-dimensional harmonic waves are expressed in Landau expansions. The numerical techniques developed in this paper, including the higher order spectrum method and the more effective algebraic mapping for dealing with the problem of an infinite region, increase the numerical accuracy and the rate of convergence greatly. With the predictor-corrector approach in the marching procedure, the normalization, which is very important for PSE method, is satisfied and the stability of the numerical calculation can be assured. The effects of different pressure gradients, including the favorable and adverse pressure gradients of the basic flow, on the “H-type“ evolution are studied in detail. The results of the three-dimensional nonlinear “H-type“ evolution are given accurately and show good agreement with the data of the experiment and the results of the DNS from the curves of the amplitude variation, disturbance velocity profile and the evolution of velocity.     

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.     

A new method for computing turbulence in compressible laminar-turbulent transition and boundary layers--PSE＋DNS

A new method for computing laminar-turbulent transition and turbulence in compressible boundary layers is proposed. It is especially useful for computation of laminar-turbulent transition and turbulence starting from small-amplitude disturbances. The laminar stage, up to the beginning of the breakdown in laminar-turbulent transition, is computed by parabolized stability equations (PSE). The direct numerical simulation (DNS) method is used to compute the transition process and turbulent flow, for which the inflow condition is provided by using the disturbances obtained by PSE method up to that stage. In the two test cases incfuding a subsonic and a supersonic boundary layer, the transition locations and the turbulent flow obtained with this method agree well with those obtained by using only DNS method for the whole process. The computational cost of the proposed method is much less than using only DNS method.     

PSE as applied to problems of transition in compressible boundary layers

A new idea of using the parabolized stability equation (PSE) method to predict laminar-turbulent transition is proposed. It is tested in the prediction of the location of transition for compressible boundary layers on flat plates, and the results are compared with those obtained by direct numerical simulations (DNS). The agreement is satisfactory, and the reason for this is that the PSE method faithfully reproduces the mechanism leading to the breakdown process in laminar-turbulent transition, i. e., the modification of mean flow profile leads to a remarkable change in its stability characteristics.     

Boundary-layer receptivity to external disturbances using multiple scales

Non-homogeneous multiple scales are introduced to solve the resonant problem of non-parallel boundary-layer receptivity originating from the quadratic mixing of environmental disturbances. The resulting algorithm is computationally inexpensive and can be efficiently included in industrial codes for transition prediction. The mutual interactions between acoustic wave, vorticity wave, wall vibration and wall roughness are discussed in detail and the receptivity coefficient, which relates the amplitude of the excited wave to the amplitude of the exciting sources, is computed. The largest effect is found for the interaction between acoustic waves and wall roughness perturbations. Other coupling mechanisms are less effective. By comparing parallel and non-parallel results, it is found that flow non-parallelism can play a non-negligible role even in Blasius’ boundary layer, although the largest effects are evident for the three-dimensional boundary layer over an infinite swept wing. For the particular case of wall roughness—wall vibration mixing, the velocity disturbance is shown to be exactly equal to the velocity perturbation induced by wall roughness alone on a wall vibrating in the normal direction.     

Accuracy analysis of predicted velocity profiles of laminar duct flow with entropy generation method

The objective of this work is to estimate the accuracy of a predicted velocity profile which can be gained from experimental results, in comparison with the exact ones by the methodology of entropy generation. The analysis is concerned with the entropy generation rate in hydrodynamic, steady, laminar, and incompressible flow for Newtonian fluids in the insulated channels of arbitrary cross section. The entropy generation can be calculated from two local and overall techniques. Adaptation of the results of these techniques depends on the used velocity profile. Results express that in experimental works, whatever the values of local and overall entropy generation rates are close to each other, the results are more accuracy. In order to extent the subject, different geometries have been investigated. Also, the influence studied, and the distribution of volumetric geometries is drawn. of geometry on the entropy generation rate is local entropy generation rate for the selected geometries is drawn.