飞机起飞性能计算中发动机推力确定方法的改进

为准确计算飞机起飞滑跑距离和起飞航迹,应用数值插值方法,对起飞性能计算中发动机推力的确定方法进行了改进,提出以发动机推力曲线图所给推力为基础,利用分片双二次插值函数法确定发动机推力。基于改进算法对某型飞机发动机推力和起飞航迹进行了实际计算,计算结果表明,改进后的方法比原有方法更准确。     

飞机起飞滑跑发动机推力数值确定方法

利用数值理论对飞机起飞滑跑距离计算中发动机瞬时推力的确定方法进行了研究,提出以发动机推力曲线图中已知推力曲线为基础,通过拉格朗日插值确定未知推力的算法。基于所提算法,对某型飞机在2 000 m压力高度上的推力进行计算,并与从推力曲线查到的实际数据比较,二者的平均相对误差为3.2‰。用所提算法确定推力,对该型飞机在5种条件下的起飞滑跑距离进行计算,计算结果与实际滑跑距离的平均相对误差为3.0%。计算结果表明,算法是准确可行的。     

非结构嵌套网格中的一种改进型径向基函数插值方法

非结构混合网格拓扑关系相对复杂,在嵌套边界区域进行流场信息插值传递容易有精度的损失。利用已知的解析函数作为插值样本,对线性插值(LINE)、距离倒数权重插值(WA)、逆二次径向基函数插值(IQ)、Wendland暞sC2径向基函数插值四种方法的精度进行对比,筛选出精度较高的C2径向基函数;提出一种改进的C2径向基函数插值方法,通过调整径向基函数的作用半径来控制插值矩阵条件数,进而消除使用径向基函数插值在嵌套边界区域出现的数值奇性问题;选用MD30P/30N 多段翼与三维AEDC典型外挂物分离模型对该插值方法进行验证。结果表明:改进后的径向基函数插值方法能够有效消除嵌套边界区域数值奇性,同时计算收敛速率更快,插值精度更高。     

模型两相脉冲爆震发动机推力的测试与研究

应用推力间接测量法,对模型两相脉冲爆震发动机的瞬时推力进行了详尽的测试,并依此计算了相应的平均推力。分析比较所得结果表明,该模型能产生很大的瞬时推力。      

基于本征正交分解的流场快速预测方法研究

对于飞行器气动设计,通过实验或CFD计算获得所有设计状态处流场信息成本高、耗时长,难以满足工程应用要求。利用基于本征正交分解(POD)降阶模型结合相应的插值方法,只需要通过实验或CFD计算获取有限个数的采样流场,可以快速预测出具有满足精度要求的设计状态处流场信息。首先,证明对POD基系数插值获取预测流场与对采样流场直接插值获取预测流场的等价性,并分析采样流场分布的要求。然后,以三段翼型流场作为研究对象,POD降阶模型分别与三次样条插值和径向基插值结合,实现设计状态处流场信息快速预测。最后,通过将流场预测结果与相同条件下的流场CFD计算结果对比,表明了POD降阶模型结合两种插值方法的各自特性。     

径向基函数在高超声速热流插值中支撑点选择的判据研究

在使用径向基函数进行高超声速热流插值时,为了避免传统单一绝对误差判据选择得到的支撑点在热流较小处插值不准确的问题,提出一种双重误差判据下的径向基函数插值过程。这一过程首先使用绝对误差判据选择一定数量的支撑点,再使用相对误差判据标准选择另一部分点。通过数值实验验证同时采用以绝对误差和相对误差为标准选取支撑点的方法,结果表明：采用双重误差判据可以同时保证了热流较大与较小处的插值精度,尤其能避免传统单一绝对误差判据导致插值结果中容易出现负热流的问题。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

涡流阀固体变推力发动机控制流参数影响规律研究

为快速求解涡流阀固体变推力发动机的工作参数、研究控制流参数对涡流阀固体变推力发动机的影响规律,改进了求解涡流阀固体变推力发动机工作参数的理论计算模型,并利用该模型研究了控制流参数对涡流阀固体变推力发动机推力调节性能的影响规律。结果表明:该模型最大相对误差为10.83%,可用于涡流阀固体变推力发动机的计算分析;在控制流压强与控制流流量相同及分子量相当的情况下,提高控制流温度可以提高推力调节比;控制流流量是决定涡流阀固体变推力发动机推力调节比最重要的参数;增加燃气发生器的燃面,不仅可以提高涡流阀固体变推力发动机的推力调节比,同时可以减小发动机的比冲损失。     

基于RBF插值边界条件的笛卡尔网格方法

介绍了一种基于RBF(Radial Basis Function,径向基函数)插值处理物面边界条件的笛卡尔网格方法.整个流场采用有限体积方法计算,但是对于物面附近切割的网格计算域不完整,因此通过RBF插值来得到物体内部实单元上的值,从而使得流场计算顺利的进行推进.相比于其他插值方法来说,RBF插值特别适合于多维离散点的插值,不会因为插值点的位置关系而产生病态的系数矩阵,而且使用比较少的网格点就可以得到同最小二乘方法相当的精度.将这种方法应用于圆柱和翼型绕流中,得到了满意的结果,证明了方法的有效性,还对边界条件的各种不同形式对结果的影响作了对比.     

两种高精度高效率静气弹数据传递函数

遵循静气弹问题中气动/结构网格之间的数据传递要求,引入了基于径向基函数的流固耦合界面数据传递方法,确定了紧支半径的选取原则,并编制相应的计算程序。以M6机翼为例,将不同径向基函数的插值精度及计算效率进行比较,结果表明体样条（Volume Spline）函数和C2函数具有良好的插值精度和计算效率。这两种流固耦合界面数据传递方法非常适用于工程应用。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

基于物理量梯度修正的RBF数据传递方法

径向基函数是复杂外形飞行器多场耦合计算中重要的数据传递方法。复杂飞行器物理量分布通常与局部流动特征密切相关且在不同区域变化剧烈,如何在具备各向同性特征的径向基函数中考虑实际物理量各向异性分布因素的影响成为制约其精度进一步提升的关键。针对径向基函数的问题,提出了一种基于局部物理量梯度对径向基三方向权重进行自适应修正的数据传递新方法。采用高超声速舵面、翼身组合体算例对该方法在单次数据传递中的可行性及效果进行了验证,同时基于该方法对压缩拐角外形开展了气动力/热/结构多场耦合数值模拟,并与风洞试验结果进行了对比分析,以验证该方法在实际耦合问题计算长时间高频次数据传递中的适用性及可靠性。研究结果表明该方法能够很好地提升如薄板样条插值（TPS）、多重二次函数插值（MQ）等全域基函数的单次插值效果,并且有效地改进了MQ方法对形态参数的鲁棒性,而对于梯度修正后的紧支基函数,其改进效果更佳,以较少的点达到了TPS、MQ等全域基函数的插值效果,实现了数据传递效率和精度的兼顾。采用该方法获得的耦合计算结果与试验对比一致,表明了该方法的可用性以及在复杂外形多场耦合问题中的良好工程应用前景。     

航空发动机压气机径向间隙设计方法研究

简述了航空发动机压气机径向间隙控制的意义、作用,以及径向间隙设计选择的依据。对发动机压气机径向间隙的计算方法进行了研究,提出了更加完善的径向间隙设计方法。采用该方法对某发动机压气机径向间隙进行了设计,具体分析了影响径向间隙的因素及其变化规律。试验验证表明,使用该方法能获得安全可靠的压气机径向间隙。     

容积与转子状态量在变循环发动机数值仿真中的应用

使用容积动态与转子动态来计算容积与转子的状态量,为压气机及涡轮的特性图插值与计算过程提供输入数据,使得变循环发动机数值仿真的运算过程无需猜值与迭代;并给出了对变循环发动机的可调几何的仿真方法。使用此模型仿真了从双外涵模式向单外涵模式切换的过程以及这两种模式下变循环发动机的高度、速度、节流特性。在高度0km且马赫数为0的大气环境下模式转换的仿真结果显示,变循环发动机涵道比由0.700下降到0.307,单外涵模式与双外涵模式相比推力提高了16.823%,单位燃油消耗率增加了15.163%。     

某型大涵道比涡扇发动机飞行推力确定方法研究

以某型大涵道比涡扇发动机为研究对象,结合动力装置适航取证＂飞行推力确定＂试验科目,对该型发动机飞行推力确定方法进行了分析,并开发了该型发动机飞行推力确定计算分析软件。结合发动机设计单位给出的设计值,对该计算方法进行了验证。结果表明,采用飞行推力确定计算方法得到的计算值与设计值间吻合程度较高。     

斜切喷管火箭发动机推力的近似计算

本文根据火箭发动机推力等于作用在火箭发动机内外表面上压力的合力这一基本公式,引入在一定条件下喷管扩张段内的燃气压力与喷管截面半径成线性变化的关系,导出大料切角时斜切喷管火箭发动机的轴向推力和侧向力理论计算的近似公式。计算结果与六分力推力试验的测试结果相比颇为接近,说明该公式具有较好的计算精度,可以用于涡轮式火箭弹推力和转矩计算以及某些反推发动机的推力计算。     

两相脉冲爆震火箭发动机性能实验

为了获得脉冲爆震火箭发动机（PDRE）的性能参数,采用液态煤油为燃料、氧气为氧化剂、压缩氮气为隔离气体,进行了一系列多循环爆震实验。使用孔板流量计测量煤油流量,使用集气法测量氧气流量,使用动态压电式压力传感器测量了爆震室轴向的沿程压力,使用火焰温度及水蒸气浓度红外光谱测量仪测量爆震管出口平面的尾焰温度,使用动态压电式推力传感器测量PDRE所产生的瞬时推力。实验获得PDRE不同频率下的平均推力和比冲。实验结果表明：爆震压力和温度随着工作频率的变化而有所变化,填充系数对于PDRE比冲大小有着显著影响。采用爆震室部分填充的策略,可以显著地提高发动机比冲。     

一种非结构网格上基于径向基函数重构的ENO格式

基于二维Euler方程,对非结构三角形网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个三角形单元构造插值径向基函数,而在计算交界面的流通量采用两点高斯积分公式以保证格式的整体精度。时间离散采用三阶TVD Runge—Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

固液火箭发动机试验瞬时燃速分析方法

介绍了一种用于固液火箭发动机试验的瞬时燃速分析方法,运用该方法对H2O2/HTPB固液火箭发动机单室双推力长时间热试车试验进行燃速分析,拟合得到该工况下的燃速公式为r=1.847×10-2G0o.7304。根据燃速公式等计算结果对一次发动机试验进行了预估,计算得到的预估内弹道性能曲线与试验结果吻合较好,验证了该瞬时燃速分析方法的可行性,为发动机工作时间较长、氧化剂流率变化较大时的燃速分析提供了一种途径。