二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

潜入喷管背壁区流场影响因素实验研究

应用相位多普勒粒子分析仪 (PDPA)对含潜入喷管、通道为矩形的固体火箭发动机冷流模型内流场进行测量 ,研究了雷诺数和背壁空腔形状对背壁区流场的影响。结果表明 ,随着雷诺数增大气流分离点向下游移动 ,回流速度增大 ,椭圆形后封头、直线形喷管背壁构成的背壁区回流速度比圆形后封头、台阶形喷管背壁大。对于潜入深度大的喷管 ,气流在燃烧室壁面分离 ,背壁区形成强旋涡 ;而潜入深度小的喷管 ,分离点会后移至后封头上 ,在背壁区后部形成弱旋涡。     

Self-sustained oscillation for compressible cylindrical cavity flows

The presence of a cavity changes the mean and fluctuating pressure distributions. Similarities are observed between a cylindrical cavity and a rectangular cavity for a compressible flow.The type of cavity flow field depends on the diameter-to-depth ratio and the length-to-depth ratio.The feedback loop is responsible for the generation of discrete acoustic tones. In this study, the selfsustained oscillation for a compressible cylindrical cavity flow was investigated experimentally. For open-type cavities, the power spectra show that the strength of resonance depends on the diameterto-depth ratio(4.43–43.0) and the incoming boundary layer thickness-to-depth ratio(0.72–7.0). The effective streamwise length is used as the characteristic length to estimate the Strouhal number. At higher modes, there is a large deviation from Rossiter's formula for rectangular cavities. The gradient-based searching method was used to evaluate the values of the empirical parameters. Less phase lag and a lower convection velocity are observed.     

湍流参数对复杂形状涡轮盘腔流场的数值研究

以罗罗公司未来某型双轴发动机的高压涡轮旋转盘腔为计算模型，目的是要对复杂的内部旋转空腔流场有一个详细的了解，通过CFD(FLUENT6．0．12)的研究为即将进行的放大比例的旋转涡轮盘腔试验提供指导。计算中采用了Spalart-Allmaras湍流模型，结果表明湍流参数相同的流动，其流体结构是非常相似的；流场主要由湍流参数λT决定涡旋的大小和数量；对于湍流参数小于0．199的旋转流动，盘腔内存在着自由涡流。     

超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

复杂旋转盘轴腔两相流动与传热数值模拟

研究了航空发动机空气系统和润滑系统中复杂的旋转盘腔和旋转轴腔内油气两相流动与传热的数值计算方法,分析了其流动传热特性.以典型小型涡扇发动机的风扇轮盘前腔、风扇轮盘后腔、轴流轮盘前腔、前轴承腔、相关连接气路等组成的多进口、多出口的旋转盘轴腔为对象,研究了用Mixture模型和Eulerian模型计算该系统的速度场、压力场、温度场的方法.结果表明:在相同的计算条件下,两种模型计算的速度场基本一致;两者计算的压力场只在轴流轮盘前腔略有差别,Eulerian模型计算的该腔压力约为Mixture模型计算值的93%;Mixture模型得到温度场较高,由Eulerian模型计算的前、后轴承温度分别约为Mixture模型计算值的93%和94%;Mixture模型计算经济性较好,其迭代一步所需时间约为Eulerian模型的63%.      

Local fracture of the burning solid propellant

The problem of fracture mechanics for the solid propellant deteriorated by a crack-shaped cavity, the surface of which burns is considered. It is assumed that the propellant reagents are uniformly distributed in the solid phase while the combustion products are gaseous. The sufficient condition for the “combustion -fracture” stability regime is obtained.     

三维湍流室内气流数值解

本文基于有限差分法,k-ε模型和p-v法提出了计算三维湍流室内气流问题的数学模型,且用它求解了某水电厂厂房的流场。同时还做了相应的模型试验,并把二者进行了比较,获得了满意的结果。     

热力学效应对液氢非定常空化流动的影响

为研究热力学效应对低温流体非定常空化流动的影响,以液氢为研究对象,采用大涡模拟计算了热力学和等温条件下液氢绕回转体的非定常空化流动。计算结果表明:热力学效应延长了液氢的空化周期,增强了空化的非定常特性,抑制了空化的发展。热力学条件下的空泡团内部包含更细、更小的气泡,呈现出多孔质特征,且使整个流域存在1.5K左右的温度波动。通过旋涡运动分析发现,热力学效应使旋涡结构由空穴交界面向空穴内部移动。基于涡量传输方程分析了热力学和等温条件下空穴和旋涡的交互作用,发现在热力学条件下,旋涡伸长项的作用位置主要位于空穴前端和尾端的交界面处,旋涡扩张项和斜压扭矩项主要位于空穴内部;在等温条件下,旋涡伸长项和旋涡扩张项主要位于空穴上方交界面处,斜压扭矩项主要位于空穴尾端。      

数值模拟有外挂物的空腔流动

数值模拟了有外挂物存在的开式空腔的超音速流动。数值计算采用全隐式无分裂方法,湍流模型采用修正的B-L代数模型。研究了有外挂物存在的空腔流动的剪切层自持性振荡,并与无外挂物存在时的现象进行类比。压强与时间的变化曲线、空腔底部流向的静压系数分布以及空腔后壁面法向的静压系数分布表明,空腔内的压力脉动与外挂物所处位置是有关的。对于Rossiter振荡模态,数值计算结果能够与实验结果相吻合。     

反旋进气盘腔内流动与换热的数值模拟

应用RNGk-ε湍流模型对围屏上带反旋喷嘴的径向内流旋转涡轮盘腔内的流动与换热进行了数值模拟,揭示了盘腔内的压力损失及冷气流量、旋转雷诺数、湍流参数等因素对盘腔内流动换热的影响.计算结果表明:盘腔内的流动结构主要由湍流参数决定;在某一旋转雷诺数下盘腔内压力损失随冷气流量的增大而呈现S型变化;反旋喷嘴的应用能有效地降低盘腔内的压力损失;转盘附近的努赛尔数随冷气流量及旋转雷诺数的增大而增大.      

直接供气预旋转静系流动和换热数值模拟

采用CFD商业软件FLUENT分别从维度模型和湍流模型两个角度对接近真实的高压涡轮旋转盘腔进行数值模拟,研究了一种直接供气预旋转静系旋转盘腔模型内的流动和换热特性.通过计算研究发现:低位腔内的流动仅由中心入流控制,二维计算不能反映高位腔内中心入流与预旋入流的掺混;二维和三维计算所得的腔内压力分布相似,实验值介于两者之间;主盘面的换热主要受中心入流控制,二维和三维计算差别不大,Realizable k-ε湍流模型在定性和定量上均能更好地反映换热实验结果.      

有限平板绕流Monte—Carlo方法仿真

本文采用直接仿真Monte-Carlo方法求解有限长平板绕流等问题。本方法是通过计算机跟踪仿真分子的运动来实现数值模拟的。仿真分子间的碰撞计算由统计抽样确定。碰撞模型分别选用了硬球分子模型和负幂律分子模型。仿真分子与固壁作用采用由完全扩散反射和镜面反射按比例混合组成的模型。为了检验方法的可靠性,还计算了激波结构和Rayleigh问题等一维流动,二维计算采用同步并行程序。数值结果表明直接仿真Monte-Carlo方法能够较好地模拟稀薄气体力学中的一些问题。对于二维计算,所花费的机时和所需的内存均在国内机器所允许的范围内。     

Aerodynamic characteristics of store during lateral jet assisted separation from cavity using free drop technique

Aimed at the problem of store separation from internal cavity, this paper innovatively puts forward a separation scheme of using lateral jet to assist store safe separation. The jet ensures that the store is continuously subjected to down head moment during separation, so as to ensure safe separation. The wind tunnel free drop test technique with lateral jet is established, which can ensure that there is no support interference in the motion process of the store and more truly simulate the motio...     

离心叶轮内稀释固液两相流动的数值模拟

 对固液两相流动进行了数值模拟。液相和固相分别采用 k- ε湍流模型和代数颗粒流( AP)湍流模型并考虑了液固两相间的相互作用。用此模型结合 SIMPLE- C算法对离心叶轮内的稀释固液两相流动进行了计算并研究了颗粒对流场的影响。计算结果同试验结果基本一致     

空腔流激振荡的简化模型分析及振荡频率预估

本文由矩形空腔后缘声振动反馈激励腔口自由剪切层不稳定波的观点出发,把空腔流动振荡问题构造成一个关于空腔内外扰动波的非齐次边值问题,应用Ｆｏｕｒｉｅｒ－Ｌａｐｌａｃｅ变换技术求解,可以预估振荡频率,对Ｌ／Ａ＝Ａ亚声速空腔流动的预估结果与实验符合很好,揭示了此类空腔流动振荡的机理。     

含潜入喷管发动机尾部流场冷流模拟

为了揭示含潜入喷管的固体火箭发动机尾部流动特征,按照几何相近和气动相似原则设计了通道为矩形的二维冷流实验模型,利用相位多普勒粒子分析仪（ＰＤＰＡ）对燃烧室尾部气流的时均速度和湍流脉动速度进行了测量。实验结果表明气流在潜入喷管入口上发生分离,再附点位于喷管前端部外侧,背部空腔内形成一个较为稳定的回流区,流场的轴向和横向湍流强度都比较大。     

有预旋进气转静盘腔中的流动和换热特性数值研究

运用RNG湍流模型对有静盘外缘预旋进气和轴向中心进气、转盘外缘轴向出气和外围屏径向出气的转-静盘腔中流动和传热过程进行了数值研究。在盘腔结构中,静止内隔片和旋转内隔片将盘腔分为内转-静盘腔室和外预旋腔室。由于盘腔流动结构具有周向周期性,取盘腔的1/30作为计算域。研究发现:旋转雷诺数Re_θ、预旋喷嘴进气无量纲流动速度C_w,_p、静盘中心轴向进气的无量纲流动速度C_w,_d和预旋比β_p都不同程度地影响盘腔内流动和换热。预旋和隔片的协同作用可以很好地起到阻止热燃气入侵转-静盘腔和提供低静温冷却气体流入涡轮叶片冷却通道。      

模型三组元喷嘴流量特性实验

三组元喷嘴内混式外喷嘴与内喷嘴组合形成的内混腔，气氢和煤油在此混合形成两相流，通过对大量模型喷嘴进行实验，详细研究了混合腔的结构参数对于喷射进行其内的气体和液体流量特性的影响规律，以及气体，流体流量特性之间的互相影响规律，实验结果表明：气体和液体流量特性互相影响强烈，多数结构参数对于气体和液体的喷入流量系数影响不明显。     

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。