带反馈控制的基本带翼外形的耦合气动/运动数值模拟研究

介绍了一种先进的、时间精度的多学科数值模拟平台的开发和应用,该平台可用于预测飞行器机动过程中的非定常气动特性和飞行力学特性.通过耦合求解非定常Navier-Stokes方程、刚体动力学方程及控制律,对基本带翼外形进行了虚拟飞行数值模拟研究.求解器采用基于结构网格的有限体积法,时间方法采用双时间步方法.使用反馈控制器获得舵偏控制指令,重叠网格方法用来模拟物体间的相对运动.模拟结果显示控制器的增益必须正确选择才能获得正确的响应.仿真结果表明本文发展的虚拟飞行数值模拟方法能够预测包含动边界的复杂非定常流动,具有较强的工程适用性.     

基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

基于流-固耦合传热的热气防冰系统干空气飞行蒙皮温度场计算研究

以机翼热气防冰系统为研究对象,建立了包含热气防冰系统防冰腔内外流场对流换热和固体结构导热的三维稳态流-固耦合传热物理模型,对整个计算区域生成混合网格,边界条件为第三类边界条件,采用计算流体力学方法以 FLUENT 软件为工具,对干空气飞行状态下流-固耦合传热模型进行了求解,获得防冰腔蒙皮内外表面对流换热系数分布和温度场结果,并对计算结果进行了分析。结果表明：防冰腔铝合金蒙皮沿展向和厚度方向导热显著,温度分布较均匀,防冰引气温度为200℃时,防冰腔蒙皮内外表面上最高温度为101℃,最低温度为21℃,3 mm厚的蒙皮同一点处内外表面最大温差仅为4℃,防冰腔排气口处气体的平均温度为63℃。热气防冰系统蒙皮温度场计算方法和计算结果,能够为热气防冰系统干空气飞行试验设计和测试中温度传感器的选型与布置提供依据。     

煤油两相连续旋转爆震燃烧室工作特性试验研究

为了深入分析煤油燃料两相连续旋转爆震燃烧室的工作特性,采用富氧空气或氧气为氧化剂,通过试验得到爆震波的时域、频域特征,对两相连续旋转爆震燃烧室中爆震波的起爆过程和稳定后的传播过程进行研究。利用基于激光散射相位多普勒分析(PDA)技术对雾化流场进行了测量,得到喷注器出口不同平面处煤油液滴速度与直径的统计分布。试验结果表明当煤油流量为78g/s,氧气流量为224.0g/s,空气流量为72.5g/s,当量比为1.083时,燃烧室在单波模态下工作,爆震波传播频率为0.904k Hz,平均转速为649m/s。使用氧气作为氧化剂,当煤油流量为81.8g/s,氧气流量为231.8g/s,当量比为1.222时,燃烧室在双波模态下工作,爆震波传播频率为5.882k Hz,平均转速为1848m/s,传播过程中表现出很强的非定常性。在当量比为0.805~0.908的富氧工况下,随着氧化剂中含氧量的增加,爆震波的速度逐步增大,最终达到2440m/s;在当量比为1.057~1.220的富燃工况下,随含氧量的增加爆震波速度呈现线性增长的特征。     

α，β相分布形态对TC11伏合金高温循环蠕变速率的影响

对TC11钛合金施行了不同的热处理制度,得到了3种组织形态,研究了不同的组织形态的高温循环蠕变,结果表明,晶间β组织具有低的循环蠕变速率,而随最大应力的增加,蠕变速率增加快,而等轴组织 片状组织循环蠕变速率较高,但随应增加,蠕变速率增加较慢,分析认为,晶界滑移行为对高温循环蠕变有重要影响。     

基于气动外形优化和主动流动控制的减阻技术

增升和减阻是民用飞机设计的永恒目标。本文对民用飞机减阻技术最新进展及未来发展进行综述,在此基础上,报道了计算所团队通过数值优化设计及主动流动控制实现减阻的研究进展。为此,发展了基于壁面模化大涡模拟(WMLES)准确模拟湍流边界层的先进数值方法,以及利用伴随方程方法实现气动外形快速数值优化的技术。在具体实践中,通过组合基于自由变形(FFD)的几何参数化模块、基于线弹性体法的网格变形模块,以及基于伴随方程方法和序列二次规划(SQP)的优化算法模块,搭建了面向工程、适于超大规模变量优化设计的整机气动优化设计平台,在NASACRM模型上实现了设计变量超过600的大变量、跨声速气动数值优化设计,在合理的工程约束条件下,有效削弱了机翼表面激波,总体减阻率达到2%以上;通过对湍流平板边界层外层结构进行射流非定常控制,在中等雷诺数Reτ=4700条件下,实现了5%～6%的当地减阻率。研究证实,数值优化设计和基于外层射流控制的主动流动控制技术将是大型民用飞机减阻中非常有前景和值得重视的两种方法。     

用于级间分离研究的TBCC动力TSTO气动布局概念设计

基于涡轮/冲压组合动力的水平起降两级入轨飞行模式是重复使用飞行器降低发射成本、缩短发射周期的重要途径之一,气动布局设计需要在飞行器总体规模尺度约束下满足全速域全空域的气动特性、操稳特性及防热特性需求。为研究级间分离特性,讨论了任务使命、动力配置、飞行模式及总体规模限制下的两级入轨(TSTO)重复使用飞行器的气动布局设计,针对二级及一级面临的飞行任务需求和气动特性需求,分别提出了多种气动布局方案。对两种改进方案进行了初步的气动计算,并进行了上升段升重平衡下的飞行剖面重构。为提高级间分离的安全性并提高超声速/高超声速的升阻效率和航向稳定性,对双垂尾及可下折翼梢进行了适当修改,形成了TSTO系统新一轮研究方案,并在此基础上规划了后续研究工作。     

考虑空腔的高超声速多流动区域同步数值模拟

先进的高超声速飞行器具有薄壁空腔结构,在飞行过程中受热会产生空腔内气体流动现象,从而影响流场和结构的温度分布。采用数值方法准确模拟高超声速流场、结构温度场和空腔内流动对热结构分析是很有必要的。以研究空腔流动对结构温度分布影响为目的,发展了一种适用于多流动区域流场/结构温度场耦合问题的同步计算方法,并以高超声速带空腔结构物体为例,数值研究了其外部气动热/结构热传导引起的空腔热对流问题。以已发展的高超声速外流场/结构温度场同步计算方法为基础,为了进一步考虑空腔内低速流场,采用了预处理矩阵方法。在流场与结构温度场的交界面两侧分别引入虚拟单元,从而高效地实现相邻场之间物理信息交换。首先通过标准算例验证了方法在求解单独气动热/结构热传导问题以及空腔自然对流问题中的准确性。进而对封闭和带有开孔的两种高超声速运动圆环分别进行多流动区域同步数值模拟。计算结果表明,由于结构温度不均匀引起的空腔内热对流反之也会对结构温度场分布产生轻微的影响。在空腔内气体流动的影响下,封闭圆环的前缘温度在35s内最多下降0.8%左右。对于带开孔空腔的圆环,其孔壁周边温度在0.5s内能够超过外流前缘驻点温度。     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。