展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

流场可压缩性对涡相互作用影响的数值研究

涡相互作用作为冲压发动机喷注装置的典型抽象流动现象,研究其可压缩性影响对于认识包含化学反应的真实燃料喷射场具有一定基础理论价值。基于经过数值验证的可压缩Navier-Stokes算法与压力泊松方程初始条件设置方法研究流场可压缩性对涡相互作用演化过程的影响。结果表明,以较高涡旋马赫数表征的可压缩性在改变涡对形态的同时,具有延缓涡旋相互靠近,迟滞融合进程的作用。关于涡对系统开始融合的临界条件,可压缩相互作用开始的临界展弦比与无量纲时间相比于低速涡对明显提高。为在无量纲时间意义下统一不同马赫数涡对相互作用的进程,从涡心密度随时间变化规律出发在不可压涡对特征时间的基础上,初步构建了考虑可压缩性的时间尺度修正关系。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

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.     

Advances in critical technologies for hypersonic and high-enthalpy wind tunnel

Hypersonic and high-enthalpy wind tunnels and their measurement techniques are the cornerstone of the hypersonic flight era that is a dream for human beings to fly faster, higher and further. The great progress has been achieved during the recent years and their critical technologies are still in an urgent need for further development. There are at least four kinds of hypersonic and high-enthalpy wind tunnels that are widely applied over the world and can be classified according to their operation modes. These wind tunnels are named as air-directly-heated hypersonic wind tunnel, light-gas-heated shock tunnel, free-piston-driven shock tunnel and detonation-driven shock tunnel, respectively. The critical technologies for developing the wind tunnels are introduced in this paper, and their merits and weakness are discussed based on wind tunnel performance evaluation. Measurement techniques especially developed for high-enthalpy flows are a part of the hypersonic wind tunnel technology because the flow is a chemically reacting gas motion and its diagnosis needs specially designed instruments. Three kinds of the measurement techniques considered to be of primary importance are introduced here, including the heat flux sensor, the aerodynamic balance, and optical diagnosis techniques. The techniques are developed usually for conventional wind tunnels, but further improved for hypersonic and high-enthalpy tunnels. The hypersonic ground test facilities have provided us with most of valuable experimental data on high-enthalpy flows and will play a more important role in hypersonic research area in the future. Therefore, several prospects for developing hypersonic and high-enthalpy wind tunnels are presented from our point of view.     

基于多Agent的机群编队导航系统设计与研究

智能化建模在装备作战仿真领域中独居优势.将多Agent技术程序设计方法引入机群编队导航系统的设计中,对模型进行了顶层设计并分析了Agent子系统之间的信息交互过程,为机群编队导航系统的研发和设计提供了新的思路.     

水气动量交换的实验研究

在风浪槽中进行水气动量交换实验，测量了不同风速和风区下的波浪及空气湍流流场。测量数据表明水面上风速分布符合对数律，水气界面附近在小范围的动量交换常通量层。阻力系数ＣＤ随风速增加线性增长，随波龄增加而减小。这一趋向与理论分析及外海帝测结果一致。对风浪发展过程中波浪和空气湍流谱结构的分析证明了水气运动的耦合关系。     

Global aerodynamic design optimization based on data dimensionality reduction

In aerodynamic optimization, global optimization methods such as genetic algorithms are preferred in many cases because of their advantage on reaching global optimum. However, for complex problems in which large number of design variables are needed, the computational cost becomes prohibitive, and thus original global optimization strategies are required. To address this need, data dimensionality reduction method is combined with global optimization methods, thus forming a new global optimization system, aiming to improve the efficiency of conventional global optimization. The new optimization system involves applying Proper Orthogonal Decomposition (POD) in dimensionality reduction of design space while maintaining the generality of original design space. Besides, an acceleration approach for samples calculation in surrogate modeling is applied to reduce the computational time while providing sufficient accuracy. The optimizations of a transonic airfoil RAE2822 and the transonic wing ONERA M6 are performed to demonstrate the effectiveness of the proposed new optimization system. In both cases, we manage to reduce the number of design variables from 20 to 10 and from 42 to 20 respectively. The new design optimization system converges faster and it takes 1/3 of the total time of traditional optimization to converge to a better design, thus significantly reducing the overall optimization time and improving the efficiency of conventional global design optimization method.     

中国核电厂操作员反应时认知模型的应用(英文)

讨论了一些有关人员行为与时间有依赖关系的可靠性模型；给出了中国核电厂操作员在清华大学培训模型机上的一些可靠性实验结果。此外，还根据已获得的实验数据计算了在ＡＴＷＳ事故条件下人的失误概率。