战术导弹绕流流场并行计算

本文采用求解NS方程的方法对战术导弹的流场进行了数值模拟.应用区域分割技术和拼接网格的并行策略,发展了一种适合于分布式存贮并行环境的隐式有限体积并行算法,在PVM并行环境下,对战术导弹绕流流场实现了单块网格多机并行计算和多块网格多机并行计算,通过负载平衡等方法得到了较高的加速比.计算结果表明并行计算结果与串行计算结果完全一致.     

多块多重网格法及其跨声速转子内流并行模拟

为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。     

三维分区对接网格Euler方程的并行数值模拟

本文根据CFD中分区对接网格的特点，应用基于工作站集群系统的MPI网络并行环境，实现了典型双三角翼模型亚声速绕流Euler方程的并行计算。文中对计算结果的收敛情况、耗时、加速比及效率等并行性能指标进行了讨论。结果表明，多域多机与单域单机的计算结果基本吻合；2节点并行效率达到88．21％，4节点并行效率达到78．99％。     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

基于图剖分的多块结构网格负载平衡方法

负载平衡是影响并行计算性能的重要因素。针对多块结构网格,给出了一种改进的多层次图剖分负载平衡方法。该方法设计了新的网格剖分算法,采用改进的子块分裂方法与图剖分算法的循环调用实现结构对接网格剖分,并通过建立不同物体重叠网格间的连接关系,实现了结构重叠网格的负载平衡。采用2个典型算例对方法进行了对比验证,数值结果表明,子块分裂方法对剖分结果具有重要影响,采用循环调用算法及改进的子块分裂方法能有效地实现计算负载均衡及通信量优化,同时显著减少了网格块数及因虚网格导致的内存需求,有利于提高并行效率。该负载平衡方法与网格拓扑无关,适用于多块结构对接网格及重叠网格,且整体型剖分方式对于多块结构重叠网格具有更好的剖分效果。     

GPU平台上的叶轮机械CFD加速计算

通过数据并行的方式对一个成熟的叶轮机多块网格气动计算程序(MAP)进行了并行化处理,利用计算统一设备架构(CUDA)技术实现了在图形处理单元(GPU)上的并行计算.保留了原程序中的2阶空间迎风格式和隐式时间离散格式,并采用了隐式迭代对线性系统进行求解.经过2个叶轮机械算例的测试,与在传统的中央处理器(CPU)上运行的原程序相比,在计算结果完全一致的前提下,单GPU的计算速度最高可达单CPU计算速度的8.89倍,与四核并行的CPU计算相比可以得到2.39倍的加速.     

非结构网格多区域动态加密及分布式并行计算

针对分布式并行系统,提出一种二维非结构网格多区域动态并行策略,可在计算过程中对各计算区域的网格进行局部自适应网格加密,可在不过多消耗计算机资源的情况下,快速,有效地实现多台同构计算机上的计算网格的动态负载平衡.并针对二维Euler方程,利用已有的高精度MmB有限体积格式,应用提出的并行计算策略,进行了多个数值试验,得到令人满意的效果.     

多GPU并行可压缩流求解器及其性能分析

为实现可压缩流问题的大规模高效数值求解,开展基于图形处理单元（GPU）的并行计算研究。在NVIDIA GTX 1070上建立了基于消息传递接口+统一计算设备架构（MPI+CUDA）的多GPU并行可压缩流求解器,该求解器基于结构网格有限体积法,空间离散采用AUSM⁺UP格式。采用一维区域分解法对计算网格进行划分,使得各GPU之间达到负载平衡。针对超声速进气道算例,对算法单GPU并行性能和多GPU可扩展性能进行分析。数值结果显示,单GPU并行计算可以获得37~46倍的加速比,极大地提高了计算效率;4块GPU并行计算加速比从47倍增加到143倍,并行效率维持在70%以上,说明并行算法具有良好的可扩展性。     

网络连接机群上CFD计算的一种负载平衡方法

讨论在网络连接的分布式存储并行机群上进行CFD计算的一种负载平衡方法和相关问题。当数值计算的方法是区域分解法时,各节点机上的负载与其处理的子区域的网格节点数成正比,所以同样大小的子区域可保证负载平衡。采用负载再分配即区域再剖分的方法来达到动态负载平衡,构造了一种对多块网格(每块内为结构网格但各块间的连接关系可以是非结构的)进行近似平均分划的方法,并用一种数据结构和算法保证了剖分出的多个子区之间正确的互联和边界信息传递关系。给出了用这种剖分算法对常见多块网格进行分划的例子,并对三维流场进行了实际计算。结果证明本文方法是正确的,能有效地均衡计算负载,达到较理想的并行效率。     

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

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

短舱及离散精度对螺旋桨桨叶载荷分布的影响

采用多块结构网格生成技术,分别生成螺旋桨加中心体和螺旋桨加短舱两种构型的动态面搭接网格。在相同计算状态下,分别对螺旋桨加中心体和螺旋桨加短舱流场进行不同时空离散精度的数值模拟。与螺旋桨加中心体构型的计算结果进行对比,发现短舱的出现较大程度改变了螺旋桨桨叶载荷分布,并且采用不同时空离散精度,载荷分布的改变程度不同。     

弹道并行计算的设计与实现

多条弹道并行计算是目前作战仿真，多目标探测等领域所急需解决的问题，利用面向对象的软件设计方法，将弹道计算过程抽象化为多个对象，分别对各种对象进行封装，有效地解决了并行弹道计算中内存混乱，可靠性差等问题，实现了多条弹道并行计算。     

ITERATIVE COMPUTATION METHOD FOR NONLINEAR FORCED VIBRATION EQUATION ON HELICOPTER MAIN ROTOR LOAD

ＩＴＥＲＡＴＩＶＥＣＯＭＰＵＴＡＴＩＯＮＭＥＴＨＯＤＦＯＲＮＯＮＬＩＮＥＡＲＦＯＲＣＥＤＶＩＢＲＡＴＩＯＮＥＱＵＡＴＩＯＮＯＮＨＥＬＩＣＯＰＴＥＲＭＡＩＮＲＯＴＯＲＬＯＡＤＬｉｕＸｉａｎｇｊｉａｎ；ＳｈｅｎＸｉｎｋａｎｇ；ＸｕｅＺｈｅｎｇｚｈｏｎｇ（...     

直升机旋翼载荷中的非线性强迫振动方程迭代算法

介绍用迭代方法求解直升机桨叶气动载荷中的非线性强迫振动方程。计算时,先假定诱导下洗已确定,求出桨叶的初始载荷;再通过非线性强迫振动方程求出弯－弯－扭弹性变形。以某一直升机为例,求出了桨叶的各阶谐波载荷和弹性变形。计算结果与实验数据吻合。     

飞机飞行载荷测量载荷方程的优选

介绍用应变计方法测量飞机飞行载荷的载荷方程,应用方差分析和下检验技术进行优选,逐次剔除无应电桥和多余电桥,经自动计算,得到仅会对某种载荷有显著贡献的少量电桥的最优载荷方程,供飞行试验测量飞机飞行载荷使用。算例结果表明,此方法获得最优载荷方程有效。     

学科间载荷参数空间插值传递方法

针对多学科优化分析中面面耦合学科网格间载荷传递的问题, 提出了在参数空间中的三次函数插值传递方法.该方法将三维空间的源学科载荷点投影到二维参数空间平面上, 用三次插值函数对目标点载荷进行计算.此方法解决了叶片表面弯曲对载荷传递精度的影响.某涡轮叶片温度传递的结果表明本文的方法具有很高的精度.      

网格生成和数值模拟的讨论

讨论了多种网格生成方法及其数值模拟。给出一种生成多部件复杂外形飞行器分区网格和Ｅｕｌｅｒ方程分区解。讨论了前缘分离涡为主体的旋涡流场。作了跨声速流分区的并行计算的初步尝试。给出一种二维非结构网格的生成方法。     

基于模糊变结构的空间飞行器姿态控制

根据空间飞行器姿态系统具有非线性、强耦合、多输入多输出的特点,对飞行器姿态模型的非线性和不确定性,提出分散模糊变结构控制方法。利用模糊系统对不确定性函数进行逼近,将获得的模糊系统函数作为系统不确定性界函数,对模糊逼近所带来的误差以及外部干扰项,采用控制补偿方法。理论分析和仿真研究表明,提出的控制方法具有姿态跟踪精度高,实时计算量小,便于工程实现等优点。     

压气机叶片自动优化设计

采用三次多项式和多圆弧方法生成叶型中弧线 ,三次多项式分布叶型厚度 ,对叶型进行参数化。用此参数化方法逼近一低速叶型和一跨音叶型 ,逼近程度令人满意。采用 N-S方程正问题数值计算结合单纯形法寻优 ,构成叶轮机械叶型自动优化设计方法。以消除叶片吸力面分离为目标 ,采用上述方法得到了满意的结果     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。