燃料分级对旋流火焰的动态特性影响研究

为深入研究分级旋流火焰特性,以分级旋流模型燃烧室为研究对象,对四个不同燃料分级比(Rf)条件下的分级旋流火焰进行了数值研究,在时均燃烧场特性分析的基础上进一步对燃料分级比为1和3两个工况进行了基于壁面建模的大涡模拟(WMLES)研究。结果表明:燃料分级比的改变会影响中心回流区(CRZ)的长度和宽度。燃烧室中截面的散点分布图能够显示出不同燃料分级比条件下的燃烧特征。燃料分级比为1时,燃烧室剪切层仅存在零散的涡破碎区;而燃料分级比为3时,伴随涡破碎区还出现了单螺旋分支进动涡核(PVC)。通过FFT变换获得的燃烧室内剪切层速度能谱主频与进动涡核的旋转频率相同,表明内剪切层速度脉动的产生与进动涡核有关。另外进动涡核会使流场内的燃料分布和燃烧模式发生周期性的变化,进而影响燃烧过程。调整燃料分级比在1附近,能够使分级火焰达到稳定燃烧降低排放的目标。     

主流攻角对双射流孔气膜冷却特性的影响

应用压力敏感漆技术,在平板上测量了不同主流攻角(i=-30°,-20°,-10°,0°,10°,20°,30°)下双射流孔的气膜冷却效率,并利用计算流体动力学(CFD)计算得到的流场对气膜冷却效率的规律进行了解析。所研究的双射流孔结构的孔间无量纲横向距离为0.5,孔间无量纲流向距离为3;射流与主流密度比为1.0,吹风比分别为0.5、1.0、1.5、2.0。结果表明小的主流攻角(i=-10°,0°,10°)下,流场中存在反肾型涡对或挤压作用,气膜层与壁面贴附良好,气膜冷却效率最高;大正值攻角(i=20°,30°)下,虽然气膜覆盖面积大,但反肾型涡对退化,气膜冷却效率下降;大负值攻角(i=-20°,-30°)下,流场中有肾型涡对,且气膜横向覆盖受限,气膜冷却效率最低。     

转子系统挤压油膜阻尼器设计与动力特性

鉴于挤压油膜阻尼器(SFD)设计必须要同时考虑转子系统的动力学特性,提出了一种转子系统与挤压油膜阻尼器耦合设计方法,给出了详细的设计流程。对所设计的阻尼器进行了CFD数值模拟、油膜压力测量以及减振效果实验,结果证明所提出的设计方法是有效的。相比未采用阻尼器,采用阻尼器后转子系统两个转盘的振幅分别下降了46%和39%。通过实验还研究了不平衡量、支承刚度、供油压力和滑油温度对挤压油膜阻尼器减振效果的影响。结果表明,相比于不平衡量和支承刚度对减振效果的影响,供油压力和滑油温度的影响并不显著。进行挤压油膜阻尼器设计时,重点应该关注转子上的不平衡量大小和支承刚度。     

Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian (ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces. In addition, the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method. Based on this method, fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained. Numerical results such as maximum Root Mean Square (RMS) drag, general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results. This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.     

考虑滑流影响的埋入式进气道气动特性分析与流动控制研究

研究滑流对埋入式进气道性能的影响机理对于发动机的稳定性具有重要意义。针对某型涡桨发动机 滑油散热器进排气道系统,基于 CFD技术建立系统的螺旋桨滑流与滑油散热器内流一体化数值模拟方法,进 行考虑滑流影响的滑油散热系统埋入式进气道气动设计仿真分析;针对埋入式进气道由于吸入边界层低能气 流导致通过散热器流量低的问题,设计4组涡流发生器进行主动流动控制。结果表明:滑流使得进气道中的气 流偏离,难以形成稳定漩涡,容易发生流动分离;而合理设计的涡流发生器可以有效改善埋入式进气道内流场 特性,并使通过散热器的流量提高12%。     

Numerical Simulation of Rotor Flow Field Based on Overset Grids and Several Spatial and Temporal Discretization Schemes

A numerical method based on solutions of Euler/Navier-Stokes (N-S) equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method (AUSM) scheme, upwind AUSM/van Leer scheme, AUSM+ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points (IGBPs). Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.     

RELATIONSHIPS FOR MOTION DESCRIBED IN ANY TWO ROTATIONAL REFERENCE FRAMES

In this paper the definitions of the relative derivatives of a vector with respect to time in any reference frame by using the notation of vectrix are proposed. It is further shown that with the help of vectrix one can treat the relative derivatives more clearly and more rigorously when multiple reference frames are concerned. Very useful relationships between the first order relative derivatives described in any two reference frames are derived. Relationships for the second order relative derivatives are also presented in detail. A new concept called "rotation induced acceleration" which is caused by the rotation of two reference frames when motion is necessarily described in them is presented. The result obtained in this paper is the generalization of the acceleration theorem given in conventional theoretical mechanics.     

Structure and kinetic properties of slow-mode shocks associated with magnetic reconnection in the near-Earth magnetotail

We study the structure and kinetic properties of slow-mode shocks near the plasma sheet boundary layer (PSBL) associated with magnetic reconnection by Cluster observation. The presence of slow-mode shocks is confirmed by traditional Rankine–Hugoniot (RH) analysis and Monte-Carlo shock fitting method. The Walén analysis, applied to the tailward flow associated with slow-mode shocks, also supports that plasma was accelerated across a Petschek-type slow-mode shock connected to the diffusion region. Back-streaming ions were observed on the shock layer, and cold ions were accelerated and heated by slow-mode shocks. In addition, whistler and electrostatic solitary waves were observed around the slow-mode shocks. These waves might be excited by the observed field-aligned electron beams near the shocks.     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

二维三角晶格的Peierls相变前后的声子色散曲面的研究

利用晶格动力学理论计算了二维三角晶格的Peierls相变前后的声子能谱。并且在第一布里渊区绘制了声子色散曲面。通过对Peierls相变前后的声子能谱比较可知，随其对称性降低，声子能谱的简并度被破缺。二维三角晶格的Peierls相变前只有两支声频支声子能谱，Peierls相变后有两支声频支声子能谱与两支光频支声子能谱，并且声频支声子能谱与光频支声子能谱不相交，两支光频支声子沿ky轴平缓地变化。