新的多重网格法求解离心压气机内部流场

受英国剑桥大学 Denton教授对叶轮机械数值计算的启发, 本文提出了一种新的多重网格法, 并将之与显式时间推进法、有限体积差分格式、粘性体积力法配合使用来求解离心压气机内部流场。数值计算表明, 作者多重网格法的计算结果与试验结果吻合很好, 同时该多重网格法的收敛速度比不用多重网格法加快了将近两倍, 也比 Denton多重网格法的收敛速度加快了将近 1倍。      

多重网格法求解离心压气机内部流场

本文提出了一种新的多重网格法,并将之与显式时间推进法、有限体积差分格式、粘性体积力法配合使用求解离心压气机内部流场。数值计算表明,作者多重网格法的计算结果与试验结果吻合比较好,同时该多重网格法的收敛速度比不用多重网格法加快了将近两倍。     

非结构化多重网格流场数值计算

本文对二维非结构网格无粘流场计算的多重网格方法作了探讨。研究了非结构化多重网格的生成方法,插值算子的构造以及流场计算．以绕流ＮＡＣＡ００１２翼型无粘流场为算例,研究了多重网格法在非结构网格上的实现方法及其加速效果表明本文的方法有显著的加速收敛效果。     

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

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

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

两种多重网格法求解径流叶轮机械可压流动

为了加快数值求解叶轮机械内部流场的收敛速度,并改善计算结果,在数值求解径流叶轮机械可压流动的计算中,采用了由微元、二重和三重网格逐级转入和同时转入的两种多重网格法。对德国宇航局Ｋｒａｉｎ Ｈ设计的离心压气机内部流场的数值计算表明,两种多重网格法的计算结果差别不大,尽管在收敛速度方面存在差异,但均与测量数据吻合比较好。     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。     

基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

用多重网格法模拟三维流动

采用多重网格法模拟三维流动。控制方程在三维直角或圆柱坐标系中通过有限体积方法进行离散 ,SIMPLE算法对速度与压力进行耦合 ,亚松弛策略提高离散方程的收敛速度。在压力方程的求解过程中 ,界面速度应用由Rhie、Chow和Majumdar等人提出的动量方程插值法来剔除不合理的压力波动。三维流动算例表明采用多重网格法模拟三维流动可以得到良好的计算精度和收敛性。     

Improved aerodynamic design of turbomachinery bladings by numerical optimization

The aerodynamic optimization of a transonic compressor is reported in this paper. The Q3D Navier–Stokes solver COLIBRI is coupled to a gradient-based method (CONMIN) and to a genetic algorithm (GADO). The suction side of a 2D blade is optimized by using both optimization methods with a significant efficiency improvement. In 3D, the performance improvement is obtained by modifying the suction surface of a transonic compressor with a Bézier surface and by using the CANARI solver coupled to the gradient method (CONMIN).     

篦齿封严结构气弹稳定性数值分析

针对航空发动机中篦齿封严装置的气弹稳定性问题,用所发展的气弹稳定性分析的能量法,考查了两种不同的篦齿封严结构模型的颤振特性,对比分析了在不同节径振动下的气动功分布和模态气动阻尼比。结果表明:能量法可以用于预测篦齿封严结构的气弹稳定性,通过数值计算给出了改进前后篦齿封严结构低阶振动模态的气动功和模态气动阻尼比,改进后的篦齿结构在2节径振型下不再发生颤振,说明改进后的结构比改进前的气弹稳定性更好。      

Blewitt周跳探测方法的改进

在精密单点定位中,能否准确可靠地探测载波相位观测量中出现的周跳,直接影响到定位的精度及可靠性。本文概述了Blewitt周跳探测方法的基础原理,并对周跳的探测算法进行了改进,将电离层残差组合观测量一次差的变化是否超限作为周跳存在与否的判断依据。最后利用一段人为加入了几个典型周跳组合的双频观测数据进行验证,结果表明,改进后的算法能准确探测出两个载波上大于1周的任意周跳,提高了周跳探测的准确性。     

改善离子推力器束流均匀性的方法

离子推力器的束流分布,直接影响离子推力器离子光学系统（亦称栅极组件）的性能和寿命。通过对离子推力器离子光学系统的改进和放电室磁场的优化,使束流均匀性系数R值达到0.7左右。离子推力器束流分布均匀性的有效提高,有助于改善离子光学系统的受热分布,降低离子光学系统的温度,并能减小其温度差,使离子光学系统的热应力和热变形降低,进而延长离子推力器的寿命。     

Optimum detection and location estimation of target lines in the range-time space of a search radar

The average likelihood ratio detector is derived as the optimum detector for detecting a target line with unknown normal parameters in the range-time data space of a search radar, which is corrupted by Gaussian noise. The receiver operation characteristics of this optimum detector is derived to evaluate its performance improvement in comparison with the Hough detector, which uses the return signal of several successive scans to achieve a non-coherent integration improvement and get a better performance than the conventional detector. This comparison, which is done through analytic derivations and also through simulation results, shows that the average likelihood ratio detector has a better performance for different SNR values. This result is justified by showing the disadvantages of the Hough method, which are eliminated by the optimum detector. To have an estimate for the location of the detected target line in the optimum detection method as the Hough method, which detects and localizes the target lines simultaneously, we present the maximum a posteriori probability estimator. The estimation performance of the two methods is then compared and it is shown that the maximum a posteriori probability estimator localizes the detected target lines with a better performance in comparison with the Hough method.     

非线性偏最小二乘回归法在均匀设计响应面法中的应用

 针对目前多维变量可靠度分析中广泛应用的均匀设计响应面法(RSM),指出了使用最小二乘(LS)法拟合拟线性回归模型时存在的局限性,并提出采用拟线性偏最小二乘(PLS)法来回归响应面系数。由于拟线性回归法限制了模型的形式,精度提高有限,结果也很不稳定,因此又提出用基于样条变换的偏最小二乘回归模型代替拟线性回归模型并用于结构失效概率的计算,既能处理最小二乘法无法解决的变量间多重相关性的问题,又能避开拟线性回归中预先对模型形式的假定。通过算例验证了基于样条变换的偏最小二乘法的适用性和有效性,尤其对于多维变量非线性程度较高的可靠度分析,与普通最小二乘法拟合的响应面相比,其模型更加精确,失效概率精度更高。     

飞机性能参数预测的不确定性处理

 利用飞机的性能参数对飞机进行故障预报和状态监控是非常重要的。飞机的性能参数不仅具有非线性而且往往包含噪声,使得故障预测结果具有不确定性。针对这些问题,研究了利用非线性支持向量机处理飞机性能参数的预测问题,通过增加线性约束的方式解决了噪声带来的不确定性问题。此种方法不仅提高了预测的精度,而且模型可以利用适用于处理大规模二次规划的序列最小最优化算法进行求解,使得其可以解决大数据量的预测问题。利用仿真数据以及实际飞机性能参数对该方法进行了实验分析,实验结果表明此方法在精度上较不考虑噪声影响的模型有所提高,对于进一步提高飞机故障预测的精度,从而提高飞机的安全性具有重要意义。     

高精度圆锥倒角刀后刀面修磨方法改进

采用试切法对高精度圆锥倒角刀后刀面修磨,一次修磨合格率低于80%。通过针对使用试切法对倒角刀后刀面进行修磨的弊端进行分析,提出了改进方法。该方法通过理论计算获得万能分度头三个转角和刀具角度的关系,并具有规范的操作流程。对该方法进行了误差分析并进行了实验验证,证明通过该方法修磨的圆锥倒角刀锥角精度在±15′以内。相对于试切法,改进方法的锥角误差降低,稳定性大幅提高,操作时间大幅缩短。     

Frequency Agility Radar Range Calculation Using Number of Independent Pulses

A method for accurately calculating the radar range improvement introduced by frequency agility is presented. This calculation is made by using the concept of independent pulses. The number of independent target or clutter pulses is calculated by a simple method based on such basic combinatorial probability theory as drawing differently colored balls from an urn.     

Elimination of Blind Velocities of MTI Radar by Modulating the Interpulse Period

The performance of a coherent MTI is determined by two sets of parameters: the interpulse periods and the weighting coefficients. The latter are expressed in terms of the former so as to maximize clutter attenuation. Interpulse periods are optimized so as to remove blind velocities. The many local optima of the interpulse periods have improvement factors that show differences of up to 6 dB for a 5-pulse canceler. A computer program has been developed that combines random search and the gradient method to produce points that combine a high value of the improvement factor with good blind speeds performance.