阵风与叶栅干涉噪声的数值模拟

采用计算气动声学方法求解二维线化欧拉方程对阵风与叶栅的干涉噪声进行了数值模拟。为减少频散和耗散误差,采用七点四阶频散相关保持有限差分格式进行空间离散;时间推进采用优化的4/6低频散低耗散龙格库塔法;在各边界处均选择或建立了适当的边界条件。首先通过模拟阵风在自由空间中的传播验证了无反射边界条件的可靠性;然后对阵风与平板叶栅的干涉问题进行了模拟,数值结果与精确解符合得很好;最后给出了阵风与NACA-0012翼型叶栅干涉的数值模拟结果,显示了数值求解方法的可行性。      

气动声学计算中一种声扰动方程的改进

研究了气动声学计算中的一种改进声扰动方程(IAPE)。采用三种不同的声传播模型,数值分析了声在均匀流以及剪切流中的传播问题。与线化欧拉方程(LEE)相比,除了个别位置预测的压力峰值较小以外,声扰动方程(APE)几乎能够预测到类似的物理特性;另外,APE方程具有无需求解密度方程,计算成本相对较少等优点。为弥补APE方程预测压力峰值较小的缺点,通过分析APE方程,这类方程中需要加入剪切流与声波相互作用的源项。为此,对APE方程进行改进,加入了一种声源项,计算结果表明,改进的APE方程(IAPE)能够和线化欧拉方程的计算结果保持一致。     

宽频声传播计算中多自由度时域阻抗边界条件的适定性

旨在验证一种多自由度时域阻抗边界条件在宽频声传播问题中的适用性与精度.选用NASA(美国国家航空航天局)流管实验数据进行数值校核;采用高精度计算气动声学方法数值求解线化欧拉方程.计算结果与实验数据及单频结果均吻合很好,表明了该多自由度时域阻抗边界条件具备精确处理宽频声传播问题的能力.      

横向喷流与超音速主流干扰流场的数值研究

 采用欧拉方程与可压缩涡对模型相结合的数值方法,研究了超音速主流与横向喷流干扰流场的气动特性。喷流横截面用一对反向旋转的可压缩涡来模拟。采用有限元Osher格式及非结构网格数值求解二维欧拉方程以得到非等熵外流场的气动特性及弓波位置。将欧拉方程数值解结果与可压缩涡对流场解结果相结合以研究喷流的轨迹及喷流横截面的变化。最后给出了数值模拟结果并与工程估算结果进行了比较,表明本文所述方法是可行的。     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

基于CFD/CSD耦合的跨声速非线性气动弹性研究

跨声速气动弹性数值数值模拟中,流场非线性特征明显,传统线化理论不再适用,需要采用时域方法,计算任意时刻的流场信息,进而得到正确的气动弹性现象。为AVICFD-X软件开发了气动弹性模块,该模块基于CFD/CSD(计算流体力学/计算结构力学)耦合模拟方法,包含流体、固体求解器以及交界面信息传输程序,在模块中分别求解流体控制方程和固体结果方程,通过界面进行流、固体之间的气动载荷以及运动状态等信息交换。AVICFD-X软件气动弹性模块数值模拟强迫振荡、静气弹、颤振等问题,验证气动弹性模块是否正确。数值模拟结果满足预期效果,与实验结果相吻合,充分说明了AVICFD-X软件气动弹性模块适用于跨声速非线性气动弹性研究。     

超燃冲压发动机凹腔稳焰的数值模拟

采用迎风格式数值模拟带凹腔(L/D=4)的超燃冲压发动机燃烧室内的高速可压缩流动,首先采用二阶NND格式对欧拉方程进行冷态数值模拟,以更清楚的捕捉、分析激波的特点;然后采用迎风三阶精度MUSCL格式求解二维N-S方程,湍流模型采用剪切修正的k-ω模型,对喷氢燃烧工况进行数值模拟。结果证明,凹腔可以达到掺混和稳定燃烧的目的,冷态时,总压损失较大,但喷氢燃烧降低了马赫数,减小了激波强度,实际工作状况凹腔不会引起很大的总压损失。     

Roe格式在水滴欧拉方程数值求解中的应用

目前普遍采用欧拉法求解水滴运动控制方程以模拟结冰表面水滴撞击特性。水滴欧拉方程近似没有压力项的空气欧拉方程,但已应用于空气欧拉方程的空间离散格式多数不能直接应用于水滴欧拉方程,有研究者将经典Roe格式应用于水滴欧拉方程时出现了奇点问题。通过对水滴欧拉方程进行形式改造,基于中心格式验证改造方法的有效性,并应用Roe格式结合线性重构与限制器函数实现对水滴欧拉方程二阶精度的求解。将NACA23012翼型计算结果与FENSAP-ICE软件进行对比,两者水滴收集系数分布基本一致。     

机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

绕侧滑锥体,三维机翼与机身的超音速非线性全位势流计算

本文在Grossman的非守恒全位势方程的方法中,提出以固定网格为基础的头激波安装法,节省计算机时。计算结果与实验和欧拉方程解符合很好,与线化理论结果的比较,显示了线化理论的不足。     

二维平行剪切层声波产生和辐射的数值模拟

 在二维线化Euler 方程的基础上解决了第3 届计算气动声学专题讨论会的第5 类标准问题: 二维平行剪切层声波的产生和辐射。频散相关保持有限差分格式用于空间离散, 低频散低耗散的龙格库塔法用于时间积分。重点使用了远场无反射边界条件。计算结果与NASA GLENN 研究中心的解析解符合得很好。     

亚、超音喷管中声传播的数值模拟

采用计算气动声学方法研究了管道近音速区的声传播和声波与激波的相互作用问题。空间离散采用高阶DRP格式,时间离散采用低耗散低频散的龙格库塔方法,边界条件采用Giles的特征变量一维非定常无反射边界条件。计算结果和解析解的精确吻合表明,计算中应用的这种高阶的DRP格式有很高的分辨率和很好的稳定性,适合声传播的计算。      

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

Investigation of gas dynamics numerical schemes based on the solution of one-dimensional test problems

The application of the optimized dispersion relation preserving scheme (DRP-scheme) in combination with the explicit optimized two-layer Runge-Kutta scheme is presented to solve a system of one-dimensional and quasi-one-dimensional Euler equations using as an example the solution of four test problems, namely, discontinuity disintegration in a tube (Sod’s problem); transfer of the lowamplitude Gaussian pulse; acoustic wave propagation through the transonic nozzle; acoustic wave-shock interaction. Also given are the comparison of the calculation results using different schemes: DRP, CABARET, CE-SE and the standard Lax-Wendroff schemes as well as the solutions obtained with the use of software packages.     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

翼型外形高气动效率/低可探测性的优化

 采用ｖａｎＬｅｅｒ矢通量分裂格式求解Ｅｕｌｅｒ方程的方法计算了绕翼型的气动特性;采用矢通量分裂方法计算了绕翼型的时域电磁散射场特性及雷达散射截面积（ＲＣＳ）;采用一种简单而有效的数值优化方法对流场解和电磁场解进行了翼型外形高气动效率／低可探测性的优化计算。算例结果表明,本方法提供了一种对翼型既可作气动优化设计亦可进行多学科综合设计的有效工具。     

WNND格式在稳态问题中的应用

基于二阶NND格式,引入Jiang和Shu的空间模板加权思想,构造了一种空间守恒变量型三阶精度的Weighted NND格式(WNND格式),通过组合波、球头以及Stardust返回舱绕流等多个算例对WNND格式进行了验证,数值结果表明,WNND格式引起的耗散和波动较小,并且能够高精度地分辨间断.在返回舱数值模拟中,引入了ADI隐式时间离散方法,加快了收敛速度,提高了计算效率.     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.