化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

叶轮行波振动下的气动阻尼计算

为展开运动边界下离心叶轮流场的数值分析,独立开发了网格变形程序和非定常流动模拟程序,实现了流场中振动离心叶轮的气动阻尼计算。采用紧支径向基函数法进行结构到气动表面变形的数据传递,应用二叉树技术进行壁面距离的计算,大幅提高了网格变形和流场分析中距离搜索的计算效率。通过振动叶栅和径向叶轮的算例,验证了程序应用于运动边界流场计算和叶轮流场模拟的正确性。以半开式叶轮为对象,展开其在行波振动条件下气动阻尼的计算分析。结果表明,叶轮在前、后行波振动下的气动阻尼均为正,后行波振动引起的气动推力和气动阻尼较大,叶轮表面的平均气动功率密度呈现循环对称的分布特征。     

蠕变时效成形回弹预测及工装型面优化

蠕变时效成形是一种将蠕变成形与时效强化同步进行的成形工艺.由于卸载后工件的回弹量较大,因此回弹的预测和成形工装型面的优化很重要.以LY12的蠕变时效成形为研究对象,对经典蠕变本构方程进行了修正并将其开发为ABAQUS/CREEP子程序.对板材的蠕变时效成形特点进行了模拟研究,建立了LY12的回弹预测函数.试验结果表明:其预测准确性较高.将变偏差调节算法与ABAQUS/CREEP相结合用于某蒙皮的成形工装型面的优化,通过7次迭代,型面精度达到97.9％,结果表明该组合算法优化速度快、精度高.     

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

基于人机闭环稳定性的舵机速率限制边界

舵机速率限制是造成电传操纵飞机人机耦合的主要原因。利用描述函数法对舵机速率限制非线性进行建模,分析舵机速率限制非线性、人机耦合发生频率及人机闭环稳定性间的关系,并基于此提出舵机速率边界的确定方法;以典型放宽静稳定性飞机为例,基于最优McRuer驾驶员模型,确定人机闭环稳定性所需的最小舵机偏转速率;基于开环起始点(OLOP)准则对所确定的速率限制边界进行验证。结果表明:本文提出的舵机速率边界的确定方法最小成本地避免了人机耦合;所确定的舵机速率限制边界与OLOP准则边界对应的舵机速率基本吻合,即所建立的舵机速率限制边界确定方法合理。     

基于不可逆热力学的宏细观粘弹性损伤本构模型

为建立复合固体推进剂的损伤本构模型,基于不可逆热力学和叠加原理,通过引入宏观损伤效应张量,推导出一个通过有效应力表征损伤的蠕变型损伤本构方程。假设材料为初始各向同性,进一步引入细观标量损伤效应函数表征材料对称性的改变,进而得到由细观损伤效应函数表征宏观损伤效应张量的一般表达式。通过选取合适的细观损伤效应函数,文中建立的本构方程可以用于表征材料的正交各向异性损伤、横观各向同性损伤和各向同性损伤。最后,基于Schapery粘弹性微裂纹扩展模型,选取相对微裂纹密度为损伤内变量,建立了一个由损伤热力学对偶力表征的幂律型损伤演化方程。数值结果表明,建立的模型能够较好地反映材料损伤的率相关性和温度依赖性,具有良好的预测精度。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

空心风扇叶片结构优化设计方法及程序实现

在空心风扇叶片结构强度分析的基础上,以叶片最大等效应力水平为约束条件,建立了叶片质量和径向位移最小化的双目标优化模型.采用正交试验设计的方法,分析了各设计变量对约束函数和目标函数的影响,减少了设计变量的数量.为了提高优化设计效率,采用径向基函数插值的方法,构造了约束函数和目标函数的响应面替代模型,从而避免了优化设计过程中大量的结构有限元分析求解.针对空心叶片结构强度优化设计的示例,探讨了替代模型和带精英策略的Pareto排序遗传算法的具体应用,得到了分布均匀的Pareto最优解,给出了空心叶片示例的具体优化设计结果.      

用Tam-Ganesan湍流模型研究排气系统红外辐射特性

通过三维流场、介质的质量浓度场与辐射场耦合求解方法,计算了1台有中心锥的轴对称喷管的空腔-喷流组合红外辐射特性.利用专门为自由剪切流动发展的Tam-Ganesan湍流模型,有效提高排气系统红外辐射模拟精度.采用有限体积法(FVM)结合窄带模型(NBM)计算了排气系统在3~5μm波段的辐射传输问题,并考虑大气衰减作用.对比了采用标准 k - ε 模型和Tam-Ganesan模型对排气系统核心区长度的影响,以及由此引起的排气系统红外辐射特性的变化.通过与试验数据对比证明了该结果的准确性,结果表明选择合适的湍流模型对计算排气系统红外辐射特性有重要的意义.      

基于RBF神经网络的开关磁阻电机在线建模 及其实验验证

 为了获取开关磁阻电机(SRM)的精确模型,提出了一种基于径向基函数(RBF)神经网络对SRM进行建模的方法,主要包括离线建模和在线建模两部分。离线建模通过实验测量得到SRM的磁链特性曲线,并利用该数据训练RBF神经网络,实现SRM磁链的离线建模;在线建模是指当SRM的运行状态发生变化时,离线模型的估计磁链与实际磁链会产生误差,通过对神经网络的输出权值进行在线调节,实现具有在线动态调节功能的SRM在线模型。为了验证该方法的可行性,针对一台12/8结构的SRM样机进行仿真和实验,结果表明SRM的离线模型和在线模型在仿真和实验条件下均能正确地估计SRM的磁链特性,而且在线模型的估计精度高于离线模型,验证了本文的研究方法合理有效。