非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

NND格式的推广及在粘流计算中的应用

本文把ＮＮＤ格式推广到了特征变量，守恒变量和原始变量形式。通过引入通量差分分裂方法，由守恒变量形式的ＮＮＤ格式得到了Ｒｏｅ的近拟Ｒｉｅｍａｎｎ解法，并讨论了限制器的作用。在隐式ＮＮＤ格式中引入了一种分解算法，提高了计算效率。本文的数值计算表明，ＮＮＤ格式对复杂波系干扰，边界层分离与再附等流动现象均有较高的分辨率。     

计算气动声学中的紧致滤波格式

从理论上分析了数值伪波产生的原因,要消除数值伪波可用滤波方法。为了使滤波器具有理想的截断特性,引入了修正的高斯函数。通过让滤波器的频响函数与修正的高斯函数逼近,利用序列二次规划(SQP)方法优化了五对角紧致滤波格式。优化的滤波格式提高了计算精度和效率,增强了数值稳定性,更易于实施。一维和二维算例体现了优化滤波格式的性能改进。     

结构动力学方程的显式与隐式数值计算

从理论上详细推导了动力学方程的显式隐式时间积分数值计算方法,针对经典中心差分格式在求解阻尼结构计算效率降低的问题,研究提出了一种新的显式积分格式,理论证明该数值积分格式能够处理各种复杂阻尼状态而不降低计算效率。对于不同的时间步长通过设定调节加速度参数,能得到比经典显式计算格式和NEWMARK逐步时间积分更高的精度,最后通过数值算例对三种计算格式的理论分析进行了验证。     

一种改进的紧致WENO混合格式

给出了一种用于激波捕捉计算的守恒型混合紧致WENO格式.混合格式可视为五阶WENO格式和五阶守恒紧致格式的加权平均.在格式构造中,本文采用了间断分辨率更好的WENO权因子计算方法以及新的子格式权因子计算方法.对于Euler方程,仅对混合格式中WENO格式部分做特征投影处理,不仅获得了特征型WENO格式良好的间断分辨率,同时又保证了计算效率.数值实验验证了本文混合格式的高时间效率和高流场分辨率.     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

某型7级轴流压气机中间级三维气动数值优化

针对多级压气机三维气动数值优化中存在的级间匹配与计算成本难以兼顾的问题,采用了只对待优化级进行单级三维数值求解,而在边界条件中考虑级环境影响的优化策略,对某型7级高负荷轴流压气机存在叶背分离的第4级静子进行了数值优化.在有限的计算时间内,多级压气机绝热效率相比原始设计提高了近0.8%,压比增加了近4%,压气机的整体性能得到提高.应用表明:该优化策略较好地平衡了级间匹配与计算成本的矛盾.      

采用二阶迎风格式的 PISO 算法研究

在ＰＩＳＯ算法基础上,发展了一种新的二阶迎风格式以提高计算精度又保证算法的稳定性。新格式应用于液体火箭发动机燃烧流动过程计算,其数值稳定性很好。     

一种空心杯直流测速发电机的设计分析与方案改进

空心杯绕组是空心杯直流测速发电机的核心组件,由于原样机绕组匝数多、线径细,现有绕线设备无法进行绕制,为了从工艺上实现绕线设备的利用,加速研制进度,在保证空心杯直流测速发电机输出特性不变的前提下,通过建立二维场路耦合有限元模型,对原方案进行结构尺寸和绕组参数改进,计算并比较了原方案和改进方案的磁场分布、输出电压和输出斜率、峰峰值纹波系数,分析了电刷安装偏移的影响,通过对输出特性进行测试,证明了有限元计算方法的正确性,改进方案的性能指标与原方案一致,满足设计指标的要求.     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

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

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

预估-校正LU-SGS的隐式算法

为了进一步提高Euler方程求解效率,在传统LU-SGS格式基础上提出了预估-校正LU-SGS。结合高分辨率迎风通量格式,发展了一套高效的Euler方程求解器。以NACA0012翼型和RAE2822翼型跨音速无粘流动作为算例,对比分析了预估-校正LU-SGS、传统LU-SGS和GMRES格式的计算效率。结果表明:预估-校正LU-SGS算法能显著提高求解效率,其效率接近GMRES算法,比传统LU-SGS方法提高了近三倍。     

高超弯曲激波二维进气道初步研究

研究了一种压升规律的曲面压缩面,设计了高超弯曲激波二维进气道,并用数值模拟手段对该进气道和同等条件下的常规高超二维三楔进气道、楔+等熵进气道的性能进行了比较.数值模拟表明:通过给定合理压缩面压升规律来设计压缩面并改善压缩面附面层稳定性是可行的,弯曲激波二维进气道的长度比同等条件下的常规二维三楔、楔+等熵进气道分别缩短12%和10%,并且对来流 Ma 变化不敏感,综合性能优势明显、应用前景大.      

运动激波绕射的数值模拟

用高分辨率的ＴＶＤ格式对复杂的激波反射过程做数值模拟。激波斜壁反射计算结果与实验干涉照片、压力分布、密度分布均符合得很好。通过对运动激波绕射凸起物的数值模拟,得到非定常反射流场发展过程及复杂流场间断的清晰图谱。     

航空发动机整机结构效率评估参数与计算方法

根据航空发动机整机结构系统功能和设计要求,对整机结构效率内涵进行了详细阐述,基于整机载荷特征,建立了从抗变形能力和力学环境适应能力两方面对整机进行结构效率评估的参数体系.通过对转静间隙变化量、应变能分布系数、敏感度系数等评估参数归一化方法及结构效率系数运算法则的建立,明确了结构效率的分析流程.以某小涵道比涡扇发动机为对象进行结构效率评估,详细给出了结构效率的计算过程.针对该小涵道比涡扇发动机的另一种结构方案2,基于结构效率评估方法进行评估并与原结构方案进行了对比分析,结果表明:原结构方案在抗变形能力和力学环境适应能力两方面均优于结构方案2,结构效率系数比结构方案2高出0.11,初步验证了结构效率评估方法的工程实用性.      

直升机鲁棒控制器的设计研究

提出一种改进直升机悬停性能的控制方案。在直升机悬停控制器原速度回路的基础上增加鲁棒控制器,可提高直升机悬停参数变化的鲁棒性能。详细叙述了鲁棒控制器结合悬停控制器的设计方法,并用仿真结果说明了鲁棒控制器在本算例中所起的作用     

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。     

大功率九相变频异步电机通风系统热流耦合仿真分析

为了有效解决大功率低速异步电机通风散热难题,基于热流耦合仿真分析计算方法,以1台2 800 kW、8极九相变频异步电机为例,系统分析了顶吸式强迫风冷电机内定子齿压板宽度距离对电机内风量分配、绕组温升分布以及内部风阻特性的影响,并筛选了匹配方案;同时,对比分析了九相减薄绝缘结构与三相高压绝缘结构的换热性能。最后,进行了电机的温升试验测试,温度场计算值与试验结果基本吻合,验证了基于热流耦合仿真分析计算的准确性和有效性。