基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

多段翼型大迎角分离流动的DELAYED RANS/LES混合算法

多段翼型的大迎角绕流发生大范围附面层分离,具有明显的三维与非定常流动特性。RANS/LES混合算法继承了LES对流动分离区大尺度漩涡准确模拟的优点,避免了纯LES算法需求网格量巨大与亚格子模型壁面函数不成熟等问题,对分离流动的模拟效果优于RANS算法。以S-A湍流模型与Smagorinsky亚格子模型为基础,借鉴DDES的附面层延迟控制思想,构造了可用于对接网格、重叠网格的DELAYED RANS/LES混合算法。研究了GA（W）-1多段翼型的大迎角分离流动及其气动特性。     

"回收/调节"方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离...     

DNS/LES方法在剪切湍流模拟中的应用

提出采用DNS/LES混合方法进行自由剪切类混合流动湍流研究,该方法兼顾了直接数值模拟(DNS)方法在转捩线性范围内的高精确度模拟和大涡模拟(LES)方法在转捩非线性及全湍流区的高效模拟两方面优点,通过对低对流Mach数流场计算表明该方法可以较好地模拟流场的转捩过程.同时,文中还对比实验分析了湍流场对光学传输效应的影响,给出了与实验一致的结论.     

一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。     

基于当地流动结构的RANS/LES混合模型

针对传统雷诺平均Navier-Stokes/大涡模拟(RANS/LES)混合模型固有的RANS和LES区域衔接难题,借鉴Vreman提出的亚格子(SGS)模型,提取其中描述流动状态的判断函数,利用不可压缩槽道湍流标定模型系数,构造了一个新的RANS/LES混合模型。评估算例包括不可压缩槽道湍流、稳态超声速平板边界层、NACA4412翼型绕流和亚临界雷诺数圆柱绕流。该模型作为壁面模化LES,解决了传统RANS/LES方法中速度型偏离对数律(LLM)的缺陷,且拓宽了原Vreman模型在粗糙网格上的适用性,具有较高精度。该模型同时解决了模化应力折损问题,且对大尺度分离流动的模拟精度相较传统的RANS/LES混合模型如脱体涡模拟(DES)有进一步提升。     

复杂湍流流动的混合RANS/LES方法研究

混合RANS/LES方法是近年来复杂湍流模拟的重要方法。简要回顾了混合RANS/LES模拟方法的发展历程,着重总结和分析了分离涡混合模型、类Menter SST混合模型以及湍流能量谱一致混合方法三类混合模型的构造方法和发展历程,对这三种方法的应用以及优缺点进行了简要评述,最后指出了混合RANS/LES方法发展应该兼顾的问题,为下一步的混合RANS/LES模拟提供了参考。     

后向台阶超声速流动的动力混合RANS／LES模拟

针对可压缩湍流流动,构建了1种动力混合RANS/LES模式。运用该模式模拟了超声速后向台阶流动,结果再现了二维后向台阶剪切层的拟序结构,计算得到了统计平均的压力分布,回流区长度与实验结果符合较好。与定系数的一般混合RNAS/LES模拟结果的对比表明,动力混合RANS/LES模式在对流场拟序结构的捕捉、流动时均压力分布方面表现较优。     

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架，设计开发了LES(Large Eddy Simulation)湍流模型，主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法，给出设计的软件框架和计算流程，尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例，验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型，具备高精度数值格式、亚格子模型、湍流统计等通用模块，可为国内学者提供一个LES湍流模拟研究的开放平台。     

一种新型的动态RANS/LES混合方法

为解决RANS/LES混合方法过渡依赖网格疏密的问题,使用熵函数fs作为交界面函数构造一种新型的动态RANS/LES混合方法。同时采用脱体涡模拟方法 (DES97),对平板边界层流动、超声速凹腔-压缩拐角流动以及超声速底部流动进行数值模拟,并与实验数据进行对比。结果表明,无论对于附体流动区域还是分离区域,fs都能够动态地指示出边界层的范围。构造的动态RANS/LES混合方法能够显著缓解DES97由于网格因素导致的雷诺应力模化不足的现象,同时保证分离区域涡结构的求解精度。另外,熵函数fs是从能量耗散的角度进行构造,并不依赖于特定的湍流变量,因此更具有普适性。     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

μ理论在电液负载模拟器中的应用

 研究μ理论在电液负载模拟器中的应用，提出电液负载模拟器的鲁棒控制策略。电液负载模拟器是一个复杂的机-电-液复合系统，且是一个强耦合、时变受控对象。综合考虑参数变化、模型变动和干扰等不确定性的影响，利用μ综合控制理论设计电液负载模拟器的鲁棒力控制器，并通过μ分析使用鲁棒力控制器时系统的鲁棒稳定性和鲁棒性能。给出使用鲁棒力控制器和经典力控制器时的实验结果，实验结果表明所设计鲁棒力控制器的有效性和优越性。     

有翼高超声速再入飞行器气动设计难点问题

有翼高超声速再入飞行器是近年来的研究热点,气动设计是飞行器设计的关键。为了更清楚地认识有翼高超声速再入飞行器气动设计的难点问题,对有翼高超声速再入飞行器的发展、优势及总体任务剖面进行了介绍,从5个方面详细介绍了该类飞行器气动设计的难点问题,包括多约束复杂面对称气动布局设计、高温真实气体效应对气动特性影响、天地差异与天地换算方法、反作用控制系统(RCS)喷流干扰对气动特性的影响以及气动数据不确定度等,简要阐明了这些难点问题对总体设计的重要性以及初步的解决思路,为有翼高超声速再入飞行器气动设计提供了一些参考。     

基于圆型限制性三体问题模型的行星卫星逃逸能量研究

 基于由飞行器、行星及其卫星组成圆型限制性三体问题模型，通过庞加莱映射的方法，研究了飞行器从行星卫星附近逃逸的问题。在Jacobi常数确定的前提下，通过逆向积分，飞行器从L1或L2点附近返回近月点，得到近月点速度出发速度。研究结果表明绕飞L1点和L2点逃逸行星卫星需要的最低能量是不同的，从月球表面逃逸所需的速度脉冲分别比开普勒算法节省46.5m/s和42.3m/s，且均小于Villac等人在Hill模型下得到38.9m/s，从而改进了Villac等人的相关工作，同时也给出了从太阳系主要行星卫星表面逃逸所需的最小能量。     

Model Development and Adaptive Imbalance Vibration Control of Magnetic Suspended System

对星载磁悬浮单框架控制力矩陀螺,建立其转子的平动与转动的动力学方程。所建模型反映了转子的动、静不平衡特性,以及转子运动与框架、星体运动的耦合关系。模型是主动磁悬浮控制与动框架效应分析的基础。对动、静不平衡量未知的转子,设计了自适应对中控制器。在分散PD、比例交叉反馈控制器与高低通滤波器的基础上,通过动、静不平衡量的辨识与补偿控制,使转子绕其惯量主轴旋转,衰减角动量交换执行机构高速转子因不平衡而对基座产生的干扰力,消除星上主要颤振源,达到降低卫星姿态抖动的目的。     

Love Waves in Layered Graded Composite Structures with Imperfectly Bonded Interface

Theoretical analysis and numerical calculations of Love wave propagation in layered graded composites with imperfectly bonded interface are described in this paper. On the basis of WKB method, the approximate analytic solutions for Love waves are obtained. By the interface shear spring model, the dispersion relations for Love waves in layered graded composite structures with rigid, slip, and imperfectly bonded interfaces are given, and the effects of the interface conditions on the phase velocities of Love waves in SiC/Al lay- ered graded composites are discussed. Numerical analysis shows that the phase velocity decreases when the defined flexibility parameter is greater. For the general imperfectly bonded interface, the phase velocity changes in the range of the velocities for the rigid and slip interface conditions.     

Experimental Investigation on Performance of Pulse Detonation Rocket Engine Model

The PDRE test model used in these experiments utilized kerosene as the fuel, oxygen as oxidizer, and nitrogen as purge gas. The solenoid valves were employed to control intermittent supplies of kerosene, oxygen and purge gas. PDRE test model was 50 mm in inner diameter by 1.2 m long. The DDT （deflagration to detonation transition） enhancement device Shchelkin spiral was used in the test model. The effects of detonation frequency on its time-averaged thrust and specific impulse were experimentally investigated. The obtained results showes that the time-averaged thrust of PDRE test model was approximately proportional to the detonation frequency. For the detonation frequency 20 Hz, the time-averaged thrust was around 107 N, and the specific impulse was around 125 s. The nozzle experiments were conducted using PDRE test model with three traditional nozzles. The experimental results obtained demonstrated that all of those nozzles could augment the thrust and specific impulse. Among those three nozzles, the convergent nozzle had the largest increased augmentation, which was approximately 18%, under the specific condition of the experiment.     

Compression Creep Behavior of High Volume Fraction of SiC Particles Reinforced Al Composite Fabricated by Pressureless Infiltration

The compression creep deformation of the high volume fraction of SiC particles reinforced Al-Mg-Si composite fabricated by pressure-less infiltration was investigated. The experimental results show that the creep stress exponents are very high at temperatures of 673 K, 723 K and 773 K, and if taking the threshold stress into account, the true stress exponent of minimum creep strain rate is still approximately 5, although the volume fraction of reinforcements is very high. The creep strain rate in the high volume fraction rein- forced aluminum alloy matrix composites is controlled by matrix lattice diffusion. It is found that the creep-strengthening effect of high volume fraction of silicon carbide particles is significant, although the particles do not form effective obstacles to dislocation motion.     

Simulation and Analysis of Crashworthiness of Fuel Tank for Helicopters

Crashworthiness requirement of fuel tanks is one of the important requirements in helicopter designs. The relations among the protection frame, textile layer and rubber layer of the fuel tank are introduced. Two appropriate FE models are established, one is for an uncovered helicopter fuel tank without protection frame, and the other is for fuel tank with protection frame. The dynamic responses of the two types of fuel tanks impinging on the ground with velocities of 17.3 m/s are numerically simulated for the purpose of analyzing energy-absorbing capabilities of the textile layer and protection frame. The feasibility of the current crashworthiness design of the fuel tank is examined though comparing the dynamic response behaviors of the two fuel tanks.     

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。