三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

An adaptive integration surface for predicting transonic rotor noise in hovering and forward flights

In this paper, a new permeable adaptive integration surface is developed in order to evaluate transonic rotor noise in accordance with FW-H_pds equations(Ffowcs Williams-Hawkings equations with penetrable data surface). Firstly, a nonlinear near-field solution is computed on the basis of Navier-Stokes equations, which is developed on moving-embedded grid methodology.The solution calculated through the present CFD method is used as the input for acoustic calculations by FW-H_pds equations. Then, two criteria for constructing integration surfaces are established based on the analysis of the quadrupole source strength and the nonlinear characteristic.A new surface is determined adaptively by the pressure gradient or density in a given flowfield,eschewing the uncertainties associated with determining cylinder-shaped integration surfaces. For varying hover cases, transonic noises are simulated with new integration surfaces for a UH-1 model rotor. Furthermore, numerical results of the new integration surface derived from the density perturbation value conform better to experimental data than results derived from the pressure gradient.Finally, the integration surface given by jrqj being 0.1, which is an applicable criterion obtained from hover cases, is used to predict transonic rotor noise in forward flight. The computational accuracy of the new integration surface method has been validated in predicting transonic rotor noise of an AH-1 model rotor at different advance ratios.     

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

分布式波浪前缘静子叶片对单级轴流风扇单音噪声影响的数值研究

为了适应逐年严苛的适航噪声标准，改善飞机噪声环境，突破航空发动机降噪技术的瓶颈，有必要开发新的降噪手段以指导叶轮机械降噪设计。本文采用URANS与FW-H方程混合方法，通过将波浪前缘构型运用在叶片流动损失较大、声源强度较强的部分后，进一步观察其气动性能和降噪效果。研究表明：与基准叶片相比，分布式波浪前缘静子叶片可以在1BPF（Blade Passing Frequency）时降低风扇入口声功率级0.67~1.9dB，2BPF时降低风扇入口声功率级1.87~4.18dB，3BPF时降低风扇入口声功率级2.4~6.8dB，同时，总压比最高提升0.027%，等熵效率最高提升0.28%。因此，分布式波浪前缘静子叶片在叶轮机械降噪方面有很好的运用前景。     

线圈炮——电磁同轴发射器及其系统分析

线圈炮在宇航与国防中有重要应用前景。本文先简述其分析模型与方程。方程为非线性变系数常微分方程的初值问题，有时呈刚性。用Ｔｒｅａｎｏｒ方法求解，所得结果与国外已有装置的试验数据符合良好。并对一个小炮做了初步的原理设计。     

利用透视投影基本方程选定透视形象的参数分析

在工程设计中，有些内容往往是工程设计人员根据一些不完整的历史资料萌发创作意念并概略地草构出透视表现图，虽然，在绘制标准的透视图时，这些资料所提供的参数是不够的，本文通过对物体三个主方向及其透视图中几何关系和矢量关系的分析，利用透视投影本方程由一已知物体的部分透视，确定其进行透视投影的全部参数，为透视变换和利用计算机绘制透视图提供了有效的基础条件。     

边界层稳定性的非平行性问题研究

采用抛物化稳定性方程的新方法研究稳定性的非平行问题，其控制方程从涡量形式的N－S方程导出，为提高计算精度和收敛速度，在法向采用高精度的高阶谱方法，结合边界层特征安排边界条件的方程，以及对无限区域数值问题作高效的代数变换处理，通过对预估校正迭代过程中推进步长有效的控制，以满足正规化条件，研究了不同类型压力梯度对稳定性的作用，由系列算例得到的增长率变化曲线和扰动中性曲线等，精确地给出了非平行性对流动稳定性的影响，计算的结果与相关数据符合甚好。     

波瓣混合结构三维流场数值计算

在贴体曲线坐标系上，采用有限控制体的非交错网格设置，对波瓣混合结构的三维混合流场进行了数值研究。经计算得到了波瓣混合结构气流混合的不同轴向截面上的速度分布，可以看出，由于波瓣尾缘处轴向旋涡的作用，主流和次流形成高效率的对流型混合，其轴向速度分布型面发生了明显的变化，气流间梯度较大的参数分布趋于均匀，这些特征与实验结果是一致的。     

求解大型非对称线性方程组的UNSYMMLQ方法

本文给出了求解大型非对称线性方程组的Lanczos方法的一个判据,提出了求解非对称方程组Ax=b的UNSYMMLQ方法,它是Paige和Saunders求解对称线性方程组的SYMMLQ方法的推广。文中描述并讨论了一些数值试验。     

不可压缩Navier－Stokes方程的基本解法

利用不同阶次的Ｈｅｌｍｈｏｔｚ算子方程的基本解，建立了时间离散不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的完全边界积分表示式。对俯仰谐振的ＮＡＣＡ００１２翼型所计算得的升力和阻力的滞后现象，与实验结果相一致