一种新的动网格方法及其应用

 介绍了一种新型的动网格方法。该方法以Delauney图为基础进行整个网格的移动变形：首先是生成满足Delauney特性的Delauney图;然后把对网格移动的特性映射到Delauney图的移动上;最后把网格点在Delauney图以映射的方式进行重新定位而得到新的网格点坐标,从而得到移动后的网格。该方法不需要迭代,因而效率高,并且能适用于任意拓扑结构的网格。以翼身组合体的机翼颤振、翼身融合体的变体过程的非定常流动等为例来说明该方法的高效性及适用面广等特点。     

基于Reynolds平均的高超声速二方程湍流模型

 在Reynolds平均的框架下推导了可压缩湍流Reynolds应力方程和湍动能方程。根据一定的假设和尺度分析简化并封闭了所推导的湍动能方程。在均匀湍流假设下,湍动能耗散率可分解成可压缩性耗散和旋度耗散,并对其中的可压缩性耗散进行了封闭;同时认为旋度耗散不受可压缩性影响,直接引用传统的Reynolds平均不可压缩湍动能耗散率模型方程。由此构造了适用于高马赫数的二方程可压缩湍流模型。应用所发展的模型计算了高超声速平板绕流,并与若干现有模型的计算、实验与半经验公式的计算结果进行了对比,验证了所发展的模型。在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对所发展的湍流模型,以及若干现有模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了湍流模型可压缩修正的必要性,得到了所发展模型的计算结果最接近实验结果的结论。     

气流速度振荡场中幂律液膜不稳定性分析

采用线性稳定性分析研究了处于气流速度振荡场中幂律液膜时间模式的不稳定性。振荡的气流速度导致动量方程为含有时间周期系数的希尔方程，采用Floquet理论进行求解。详细研究了不同振荡幅值和振荡频率下表观雷诺数、幂律指数及无量纲速度因子对各不稳定区间的影响。结果表明：振荡幅值的增加或振荡频率的减小会使液膜不稳定区域的个数增加，且Kelvin-Helmholtz（K-H）不稳定区域的最大增长率、主导波数和截止波数随振荡幅值和振荡频率的增加而增加；表观雷诺数、幂律指数和无量纲速度因子的增加增强了K-H不稳定区域内的不稳定性，使参数不稳定区域内的增长率先减小后增加；振荡幅值的变化不改变最大增长率发生转折时对应的流变参数，而当振荡频率较小时，幂律指数和无量纲速度因子的增加却使最大增长率单调增加。     

Influence of non-equilibrium reactions on the optimization of aerothrust aeroassisted maneuver with orbital change

The spaceplane is perspective vehicle due to wide maneuverability in comparison with a space capsule. Its maneuverability is expressed by the larger flight range and also by a possibility to rotate orbital inclination in the atmosphere by the aerodynamic and thrust forces. Orbital plane atmospheric rotation maneuvers can significantly reduce fuel costs compared to rocket-dynamic non-coplanar maneuver. However, this maneuver occurs at Mach numbers about 25, and such velocities lead to non-equilibrium chemical reactions in the shock wave. Such reactions change a physicochemical air property, and it affects aerodynamic coefficients. This paper investigates the influence of non-equilibrium reactions on the aerothrust aeroassisted maneuver with orbital change. The approach is to solve an optimization problem using the differential evolution algorithm with a temperature limitation. The spaceplane aerodynamic coefficients are determined by the numerical solution of the Reynolds-averaged Navier-Stokes equations. The aerodynamic calculations are conducted for the cases of perfect and non-equilibrium gases. A comparison of optimal trajectories, control laws, and fuel costs is made between models of perfect and non-equilibrium gases. The effect of a chemically reacting gas on the finite parameters is also evaluated using control laws obtained for a perfect gas.     

Nacelle-airframe integration design method for blended-wing-body transport with podded engines

Strong shock waves and flow separation often occur during the integration of nacelle and airframe for blended-wing-bodies with podded engines. To address this problem, this paper presents an integration method with numerical simulations. The philosophy of channeling flow and avoiding the throat effect on the nacelle and airframe is established based on the analysis of flow interference in the initial configuration. A parametric integration design method is proposed from twodimensional plane to three-dimensional space with control mechanisms and selection principles of the key parameters determined by their influences. Results show that strong shock waves and flow separation can be successfully eliminated under the influence of both the reshaped channel and decelerated inflow below the nacelle. Supersonic regions around the nacelle are effectively reduced, concentrating mainly on the lip position. Thus, a significant cruise drag reduction(8.7%) is achieved though the pressure drag of the nacelle increases.     

Recent advances in physical understanding and quantitative prediction of impinging-jet dynamics and atomization

Impinging-jet injectors are widely used in liquid propulsion applications, since their simple configuration provides reliable and efficient atomization. The flowfield involves a series of complicated spatio-temporal evolutions. Much effort has been directed toward understanding the underlying physics and developing quantitative predictions of impinging-jet atomization. This paper summarizes the recent advances in this direction, including state-of-the-art theoretical, experimental, and numerical studies, along with representative results. Finally, concluding remarks address remaining challenges and highlight modeling capabilities of high-fidelity simulations.     

升力偏置对共轴刚性旋翼前飞气动特性的影响

为了分析升力偏置对共轴刚性旋翼前飞气动特性的影响，建立了基于雷诺平均Navier-Stokes方程的计算流体力学方法进行共轴旋翼流场求解，采用嵌套网格方法模拟桨叶运动，采用双时间方法进行时间推进。针对不同升力偏置状态，采用基于"差量法"的共轴旋翼高效配平策略进行操纵量配平。通过对Harrington-1旋翼性能的计算，验证了方法的有效性。对比计算了共轴刚性旋翼在不同前进比和升力偏置量下的气动性能和流场特征，结果表明：双旋翼操纵量在小前进比状态有明显差别，在大前进比状态基本一致；在相同拉力状态，随着升力偏置量的增大，共轴旋翼升阻比先升高后降低，其阻力却不断增大，不同前进比状态的最大升阻比对应的升力偏置量不同；双旋翼相遇时桨叶拉力出现脉冲式波动，由于流场被前行桨叶所主导，因此后行桨叶拉力波动幅值更大，且波动幅值随升力偏置量的增加而增大。     

近壁区流向涡生成的数值模拟

对于在湍流边界层近壁区,作为相干结构的主要特征的流向涡的产生及发展的全过程进行了研究.采用准二维的共振三波作为湍流边界层近壁区相干结构初值,用直接数值模拟方法模拟了准二维波发展到明显的三维扰动以及其中流向涡生成的整个过程,经分析发现,在该过程中,由于非线性作用,相干结构出现了明显的三维性,并产生了马蹄涡,与实验结果相一致,而且也研究了压力梯度对于流向涡生成以及发展的影响, 在逆压梯度下,流向涡的生成得更早、更快,相干结构幅值增长较快,结构变得更加复杂.这些都说明了逆压梯度对相干结构具有激励作用.     

Delaying stall of morphing wing by periodic trailing-edge deflection

Morphing wings can improve aircraft performance during different flight phases. Recently research has focused on steady aerodynamic characteristics of the morphing wing with a flexible trailing-edge, and the unsteady aerodynamic and stall characteristics in the deflection process of the morphing wing are worthy further investigation. The effects of the angle of attack and deflection rate on aerodynamic characteristics were examined, and based on the aerodynamic characteristics of the morphing wing, a method was developed to delay stall by using the flexible periodic trailing-edge deflection. The numerical results show that the lift coefficients in the deflection process are smaller than those in the static situation at small angles of attack, and that the higher the deflection rate is, the smaller the lift coefficients will be. On the contrary, at large angles of attack, the lift coefficients are higher than those in the static case, and they become larger with the increase of the deflection rate. Further, the periodic deflection of the flexible trailing-edge with a small deflection amplitude and high deflection rate can increase lift coefficients at the critical stall angle.     

Experimental Study of Corner Stall in a Linear Compressor Cascade

In order to gain a better knowledge of the mechanisms and to calibrate computational fluid dynamics (CFD) tools including both Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES),a detailed and accurate experimental study of corner stall in a linear compressor cascade has been carried out.Data are taken at a Reynolds number of 382 000 based on blade chord and inlet velocity.At first,inlet flow boundary layer is surveyed using hot-wire anemometry.Then in order to investigate the effects of incidence,measurements are acquired at five incidences,including static pressures on both blade and endwall surfaces measured by pressure taps and the total pressure losses of outlet flow measured by a five-hole pressure probe.The maximum losses as well as the extent of losses of the corner stall are presented as a function of the investigated incidences.