一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

基于纵向滤波的星敏感器低频误差在线估计

多敏感器数据融合是获得更高精度姿态测量的有效方法,敏感器数据融合前必须先修正低频误差。首先,介绍了星敏感器低频误差（LFE）的产生机理及对其在线估计的必要性。其次,针对传统算法的不足,提出了基于纵向滤波的低频误差在线估计算法,该算法将传统低频误差估计问题转化为若干个常值误差估计问题,提高了估计精度。最后,给出了该算法具体实施方式,说明相关参数物理意义及选取原则。通过理论分析及仿真,算法误差可忽略不计。通过在轨数据仿真,星敏感器轨道周期低频误差可被消除。     

流-固-推进耦合下的机体/推进一体化性能分析

吸气式空天飞行器的一体化性能随扰动变化的敏感性高,在高马赫数飞行条件下,有必要开展流-固-推进耦合性能分析。针对机体/推进一体化布局的吸气式飞行器,明确一体化部件之间的耦合关系和耦合问题,利用CFD、有限元和准一维流方法,结合本征正交分解(Proper Orthogonal Decomposition,POD)降阶手段,建立吸气式空天飞行器流-固-推进多物理场耦合快速分析方法,并开展多场耦合特性分析。结果表明:(1)进气道压缩面的流-固耦合导致出口静压的最大振荡振幅约为平均静压的21.6%,而出口马赫数的最大振幅约为平均马赫数的8.45%。(2)进气道出口性能的振荡会影响发动机的推力性能,导致推力振荡幅值可达平均值的31%,且随着时间的推移,会在进气道外压缩流场产生大量的气动涡,涡结构进入进气道后会导致进气道出口性能的持续下降,进一步削弱了发动机的平均推力性能。     

阵列式表面电弧等离子体气动激励控制三角翼流动分离的实验

为探索多路阵列式微秒脉冲表面电弧放电(μs-SAD,Microsecond pulse surface arc discharge)对尖前缘小后掠三角翼流动分离的控制效果和作用机理,首先通过放电测试和纹影测试对多路阵列式μs-SAD的激励特性进行研究,揭示其对流场的作用原理,进一步将多路阵列式μs-SAD用于三角翼流动控制,开展了小后掠三角翼流动分离控制低速风洞实验,研究了来流速度、激励电压和激励频率等参数对控制效果的影响规律。结果表明:多路阵列式μs-SAD能够快速放热,单路瞬间放电能量可达68mJ,在流场局部可诱导产生冲击波;机翼前缘多路阵列式μs-SAD能有效改善三角翼大迎角气动特性,当来流速度为30m/s时,使最大升力系数提高27.2%,失速迎角推迟4°;来流速度增大到40m/s时,流动控制效果减弱,使最大升力系数提高15.5%;存在最佳激励频率使无量纲频率F+=1时,控制效果最好;激励电压存在阈值,其随来流速度的增加而增大,当激励电压超过阈值电压继续增大时,流动控制效果不再增强。     

背压对撞击式喷嘴雾化特性影响研究

为研究背压对撞击式喷嘴雾化特性的影响,将压力的变化等效为气体密度的变化,基于一种树形自适应加密算法,通过直接数值求解不可压Navier-Stokes方程组实现了不同背压条件下射流撞击雾化的数值模拟。首先将数值模拟结果与试验数据进行对比,验证数值模拟的有效性,在此基础上开展了高背压条件下雾化过程的数值模拟。结果表明,随着背压的提高,气动力相应增强,液膜的破碎更加剧烈,一次雾化区域的液滴数密度增大,雾场由稀疏向稠密发展;液膜在向下游运动过程中波动速度的幅值逐渐增大,并且随着背压的提高,液膜的波动由线性向非线性转变;背压增大导致破碎长度减小,经过参数修正得到了液膜破碎长度的经验公式,并与试验数据进行了对比;背压对液滴尺寸分布规律没有显著影响,但随着背压的提高,同一时刻大液滴所占的比例提高,整个雾场的Sauter平均直径有增加的趋势,当背压从0.1MPa增大到1MPa时,雾场的Sauter平均直径由155.5μm增大到166.9μm;背压增大,液滴粒径分布的均匀度指数减小,液滴尺寸分布更加不均匀。     

翼身融合运输机分布式电推进系统设计及油耗评估

针对翼身融合运输机开展了分布式电推进系统的总体设计与油耗评估。通过数值计算完成了70t载质量翼身融合飞机的气动设计与优化。在巡航马赫数为0.80和10km高度的设计点,最大升阻比达到了24。通过求解积分边界层方程组,完成了电推进系统的总体设计。电推进系统包含10个推进风扇,风扇直径为1.45m,压比为1.35,巡航功率为2.94MW。建立了考虑燃烧过程的发动机一维性能模型,对发动机油耗进行了评估,获得了不同发动机循环参数下燃油消耗。建模结果表明,基于翼身融合布局和分布式电推进技术,可使运输机的油耗较C-17节省近50%。      

反推装置对翼身融合分布式推进系统气动性能影响数值模拟

为了研究新型反推装置对翼身融合分布式推进系统气动性能的影响,将机身简化为二维模型,利用数值模拟的手段来研究攻角、反推开度,以及风扇打开或关闭对飞行参数的影响,为翼身融合分布式推进系统的反推装置设计提供初步的意见和建议。结果显示:反推闭合和风扇打开之后该机型的升力系数和俯仰力矩系数都会有较大增加;随着反推开度的增加,升力系数和俯仰力矩系数也随之增加;来流攻角对翼身融合分布式推进系统的气动性能也有较大影响。      

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符;反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。      

表贴式永磁同步电机直接转矩控制变幅值预测控制研究

基于以电压矢量幅值为变量的表贴式永磁同步电机(SPMSM)磁链和转矩方程,采用预测控制计算得出施加不同电压矢量幅值的下一时刻磁链和转矩值,选择使磁链和转矩误差目标函数最小的电压矢量幅值作为最优值,从而确定下一时刻施加的电压矢量。仿真结果表明:在电压矢量变幅值预测下,SPMSM直接转矩控制(DTC)系统运行良好,定子磁链轨迹为理想圆,磁链和转矩均符合控制要求,转速跟踪良好,定子电流波形正弦。预测控制需要从备选电压矢量幅值集合中选择最优幅值。理论上,备选电压矢量幅值个数越多,系统优化效果越好,但也带来更大的计算负担。研究了9种等分电压矢量情况下变幅值预测控制系统的控制效果。根据控制效果和计算负担,提出变幅值预测控制将电压矢量幅值三等分较为理想。     

翼身融合布局客机总体参数分析与优化

在翼身融合布局客机总体设计阶段,为评估设计方案的总体性能,建立了翼身融合布局客机总体参数综合分析与优化平台,该平台以翼身融合布局客机的几何参数为输入,完成动力、几何、重量、气动、性能和经济性等模块分析,并以此为基础建立优化设计模型。为快速评估设计方案的性能及优化设计效果,动力分析模块采用了部件级分析模型,重量分析模块采用半经验估算方法,气动分析模块采用面元法结合工程估算方法,性能分析模块采用简化运动学方法,优化模型采用可并行计算的子集模拟优化算法。以某555座级翼身融合布局客机方案为例,应用开发的分析与优化平台,完成了总体参数分析,结果表明分析模型合理。在此方案的基础上,以客机的外形参数和发动机海平面最大推力为设计变量,分别建立了以最大起飞重量最小为目标的单目标优化,以及同时以直接使用成本和进场速度最小为目标的多目标优化,单目标优化结果最大起飞重量降低了约7.17%,多目标优化结果表明直接使用成本降低8.77%的同时进场速度会增加3.32%。