跨音松弛法在叶栅绕流中的应用

文章介绍跨音松弛法应用于叶栅绕流的定常与非定常流场计算。在定常流中,采用弦线法向小扰动假定;在非定常流中,采用线化非定常分量小扰动假定,从而简化了数值求解过程。由于松弛法大大地节省机时和内存,便于工程上推广应用。     

对旋轴流风机非定常流场数值模拟

利用全隐格式的非稳态压力修正求解方法,分析和研究了对旋叶轮级间流场非定常流动情况,上游尾迹与下游叶片位势作用对流场的影响,并与动态实验结果进行了有益的对比.结果表明:非定常流场的计算更能反映叶轮机实际运行性能,如果能很好地组织对旋叶轮的非定常流动,还可以进一步提高其压比和效率.     

局部振动对火星环境下薄翼型气动性能的影响

火星的稀薄大气环境迫使无人机在亚临界雷诺数范围工作，低雷诺数层流分离问题给无人机气动性能带来极其不利的影响。同时，火星大气的声速较低，使无人机运行的马赫数更高，压缩效应增强并可能产生激波。为研究火星环境下翼型局部振动的流动控制作用，采用基于动网格的数值方法对非定常流场进行模拟。选取NACA5605低雷诺数薄翼型，雷诺数为1.5×10⁴，马赫数为0.43和0.63。时均流场和时均气动力系数结果显示：翼型局部振动能够明显减少时均分离区的大小，起到增升减阻的作用。非定常流场表明流动控制机理在于振动产生的涡流运动抑制了翼型尾缘附近的层流分离。研究了不同振幅、频率和振动位置下的流动控制效果。最佳参数下，马赫数为0.43时升阻比最多提高24.7%，马赫数为0.63时升阻比最多提高52%。     

求解跨音位流问题的改进SIP方法

本文研究Sankar提出的,用全位势方程求解跨音流场的强隐式方法(Strongly Implicit Procedure)。针对强隐式方法在求解跨音流场中的问题作改进,即采用了新的人工压缩密度计算方法,计算绕孤立机翼、叶栅通道的跨音流场,降低了在跨音位流计算中激波前后的数值振荡,并提高了算法的稳定性,     

浮空气囊定常及非定常气动特性计算

对于0~20 km范围内不同外形的浮空气囊,采用虚拟压缩方法求解不可压缩非定常N-S(Navier-Stokes)方程,数值模拟了5个高度下气囊的定常和非定常流场.定常状态时,随高度及气囊外形的变化,模拟了气囊气动力及背风区复杂流动分离.结果表明,当气囊厚度率逐渐增大时,其背风区存在复杂的空间流动分离.运用双时间步方法及动网格技术,成功实现了气囊平移或俯仰运动状态下非定常流场的数值模拟.考察了不同的平移方向及速度对气囊非定常气动特性的影响,得到了气囊平移过程中气动力及流场的变化情况.对于气囊的俯仰运动,获得了不同俯仰速度及不同最大俯仰角时气囊的气动力及流场.计算结果表明气囊的非定常运动对其流场和气动力影响显著.     

高超声速飞行器滑模控制参数整定方法设计

本文针对高超声速飞行器滑模控制人工试凑的参数整定方法较为繁琐、效率低的问题,改进了连续动作学习自动机（Continuous action reinforcement learning automata, CARLA）算法并将其应用于滑模控制参数整定问题中,改进后的算法通过对回报函数的设计有效克服了常规CARLA算法收敛速度慢、易受干扰、求解效率低的问题。该算法引入控制性能指标评价函数,在迭代中学习阶跃响应的经验数据,实现了控制参数自整定。仿真表明,对于阶跃响应问题,本文提出的算法能够在100次迭代中整定出一组高品质的控制参数,完成对给定指令的快速准确跟踪,与遗传算法、模拟退火算法相比,在求解速度上具有显著优势。由于该方法不依赖于模型,除了滑模控制器参数整定外对其他控制方法的控制参数整定问题也有一定适用性,具有推广应用价值。     

翼伞充气过程的流固耦合方法数值仿真

为研究冲压式翼伞折叠充气过程的流固耦合动力学特性，基于自由曲面变形理论建立了多气室冲压翼伞的展向折叠模型。流体域通过时步更新技术实现了随伞载系统运动，采用任意拉格朗日-欧拉（ALE）方法开展了翼伞非定常充气展开过程的非线性动力学数值计算，数值计算结果与空投试验结果具有较好的一致性。深入分析了翼伞充气过程中的三维外形及非定常流场分布情况，表明翼伞充气过程由于翼尖涡绕流，存在“翼尖上翘，中部凹陷”的翼伞尾流再附现象；各气室的充气规律关于中央气室对称；分析了翼伞气动特性的动态变化规律，充满后翼伞滑翔比稳定在2.24。上述研究为翼伞设计及开伞性能预测提供了一定的理论依据。      

低反力度对高负荷跨声风扇进口级的影响

认为低反力度是提高跨声速风扇/压气机进口级气动负荷的必然结果.通过气动负荷系数和反力度的定义,得到了反力度随气动负荷系数的变化规律,并给出了气动负荷系数的两个临界值.以成功设计并实验验证的某单级高负荷风扇为基准,通过数值模拟对4组跨声速风扇转子方案进行比较研究.发现随着气动负荷系数的提高,转子轮毂区的气动参数分布发生了显著地变化.认为静子的展向跨声速来流是高气动负荷跨声速风扇进口级设计面临的一个问题.通过讨论,认为静子跨声速来流问题可能通过已有的设计手段解决.从研究方案中选取了一个最优的单级风扇方案,其气动负荷系数达到0.53,并分析了其三维粘性流场.      

圆弧翼型跨声速流动的动态模态分析

跨声速翼型的激波周期性自激振荡会给机翼结构带来附加的脉动载荷,从而加剧飞行器表面结构的疲劳损伤。使用动态模态分解（DMD）方法研究了跨声速下绕厚度18%的对称双圆弧翼型的压力脉动场,分析了DMD提取的各阶主模态的频率特征、压力脉动的空间分布以及压力脉动随激波振荡的时间演化过程,并使用DMD模态进行流场重构。结果表明,DMD方法能准确捕捉流场各特征频率的模态,第1阶模态是激波抖振的主频,在激波的自激振荡过程中占主导作用,前7阶模态重构的流场损失函数降低至4%以内,误差主要分布于激波间断处。      

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

航天数学中积分方程法的作用

文章通过实例阐明用积分方程法解航天科学技术中的定解问题,利用线性化积分理论求得显式解;利用卷积积分简化解题并得明确的结论,还可使解域扩大等优点。     

二维非结构网格的可压缩Euler方程流场解

讨论了非结构网格的生成方法,用Delaunay三角化方法和推进阵面法相结合的方法生成了非结构网格,提高了生成非结构网格的自动化程度和网格的质量.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼形和两段翼型的算例,并与实验结果进行了比较和分析.     

一种面向模块化可重构机翼的分步补偿优化方法

针对模块化可重构机翼结构的优化设计问题,以沿展向分布的3个翼段模块为研究对象,研究了不同翼展方案机翼的载荷相关性,通过在优化迭代过程中自动调整设计空间,解决了模块化可重构机翼优化设计时复杂的变量-约束耦合影响问题,建立了适用于模块化可重构机翼结构的分步补偿优化方法。以某无人机模块化可重构机翼结构优化设计问题为例,建立了优化模型,并分别采用所提分步补偿优化方法和传统单方案优化方法进行了优化设计。结果表明:所提方法能够稳定收敛,与单方案优化结果相比较,所提方法以较小的重量代价满足了3种重构方案的所有设计要求,且优化结果具有较好的工程实用性。      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

跨音速机翼设计方法的改进

以Takanashi提出的"正反迭代、余量修正"的原理为基础,针对其当机翼上表面存在较强的超音速流动时,解的收敛性不好的问题,对Takanashi方法进行改进.具体方法就是在控制方程的双曲部分引入基于迎风格式的附加项.通过对某干线飞机机翼的设计,证明了改进后的方法是成功的,设计结果在超音区有很好的收敛性.      

跨音速翼型的气动最优化设计方法

将计算流体动力学(CFD)与最优化技术相结合,通过数值求解欧拉方程,对翼型绕流流场作出了数值模拟.再利用几何、流动和最优化控制方程,反复迭代求得在一定约束条件下气动性能最优的翼型.本文以NACA0012为原始翼型,选取两种设计工况,都取得了满意的结果.     

一种利用正交流线坐标系求解跨音速管流的新方法

本文以正交流线坐标系为工具,采用位流模型,导出了正交流线坐标系下平面或轴对称流动的控制方程组,利用流函数与位函数的存在性,进而导出了两个首次积分。设计了一套求解匀熵管流的分凝迭代方案,模拟了管内无激波的平面或轴对称喷管的跨音速流动。数值实验表明该方法能够处理小喉道曲率半径平面或轴对称喷管跨音速流动,与已有的数值计算结果及实验结果符合较好。