三维不可压边界层抛物化稳定性方程的椭圆特性研究

针对三维平行不可压边界层流动的抛物化稳定性方程,研究其剩余的椭圆特性,提出了消除该特性的具体方法,并对抛物化稳定性方程解与全椭圆问题解的一致性进行了证明。     

人工释能法求解线性代数方程组

本文建立求解线性代数方程组的人工释能法迭代格式,并给出该迭代格式收敛的充分必要条件和最佳释能时机,同时给出实现人工释能法的数值计算方法及其数值稳定性条件。     

大型军用运输机的飞行航迹优化

军用运输机通常在紧急情况下执行空运任务,应以时间成本作为此类飞机垂直剖面飞行航迹的优化指标.将飞行时间按上升、巡航和下滑3个阶段进行分解,采用能量状态法简化飞机质点运动方程,将优化指标转化为等效能量形式.优化巡航段参数,采用最小值原理分别完成纵向最优剖面上升和下滑飞行航迹的优化.仿真结果表明:军用运输机在上升段采用各高度的最大可用速度上升,在低于升限的某高度以最大飞行速度巡航,在下滑段采用各高度的最大可用速度下降时,全剖面飞行时间最短.研究结果对提高大型军用运输机在可能存在威胁的前线环境下的运输效率和飞行安全性等均具有一定的实用参考价值.     

基于速度剖面拟合的航空器场面滑行4D轨迹预测

为提高机场的整体运行效率和综合保障能力,提出了一种基于速度剖面拟合的航空器场面滑行4D轨迹预测方法。构建了航空器的滑行动态模型,进而引入标称速度剖面的概念。将DSW算法应用到速度剖面的拟合中,得出一种标称速度剖面的生成方法,并通过建立平均速度修正参数修正不同机型对滑行速度的影响,基于动力学平衡方程,构造了标称速度到瞬时速度的映射,结合BADA数据实现对瞬时速度的修正。在上述分析的基础上得到航空器场面滑行4D轨迹预测模型。案例表明,与基于动力学的方法相比计算结果更加准确,使平均误差降低47．3％,能够有效地预测航空器的4 D轨迹。     

旋转弹气动力建模与飞行轨迹仿真

运用一种自回归滑动平均(ARMA)的时域气动力建模方法,以计算流体力学与刚体动力学(CFD/RBD)耦合仿真的输出结果为样本,对旋转弹的非定常气动力进行建模。利用建立的气动力模型与刚体动力学方程求解模块耦合,实现了旋转弹轨迹的快速仿真,并讨论了不同的建模方式对仿真精度的影响。算例结果表明:采用气动力模型与刚体动力学方程耦合仿真技术可以在不同初始发射条件下进行旋转弹飞行姿态与运动轨迹预测,且与CFD/RBD仿真结果吻合较好,证明ARMA气动力建模方法可以在保证旋转弹轨迹预测精度的同时大幅缩短仿真时间,节省计算资源。     

弹射过程座舱盖抛放轨迹的预测

通过分析座舱盖在各个不同抛放阶段的受力关系,并着重考虑了抛放过程中气动系数的修正和铆钉剪切的影响,最终建立了座舱盖抛放轨迹的预测模型.该模型可广泛用于预测不同型号座舱盖在各种飞行状态下的抛放轨迹.算例表明计算结果与实验数据一致,证明了模型的实用性和有效性.     

混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

基于粒子轨迹的结冰风洞收缩段优化设计数值模拟

为了研究收缩段型面和液滴轨迹特性对结冰风洞流场品质与结冰试验效果的重要影响,获得收缩段的优化设计参数,采用数值模拟方法,使用二维简化模型对3m×2m结冰风洞收缩段的流场特性进行了研究,对几种典型收缩曲线的流场进行了计算,对比了收缩曲线对沿轴向压力分布、出口截面速度分布、不同粒径云雾粒子的极限释放距离及试验段入口云雾粒子分布的影响,结果表明:三次曲线与x_m为0.4的五次曲线不能满足流场合格指标,其他曲线流场特性差别不大,维氏曲线具有最小极限释放距离,在收缩段出口,三次曲线和五次曲线粒子包络面积最大,综合比较,x_m为0.5的五次曲线能有效兼顾流场均匀性、极限释放距离和收缩段出口粒子包络面积比.      

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

悬停状态直升机桨叶扭转分布的优化数值计算

建立了一套基于高精度计算流体力学(CFD)技术和代理模型优化算法的旋翼气动外形设计方法。在该方法中,旋翼流场气动性能的计算采用了基于Navier-Stokes/Euler方程的CFD方法,并根据流场特点、精度和效率的要求采用Baldwin-Lomax(B-L)湍流模型,通量计算采用Roe-MUSCL格式进行。为提高网格生成质量和便于流场控制方程的求解,将流场分成两个区域,即围绕旋翼的黏性区和无黏的背景网格区。其中,桨叶网格使用了基于二维翼型网格的参数化方法生成,数值计算结果表明该方法有效地提高了网格生成质量及效率。在参考旋翼流场及桨叶细节流动分析的基础上给出设计变量及范围,有效减小了优化问题的规模;为满足优化和机理分析的需要,将基于置换遗传算法优化的拉丁超立方(PermGA LHS)方法和径向基函数(RBF)的代理模型优化方法引入到桨叶外形的优化设计中。首先以Helishape 7A旋翼为算例,检验了数值模拟方法的准确性。然后,应用所建立的优化方法针对旋翼负扭转分布进行了优化计算,结果表明优化后的旋翼悬停气动性能比优化前有了明显提高。