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

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

热空腔-喷气流的组合辐射

讨论了热空腔(喷气发动机的加力筒)和喷气流的组合辐射问题, 推演了组合辐射的辐射传输微分方程, 导出了有边界辐射条件下的解。讨论了由多个温度不同的壁面组成且腔内充满高温燃气的热空腔辐射传输方程解的具体形式, 导出了它向后半球的光谱辐射强度的表达式;同时还推导了热空腔-喷气流组合辐射的光谱辐射强度的表达式。最后给出算例, 算例结果表明, 热空腔-喷气流组合辐射的空间分布具有与实测分布较一致的梨形包络, 其光谱分布也与实测的光谱分布的一致性较好。      

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。     

结冰飞机着陆阶段飞行安全包线确定及操纵应对策略

结冰会导致飞行安全包线收缩、严重威胁飞行安全,研究结冰后安全包线变化对于操纵应对策略设计及提高飞行安全具有重要意义。以美国国家航空航天局（NASA）的项目飞机GTM（Generic Transport Model）为研究对象,通过对飞机的气动参数进行多项式拟合,建立了结冰飞机纵向通道的动力学模型。为了得到能随结冰程度变化的安全包线,将可达性分析理论引入到对结冰飞机着陆过程的安全性分析。提出将正向可达集与反向可达集的交集作为飞行安全包线,其中可达集的确定是基于水平集方法求解哈密尔顿-雅克比方程的最优解。最后针对不同程度的结冰条件进行了操纵时域验证,并提出了相应的操纵控制策略。研究结果表明,轻度结冰对安全包线影响较小,整个着陆过程的飞行状态始终能在最优控制指导下保持在安全包线以内;但对于重度结冰,飞行安全包线收缩严重,常规操纵已经很难使飞机达到着陆要求,需要进行飞行状态改出处理。研究结果为指导飞行操纵及实时包线保护打下基础。     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

径向磁气组合轴承系统的动态性能

为探究磁轴承失效时能否减轻转子跌落带来的损害,将人字槽径向动压气体轴承引入磁轴承转子系统,研究磁气组合轴承对系统动态性能的影响和动压气体轴承的支承特性。采用有限差分法和小扰动法求解气膜厚度方程和雷诺方程,研究动压气体轴承的静动态特性。对磁轴承电磁力进行了分析,采用磁轴承不完全微分PID（比例-积分-微分）控制策略,对系统进行了理论模态分析和试验模态分析,完成了系统高速旋转试验,测试了动压气体轴承在不同偏心率和转速下的承载力。结果表明,引入动压气体轴承可提高系统的动态性能,在磁轴承失效造成高速转子跌落瞬间,偏心率趋近于1,两个径向动压气体轴承能够产生较大承载力,减轻转子跌落造成的损害。      

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.