空间发展圆管转捩的直接数值模拟

利用直接数值模拟研究圆管流动中由局部壁面引入的周期性吹吸(PSB)扰动沿流向的空间发展,流动的雷诺数Re选定为3000.在临界幅值的PSB扰动下,在较短的圆管内,圆管中的扰动沿流向快速增长,在足够长的圆管内,扰动沿流向持续增长发生转捩,流动发展到湍流阶段.     

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.     

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.     

Spatial direct numerical simulation of high-speed boundary-layer flows part I: Algorithmic considerations and validation

A highly accurate algorithm for the direct numerical simulation (DNS) of spatially evolving high-speed boundary-layer flows is described in detail and is carefully validated. To represent the evolution of instability waves faithfully, the fully explicit scheme relies on non-dissipative high-order compact-difference and spectral collocation methods. Several physical, mathematical, and practical issues relevant to the simulation of high-speed transitional flows are discussed. In particular, careful attention is paid to the implementation of inflow, outflow, and far-field boundary conditions. Four validation cases are presented, in which comparisons are made between DNS results and results obtained from either compressible linear stability theory or from the parabolized stability equation (PSE) method, the latter of which is valid for nonparallel flows and moderately nonlinear disturbance amplitudes. The first three test cases consider the propagation of two-dimensional second-mode disturbances in Mach 4.5 flat-plate boundary-layer flows. The final test case considers the evolution of a pair of oblique second-mode disturbances in a Mach 6.8 flow along a sharp cone. The agreement between the fundamentally different PSE and DNS approaches is remarkable for the test cases presented.     

Closed-loop control of a free shear flow: a framework using the parabolized stability equations

In this study the parabolized stability equations (PSE) are used to build reduced-order-models (ROMs) given in terms of frequency and time-domain transfer functions (TFs) for application in closed-loop control. The control law is defined in two steps; first it is necessary to estimate the open-loop behaviour of the system from measurements, and subsequently the response of the flow to an actuation signal is determined. The theoretically derived PSE TFs are used to account for both of these effects. Besides its capability to derive simplified models of the flow dynamics, we explore the use of the TFs to provide an a priori determination of adequate positions for efficiently forcing along the direction transverse to the mean flow. The PSE TFs are also used to account for the relative position between sensors and actuators which defines two schemes, feedback and feedforward, the former presenting a lower effectiveness. Differences are understood in terms of the evaluation of the causality of the resulting gain, which is made without the need to perform computationally demanding simulations for each configuration. The ROMs are applied to a direct numerical simulation of a convectively unstable 2D mixing layer. The derived feedforward control law is shown to lead to a reduction in the mean square values of the objective fluctuation of more than one order of magnitude, at the output position, in the nonlinear simulation, which is accompanied by a significant delay in the vortex pairing and roll-up. A study of the robustness of the control law demonstrates that it is fairly insensitive to the amplitude of inflow perturbations and model uncertainties given in terms of Reynolds number variations.     

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.     

NUMERICAL INVESTIGATION OF EVOLUTION OF DISTURBANCES IN SUPERSONIC SHARP CONE BOUNDARY LAYERS

The spatial evolution of 2-D disturbances in supersonic sharp cone boundary layers was investigated by direct numerical simulation (DNS) in high order compact difference scheme. The results suggested that, although the normal velocity in the sharp cone boundary layer was not small, the evolution of amplitude and phase for small amplitude disturbances would be well in accordance with the results obtained by the linear stability theory (LST) which supposes the flow was parallel. The evolution of some finite amplitude disturbances was also investigated, and the characteristic of the evolution was shown. Shocklets were also found when the amplitude of disturbances increased over some value.     

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.     

Nonlinear evolution of Klebanoff type second mode disturbances in supersonic flat-plate boundary layer简

Studying the evolution of 3D disturbances is of crucial theoretical importance for understanding the transition process. The present study concerns the nonlinear evolution of second mode unstable disturbances in a supersonic boundary layer by the numerical simulation, and discusses the selectivity of 3D disturbances and possibility to transition. The results indicate that a Klebanoff type nonlinear interaction between 2D and 3D disturbances with the same frequency may amplify a band of 3D disturbances centered at a finite spanwise wavenumber. That is, certain 3D disturbances can be selectively and rapidly amplified by the unstable 2D disturbances, and certain small-scale 3D structures will appear.     

Mechanism of transition in a hypersonic sharp cone boundary layer with zero angle of attack

Firstly, the steady laminar flow field of a hypersonic sharp cone boundary layer with zero angle of attack was computed.Then,two groups of finite amplitude T-S wave disturbances were introduced at the entrance of the computational field,and the spatial mode transition process was studied by direct numerical simulation (DNS) method. The mechanism of the transition process was analyzed.It was found that the change of the stability characteristics of the mean flow profile was the key issue.Furthermore,the characteristics of evolution for the disturbances of different modes in the hypersonic sharp cone boundary layer were discussed.     

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.     

Spatial direct numerical simulation of high-speed boundary-layer flows Part II: Transition on a cone in Mach 8 flow

The laminar breakdown of the boundary-layer flow of an axisymmetric sharp cone in a Mach 8 flow is simulated by a synergistic approach that combines the parabolized stability equation (PSE) method and spatial direct numerical simulation (DNS). The transitional state is triggered by a symmetric pair of oblique second-mode disturbances whose nonlinear interactions generate strong streamwise vorticity, which leads in turn to severe spanwise variations in the flow and eventual laminar breakdown. The PSE method is used to compute the weakly and moderately nonlinear initial stages of the transition process and, thereby, to derive a harmonically rich inflow condition for the DNS. The strongly nonlinear and laminar-breakdown stages of transition are then computed by well-resolved DNS, with a highly accurate algorithm that exploits spectral collocation and high-order compact-difference methods. Evolution of the flow is presented in terms of modal energies, mean quantities (e.g., skin friction), Reynolds stresses, turbulent kinetic energy, and flow visualization. The numerical test case is an approximate computational analog of one of the few stability experiments performed for hypersonic boundary-layer flows. Comparisons and contrasts are drawn between the experimental and the computational results. Rope-like waves similar to those observed in schlieren images of high-speed transitional flows are also observed in the numerical experiment and are shown to be visual manifestations of second-mode instability waves.This research was supported under NASA Contracts NAS1-19831 and NAS1-20059 for the first and second authors, respectively.     

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.     

波包在后掠翼三维边界层中的演化特征

研究了点源产生的孤立波包在后掠平板三维边界层中的演化特征．理论计算的波包增长路径和ＮｉｔＳＣｈｋｅ－Ｋｏｗｓｔｙ所做流场显示得到的条纹结构一致,证明条纹是波包的等相位线;波包方程的渐近解表明在三维边界层中ｅＮ方法应沿着为实数的方向积分;ｅＮ方法过高预测了扰动幅值的增长率,与Ｄｅｒｈｌｅ的实验结论一致．     

Selective enhancement of oblique waves caused by finite amplitude second mode in supersonic boundary layer

Nonlinear interactions of the two-dimensional(2D) second mode with oblique modes are studied numerically in a Mach 6.0 flat-plate boundary layer, focusing on its selective enhancement effect on amplification of different oblique waves. Evolution of oblique modes with various frequencies and spanwise wavenumbers in the presence of 2D second mode is simulated successively, using a modified parabolized stability equation(PSE) method, which is able to simulate interaction of two modes with different frequencies efficiently. Numerical results show that oblique modes in a broad band of frequencies and spanwise wavenumbers can be enhanced by the finite amplitude 2D second mode instability wave. The enhancement effect is accomplished by interaction of the 2D second mode, the oblique mode, and a forced mode with difference frequency. Two types of oblique modes are found to be more amplified, i.e., oblique modes with frequency close to that of the 2D second mode and low-frequency first mode oblique waves. Each of them may correspond to one type of transition routes found in transition experiments. The spanwise wavenumber of the oblique wave preferred by the nonlinear interaction is also determined by numerical simulations.     

三维扰动波的非平行边界层稳定性研究

导出了三维扰动波的原始变量形式的抛物化稳定性方程（PSE）,研究了三维空间模态TS波的非平行边界层稳定性问题.采用了法向四阶紧致格式,以提高计算精度.通过给出不会导致奇性的坐标变换、修改外边界条件以及克服平行流初始值的瞬态影响和推进步长的限制,保证了计算的数值稳定.用补全元素带状矩阵法求解块三对角矩阵,大大提高了速度.计算结果清楚地显示了三维扰动波的演化过程和非平行性对边界层稳定性的影响,特别是,观察到非平行性对三维扰动波的影响,有时会使其稳定性出现逆转的现象.还研究了逆压梯度的作用.算例的结果与其他结果符合良好.     

Formation and evolution of the streamwise vortices in mixing layers

The evolution of the three-dimensional time-developing mixing layer is simulated numerically using the pseudo-spectral method. The initial perturbations used in this study consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances is made. In one case the results are similar to that of Rogers & Moser^[1], whereas a different way in which the quadrupole forms and sudden expansion of the rib are observed in another case. The simulation also confirms that the stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results are presented. The project supported by the Zhejiang Province Natural Science Special Fund for Youth Scientistis' Cultivation.     

On the Secondary Instability of Three-Dimensional Boundary Layers

One of the possible transition scenarios in three-dimensional boundary layers, the saturation of stationary crossflow vortices and their secondary instability to high-frequency disturbances, is studied using the Parabolized Stability Equations (PSE) and Floquet theory. Starting from nonlinear PSE solutions, we investigate the region where a purely stationary crossflow disturbance saturates for its secondary instability characteristics utilizing global and local eigenvalue solvers that are based on the Implicitly Restarted Arnoldi Method and a Newton–Raphson technique, respectively. Results are presented for swept Hiemenz flow and the DLR swept flat plate experiment. The main focuses of this study are on the existence of multiple roots in the eigenvalue spectrum that could explain experimental observations of time-dependent occurrences of an explosive growth of traveling disturbances, on the origin of high-frequency disturbances, as well as on gaining more information about threshold amplitudes of primary disturbances necessary for the growth of secondary disturbances. Received 13 July 1998 and accepted 7 July 2000     

Spatial direct numerical simulation of boundary-layer transition mechanisms: Validation of PSE theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to evaluate critically the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave breakdown. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the two-dimensional test problem. This discrepancy is attributed to far-field boundary-condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.     

The 3D Transition to Turbulence in Wake Flows by Means of Direct Numerical Simulation

The three-dimensional transition to turbulence in flows around bodies of non-rectangular configuration has been analysed physically by performing direct numerical simulation to solve the system of Navier-Stokes equations. The successive stages of 3D transition, beyond the first bifurcation, have been detected first in the incompressible regime, for a circular cylinder configuration. The generation of streamwise vorticity, organised according to spanwise periodic cells has been associated with the development of large-scale coherent spanwise undulations of the originally rectilinear (nominally 2D) alternating vortex rows. The wavelengths of these undulations have been determined as a function of Reynolds number. As this parameter increases, a further inherent change of the flow transition is obtained and analysed, the natural vortex dislocations pattern. Beyond this change, the increase of Reynolds number yields an abrupt shortening of the spanwise wavelength and the flow undergoes another transition step, whose critical Reynolds number is evaluated by the present DNS approach in association with the Ginzburg-Landau model. Therefore, the linear and non-linear parts of the flow transition have been quantified by means of the amplitude evolution versus time obtained by the present DNS, in conjunction with the mentioned global oscillator model. This revised version was published online in July 2006 with corrections to the Cover Date.