高频不稳定燃烧的声学数值仿真

高频不稳定燃烧一直是液体火箭发动机研制过程中所要面临的重大难题之一.采用具有低色散低耗散特点的计算气动声学方法,对自燃推进剂变轨发动机(OME)的高频不稳定燃烧进行时域下的数值仿真.由Crocco的压力时滞模型对燃烧热释放和声波之间的耦合进行模拟,并对不同的时滞模型参数对稳定性结果的影响进行了分析,给出发动机的稳定性极限图,确定一阶切向及一阶径向振型为主的不稳定振型,与地面试车实验捕捉到的不稳定振型相一致.结果表明:采用计算气动声学方法对带有Crocco压力时滞模型的声波扰动方程进行时域下的数值求解,可以对发动机的高频不稳定燃烧进行成功地预测.      

高温燃气进入舵机舱过程仿真与流动机理分析

采用计算流体力学(CFD)方法研究了火箭发动机工作拖尾段高温发动机燃气进入舵机舱的物理现象。结合导弹实际飞行弹道参数变化特点和超声速流场扰动不向前传递的空气动力学理论,提出了简化而不失真的非定常流场仿真方案,显著缩短了仿真周期;复现了某型导弹实际飞行时舵机舱先被"抽气"再进高温燃气的动态过程,并分析了高温发动机燃气进入舵机舱的流动机理,即在发动机工作段,导弹底端面压强低于舵机舱内压强,舵机舱被"抽气",在拖尾段随着燃烧室总压降低,喷口附近的马赫盘向导弹底端面移动,使导弹底端面压强增大且高于舵机舱内压强,高温燃气进入舵机舱烧毁电路致使导弹折断;明确了某型导弹折断故障产生的诱因,提出了改进措施和检测方法,并得到了大量飞行靶试的验证,解决了舵机舱热防护结构可靠性问题。     

CFD数学模型的线性化方法及其应用

计算流体力学(CFD)方法不仅仅起到数值模拟的作用,它本身是一个复杂的非线性系统。在流动稳定性分析、气动弹性分析、优化设计以及流动控制等领域,从系统的角度出发,对CFD数学模型线性化后,可以对模型的系统矩阵进行定量分析获得更多的系统特性。但是CFD数学模型往往非常复杂且阶数很高,因此其线性化系统矩阵的获得比较困难。鉴于此,采用人工编程和自动微分相结合,构造有限体积法并行CFD模型的线性化系统矩阵。其中自动微分只被用来得到每个界面通量的局部雅可比矩阵,而采用人工编程方法来实现并行环境下的稀疏雅可比矩阵的组装。线性化系统的并行求解采用了块雅可比预处理的广义最小残量法,每个并行进程内部则采用零填充不完全LU分解预处理。为了验证这种线性化方法,上述方法被用于:① NACA 0012翼型的非定常绕流线性系统构造与求解;② NACA 0012翼型稳态流动的伴随方程构造与求解;③ AGARD wing 445.6机翼颤振问题降阶建模。上述三个算例的结果与CFD模拟的吻合一致。     

旋翼VPM/CFD耦合模型的建立及其在小前进比气动分析中的应用

为了提高旋翼计算流体力学（CFD）流场计算效率及克服其在尾迹捕捉上的不足,将旋翼黏性涡粒子方法（VPM）与CFD分析相结合,建立了一个新的旋翼VPM/CFD耦合气动分析模型。在该模型中,采用VPM分析以实现对旋翼尾迹中黏性涡的高效捕捉而不引入数值耗散,而采用CFD分析用于精确地模拟旋翼桨叶近体区域内复杂的流动现象,同时也为VPM分析提供一个较高精度的涡源模型。至于两者信息交换,则使用集中涡源法将CFD信息传递至VPM分析,而VPM计算得到的尾迹信息则通过边界修正施加至CFD域的远场边界上,从而可以鲁棒地实现CFD域与VPM域的耦合计算。在此基础上,对“Helishape 7A旋翼”小前进比前飞桨-涡干扰（BVI）状态进行了较为深入的气动分析,计算结果表明：与全CFD计算比较,建立的VPM/CFD耦合分析模型可以有效地避免旋翼尾迹区桨尖涡的数值耗散,从而更加可靠地捕捉桨-涡干扰状态下的桨叶非定常气动载荷脉动,同时对于本文算例,计算效率可以提高30%以上。     

吸气式高超声速飞行器冷流试验设计及验证

对于吸气式飞行器而言,地面冷流试验是检验其进气道性能及气动特性的一项重要手段.以二元混压式进气道、机体/推进系统耦合为基本特征,设计了采用超燃冲压发动机为推进系统的内外流一体化巡航飞行器,针对其高超声速特性开展了冷流风洞试验,来流速度范围Ma=5.0~7.0,攻角范围α=-4°~8°.测压试验结果表明,随着来流马赫数的增大,进气道的总压恢复系数下降;而流量系数先上升,在设计点达到最大值;在一定攻角范围内,进气道的总压恢复系数和流量系数提高,但当攻角增大至巡航攻角时,随着攻角的增大,进气道的总压恢复系数和流量系数逐渐下降.测力试验验证了数值算法的有效性,除轴向力系数以外,其余气动特性系数的发展规律及数值基本吻合,可通过修正试验值的方式外推出飞行器的气动特性数据.     

基于风洞试验的非定常气动力微分方程建模方法

为准确描述大攻角非定常运动中非定常气动力的气动特性,研究了以微分方程为基础的非定常气动建模方法。在某战斗机大振幅偏航、滚转单自由度及偏航-滚转耦合运动风洞试验的基础上,对偏航力矩、滚转力矩建立微分方程模型,分析研究不同频率、同一侧滑角情况下,气动力的迟滞时间;然后对高频运动下的气动力进行拟合,进而建立高频运动下的气动力微分方程模型。研究表明:模型可精确预测不同机动下非定常气动特性;建模方法具有较强的工程可行性。     

某型串列静子的三维气动数值优化设计

针对某型高负荷轴流压气机串列静子近轮毂处存在的分离,利用自行开发的压气机三维气动数值优化设计平台,选择复合形优化方法,通过修改串列静子近轮毂处的三维积叠方式对串列静子进行了优化设计.结果表明,优化设计取得了良好效果,在消除分离的同时,有效提高了压气机的整体性能.      

基于GPU和隐式格式的CFD并行计算方法

从图形处理器(GPU)架构特点出发,提出了基于数据并行的隐式计算流体力学(CFD)求解方法,空间离散格式采用迎风Roe格式,计算网格适用于结构和非结构网格。采用统一计算设备架构(CUDA)技术实现了GPU上的隐式CFD并行计算。分别在Intel Core2 Quad3.0GHzCPU和NVIDIAGTX280 GPU上进行了计算,结果表明隐式格式计算速度是显式格式6倍以上,采用显式格式的计算加速比达到28倍,采用隐式格式计算加速比达到了28.7倍,同时计算加速比随计算规模的增加而增加。计算结果和实验结果较为吻合。     

重叠网格洞面优化技术的改进与应用

针对传统割补法在壁面狭缝处网格重叠失败的问题,提出了以单元顶点状态进行洞面优化的vertex方法及以单元大小进行洞面优化的volume方法的两种改进方法。这两种方法将传统割补法中切割和填补两步搜索方法改为一步搜索方法进行重叠网格的洞面优化。通过数值试验证明:两种洞面优化方法都能很好地解决壁面狭缝网格重叠问题,网格重叠区域流场数值传递正确,计算结果与试验值吻合,优化效率高,收敛速度快。volume优化方法与vertex优化方法相比可以减少由网格大小不对等引起的数值误差,计算结果和计算精度均优于vertex优化方法。     

严格随机正则(3,s)-SAT模型及其相变现象

研究变元和文字出现次数受限制的规则3-SAT问题,提出了一种严格随机正则（3,s）-SAT问题,并给出了该问题的实例产生模型--SRR模型。结合一阶矩方法和生成函数展开项系数的渐近近似技术,证明了严格随机正则（3,s）-SAT问题相变点的上界,即当变元规模N较大且变元出现次数s>11时,严格随机正则（3,s）-SAT实例是高概率不可满足的。实验结果表明：由SRR模型所生成的随机实例中,当N>60且s>11时,所有的（3,s）-SAT实例均是不可满足的,而当N>150且s<11时,所有的（3,s）-SAT实例均是可满足的,即严格随机正则（3,s）-SAT实例的相变点位于s=11处,且在s=11处（子句变元比为11/3）的严格随机正则（3,s）-SAT实例,比在相变点（子句变元比）4.267处同规模的均匀随机3-SAT实例更难求解,因此,SRR模型可以很方便地在s=11处构造难解的随机3-SAT实例。