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

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

一种适用于压气机非定常流动的相位滞后算法

应用相位滞后边界条件,以1∶1的转、静叶计算域,在不改变转子尾迹及静子势干扰的频率的条件下,计算叶轮机非定常流场.以NASAstage₆₇为算例,计算结果与实验数据相吻合,表明采用的方法具有较大的工程实用性.      

陶瓷基复合材料高精度宏细观统一本构模型研究

基于高精度通用单胞模型(HFGMC)建立了陶瓷基复合材料(CMCs)的宏细观统一本构模型.首先基于CMCs的细观结构和损伤特点建立代表体元(RVE)并定义细观损伤变量.然后采用HFGMC模型计算细观损伤变量的演化,由此获得宏观应力应变曲线,从而建立了CMCs的宏细观统一本构模型.采用该本构模型预测了CMCs的宏观应力应变曲线,与实验结果吻合良好.结果表明该本构模型能够模拟CMCs的拉伸与损伤行文,为CMCs强度分析提供了有效手段.      

3种铝合金材料动态性能及其温度相关性对比研究

 对3种铝合金2024-T351、7050-T7451和LY12-CZ进行了温度在77～573 K的静、动态压缩（应变率10^-3～6 000/s）和拉伸（应变率10^-3～3 000/s）试验,分别得到了3种铝合金材料的应力 应变关系和失效应变。对比研究结果表明,随温度升高3种材料的塑性流动应力降低,应变率敏感性增加。最后基于Johnson-Cook模型,拟合了用以预测铝合金材料塑性流动应力的模型参数,其预测结果与试验结果吻合较好。     

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

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%。因此，分布式波浪前缘静子叶片在叶轮机械降噪方面有很好的运用前景。     

晶须增韧陶瓷的微裂纹演化及本构描述

本文利用作者改进的等效夹杂理论研究了含随机分布晶须和微裂纹陶瓷的微裂纹演化及由其引起的材料本构非线性，定量计算了热残余应变、晶须含量及尺寸比的影响，对材料的微裂纹演化及破坏过程进行了数值模拟，得到了一些有意义的结论     

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

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