一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面,设计了局部化的数据组织方式,节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法,该算法仅需少量的逻辑运算与代数运算,不仅可以快速跟踪模拟分子在网格之间的迁移,而且可以准确判别分子与物面是否相互作用,搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置,避免了重新计算。引入碰撞距离的思想,既减少存储又保证正确的模拟结果。在程序编制过程中,我们充分展现了Fortran90高级语言的主要特性,引入动态数组、指针、链表以及派生类型数据,编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

动态负载平衡的二维非结构网格DSMC并行算法研究

基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,研究了二维非结构网格DSMC并行算法。提出一类非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构造了符合DSMC并行原理的单步通讯法。采用单控制多数据流(SPMD)以及Master/Slave并行模式,设计了二维非结构网格DSMC整体并行算法。在程序的编制过程中,充分展现了Fortran90高级语言的主要特性,引人动态数组、指针、链表及派生类型数据。最后对过渡流域高超声速绕流进行了并行计算,数值试验的结果表明本文设计的并行算法可以取得较为理想的加速比。     

高超声速飞行器DSMC／EPSM自适应当地时间步长混合算法

采用优化的DSMC/EPSM混合算法,研究了高超声速飞行器三维复杂流场在近连续区的气动力热特性.发展了一种DSMC/EPSM自适应当地时间步长混合算法,根据网格内仿真分子的碰撞数,自动调整网格中进行仿真分子运动和碰撞计算的时间步长,有效地提高了程序的计算效率和计算精度.用DSMC方法、DSMC/EPSM混合算法和DSMC/EPSM自适应当地时间步长混合算法对比计算了三维复杂飞行器流场压力、热流分布量,飞行器表面气动力、热分布参数,证明了DSMC/EPSM自适应当地时间步长混合算法大幅度缩短了流场稳定所需的CPU时间,并且不降低流场模拟结果的精度.     

高超声速化学非平衡流动Navier-Stokes/DSMC耦合算法

基于相同的化学反应模型,在已有计算流体力学（CFD）和直接仿真蒙特卡罗（DSMC）方法及程序的基础上,采用Modular Particle-Continuum（MPC）耦合技术,建立了包含化学非平衡Navier-Stokes/DSMC耦合算法。算法结构中DSMC计算区域在CFD计算结果上根据当地克努森数自动选取。发展了适用于流场分区信息交换的亚松弛技术,抑制DSMC方法对CFD计算的影响。把DSMC方法和CFD的应用范围拓展到过渡流区,为复杂飞行器近连续过渡流区高超声速化学非平衡流动数值模拟研究提供了一种工程适用的预测分析手段。通过对二维圆柱高超声速化学非平衡绕流的算例与其他结果的比较研究,表明耦合算法不论在流场结构、流场非平衡现象,还是飞行器表面参数、整体气动力/热特性方面,都能够得到与全DSMC计算吻合的结果,证实了所建立的Navier-Stokes/DSMC耦合计算模型与方法的有效性和可靠性。仿真了某航天器解体碎片在过渡区的化学非平衡流动,得到碎片在过渡区的气动力/热特性,为碎片的陨落计算提供依据。     

一种基于自适应碰撞距离的DSMC虚拟子网格方法

为提高DSMC方法的计算效率及其对较高密度稀薄流问题的模拟能力,针对其碰撞对抽样环节,发展了基于自适应碰撞距离的虚拟子网格方法,并通过平板、钝锥等典型算例验证了该方法的正确性。结果表明:虚拟子网格技术放宽了DSMC方法对网格尺度的要求,可以在较大的网格尺度下得到满足计算精度要求的结果。     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术,提出适合于直升机旋翼运动的动态非结构重叠网格算法,该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格,不需任何人工干预。流动控制方程采用非定常可压缩Euler方程,空间离散采用Jameson有限体积法,时间推进采用双时间步长法,为加速收敛,在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性,开展一系列的数值实验,模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场,结果表明采用本文算法得到的结果是可信的,网格处理是自动高效的,可以应用于直升机全机流场模拟,为直升机设计提供可信指导数据。     

基于两级直角网格结构的三维DSMC算法研究

结合直角网格计算的高效率以及表面非结构网格对飞行器几何外形的精细描述,本文建立了基于两级直角坐标网格结构的DSMC数值模拟算法。一级网格为计算区域离散后的背景网格,二级网格为独立的碰撞网格和取样网格,是在背景网格的基础上根据流场性质的变化经过自适应后获得的亚网格,碰撞分子对在亚网格内选取,以保证选取的分子是"最相邻的"。计算的Orion飞船外形从自由分子流到稀薄过渡流的气动力系数,与DS3V和DAC软件计算结果的一致性非常好,验证了本文算法的可靠性和有效性。同时对过渡区的流动特征和表面气动参数的分布规律进行了分析。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

侧向喷流与飞行器三维绕流干扰流场特性的DSMC/EPSM数值模拟

采用DSMC/EPSM混合算法,研究了带有侧向喷流的复杂外形高超声速飞行器稀薄区和连续区三维混合流场的干扰特性。采用一种飞行器物面网格与DSMC计算域网格分别标识的方法,通过判断模拟分子与表面碰撞来完成飞行器物面网格与DSMC计算域网格间的信息传递和信息存贮,对于复杂外形飞行器的精确描述的物面网格不需做进一步处理,直接应用于不依赖于飞行器外形的DSMC计算的通用子程序中。引入碰撞数做为混合参数,对流动区域进行划分,并通过网格中碰撞数的统计和处理,有效地将DSMC和EPSM方法结合在一起。仿真计算了三维复杂外形飞行器带喷流流场压力、热流分布量,飞行器表面气动力、气动力矩和气动热参数,对喷流与高超声速空气绕流相互作用以及它们与物面之间的相互扰动进行了分析,证明了采用的方案和技术的有效性。     

基于非结构网格有限差分法的扎染算法

应用有限差分法时遵循贴体坐标系下离散等价方程的离散准则,格式中仅用到当地网格点的坐标变换系数,据此可以构建新的非结构网格有限差分法,数值算例表明在空间二维情况下用连接离散点的3条网格线构造出来的一阶迎风格式可以稳定收敛。把这种非结构网格有限差分法推广到三维,在4条网格线基础上进行离散计算,然后提出一种局部区域多重网格、重复计算的新算法。首先在包含物体的整个区域采用直角坐标系下的均匀网格进行计算,然后消除物体内部点对物体外部流场区域的影响,最后在物面和直角网格之间很小的局部区域填补非结构网格进行计算,计算过程类似于中国传统的扎染工艺。处理流场内存在曲面物体边界问题时,相较于常规结构网格差分算法,扎染算法尽管在计算中包含需要特殊处理的无用网格点,也存在重复计算的局部区域,但是以这些计算能力的浪费为代价,换取了编程简便、内存量小、网格生成快捷、易于扩展网格规模等优势。定性分析表明,这种扎染算法非常适合研制大规模并行计算,建立的计算流体力学应用软件可以充分发挥数十万至百万处理器核心的超级计算机的效能。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

用多重网格方法计算旋翼跨声速无粘流场

发展了一种加快悬停旋翼无粘流场计算收敛速度的多重网格方法。由于悬停旋翼流场中存在不可压区域，同时旋翼尾涡系统的发展需要较长的时间，使得旋翼流场的收敛速度远低于固定翼流场，因此研究旋翼流场的多重网格算法具有重要意义。空间离散格式采用了中心有限体积方法，时间推进应用了五步龙格-库塔法。采用3层网格的V循环，对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明：尽管多重网格方法对旋翼流场的加速收敛作用不如对固定翼流场的加速收敛效果，但是多重网格方法仍然可以显著地加快旋翼流场收敛。     

复杂外形下的无网格算法研究

研究了二维、三维复杂外形下的无网格算法。在一种布置点云方法的基础上,发展一种曲面拟合的重构方式构造流场物理量,并应用于AUSM+_up格式计算欧拉方程的数值通量;计算采用了隐式时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对某三段翼型低速流动、M6机翼跨音速流动、某全机跨音速流动进行了数值模拟,表明本文发展的无网格算法能有效地模拟复杂外形的无粘绕流。     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

Revisiting the space-time gradient method: A time-clocking perspective,high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternativ...     

首都机场终端区多进港队列排序优化策略与仿真

在现行可用空域中，通过流量管理策略与进港航班队列优化相结合可提高机场运行容量。以终端区航路航线间隔和扇区容量为约束，通过合理分配机场终端区边界的进港口流量．以总进场排队等待时间最短为目标，提出了协同调配流量和机动分流两种新策略，达到提高机场运行效率的目的。以首都机场为例，利用计算机仿真验证，结果表明：与现状运行相比较采用的优化策略排队进场可使进场延误和等待时间减少，并使离场正常率提高8％左右。     

快速声散射方法在变截面管道中的应用

 快速声散射方法(FSM)是一种基于无流动Helmholtz方程边值问题的声散射预测工具,具有快速、灵活的特点.以航空发动机消声短舱的声学设计为背景,用该方法对变截面管道声传播特性进行数值模拟研究,用直接边界元方法(DBEM)进行数值求解,避免了求解管口反射系数,有效地提高了计算速度.数值研究了刚性壁面和不同位置壁面声衬组合对变截面圆环管道形状声传播的影响,并分析了管道厚度对散射声场的影响,数值结果与声类比方法结果进行了相互验证.最后,还对一种真实转子声源进行了管道声散射的数值研究,结果表明本方法在航空发动机声学设计中具有工程应用价值.     

A New Single-blade Based Hybrid CFD Method for Hovering and Forward-flight Rotor Computation

A hybrid Euler/full potential/Lagrangian wake method,based on single-blade simulation,for predicting unsteady aerodynamic flow around helicopter rotors in hover and forward flight has been developed.In this method,an Euler solver is used to model the near wake evolution and transonic flow phenomena in the vicinity of the blade,and a full potential equation(FPE) is used to model the isentropic potential flow region far away from the rotor,while the wake effects of other blades and the far wake are incorporated into the flow solution as an induced inflow distribution using a Lagrangian based wake analysis.To further reduce the execution time,the computational fluid dynamics(CFD) solution and rotor wake analysis(including induced velocity up-date) are conducted parallelly,and a load balancing strategy is employed to account for the information exchange between two solvers.By the developed method,several hover and forward-flight cases on Caradonna-Tung and Helishape 7A rotors are per-formed.Good agreements of the loadings on blade surface with available measured data demonstrate the validation of the method.Also,the CPU time required for different computation runs is compared in the paper,and the results show that the pre-sent hybrid method is superior to conventional CFD method in time cost,and will be more efficient with the number of blades increasing.     

基于POD和DMD方法的跨声速抖振模态分析

跨声速抖振现象是由于非定常跨声速流动中激波的自激振荡而引起的结构强迫振荡,这种现象在跨声速飞行器中普遍存在,对飞机的结构强度和疲劳寿命有不利影响。基于模态分解的分析方法是进一步发展抖振控制手段的有效工具。本文通过两类典型模态分析方法（本征正交分解（POD）和动态模态分解（DMD））对OAT15A翼型的跨声速抖振现象进行分析,通过对模态频率、翼面压力分布、流场重构误差等方面的研究,将两种模态分解方法进行对比。发现基于频率特征的DMD方法能够准确捕捉抖振的临界稳定特征和抖振主频的典型模态,同时能够更准确反映流场变量在激波间断附近随时间的变化过程;而POD方法尽管在流场重构时具有较小的总体误差,但对激波附近压强随时间的变化历程拟合较差。