The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.
This article studies the nonlinear evolution of disturbance waves in supersonic nonparallel boundary layer flows by using nonlinear parabolic stability equations （NPSE）. An accurate numerical method is developed to solve the equations and march the NPSE in a stable manner. Through computation,are obtained the curves of amplitude and disturbance shape function of harmonic waves. Especially are demonstrated the physical characteristics of nonlinear stability of various harmonic waves,including instantaneous stream wise vortices,spanwise vortices and Λ structure etc,and are used to study and analyze the mechanism of the transition process. The calculated results have evidenced the effectiveness of the proposed NPSE method to research the nonlinear stability of the supersonic boundary layers.
Formation and evolution of secondary streamwise vortices in the compressible transitional boundary layers over a flat plate are studied using a direct numerical simu- lation method with high-order accuracy and highly effective non-reflecting characteristic boundary conditions. Generation and development processes of the secondary streamwise vortices in the complicated transitional boundary flow are clearly analyzed based on the of numerical results, and the effects on the formation of the ring-like vortex that is vital to the boundary layer transition are explored. A new mechanism forming the ring-like vortex through the mutual effect of the primary and secondary streamwise vortices is expressed.
At the late stage of transitional boundary layers, the nonlinear evolution of the ring-like vortices and spike structures and their effects on the surrounding flow were studied by means of direct numerical simulation with high order accuracy. A spatial transition of the flat-plate boundary layers in the compressible flow was conducted. Detailed numerical results with high resolution clearly represented the typical vortex structures, such as ring-like vortices and so on, and induced ejection and sweep events. It was verified that the formation of spike structures in transitional boundary layers had close relationship with ring-like vortices. Especially, compared to the newly observed positive spike structure in the experiments, the same structure was found in the present numerical simulations, and the mechanism was also studied and analyzed.
The parabolized stability equations （PSEs） for high speed flows, especially supersonic and hypersonic flows, are derived and used to analyze the nonparallel boundary layer stability. The proposed numerical techniques for solving PSE include the following contents： introducing the efficiently normal transformation of the boundary layer, improving the computational accuracy by using a high-order differential scheme near the wall, employing the predictor-corrector and iterative approach to satisfy the important normalization condition, and implementing the stable spatial marching. Since the second mode dominates the growth of the disturbance in high Mach number flows, it is used in the computation. The evolution and characteristics of the boundary layer stability in the high speed flow are demonstrated in the examples. The effects of the nonparallelizm, the compressibility and the cooling wall on the stability are analyzed. And computational results are in good agreement with the relevant data.
通过研究一种基于函数值的(3,2)1阶二元有理插值样条函数中诸如边界插值、极限、解析和正则等性质,指出极限曲面是双曲抛物面,揭示了参数对这种插值曲面的影响.首先引入双8次矩阵表示的凸性判别函数,推导了判定插值曲面凸性的充要条件;然后根据该条件给出数值实例,展示如何适当选取参数实现有理插值样条曲面的局部保凸性.特别发现了这种插值曲面凸性在某些点处即使型值是凸的数据也是相对刚性的,并提出了插值曲面局部保凸的必要条件.最后还讨论了文献(Zhang Y,Duan Q,Twizell E H.Convexity control of a bivariate rational interpolating spline surfaces.Computers&Graphics,2007,31(5):679-687)中存在的部分计算问题.
