跨音速差分的并行运算初探

本文就跨音速平面定常小扰动位流的差分计算问题,在“银河”计算机(YH-1)上进行并行运算的初步探讨。差分方程采用隐式近似因式分解法(AF2)迭代计算。与在同一机器上的串行运算相比,并行运算比串行运算提高计算效率约16倍。     

稀薄气体圆柱绕流的并行化DSMC方法研究

本文针对ＤＳＭＣ方法的内在并行特征构造了两种并行算法,并在ＰＶＭ环境下实现了并行化计算。计算结果表明,实现ＤＳＭＣ方法的并行化能够显著缩短计算时间,提高计算效率,不同并行算法能够达到的加速比存在差异。     

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率,对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序,该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想,验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合,说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动。为了提高计算规模及计算速度,作者对程序进行了并行化并针对微机机群系统进行了并行优化。实际计算表明本程序具有较高的并行效率。      

轴对称跨音速进气道外罩绕流的高效差分算法

 本文考虑轴对称跨音速进气道外流场的有限差分计算方法。采用守恒型位流方程、贴体坐标网格和精确边界条件,根据最佳收敛准则对轴对称情形设计出新的近似因式分解迭代格式,并将所设计的格式应用到计算皮托式进气口的轴对称跨音速流场。对几种典型进气口的计算表明,本文格式收敛快,计算结果与实验符合很好。     

跨、超音速流动的区域分解方法与并行算法

研究二维跨、超音速无粘性流动的 Euler方程区域分解方法、并行算法及其应用。通过内边界耦合条件实现相邻子区域解的光滑过渡,以得到总体流场的数值解。发展了一种多块区域之间守恒型的有效内边界耦合方法,对二维翼型跨音速流动和钝头体超音速流动等进行了分区数值求解,分区计算结果与其他单区计算结果作了比较,并讨论了多种区域分解数目的分区计算效率。并行计算采用纯结点并行编程方式和“先进先出”的同步控制等待机制,利用 PVM并行环境对二维绕翼型跨音速流动做了二区和四区分区并行计算。     

圆管内不可压粘流的分布式并行算法

用并行算法求解了圆柱通道内不可压层流流动。采用有重叠的条形区域分裂法时计算域进行分区，分别求解动量方程和压力校正方程，对压力校正方程采用块修正法来加速线代数方程的收敛速度。计算结果表明有较好的并行加速比和并行效率。     

跨音速叶轮机械叶片的完全三维设计方法

设计方法由一种平面叶栅跨音速流动反问题推广得到,其特点是既计及叶片气动性能,又能在一定程度上兼顾强度、冷却、工艺等方面的要求。由于应用了一种适合叶轮机械中三维流动特点的广义ＶｏｎＭｉｓｅｓ坐标系,克服了反问题所特有的求解域边界不定的困难;并使正问题计算和反问题设计统一在同一算法之中。典型算例的结果表明：提出的跨音速叶片完全三维设计理论模型是合理的。     

跨音速翼型非守恒全位势流AF3高效算法的应用

本文利用Baker的AF3高效迭代算法,改进了Jameson提出的计算跨音速翼型非守恒全位势流的有限差分方法,研制了任意翼型跨音速分析的高效算法程序。许多试验算例表明,本算法,收敛速度快,且计算结果与Jameson的原方法非常一致。     

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

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

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

皮托管式轴对称进气道跨音速内,外流场的AF—2方法计算

本文构造了轴向大扰动,径向小扰动轴对称速势方程的AF-2分解格式,计算了IC型皮托管式进气道的跨音速内外流场,并对这种分解格式进行了线化的稳定性分析。分析表明:AF-2方法的稳定性优于SLOR方法,其收敛速度至少比SLOR方法快10倍。     

AF-2格式的加速收敛参数研究

一、引言 Ballhaus等人提出的求解跨音速定常流动的隐式近似因式分解法（AF-2）具有收敛快速,占计算机存储少等优点。文献[1,2]分别对于完全小扰动方程和全位势方程的计算给出了很好的例证。文献[3]表明上述结论对于横向小扰动位势方程也是成立的。AF-2格式的快速收敛主要是通过加速收敛参数的循环取值来实现的。因此,合理正确地选取加速收敛参数序列对AF-2格式至关重要。     

双三角翼大迎角绕流的并行计算

本文在ＬＵ－ＳＧＳ方法求解可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的串行程序基础上,发展了适合ＭＰＰ系统类型的并行计算机的并行算法程序,并针对双三角翼大迎角绕流问题,就不同的并行计算方法进行了对比研究,探讨了优化并行策略,提高计算效率和计算加速比等问题。     

绕升力机身跨音速流动的有限元素法数值分析

 <正> 本文应用全速位方程最小压强积分的有限元素法解绕升力机身的跨音速流动。机身头部可以是尖的或带有进气道的。用头部伸出无穷长圆柱来模拟头部进气。尖头机身在元素形状处理上具有一定的复杂性,从数值实验角度考虑我们分别按尖头和头部伸出无穷长细圆柱来近似尖头头部两种方法对尖头机身进行计算。采用人工密度法捕获激波,但由于绕机身流线形状比较复杂,直接采用机翼问题中的人工粘性公式,解往往不收敛。我们推导出一种较精确的,适用于复杂流动情况的人工粘性计算公式,为了加速收敛,应用网格逐次加密技术,可使收敛速度提高2～3倍。本文还给出绕升力机身跨音速流动的远场速位解析式。     

高可扩展格子Boltzmann方法

格子Boltzmann方法(Lattice Boltzmann Method,LBM)是计算流体力学中的一种常用方法.由于LBM中的格点仅与相邻的格点间存在数据传递,因此具有良好的并行性.LBM并行算法中的数值通信部分通常采用的是格点上的微观量——分布函数.每次传递的分布函数具有多个不同的速度方向,为了进一步减少LBM并行算法的通信开销,从格子 Boltzmann方法的物理特性以及相应的串行程序为切入点,深层次地挖掘可并行的因子,设计了专门用于通信面的类,使用数量较少的宏观量进行通信,降低了通信所占的比重,缩短了通信时间,提高了加速比和效率.实验表明,在4 096个计算核上依然有良好的加速比和效率.     

多块并行计算中负载平衡策略及时间成本估算方法

研究了三维Euler方程的全隐式无分裂方法的多块分布式并行算法.为提高并行效率,给出了一种负载平衡方法,负载平衡前后的墙上时间、并行效率的比较验证了该方法的有效性.针对研究所用的大型并行机特性,基于最小二乘原理,给出了一种估算计算时间成本和通讯时间成本的方法,并行算法的并行性能统计表明该估算方法的有效性及合理性.不同流场的数值模拟表明,多处理器并行计算得到的数值结果能够和单处理器上的串行结果保持一致.     

分布存储结构并行处理机上隐式有限差分算法并行尝试

通过对传统的Beam-Warming隐式格式的并行化尝试,实现了无粘和粘性空气动力学方程的并行求解。计算中采用了SPMD并行编程模式,并将静态非重叠网格区域分裂技术用于建立计算区域和并行处理系统处理器节点之间的映射。结果表明原等价串行算法的稳定性和收敛性被基本保留,获得的加速比和并行效率较为满意。     

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

姿控发动机喷管内流场的DSMC并行算法

在分析了直接仿真Monte-carlo（DSMC）方法内在并行度的基础上，针对工作站群机并行系统（NOWS）分布主存的特点，构造了适合于DSMC仿的并行算法模型，在PVM环境下将其应用于姿控发动机喷管内流场的DSMC仿真，结合并行算法效率评价准则对如何提高并行计算效率进行研究。仿结果表明了实施DSMC并行计算的必要性与有效性。     

机翼跨声速绕流的压强极小积分有限元算法及应用

本文从全速位方程出发。利用压强极小积分,给出了机翼跨声速绕流的有限元算法,并编制出相应的FORTRAN语言程序。文中给出了适于有限元法的网格生成技术,采用了线松弛迭代法和逐步加密法进行数值求解,提供了M6机翼和一个三角翼的计算结果。