桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

涡轮电推进系统建模及控制规律设计

为了研究涡轮电推进系统的工作特性及控制规律设计方法，建立了涡轮电推进系统总体性能仿真模型。基于部件法建立了涵道风扇性能计算模型；基于电机效率图建立了发电机及电动机性能计算模型；基于面向对象的建模方法建立了电推进模块计算模型。将建立的电推进模块计算模型加入涡扇发动机部件级仿真模型中，构建了涡轮电推进系统的共同工作方程组，即其性能仿真模型。采用逆算法研究了不同几何参数控制规律下，涵道风扇推力调节时推进系统的工作特性。通过在性能仿真模型中添加推力平衡方程，并结合差分进化算法构建涡轮电推进系统优化模型，形成了涡轮电推进系统控制规律优化设计方法，应用此方法计算了推进系统节流过程控制规律和性能参数。计算结果表明：本模型设计点性能计算结果与国外发动机性能计算软件计算结果最大偏差仅为0.22%，验证了建模方法的有效性；当涵道风扇推力改变时，可能导致涡扇发动机出现超温、超转和喘振等问题，通过联合调节尾喷管喉部面积和涵道引射器面积，可以有效改善这些问题；设计出的节流过程控制规律可以保证不超温、不超转，具有足够的喘振裕度，且在节流过程中还能进一步降低耗油率；该控制规律还具有连续、单调的特点，利于工程实现。     

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

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

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

结构动力响应精细积分级数解的并行计算

讨论了精细积分法的存储需求问题,指出级数解更适合计算大型结构的动力响应分析.研究了精细积分法级数解的并行计算,精细积分法的时程积分公式包含两项,采用第一项串行、第二项并行的并行算法对精细积分级数解实施并行计算,给出了相应的实现流程.讨论了并行算法的负载分配策略,为减少处理器等待初值的时间,对时间步数实施非平均分配.算例表明,并行算法具有好的加速比.     

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

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

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

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

亚音速位势方程的并行计算

将复杂组合体亚音速气动力特性串行算法改造成并行算法，并在PAN95并行机上进行了数值试验。试验表明，并行计算不仅缩短了计时，而且可增加到18000个面元，计算结果与实验结果吻合较好。     

Linux环境下MPI并行编程与算法实现研究

并行计算实现求解一些复杂的问题,并行编程技术的进步使得并行计算逐步应用到很多领域。MPI是一种消息传递编程模型并可以被广泛使用的编写消息传递程序的标准,已扩展为一种优秀的并行程序设计语言。该文简要介绍了有关并行编程的基本概念与方法,深入讨论了MPI并行编程环境与并行编程模式,对MPI并行性进行了分析,给出了Linux环境下MPI编程环境的配置方法,最后结合一个实例详细描述了用MPI实现并行算法的方法。通过对实验程序的结果分析,在Linux机群环境下用MPI实现复杂的并行算法是简单的,并且能取得较好的性能。     

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

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

计算流体力学并行计算技术研究综述

本文综述了中国空气动力研究与发展中心计算空气动力学研究所近年来开展的并行算法研究应用工作.包括有基于工作站机群和微机机群建立了分布式网络并行计算环境;在此并行环境下开展了CFD并行算法研究和并行应用软件开发研制;与国防科技大学计算机学院合作,开展了大型系列并行机研制过程中所选择的CFD并行程序的研制,并对计算机并行性能进行了测评.     

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

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

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

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

面向流场计算的 PVM 并行程序设计研究

针对流场计算的特点,提出了分布共享变量思想和异-同步并行算法,采用SPMD的编程模式在并行虚拟机（PVM）上初步实现了流场并行计算。     

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

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

非结构混合网格自适应并行技术

计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。     

混合小推力航天器轨道保持高性能滑模控制

针对采用太阳帆、太阳电混合小推力推进的航天器，研究了其在日心悬浮轨道的保持控制问题。为解决已有控制方法中未综合考虑内部未建模动态和外部未知扰动的问题，以及进一步提高系统控制性能，设计了一种高性能滑模控制策略。首先，考虑模型不确定性，建立了混合小推力航天器在日心悬浮轨道柱面坐标系的动力学方程；其次，基于改进型条件积分滑模面和径向基（RBF）神经网络设计了控制律，结合自适应方法在线估计不确定参数；接着，将求取的虚拟控制量在推进剂最优条件下转换成实际控制量，即太阳帆姿态角和太阳电推进力；最后，数值仿真验证了上述设计方法提高了系统鲁棒性，减小了轨道位置超调，并且混合推进相比于单一太阳帆推进，在更短收敛时间内控制精度提高了4个数量级，相比于单一太阳电推进，一年可以节省约89.6%的推进剂。     

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

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

一种DSMC方法的并行策略

为提高直接模拟蒙特卡罗（DSMC）仿真模拟的并行计算效率,基于消息传递接口（MPI）的并行环境,通过对比分析主从模式及对等模式两种程序设计模式下的并行效率,探讨了对等模式下非结构网格DSMC并行程序实现的关键技术及实施途径。提出了一种非结构网格下动态负载平衡DSMC仿真模拟的并行策略,设计了基于对等模式动态负载平衡的DSMC并行算法。最后以钝锥外形的高超声速绕流问题进行仿真模拟,验证本文并行算法的有效性,结果表明,本文设计的基于对等模式动态负载平衡的DSMC并行算法能够以高效的并行效率给出合理的结果。     

巨型星座覆盖性的高效分析方法

针对巨型星座设计及建设过程中的覆盖性分析问题，提出了一种高效覆盖性分析新方法：并行墨卡托投影图叠加法（PMPS）。该算法通过对墨卡托投影图的图像叠加处理，进行覆盖域计算，突破了传统网格法（GPA）中网格数量与星座卫星数量的关联约束，降低了计算消耗对卫星数量增长的敏感度。给定计算精度下，随着卫星数量的增加，GPA与PMPS计算消耗比值逐渐提高，并最终趋于极限值，该极限值和GPA选取网格量呈正相关。仿真结果表明；全球覆盖分析中GPA网格数为10⁴量级、满足指定计算误差时，未引入并行计算的PMPS效率相对GPA最大可提高1～2个量级，当GPA网格量更大时，如10⁵量级时，未引入并行计算的PMPS效率相对GPA预计最大可提高2～3个量级。