GH2／GO2涡流冷却推力室设计与数值计算

涡流冷却是一种新型液体火箭发动机推力室冷却技术。采用该技术可以简化推力室结构、降低成本,并可提高系统可靠性。对涡流冷却推力室进行了初步设计,并采用PDF非预混燃烧模型和DO辐射模型对所设计的推力室进行了数值仿真。根据计算结果：推力室内部形成了双向涡流;推力室圆筒段壁面温度低于760K;在考虑辐射条件下,推力室圆筒段壁面温度平均升高约140K,最高温度低于900K;涡流冷却技术是可行的,但目前存在燃烧效率相对较低的问题。     

燃烧室辐射换热一维区域法模型研究

目前用于工程计算燃烧室辐射热流的Lefebvre零维辐射模型,不考虑燃气区域对附近壁面的辐射影响,过于简化,不能反映实际燃烧室火焰筒内的辐射换热过程。本文发展了一种机理严密,精度高的用于燃烧室壁温计算的区域法辐射换热模型。文中采用Lefebvre模型和区域法模型计算了一个圆筒型不冷却火焰筒结构的换热;对比分析了两种方法计算的壁面辐射热流和壁面温度分布情况。结果表明,使用Lefebvre模型计算时,由于将燃气按一维划分为定向辐射,在燃烧室前后壁面辐射热流较小,存在一定不足;而使用区域法模型计算时,辐射热流分布更趋合理。     

基于共振三波的大涡模拟壁面模型

大涡模拟的一个重要发展方向是处理工程流动问题,但高雷诺数下近壁湍流边界层对计算网格的过高要求成为制约其应用的主要瓶颈;要减少计算量,构建合理且能够正确反映近壁湍流动力学特性的壁面模型是一条有效途径.在此背景下,构造了一种计算中便于实现的基于共振三波的壁面模型.通过对充分发展槽道湍流算例的研究,完成了对该模型构造方法、方案的论证和验证计算.结果表明使用壁面模型后,在减少计算量的同时,定性和定量的结果均较为合理,能够得到重要的瞬态特征结构和正确的一阶、二阶统计量,从而部分验证了该方法的有效性.     

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。     

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

一种自动化并行重叠网格隐式装配方法

重叠网格装配是处理计算流体力学领域多体相对运动问题的关键技术之一。针对常见重叠网格隐式装配方法中几何分析过程复杂、节点无差别并行查找操作影响并行装配效率等问题,提出一种高度自动化的重叠网格隐式装配方法。首先,基于协方差分析、切割盒子等快速算法,将壁面距离计算与贡献单元存在性判断解耦,实现网格组动态重叠关系的自动化识别;其次,结合集合分析,设计出并行化的自动挖洞算法;最后,通过快速查询方法建立重叠单元与贡献单元的插值关系。针对所实现的并行重叠网格隐式装配工具库,采用某五球体部件验证了自动挖洞逻辑的准确性,并通过机翼-挂架-外挂物（WPFS）模型检验了重叠单元与贡献单元插值关系的准确性。     

本生灯预混火焰淬熄距离实验分析

淬熄距离是表征壁面淬熄效应的一项重要参数.通过直接拍摄火焰位置图片,得到甲烷本生灯预混火焰淬熄距离随来流预混气速度与当量比的变化关系.实验结果表明,本生灯火焰淬熄距离随来流速度的变化规律在不同当量比范围内是不同的.而在同一速度下,随来流当量比增大,淬熄距离贫油时减小,富油时先增大再减小,最后在某一当量比时淬熄距离保持一个定值.这一变化关系是由来流预混气对壁面的冷却作用和火焰面对壁面的加热作用共同影响所造成的,尤其是富油预混火焰的双层结构会极大地影响淬熄效应.      

CFD/CSD-based flutter prediction method for experimental models in a transonic wind tunnel with porous wall

To predict the flutter dynamic pressure of a wind tunnel model before flutter test, an accurate Computational Fluid Dynamics/Computational Structural Dynamics (CFD/CSD)-based flutter prediction method is proposed under the conditions of a 2.4 m × 2.4 m transonic wind tunnel with porous wall. From the CFD simulations of the flows through an inclined hole of this wind tunnel, the Nambu’s linear porous wall model between the flow rate and the differential pressure is extended to the porous wall with inclined holes, so that the porous wall can be conveniently modeled as a boundary condition. According to the flutter testing approach for the current wind tunnel, the steady CFD calculation is conducted to achieve the required inlet Mach number. A time-domain CFD/CSD method is then employed to evaluate the structural response of the experimental model, and the critical flutter point is obtained by increasing the dynamic pressure step by step at a fixed Mach number. The present method is applied to the flutter calculations for a vertical tail model and an aircraft model tested in the current transonic wind tunnel. For both models, the computed flutter characteristics agree well with the experimental results.     

短肢剪力墙结构体系的动力特性分析

应用三维有限元模型,通过对一幢12层短肢剪力墙结构的动力分析,研究并比较了短肢剪力墙结构体系与普通剪力墙结构体系的动力特性。结果表明,与普通剪力墙结构体系相比,短肢剪力墙结构体系抗侧刚度较小,基本周期较长,加速度反应和层间剪力都比较小,这对结构抗震是有利的。此外,本文还将短肢剪力墙在水平载荷作用下的侧移曲线反弯点的相对高度作为此类结构整体性和墙肢相对强弱的综合评价指标,设计中应使反弯点的相对高度不小于0.4。     

一种高速风洞三维洞壁干扰壁压信息修正法

给出以壁压信息法为基础的高速风洞三维模型试验洞壁干扰修正方法。该方法要求洞壁附近为亚音速流动,但允许模型附近出现超音速区及激波。方法可用于各种通气壁或实壁风洞的亚跨音速试验。最后对四个模型在三种试验段中的二十多种试验进行了修正,其结果和NASA非线性修正方法的结果吻合得很好。