混合解析/数值方法在含高频振荡点源项的欧拉方程数值计算上的应用

 通过模型方程分析看到了一个重要的现象, 如果源项涉及的时间尺度远小于对流项的时间尺度, 那么基于对流时间尺度作为步长的传统数值方法, 即使源项相对于墩流项是个小量, 也会导致平均尺度上错误的结果。为了克服这种困难, 采用时间分裂方法, 把方程分裂成含对流项部分的偏微分方程( PDE) 和包含源项的常微分方程(ODE) 基础上, PDE 使用传统的数值方法, ODE 用解析的方法求解。该混合方法在数值格式时间步长小于平均流动时间尺度时, 得到正确解, 而与点源所隐含的时间尺度无关。把这个方法应用在含源流动的Euler 方程的计算中, 计算了翼型振荡问题, 取得了理想的结果。     

带源项的守恒系统时间分裂算法定常解的收敛性研究

时间分裂算法很难获得带源项的守恒系统的定常解,本文中,考虑使用基本的隐式方法积分线性、二次和三次源项,研究了时间分裂算法不能收敛到数值定常态的原因.采用一个刚性参数的方式预测定常态的误差,刚性参数依赖于特定的源项.当源项是解的线性函数时,不存在定常态误差;当源项是解的非线形函数时,发现定常状态的误差是刚性参数的单调增加函数.定常态误差的分析将推广到标准k-epsilon双方程湍流模型计算的情况.     

航空支架类零件加工质量提升实践

针对航空支架类零件刚性差、装夹与加工过程中易变形的问题,以某型机载光电探测系统中上支架零件为例,提出了三种支架类弱刚性零件的装夹方法,提升支架类零件的加工质量,为其他支架类零件的装夹提供借鉴。     

Φ3.2 m风洞共轴刚性旋翼试验台研制

为满足共轴刚性对转旋翼气动特性及流动机理研究的风洞试验需要,在中国空气动力研究与发展中心Φ3.2 m风洞开展了共轴刚性对转旋翼试验台（以下简称共轴旋翼试验台）研制。该试验台可开展Φ2 m量级双旋翼模型的桨尖马赫数相似试验,是研究双旋翼构型高速直升机气动特性的重要试验设备。采用上、下旋翼模型独立支撑、单一电机驱动,实现了同步反向旋转、同步改变旋翼轴倾角、旋翼总距、旋翼周期变距等远程实时控制功能;通过设计计算-测试-调节的迭代,解决了试验台台体的动力学特性匹配问题;合理选择设计参数,解决了传动系统稳定性问题;利用该试验台获得的某旋翼模型风洞试验结果,拉力系数的重复性精度优于0.58%,扭矩系数的重复性精度优于0.11%。     

航天器刚性太阳翼在轨一维展开时间近似算法

针对航天器总体方案设计阶段迭代数值仿真太阳翼在轨展开时间计算复杂问题,以航天器刚性基板一维展开太阳翼为研究对象,将太阳翼的在轨展开运动简化为单自由度的刚性运动。利用动-势能守恒原理,推导了太阳翼在轨展开时间的理论公式。根据实际工程参数对太阳翼展开时间的理论公式进行了简化、拟合,建立了太阳翼展开时间的近似算法。近似算法解析计算的太阳翼在轨展开时间与数值仿真计算结果比对表明:相对误差小于10%,可以满足工程要求。用该近似算法,可在航天器总体方案设计阶段较为简单地计算并合理提出太阳翼在轨展开时间指标。     

共轴刚性旋翼桨毂减阻优化设计方法

针对共轴刚性旋翼桨毂阻力全机占比较大的特征,开展桨毂减阻优化设计方法研究。首先分析了共轴刚性旋翼桨毂外形特点,确定减阻整流罩基本型式,建立了减阻整流罩参数化模型。分析了整流罩参数对桨毂阻力的影响,进行参数样本选取。在此基础上建立代理模型,采用优化算法求解代理模型进行参数优化设计,最终建立了适用于共轴刚性旋翼桨毂的减阻优化设计方法。通过计算流体力学(Computational fluid dynamics,CFD)计算结果、试验结果以及优化结果对比验证,结果相差小于2%,说明了建立的共轴刚性旋翼桨毂减阻优化设计方法有效、可信。     

分离式共轴刚性旋翼风洞试验技术研究

为在■3.2 m风洞开展高速直升机共轴刚性旋翼风洞试验,中国空气动力研究与发展中心(China Aerodynamic Research and Development Center,CARDC)开展了分离式共轴刚性旋翼风洞试验关键技术研究。解决了双旋翼共轴反向驱动、双旋翼联动和差动变距操纵、旋翼气动力分离测量、旋翼间距调整等关键技术问题,研制了■2 m直径共轴刚性旋翼试验台和旋翼天平、旋翼操纵系统,实现了旋翼共轴等速反向旋转、双旋翼联动和差动变距操纵、旋翼升力偏置调整、旋翼气动力孤立测量、旋翼间距调整等功能。通过开展风洞试验,验证了该试验技术,表明试验技术具有技术成熟度高、数据重复性精度高、可调节参数多的优点。     

数字化协同测量在飞机弱刚性零件加工中的应用

针对飞机弱刚性蒙皮零件的高精度加工要求,通过机器人铣削系统协同标定、蒙皮边界提取与逆向建模、末端执行器定位修正等技术研究,实现飞机蒙皮、口盖、壁板类结构件边部余量检测及加工在线修正,边界测量精度0.1mm,末端执行器定位误差不大于0.2mm。经验证,本研究满足了四代机制造对数字化测量的要求,显著提高了飞机蒙皮壁板自动化铣削和数字化水平,促进了飞机制造从传统的手工装配向智能化装配模式的跨越。     

基于解耦算法的化学非平衡流并行数值模拟

化学非平衡流的数值模拟需要求解化学反应源项,而详细化学反应机理中各个基元反应的时间尺度相差显著,因而导致空间离散后形成刚性常微分方程组,其求解计算量巨大。为此,在非结构有限体积法化学非平衡流解耦算法的基础上,进一步采用基于MPI的并行计算实现大规模化学非平衡流数值模拟的加速,从而为进一步求解复杂的工程应用建立条件。由于基于解耦算法的求解器中计算量最大的化学反应源项求解是基于各个单元格心变量实现,每个时间步求解时不需要进行交互,因此便于实现分区并行计算。经典的激波诱导振荡燃烧算例的数值模拟表明,开发的并行求解器具有较高的计算精度和良好的计算效率,并且由于其较低的内存需求,具备进一步开发图形处理器(GPU)计算能力的基础。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。