液体发动机燃烧室流场模拟的并行SIMPLE算法

用串行和并行ＳＩＭＰＬＥ算法对液体火箭发动机燃烧室内的复杂流动进行了数值模拟。控制方程组用欧拉坐标系下的Ｎ－Ｓ方程组描述。在并行虚拟机环境下研究了交错网格系统中参数的传输和接受模式, 从计算结果可以看出并行计算的效率较高, 是解决此类复杂流动大规模数值模拟问题的有效手段之一。     

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

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

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

三维低速NS方程的并行计算

本文在求解低速NS方程的串行程序基础上,针对分布式存贮网络异构的并行编程环境建立了隐式并行计算软件,并行计算的方法采用区域分解法.并行软件对高速列车进行了气动特性计算,并对并行效率和并行可靠性进行了研究,结果表明并行计算结果与串行计算吻合良好,并行改造后的软件可以应用于工程实际问题.     

喷管内流场并行计算方法研究

建立了基于区域分解的并行计算方法,将已有的后台阶喷管内流场三维ＴＶＤＭａｃ－Ｃｏｒｍａｃｋ显式格式的串行计算程序改造成为并行程序,获得了较高的并行加速比。并行计算结果同串行结果进行了对比,证明计算是成功的。     

流场计算程序并行化方法的研究

利用“流场分区”的方法把流场计算程序并行化,并成功地在Transputer并行计算机上运行,并行效率可达到80％。介绍了流场划分、并行编程以及提高并行效率的方法。     

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

有限元的OpenMp并行计算技术

传统有限元串行计算技术耗时长、效率低,已远远不能满足工程实践需要。通过分析串行计算的耗时分布、程序可并行性及数据相关性,在保持串行程序不变的基础上,采用基于多线程的OpenMP技术实现单元刚度矩阵形成总体刚度矩阵的程序并行化。针对有限元计算机翼实例进行了不同线程数目下并行加速比的研究,结果表明,OpenMP技术可有效提高有限元计算效率,且存在一个最佳的线程数,使得并行加速比达到最大,此外加速比还与问题的计算量等因素有关。     

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

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

Adaptive mesh refinement computation of acoustic radiation from an engine intake

A block-structured adaptive mesh refinement (AMR) method was applied to the computational problem of acoustic radiation from an aeroengine intake. The aim is to improve the computational and storage efficiency in aeroengine noise prediction through reduction of computational cells. A parallel implementation of the adaptive mesh refinement algorithm was achieved using message passing interface. It combined a range of 2nd- and 4th-order spatial stencils, a 4th-order low-dissipation and low-dispersion Runge–Kutta scheme for time integration and several different interpolation methods. Both the parallel AMR algorithms and numerical issues were introduced briefly in this work. To solve the problem of acoustic radiation from an aeroengine intake, the code was extended to support body-fitted grid structures. The problem of acoustic radiation was solved with linearised Euler equations. The AMR results were compared with the previous results computed on a uniformly fine mesh to demonstrate the accuracy and the efficiency of the current AMR strategy. As the computational load of the whole adaptively refined mesh has to be balanced between nodes on-line, the parallel performance of the existing code deteriorates along with the increase of processors due to the expensive inter-nodes memory communication costs. The potential solution was suggested in the end.     

二维翼型微吸吹气减阻控制新技术数值研究

结合抽吸气转捩控制和微吹气湍流减阻控制的特点,探索了一种新的吸吹气减阻控制技术。使用Fluent求解器,并利用用户自定义函数(UDF)二次开发对其自带的Wilcox转捩模式进行了修正。在此基础上,数值研究了吸吹气控制对翼型阻力性能的影响。结果表明：在一定的吸气量范围内,吸气、吸吹气控制都能使翼型总阻力减小,且在同一雷诺数下,吸气控制能使翼型总阻力减小约3%,而吸吹气联合控制使翼型总阻力减小约16%。由此可见,吸吹气控制技术是一种行之有效的减阻控制技术。     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method" and ‘‘minimum distance scheme of donor elements method" are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.     

径流与斜流涡轮设计的通流反问题方法

自主发展了径流及斜流涡轮的准三维流线曲率法通流设计程序,应用该程序给出了一项500 N推力级微型涡喷发动机斜流涡轮的气动设计和叶片造型.通过计算流体力学检验和整体叶轮强度检验,该微型斜流涡轮超过了预定的设计指标.为实现径流涡轮和斜流涡轮设计反问题的统一,流线曲率法的S_2,m流面准径向平衡方程是沿流线方向m和流线的准法线方向口写出,该系统也可用于轴流涡轮的设计.采用了在所有计算站上的"全可控涡法"进行气动设计.      

一种DSMC方法的并行策略

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

An efficient SAR parallel processor

A parallel architecture especially designed for a synthetic-aperture-radar (SAR) processing algorithm based on an appropriate two-dimensional fast Fourier transform (FFT) code is presented. The algorithm is briefly summarized, and the FFT code is given for the one-dimensional case, although all results can be immediately generalized to the double FFT. The computer architecture, which consists of a toroidal net with transputers on each node, is described. Parametric expressions for the computational time of the net versus the number of nodes are derived. The architecture allows drastic reduction of the processing time, preserving elaboration accuracy and flexibility     

再入飞行器化学非平衡流动的隐式向量算法研究

在YH－2计算机上对再入飞行器化学非平衡流流场的NND隐式程序进行向量化计算，在对程序作了重构优化后，向量加速比达3．4。这是对化学非平衡流N－S方程隐式程序进行向量化计算的成功尝试，也为进一步实现程序的并行执行奠定了良好的基础。     

并行环境下外挂物动态分离过程的数值模拟

 并行工程的发展提出了计算流体力学的发展新概念,在此基础上研究开发了一个以非定常流计算为核心的多体分离数值仿真集成系统。该系统采用运动嵌套网格方法,能够处理飞行器在空间中的大尺度运动,结合六自由度运动学方程,具有实时仿真运动部件在气动力、弹射力和重力等耦合作用下运动的能力。基于外挂物投放数值模拟具有网格运动、非定常的特点,从网格重建、多重网格方法等方面探索了提高计算效率的途径和实现方法,结合并行化计算使分析周期大大缩短。以20个中央处理器(CPU)为例,并行计算和3层多重网格使计算速度提高约40倍,同时分布式存储、并发计算的特点缓解了对单节点的内存需求,解决了工程应用中受计算周期、网格规模等限制的难题。为了验证本文方法的正确性,对外挂物投放标模实验进行了数值仿真,计算得到的外挂物运动轨迹在0.5 s时间范围内与实验值基本吻合;对分离过程后半段出现的偏差,从气动力变化的时间历程上进行了说明,提出偏航力矩的误差和累积效应导致了预测出现偏离,对关键的分离初始阶段,该方法可以取得令人满意的效率和精度。     

非稳态单喷口射流控制扩压器内边界层的分离流动

对扩压器内单喷口射流产生的流向涡控制扩压器内部边界层分离流动行为进行了研究。研究的对象是一个扩张角为10°的二维扩压器,扩压器进出口面积比为1∶4.7,入口高度为15 mm。采用数值计算方法研究了当射流出口速度为非稳态时对分离流动的影响。给出了不同时刻,不同截面上的旋涡涡量分布。结果表明,采用非稳态喷射流向涡方法可以在不改变主流流量特性的情况下,有效地减小分离流动区域。     

跨音速翼型分离绕流的粘性/无粘干扰迭代计算

本文用有粘/无粘干扰迭代的概念计算了跨音速任意翼型的绕流问题。位流的速位方程用AF2格式求解,而边界层微分方程用C-S盒式法求解,逆算法的引用可以克服边界层方程在分离点处的奇性问题,对分离区湍流代数模型的修正可以得到与实验更吻合的结果。计算结果表明有粘/无粘干扰迭代概念在小分离泡的情况下也是适用的。