塞式喷管性能的数值模拟与实验验证

通过特征线法在塞式喷管中的应用，研究了塞式喷管主要结构参数对其性能的影响，同时，研究了塞式喷管的高度特性，最后还对塞式喷管合推力与发动机轴向夹角的高度特性进行了研究和分析。理论研究结果得到了实验的验证。其结果可用于塞式喷管的设计。     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

一种适用于大幅值Quasi周期轨道的数值构造方法

针对现有方法在构造大幅值quasi周期轨道时收敛性和通用性不足等问题,提出一种适用于大幅值quasi周期轨道的数值构造方法。首先针对2维不变环面的计算问题,设计了映射间隔约束的环面多步微分修正算法。在此基础上,采用自然参数延拓法和伪弧长延拓法两种不同的思路构造等映射间隔的quasi周期轨道族。接着以quasi-Vertical为对象分析了两种思路各自的特点。最后以quasi-Axial, quasi-DRO等为对象验证了方法对三角平动点及多对复特征值情况的适用性。仿真结果表明,所提出的方法能够简单有效地解决多种类型幅值比超过1.2的大幅值quasi周期轨道的数值计算问题。     

密集型空气雾化流场破碎特征数学模型研究

低速液体射流在高速湍流气体作用下的气液同轴射流雾化流场是瞬态密集型喷雾场,高速湍流中影响和控制雾化的因素很多。针对空气雾化流场液核分布与特征进行分析,开发了模拟液核的随机浸入体模型,结合大涡模拟方法对同轴射流空气雾化喷嘴下游流场进行数值模拟。模拟结果与实验结果对比表明：随机浸入体模型可快速捕捉液核的长度和位置信息。当气液动量比M为3~10000时,能够较准确地预测液核长度,当M>10时,其预测结果远优于唯象模型;该模型能够捕捉回流区、大尺度涡等流场结构。同时,可以准确地预测喷嘴附近的液滴粒径,特别是在气流速度较大（>60m/s）时,液滴平均直径预测误差<10%。     

供油瞬态喷油杆内燃油结焦机理及规律数值研究

为获得喷油杆通断循环供油工作模式下,供油瞬间燃油结焦速率陡增机理,本文建立起沿喷杆轴向的一维非稳态热-流-固耦合换热与燃油结焦计算模型。获得不同时刻壁温、油温与结焦速率等参数轴向分布。通过与连续供油模式下的结果进行对比,来分析供油瞬间灼热喷杆内壁对燃油热冲击的影响,以获得结焦速率陡增机理。结果表明通断供油模式下,供油瞬间会迅速形成一尖锐的油温和结焦前体浓度峰值,结焦固相表面反应被迅速激发,结焦速率陡增。在喷杆前端（x=100mm）与底部（x=260mm）附近各存在一结焦速率峰值区域。通断供油模式下,在x=100mm和x=260mm位置处,每个周期供油起始时刻结焦速率峰值分别在5000μg/cm2h和2000μg/cm2h左右,远大于连续供油下的值。但供油结束时刻,供油模式的影响大幅下降。此外通断供油下,供油瞬间热冲击对燃油结焦的影响沿喷杆轴向逐渐减小。减小喷杆供油流量,热冲击的影响有所下降,且随时间先升高后降低。提高进口油温,热冲击的影响显著减小,且随供油时间逐渐降低。     

基于蒙特卡罗法的航空发动机空气系统稳态算法优化

针对目前航空发动机空气系统稳态算法中收敛性依赖初值的问题,将蒙特卡罗方法与流体网络法综合应用到空气系统可压缩流体一维网络计算中,提出了一种新的计算方法Monte Carlo-Fluid Network（MC-FN）。该方法将空气系统简化为由节点和元件组成的网络,借助蒙特卡罗方法获得空气系统内各节点压力分配,再根据空气系统中各元件流阻特性和换热特性计算流量、温度。计算中通过将游动次数比较少的蒙特卡罗方法的计算结果作为流量残差法节点压力、温度的初始值,实现快速求得精确收敛解。与流量残差算法相比,MC-FN方法计算精度不变,收敛速度提升了66.5%;与线性求解法相比,MC-FN方法的计算精度提升了25.2%,收敛速度提升了43.8%。     

有声腔推力室的声学振型及其阻尼特性研究

为确定多声学模态压力振荡条件下带有1/4波长管声腔推力室的声学振型及其阻尼特性,揭示1/4波长管声腔对推力室压力振荡的抑制作用机理,对有声腔推力室和无声腔推力室进行近圆周壁面的定容弹激励仿真,激发了多模态的声学振型,给出了推力室压力分布的时空演化,并采用半带宽法定量评价每个激发声学振型的阻尼特性。结果表明：1/4波长管声腔成功抑制了目标振型（一阶切向振型）压力振荡,大幅度减小其幅值,而不是大幅度增加其半带宽,但可能增强其它非目标振型的压力振荡。声腔通过削弱目标振型波峰波谷压力差、声腔入口漩涡和增加壁面面积等三种方式来抑制目标振型压力振荡,主要以前者为主,前者旨在降低目标振型幅值,后两者是增加其半带宽。     

基于有序BDL树集的船舶电网潮流算法及执行时间分析

针对嵌入式船舶电力模拟训练系统在频繁模拟操作指令下对潮流计算的快速性需求以及嵌入式模拟训练系统极其有限的运算资源之间的矛盾,提出一种基于有序BDL树集的船舶电网潮流算法。该算法在常规回推前推法基础上对节点编号进行优化,将多电站船舶电力系统在实际开环运行模式下等效为由多个辐射状电网构成的有序BDL树集,树集中的每个有序BDL树包含单个电站中的主配电板层、配电层和负载层,并把发电机潮流整合入主配电板层,避免了常规潮流算法中节点编号突变问题。对该算法的时间复杂度分析结果及实际运行结果进行比较分析,证明该算法比常规船舶电网潮流算法更优。     

基于VBA的Maxwell二次开发在无刷交流励磁机电磁计算中的应用

针对无刷交流励磁机复杂的结构特点,将Maxwell 16.0作为开发平台,使用Visual Basic for Application(VBA)程序开发语言,对该电机的建模、添加材料、施加激励、剖分以及后处理进行二次开发。对无刷交流励磁机建模中的定子绕组、磁极、转子槽、电枢绕组实现参数化,具有一键求解、自动剖分和后处理的自动出图等功能。基于Excel的开发平台,设计可视化界面,加强人机的信息交互,减少重复性的步骤,提高设计人员的效率,缩短电机的设计周期。     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。