非轴对称扰动下含悬浮固粒旋转射流场稳定性研究

考虑含有悬浮固粒的理想不可压旋转圆射流的运动方程，借助势流理论推得扰动增长率随径向波数变化的表达式，进而得到含悬浮固粒射流稳定性的修正瑞利稳定性准则，求出了非轴对称扰动下不同扰动阶数、固粒质量密度、固－气脉动速度比值、固－气脉动速度相位差、斯托克斯数时旋转射流的增长率曲线。然后根据修正瑞利稳定性准则分析给出了关于悬浮固粒的属性对旋转射流场稳定性影响的几个重要结论，为控制旋转射流场和后续发展提供了合理的依据。     

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

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

带展向曲率的超音速混合层中扰动演化的研究

超音速混合层的流动不稳定性较之亚音速或不可压的混合层大大减弱,为了提高混合效率, 通过数值模拟的方法分别研究了展向曲率、展向速度、来流马赫数等因素对混合效率所 起的作用. 计算结果表明：在给定展向速度的情况下,带有展向曲率的三维混合层,曲率越 大三维扰动增长率越大,而且法向的卷起范围也越大. 当展向曲率不为零时,展向速度的增 大也能有效地增强混合能力. 由流场中的高频扰动波产生的涡,在向下游发展过程中会有破 碎、拉伸,低频扰动波没有发现这一现象. 对于有展向曲率和展向速度的混合层,提高来流 马赫数时,流场中最不稳定扰动的增长率仍很大. 因此,这是一种提高混合层混合效率的新 途径.     

非线性强弱对轴对称射流拟序结构影响的研究

本文采用三维离散涡方法对轴对称不可压射流进行了数值模拟。由于Re数的大小体现出流场方程中非线性项的相对强弱,所以文中研究了不同Re数对射流中拟序结构的影响。文中提出并运用了一种有效的涡增加技术以提高涡方法的分辨率和精度,描绘了射流中展向涡的空间和时间发展,细致研究了涡的卷起、配对与合并现象。结果说明,在高Re数即非线性项相对较大的范围内,轴对称射流的发展主要受大尺度拟序结构的支配,而粘性则抑制大涡的产生与发展;Re数越大,非线性越强,同一时刻流场中大涡的数量越多,发展越充分。计算结果和实验一致,对实际射流的控制有重要意义。     

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

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

入口扰动频率对矩形射流影响的大涡模拟

为研究入口扰动对矩形射流的影响,对不同扰动频率下的三维矩形射流进行了大涡模拟,用分布投影法求解动量方程,对压力泊松方程采用FACR直接解法.计算结果预报的速度在横向的单峰分布和速度半宽沿流向的近似线性增长与实验结果吻合良好.计算结果表明扰动频率为射流的本征频率fe=0.22时射流的卷吸作用最强,而高频扰动则抑制卷吸作用;激励频率接近本征频率时射流提前出现速度在展向的马鞍型双峰分布,更早趋于圆射流形态.进一步的流场二维和三维瞬态涡结构图谱表明:上述现象与f=0.1和f=0.22时射流近场更频繁卷起、更丰富的大涡结构存在密切关系,正是这些大涡运动不断地在横向卷吸入外界流体,实现射流在横向的迅速扩展.     

金属射流失稳断裂的理论分析

基于Hamilton原理,提出一个包含射流强度、剪切、应变率效应、动力学黏性、表面张力和速度梯度等多因素耦合的金属射流拉伸运动方程,具体分析了各种失稳因素,并由数值解定量给出其影响大小,以及最不稳定波长与初始应变率乘积\lambda_{m}\dot{\varepsilon}_{0} 值的变化范围;给出了射流断裂的时间判据和近似理论公式,计算得到的 t_{b}-\dot{\varepsilon}_{0}曲线与射流实验点、Chou\&Carleone拟合公式三者符合较好.      

基于Busemann双翼构型的超音速导弹减阻技术研究

针对超音速飞行器在飞行过程中要承受的强激波带来的不利波阻,本文与传统单翼进行对比,分析了Busemann超音速双翼构型的减阻原理并充分利用了双翼间激波膨胀波的有利干涉和机翼厚度减小所带来的激波减弱效应。以常规气动布局的超音速导弹为研究对象,数值计算结果表明:设计巡航条件下,来流马赫数为2.5时,采用新型双翼气动布局能够使波阻减小42%。同时,为了消除非设计马赫数下Busemann双翼构型的壅塞问题,本文探索了一种转折变形翼面技术,计算结果表明:通过控制机翼前缘入口处和最大厚度处的面积比,该方案在非设计条件下能够基本消除阻力剧增问题。此外,在Busemann双翼基础上改进的上下翼非对称的超音速双翼构型可进一步改善实际有升力飞行条件下模型的气动性能,将所计算导弹模型巡航状态的升阻比提高了22%。综合以上结果表明,本文介绍的减阻技术可以为超声速导弹的研制和发展提供新的设计思路。     

一种适用于超音速边界层的湍流转捩模式

建立一种合理反映扰动模态和可压缩性影响的新型k-ω-γ转捩模式.其主要特点为:(1)假设脉动动能k 由湍流脉动部分和非湍流脉动部分组成,且在模化后者时采用了稳定性分析的结果; (2)在间歇因子γ方程的源项中,构造了具有``自动判断转捩起始位置'功能的函数; (3)通过构造新型的物面法向长度尺度,保证了模式中所有的表达式均由当地变量构成,可以方便地应用于现代CFD程序之中. 该模式在亚音速、超音速和高超音速条件下的边界层流动中进行了验证. 计算结果表明,该模式可应用于较宽马赫数范围内的自然转捩以及旁路转捩过程,所具有的捕捉流动转捩的性能优于国际上的现有模式.      

层流射流流动的简化Navier-Stokes方程(SNSE)解

本文利用十一种简化 Navier-Stokes 方程(SNSE) 求解已知Navier-Stokes(NS)方程准确解的层射流流动,表明:多数SNSE~([1-6])的解与NS方程的准确解不一致;少数SNSE~([7,8])的解与NS方程的准确解一致,文中在射流的喉部和拐点位置,给出几种SNSE解与准确解的相对偏差,并把粘性及惯性诸项加以定量比较,强调指出:按照边界层理论量级分析为Re~(1/2)和Re~1量级的惯性项以及Re~(-1/2)量级的粘性项具有重要影响;据此从力学角度论证了简化 NS 方程时,保留全部惯性项和合理取舍粘性项的必要性。     

Linear stability analysis of supersonic axisymmetric jets

Stabilities of supersonic jets are examined with different velocities, momentum thicknesses, and core temperatures. Amplification rates of instability waves at inlet are evaluated by linear stability theory (LST). It is found that increased velocity and core temperature would increase amplification rates substantially and such influence varies for different azimuthal wavenumbers. The most unstable modes in thin momentum thickness cases usually have higher frequencies and azimuthal wavenumbers. Mode switching is observed for low azimuthal wavenumbers, but it appears merely in high velocity cases. In addition, the results provided by linear parabolized stability equations show that the mean-flow divergence affects the spatial evolution of instability waves greatly. The most amplified instability waves globally are sometimes found to be different from that given by LST.     

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.     

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.     

Existence of time periodic solutions for a damped generalized coupled nonlinear wave equations

IntroductionInRef.[1 ] ,theauthorsestablishedtheuniqueexistenceofthesmoothsolutionforthefollowingcouplednonlinearequationsut=uxxx+buux+ 2vvx, (1 )vt=2 (uv) x. (2 )Thesewereproposedtodescribetheinteractionprocessofinternallongwaves.InRef.[2 ] ,ItoM .proposedarecursionoperatorbywhichheinferredthatEqs.(1 )and (2 )possesinfinitelymanysymmetriesandconstantsofmotion .InRef.[3 ] ,P .F .HeestablishedtheexistenceofasmoothsolutiontothesystemofcouplednonlinearKdVequation[4 ]ut=a(uxxx+buux) + 2bvvx,(…     

Internal waves in a two‐layer system using fully nonlinear internal‐wave equations

In order to understand the nonlinear effect in a two‐layer system, fully nonlinear strongly dispersive internal‐wave equations, based on a variational principle, were proposed in this study. A simple iteration method was used to solve the internal‐wave equations in order to solve the equations stably. The applicability of the proposed numerical computation scheme was confirmed to agree with linear dispersion relation theoretically obtained from variational principle. The proposed computational scheme was also shown to reproduce internal waves including higher‐order nonlinear effect from the analysis of internal solitary waves in a two‐layer system. Furthermore, for the second‐order numerical analysis, the balance of nonlinearity and dispersion was found to be similar to the balance assumed in the KdV theory and the Boussinesq‐type equations. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

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.     

Nonlinear evolution analysis of T-S disturbance wave at finite amplitude in nonparallel boundary layers

ＩｎｔｒｏｄｕｃｔｉｏｎＮｏｎｌｉｎｅａｒｓｔａｂｉｌｉｔｙｐｒｏｂｌｅｍｉｎｎｏｎｐａｒａｌｌｅｌｂｏｕｎｄａｒｙｌａｙｅｒｓｎｅｅｄｓｍａｉｎｌｙｃｏｎｓｉｄｅｒｔｈｅｅｖｏｌｕｔｉｏｎｏｆｕｎｓｔａｂｉｌｉｔｙＴ＿Ｓｗａｖｅａｎｄｉｎｔｅｒａｃｔｉｏｎｂｅｔｗｅｅｎｗａｖｅｓｏｆｄｉｆｆｅｒｅｎｔｆｒｅｑｕｅｎｃｉｅｓ.Ｉｔｉｓｖｅｒｙｉｍｐｏｒｔａｎｔｔｏｓｉｍｕｌａｔｅａｃｃｕｒａｔｅｌｙｓｐａｔｉａｌｅｖｏｌｕｔｉｏｎｒｏｌｅｏｆｄｉｓｔｕｒｂａｎｃｅｉｎｂｏｕｎｄａｒｙｌａｙｅｒｓ,…     

On a new algorithm of constructing solitary wave solutions for systems of nonlinear evolution equations in mathematical physics

ＩｎｔｒｏｄｕｃｔｉｏｎＳｏｌｖｉｎｇｎｏｎｌｉｎｅａｒｅｑｕａｔｉｏｎｓｉｓａｌｗａｙｓａｖｅｒｙｉｎｔｅｒｅｓｔｉｎｇｓｕｂｊｅｃｔｆｏｒｍａｔｈｅｍａｔｉｃｉａｎａｎｄｐｈｙｓｉｃｉａｎ ,ｉｎｐａｒｔｉｃｕｌａｒ,ｓｏｌｉｔａｒｙｗａｖｅｓｏｌｕｔｉｏｎｓｆｏｒｎｏｎｌｉｎｅａｒｅｑｕａｔｉｏｎｓａｒｅｏｆｂｏｔｈｔｈｅｏｒｅｔｉｃａｌａｎｄｐｒａｃｔｉｃａｌｉｍｐｏｒｔａｎｃｅ.Ｒｅｃｅｎｔｌｙ ,Ｙａｎ[1]ｏｂｔａｉｎｅｄａｔｒａｎｓｆｏｒｍａｔｉｏｎｄｉｒｅｃｔｌｙｆｒｏｍｔｈｅｆａｍｏ…     

Gravity wave radiation from unsteady rotational flow in an f-plane shallow water system

Spontaneous gravity wave radiation from an unsteady rotational flow is investigated numerically in an f-plane shallow water system. Unlike the classical Rossby adjustment problem, where free development of an initially unbalanced state is investigated, we consider development of a barotropically unstable zonal flow which is initially balanced but maintained by zonal mean forcing. Gravity waves are continuously radiated from a nearly balanced rotational flow region even when the Froude number is so small that balance dynamics is thought to be a good approximation for the full system. The source of gravity waves is discussed by analogy with the theory of aero-acoustic sound wave radiation (the Lighthill theory). It is shown that the source regions correspond to regions of strong rotational flow. The gradual change of rotational flow causes gravity wave radiation. We propose an approximation for these strong sources on the assumption that the dominant flow in the jet region is non-divergent rotational flow. In addition, we calculate the zonally symmetric component of gravity waves far from the source regions, solving the Lighthill equation. Using scaling analyses for perturbations, these gravity waves can be calculated with only one approximated source term that is related to the latitudinal gradient of the fluid depth and the latitudinal mass flux. In spite of its simplicity, this approximation not only explains the physical cause of gravity wave radiation, but gives an amount of source close to that obtained by classical approximation derived from vortical motion.