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

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

三点五阶优化广义紧致格式

根据两种差分格式色散误差目标函数对三点五阶广义紧致格式进行了优化。理论分析及数值实验证明了优化方法的正确性与有效性。该方法具有普遍意义，可以推广用来确定其它各阶精度的广义紧致格式的系数。并且利用类似的思想，可以通过其它限制条件对紧致格式进行优化。     

自适应的高精度中心﹣WE NO混合格式

为对可压缩流动进行高精度高分辨率的数值模拟,基于模板的光滑度量函数,构造了新型的加权因子,在此基础上将7阶精度的WENO格式与8阶中心格式进行加权组合,构造了一种自适应低耗散的中心-WENO混合格式(H-WENO7-CD8),并采用Fourier方法对离散格式数学特性进行了理论分析。该自适应的高精度中心-WENO混合格式相比于七阶的WENO格式具有更小的耗散误差。通过对激波-密度干扰问题、Rayleigh-Taylor不稳定性问题和双Mach反射问题的数值模拟,并与WENO-JS格式的计算结果对比,结果表明:该格式结合了WENO格式和中心格式的特点,能更好地捕捉激波,对物理脉动也具有较高的分辨率,适合于复杂流体流动的数值模拟。     

基于5阶精度格式WCNS-E-5的p-multigrid方法研究

p-multigrid方法的基本思想是:在保证收敛结果为高阶精度的同时,利用低阶精度格式耗散大的特点,来改善高精度有限差分格式在迭代计算时收敛速度慢的弱点.本文基于5阶精度WCNS-E-5差分格式,引入1阶精度迎风格式和3阶精度加权格式,构造了p-multigrid方法,在迭代过程中采用了V循环、W循环、S循环、PreV和FMG循环等不同方式来应用这三种格式,并通过典型算例考察了这些循环方式对收敛速度的影响,初步数值试验表明,采用恰当的循环方式,本文所设计的p-multigrid方法能够加快收敛速度,并保证了最终收敛结果与5阶精度WCNS-E-5差分格式的一致性.     

二维非结构网格上的高精度有限体积WENO格式

构造了非结构网格上二维双曲型守恒律的一类新的高精度有限体积WENO格式。其主要思想是:根据格式精度的要求,按照谱体积方法对三角形单元网格进行剖分,通过选取适当的子单元组成模板,利用WENO重构方法重构二阶和三阶多项式,利用有限体积公式和高阶Runge-Kutta TVD时间离散方法,构造了非结构网格上二维双曲型守恒律的一致二阶和三阶精度的有限体积WENO格式。然后,推广到二维Euler方程组。最后,给出几个数值算例,验证了格式的稳定性、高阶精度和高分辨捕捉激波等间断的能力。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用。与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要。因此,基于两步四阶时间离散和高精度WENO重构,发展了具有四阶时间精度的气体动理学格式。该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力。目前,两步四阶格式已经成功地应用到低雷诺数湍流直接数值模拟和高雷诺数工程湍流RANS模拟中,包括低速槽道湍流、超声速均匀各向同性衰减湍流、二维亚声速翼型湍流和三维跨声速翼身湍流等。数值结果表明该格式对湍流直接数值模拟和湍流RANS模拟具有高数值精度和高数值稳定性。下一步将利用高阶气体动理学格式研究更具有挑战性的可压缩湍流问题,例如超声速湍流边界层和激波边界相互作用等。     

高精度有限差分格式的色散优化及耗散控制

本文综述了我们在高精度有限差分格式的色散优化和耗散控制方面的研究进展。首先,我们提出了半离散有限差分格式色散和耗散相互独立的充分条件,实现了差分格式色散和耗散特性的独立调节。在此基础上,提出了色散最小、耗散可控的高精度差分格式,称为MDCD格式。MDCD格式已经得到了广泛应用,取得了很好的计算效果,但其主要缺点是耗散的调节依赖于经验。为了解决这一问题,我们进一步提出耗散的自适应调节方法。具有自适应耗散特性的高精度有限差分格式的基本特征是,差分格式的耗散能够随解的局部尺度自适应调节。为了构造这类格式,我们提出了一种新型的尺度识别器,它能够以等效无量纲波数的形式来定量衡量数值解的局部长度尺度。在此基础上,设计差分格式耗散参数与尺度识别器得到的等效无量纲波数之间的关系,从而构造了一类色散最小、耗散自适应的差分格式,称为MDAD格式。为了计算含有间断的问题,同时保持在光滑区的良好耗散特性,我们利用尺度识别器对一种经典的激波探测器进行改进,提出了一种新的激波探测器,并将自适应耗散格式与对应的WENO格式相混合,得到自适应耗散混合格式。近似色散关系显示该混合格式兼具高分辨率和鲁棒性。多个含间断流场的标准算例测试结果显示,自适应耗散混合格式具有良好的分辨率和激波捕捉能力。     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.      

气动声学直接数值模拟的空间差分格式分析

流场内声波的直接数值模拟对数值格式的耗散、色散特性提出了严格要求,首先介绍了两种适于声学模拟空间差分格式——紧致差分格式与色散关系保持(DRP)格式的基本思想,然后分析了格式的色散及其与四阶龙格-库塔法结合时的耗散特性。根据Fourier误差分析和数值模拟发现:无论从格式可有效求解的波数范围还是从与龙格-库塔时间推进结合后的耗散特性来看六阶紧致格式的可有效求解范围要大于DRP格式;不过,数值计算表明紧致格式的计算量远远超过DRP格式,对计算机资源提出了较高的要求。      

Burgers方程高精度差分格式分析

以非线性扩散波动方程Burgers方程为模型,通过调节方程右端粘性项的大小,在分析了方程的性质的基础上.探讨方程左端项和右端项差分离散格式匹配问题,以及高精度格式稳定性影响和计算格式精度与粘性系数和计算网格之间的关系等问题.研究结果将对方程左右端项离散时应如何采取合适的差分格式有一定的指导意义.     

一种求解双曲守恒律方程的新方法

近年来,在计算流体力学中,WENO(WeightedEssentiallyNon-oscillatory)格式得到了很大的发展,这也促进了基于WENO格式基础上的CWENO(CentralWENO)格式的发展。WENO格式的主要思想是通过凸组合的方式进行插值重构,以得到高阶逼近,并且在间断附近防止伪震荡的产生。本文构造了一种新的五阶紧凑CWENO格式,该格式所选模板节点少,并且保留了原有中心格式的优点。最后通过一些数值算例来验证该格式的高精度,本性无震荡性质。     

基于复合模型及FSQP算法的发动机性能寻优控制试验

提出了一种在保证优化效果前提下提高发动机在线性能寻优实时性的方法,即一方面基于相似理论建立了适用于非加力状态的发动机复合稳态模型,在优化过程中替代原部件级模型以缩短优化计算时间,另一方面采用先进的可行性序列二次规划算法寻求全局最优解,以提高优化精度。最终实现了上述方案在最小油耗模式及最大推力模式下的数字仿真及半物理模拟验证,相比于传统的基于发动机部件级模型进行性能寻优控制,其优化精度相当,但优化实时性得到大幅度提高,从仿真结果可以看出上述方案的有效性及可行性。     

点火药量对双脉冲固体火箭发动机点火过程影响

为了确定第二脉冲点火药量对双脉冲发动机点火过程的影响规律,以某双脉冲固体火箭发动机为模型,采用轴对称非定常Navier-Stokes方程和标准k-e湍流模型,对第二脉冲点火过程进行了仿真,研究了点火装置工作持续时间和点火燃气流量对于点火延迟的影响.研究结果表明:点火燃气流量相同时,增加点火装置工作持续时间可以有效减小第...     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method" and ‘‘minimum distance scheme of donor elements method" are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.     

基于SVR-PSO改进算法的航空发动机稳定性控制

在多变量发动机寻优控制中,用支持向量回归算法(SVR)对粒子群优化算法(PSO)进行改进可以有效避免局部最优解的出现.将改进算法应用于航空发动机实时稳定性控制,根据发动机仿真计算程序计算出发动机在各工作点处的稳定裕度,根据控制参数的变化域进行全局寻优,寻找满足压缩系统稳定裕度最小的工作点.仿真和分析表明:该算法实时性高,收敛速度快,具有较强的全局寻优能力,能在保证发动机稳定裕度最小的同时有效降低涡轮前温度和耗油率.      

基于DMOM算法的航空发动机性能寻优控制

提出一种分散迁移优化算法(DMOM),可实现多峰值优化问题的全局最优解搜索.该算法通过随机选择参考粒子,不断迁移搜索自身所处区域峰值点,再通过分散操作排除局部最优点,重新生成新个体,可快速搜索到全局最优区域.将DMOM应用于航空发动机性能寻优控制仿真,结果表明:在最小油耗和最低涡轮温度模式下, DMOM的寻优速度相比遗传算法(GA)和粒子群算法(PSO)提高了2倍以上;同时DMOM的优化精度相比自组织迁移算法(SOMA)提高了60%以上,相比可行性序列二次规划(FSQP)算法提高了20%以上.验证了DMOM相比其他优化算法有更强的跳出局部最优的能力,在航空发动机最小油耗和最低涡轮温度这类多峰值寻优问题中具有明显的优势.      

基于云自适应粒子群算法的高压转子复杂接触有限元模型参数修正

为建立更加准确的航空发动机高压转子的有限元模型,提出一种修正有限元模型描述航空发动机复杂接触的方法.将修正问题转化为求解定义在时域的误差函数的极小值,运用云自适应方法动态调整粒子群算法的惯性权重,使得算法在接近较优解时,惯性因子分布在云低端,有利于收敛得到更优解;当问题解较差时,其惯性因子分布在云顶端,有利于跳出局部极小点,扩大搜索范围.以仿真算例和实际航空发动机高压转子模型为例,通过与相关算法的修正结果比较,证明该算法是可行且有效的.      

Entropy based inverse design of aircraft mission success space in system-of-systems confrontation

This paper presents a practical and efficient design method for aircraft Mission Success Space (MSS) based on the entropy measurement (EM). First, fundamentals regarding MSS, Inverse Design (ID) and entropy are discussed. Then, two EM schemes of entropy-based ID and the whole MSS ID procedure are given to demonstrate alternative ways of entropy quantification and MSS design. After that, Genetic Algorithm (GA) is utilized as a search algorithm to find the optimal MSS design with the minimum objective, entropy, in each EM scheme. A simulation case of aircraft penetration mission is given for the method validation where the best aircraft MSS design is obtained by GA according to the minimum entropy. Results from two schemes are compared at the end.     

Structural optimization of elliptical-tip star propellant grain using sub-problem approximation

A method is presented here for structural optimization of elliptical-tip star grain.The grain structural integrity was improved by minimizing the most critical area of inner bore hoop strain during cool down.Optimization was done by sub-problem approximation method in conjunction with finite element analysis.Both radii of the ellipse were varied during optimization to find the optimal ellipse.The optimization resulted in grain geometry having minimum level of Inner bore hoop strain without violating the preset limits of burning perimeter.The von mises strain at grain inner bore was also reduced in resultant grain.      

航线网络中基于最小费用流的航线选择问题

针对航线网络中航线选择优化问题,分析了空中交通航线网络结构和空中交通流的特点,给出了空中交通航线网络的数学描述,考虑从求解最小费用流问题的角度寻求最优航线,建立了航线网络中求解最小费用流问题的数学模型。借鉴了连续时间网络最小费用流问题中结点到达费用的概念,并结合简单算例的求解验证了航线网络中用最小费用流方法选择航线的可行性,具有良好的计算性能。