液体火箭发动机尾焰对发射平台冲击效应

为研究液体火箭发动机尾焰对发射平台的冲击效应特性,建立液体火箭发动机尾焰对发射平台冲击数值计算模型.针对液氧/煤油发动机尾焰对发射平台冲击特性,基于建立模型研究了喷管出口距离平台3,5m工况下推进剂流量和复燃对冲击特性的影响,并分析了影响差异及其产生差异的原因.结果表明:尾焰自由射流区的激波膨胀、压缩距离和壁射流区面积随推进剂流量的增大而增大;考虑复燃化学反应不仅改变了自由射流区和滞止区的形状结构,而且增大了壁射流区的面积和温度;复燃和推进剂流量均是通过影响尾焰结构对冲击特性产生影响,具体影响效果与喷管出口和发射平台间距离有关.      

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

复合材料层合板结构冲击损伤数值模拟的损伤力学模型

针对复合材料结构低速冲击损伤问题,基于连续损伤力学提出了一种动力学冲击条件下的三维损伤数值模型。模型中区分了层内损伤（纤维拉伸与压缩失效、纤维间拉伸与压缩失效）和层间分层损伤不同的失效模式。采用三维Puck失效准则与考虑压缩抑制效应的Aymerich准则对上述两类损伤进行判定,材料失效后基于连续损伤力学中线性软化模型对材料损伤进行演化。模型中考虑了复合材料层合板结构中子层的就位效应和损伤分析中的“连锁效应”。通过对Shi的冲击试验进行数值模拟,模型预测的冲击接触载荷、分层形状和尺寸与试验结果吻合较好,证明了所提出的数值模型对复合材料层合板结构低速冲击损伤预测的有效性。     

考虑低温增压真实气体效应的运输机气动特性数值模拟研究

为了给低温风洞试验数据修正提供参考,本文利用数值模拟手段研究了低温真实气体效应相比较于完全气体对飞行器气动特性的影响,以及该影响与雷诺数影响相比所占比例的大小等问题。文章应用 Aungier-Redlich-Kwongz 方程,发展了适用于模拟氮气低温高压真实气体效应的 RANS 求解软件。与 NIST 数据的对比表明,该状态方程在5倍大气压下,比热比等参数误差在0.3%以下。同时,标模测试结果表明本文软件计算精度与国外软件相当。应用本文方法研究低温增压风洞中氮气真实气体效应：以典型运输机构型 DLRF6为模型,分别计算了高速、低速状态下各种不同温度和压力工况下的流动。计算结果表明,在低温增压情况下,真实气体效应引起的气动力差异很小,升力、阻力、力矩最大相对误差均在0.3%以下,与雷诺数效应引起的偏差相比可以忽略不计。因此,以氮气为介质的低温风洞试验研究可以采用完全气体假设。     

多级叶盘结构动力特性分析的一种方法

因各级转子叶片数的不同,基于传统循环对称法的通用有限元程序无法直接用于多级叶盘结构动力特性分析.发展了一种方法,分别对每一级扇区模型施加相应的循环对称边界条件,级间则通过施加自由度耦合来保证位移场的连续.若对每一级叶盘结构选取合理的扇区数目进行建模,最大程度地考虑级间耦合效应,则精度可进一步提高.由于循环对称边界条件和级间耦合通过对位移自由度施加约束方程来实现,因此该方法易于在通用有限元程序系统中实现.对两级叶盘结构分别采用该方法和全环模型进行动力特性分析,两种方法的频率差控制在3%内.探讨了多级叶盘结构中级间耦合对模态的影响,在低节径模态两级叶盘结构耦合效应较为明显,随着节径数的增加,耦合效应降低.      

高超声速进气道快速破膜开启的流动特性

进气道处于起动状态是保证超燃冲压发动机正常工作的前提,进气道帽罩快速开启时的非定常效应可以有效提高进气道的起动能力。采用非定常数值计算深入研究了唇口帽罩不同安装位置开启时的非定常效应对进气道起动过程的影响,分析了不同帽罩安装位置开启时进气道流场的演化过程,并揭示了喉道分离泡的形成机理。研究结果表明,当帽罩上游不存在分离泡时,破膜非定常激波在压缩面反射,与上游复杂波系作用形成沿壁面的低速流,在唇口激波作用下在喉道形成分离泡。帽罩安装靠近唇口可通过缩短激波/边界层作用距离减小低速流动区范围,进气道临界起动内压比随之增大;而当帽罩上游出现大分离泡时,分离泡会先演变为低速流,之后在唇口激波作用下重新聚集形成大尺度分离,进气道临界起动内压比显著降低。     

一种三维并行自适应计算方法

基于串行网格划分软件METIS与并行化消息传递编程接口（ MPICH2）对现有串行自适应程序进行简单的并行化改造,给出了一种三维可压缩无粘流数值模拟的并行自适应方法。首先利用单个进程调用METIS,串行划分网格;然后对所有进程并行计算以获得初始网格下的流场解;再次利用单个进程对整个流场运用自适应方法进行局部网格加密并调用METIS串行划分网格;最后全部进程在流场初始解的基础上继续并行计算,以获得自适应网格下的流场解。数值模拟算例验证了此方法的可靠性与高效性。     

混合网格方法研究栅格翼亚跨超声速气动特性

采用基于结构/非结构混合网格的CFD方法对栅格翼的亚跨声速气动特性进行了研究,计算方法经过试验数据的验证,可达到工程精度。对单独栅格翼的研究表明,跨声速壅塞时,流动通过栅格前方亚声速气流减速从栅格外侧溢流而实现流量调节;超声速壅塞时,流动通过在栅格翼上产生网状脱体激波使来流减速至亚声速来进行流量调节,且计算结果与理论估算上下临界马赫数范围一致。对某栅格翼气动布局导弹的计算结果表明,在跨声速壅塞区间内,壅塞有自我调节作用,通过溢流使栅格通道内维持相似流动,从而使栅格气动性能保持基本平稳;同时由于溢流使流量下降,栅格翼升力效率下降,导弹静稳定性减弱。     

非结构网格限制器研究

研究了6种限制器形式及其参数计算方法,其中新的自适应VenkatakrishnanM限制器能够根据当地流场特征自动调节参数从而实现不同区域不同特性.数值试验表明:非结构网格限制器对邻居集合选择方法较敏感,其中共面法比共点法精度高约15%;高马赫数时压缩型限制器比耗散型限制器计算精度高约10%~40%;参变量算法中,方向梯度法得到的流场比增量法平滑;自适应VenkatakrishnanM限制器能够模拟复杂超声速流动.