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

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

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

吸气式空天飞行器的一体化性能随扰动变化的敏感性高,在高马赫数飞行条件下,有必要开展流-固-推进耦合性能分析。针对机体/推进一体化布局的吸气式飞行器,明确一体化部件之间的耦合关系和耦合问题,利用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;背压增大,液滴粒径分布的均匀度指数减小,液滴尺寸分布更加不均匀。     

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

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

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

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

表贴式永磁同步电机直接转矩控制变角度预测控制

基于以电压矢量相角为变量的表贴式永磁同步电机(SPMSM)磁链和转矩方程,采用预测控制计算得出施加不同电压矢量相角的下一时刻磁链和转矩值,选择使磁链和转矩误差目标函数最小的电压矢量相角作为最优值,从而确定下一时刻施加的电压矢量。仿真结果表明:在电压矢量变角度预测下,SPMSM直接转矩控制系统运行良好,定子磁链轨迹为理想圆,磁链和转矩均符合控制要求,转速跟踪良好,定子电流波形为正弦。研究了不同等分电压矢量相角选择区域方式下变角度预测控制系统的控制效果。综合考虑控制效果和计算负担,变角度预测控制下将相角选择区域四等分比较理想。     

轴流压气机插板式进气畸变数值仿真

以航空发动机稳定性评价体系中常用的插板实验为背景,分别进行低速、高亚声速、跨声速轴流压气机与插板畸变发生器的耦合数值仿真,研究插板畸变条件下压气机的流场特性以及不同类型压气机在插板畸变影响下的失速起始机制。研究发现:转子进口截面的周向流动是影响低速压气机稳定边界的主要因素;插板角涡对转子叶顶泄漏涡的扰动是影响高亚声速压气机稳定边界的重要原因;畸变气流影响下的叶片通道内激波强度和位置的变化是影响跨声速压气机稳定边界的关键因素。      

一种基于方钢管骨架的壁板成型模具

针对大型飞机机身复合材料壁板试验件,从模具的刚度要求、使用要求、成本要求等方面确定了模具设计边界条件;为了满足模具设计的边界条件要求,提出了一种基于方钢管骨架的壁板成型模具。在进行模具设计时,首先对比各种常用模具材料的优缺点并最终确定模具采用“INVAR”钢材料;然后介绍了模具的结构布局和结构尺寸参数,并通过有限元模拟的方法,对模具的整体刚度和起吊工况的强度进行了模拟计算。计算结果表明模具的整体刚度较好,能够满足设计要求。从模具整体刚度和模具重量两方面与框架式隔板模具对比,方钢管骨架的壁板成型模具具有更大的优势。     

热压罐成型过程中框架式模具温度场模拟与分析

为了研究复合材料构件成型模具温度场,本文以某型飞机壁板的热压罐固化工艺为例,通过模型简化,利用FLUENT等仿真软件建立模具温度场的数值模拟模型,并将模拟数据与实验数据进行对比。结果表明二者平均相对误差为7. 4%。此外通过仿真模拟一组以支撑板厚度为变量的实例,两组以U、V两个方向的支撑板厚度为变量的对照组,通过判断模具型板表面温度方差大小,研究了支撑板厚度对模具温度场分布的影响规律。结果表明,模具温度场的均匀性随着支撑板厚度的增加而逐渐降低,其中U向支撑板厚度的变化对温度场均匀性的影响比V向大。