用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

关于变形网格“几何守恒律”概念的讨论

在国内外文献中应用包含有变形的动网格进行流动数值模拟时,要满足所谓的几何守恒律条件;通过对这一概念最早提出推导过程和近年来应用发展情况分析以后发现,逻辑上存在不适当之处.从流体力学基本理论出发,对包含有变形的动网格技术的理论基础进行讨论,证实几何守恒律是流体力学控制方程的伴随方程或退化方程;通过简单模型分析有限体积方法离散过程,发现计算过程中不满足几何守恒律所引起的非物理现象本质上是目前有限体积离散处理中微元体积计算方法不符合物理定律所致;因此在变形动网格计算方法中不存在必须要满足的几何守恒律.最后根据理论提出物理背景清晰的微元体积计算方法,数值验证可行.     

高超音速飞行器前体优化设计方法研究

基于响应面法对高超音速飞行器前体进行气动优化设计,流场求解采用欧拉方程,并结合非结构动网格技术,便于适应复杂外形而且在每次外形改变后不必重新生成网格,有利于保证计算效率,响应面模型采用二次多项式来构造,数值结果表明由本文设计的前体/进气道布局结构简单,能获得更高的压缩率和总压恢复系数,以及更高的升阻比,该方法在多约束多目标条件下进行气动优化设计,设计质量较高,有一定的工程应用价值.     

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。     

空中加油系统的多体动力学和CFD仿真

介绍了作者基于Fortran程序开发的一个多体动力学求解器.依托多刚体系统理论、基于绝对节点坐标(ANCF)法描述的常质量/变质量索梁单元和计算多体动力学方程组求解技术,该求解器可用于分析飞行器拖曳结构中的索-刚体约束系统的动力学问题.本文基于并行CFD求解器PMB3D的动态重叠网格功能,构建了多体动力学和计算流体力学的非定常耦合计算程序.以软式空中加油系统为例,对真实外形飞行器拖曳系统的气动-动力学特性进行了仿真研究.探讨了战斗机构型进行软式加油对接时产生的气动干扰问题,研究了受油机的接近速度和初始偏移量对对接效果的影响.     

球轴承多体接触动力学研究

考虑钢球、套圈和保持架的动态接触关系,提出了机械系统中球轴承多体动力学分析的新方法.基于套圈滚道的三角网格模型,实现了钢球和套圈滚道的动态接触力的预测搜索算法,建立了计及润滑和Hertz接触作用的三维角接触球轴承多体接触动力学模型.运用广义-α方法计算分析了预紧力和旋转径向力作用下角接触球轴承的多体接触动力学特性,获得了球轴承的动态接触力、拖动力和运动轨迹及频谱等振动响应,并利用Gupta经典实例模型进行了实验验证.轻载中等速度下钢球的角速度以184.5rad/s-2波幅周期变化,旋滚比以0.01波幅周期变化,角加速度与动态接触载荷的频谱具有相同的56.1,112.2Hz等谐波倍频成分.中等载荷高速下保持架中心的运动轨迹呈现出以83.3Hz和200Hz双频率拟周期的平动运动.      

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

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

子母弹抛壳过程非定常流动的数值模拟

利用弹簧近似和网格重构相结合的非结构动网格技术,耦合求解Euler方程及弹道方程,模拟飞行在46千米高空、0°攻角、马赫数6.5、壳片分离初速分别为6m/s,12m/s、壳片分离角速度分别为0°/s,50°/s,150°/s的子母弹抛壳动态过程。研究表明,抛壳初期,壳片与母弹周围形成复杂流动,在壳片内表面存在激波的生成、发展、减弱等非定常过程,导致壳片受到气动力随着分离距离增加表现为很强的非线性特征;壳片分离初始角速度对壳片姿态的影响不是线性的,尽管分离初始角速度增加到150°/s,由于转动惯量小、俯仰力矩大,壳片实际飞行姿态达不到设计要求,壳片攻角较小,气动阻力小,导致壳片伴随母弹较长时间飞行;增大分离初速是实现快速分离的有效措施。     

应用维数分裂方法推广MUSCL和WENO格式的若干问题

首先对有限差分法和有限体积法的差异进行了讨论,在已有文献观点的基础上补充了二者在边界条件处理及网格需求等方面存在差异的新论据。介绍了逐维推导的MUSCL和WENO格式计算控制体界面通量的过程,认为此类格式计算界面通量的方法直接应用于高斯积分型有限体积法不够严谨,从而得到了应用维数分裂方法构造的MUSCL格式和WENO格式不属于高斯积分型有限体积法的观点,“积分格式”这一定义更能准确反映这类格式的特点。此外,讨论了MUSCL格式和WENO格式在曲线坐标系下不能保证守恒的原因,并简单介绍了消除方法。     

空间近程高速相对运动方程的状态转移矩阵

分析了空间近程高速相对运动的特点,建立了描述近程高速相对运动的线性微分方程组,并给出了状态转移矩阵与解析解。仿真结果表明:对于HEO-GEO近程高速相对运动,所建方程的相对距离计算误差仅为16 cm,不足Hill方程计算误差的15%。所给出的状态转移矩阵对于实现空间目标拦截等任务的实时高精度相对导航和制导算法具有一定的应用价值。     

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

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

An experimental method for squealer tip flow field considering relative casing motion

Squealer tip is widely used in turbines to reduce tip leakage loss. In typical turbine environment, the squealer tip leakage flow is affected by multiple factors such as the relative casing motion and the wide range of variable incidence angles. The development of experimental methods which can accurately model the real turbine environment and influencing factors is of great significance to study the squealer tip leakage flow mechanism. In the present paper, a low-speed turbine cascade test faci...     

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

模型燃烧室两相燃烧大涡模拟的并行计算

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对模型燃烧室气液两相瞬态喷雾燃烧过程进行大涡模拟并行计算研究。通过大涡模拟计算结果与雷诺平均计算结果对比表明大涡模拟方法能更好的模拟流场细观结构。同时分别采用多个处理器和动态内存分配对大涡模拟并行计算程序进行数值模拟,并行计算结果表明:多个处理器之间并行计算结果与PIV测量的瞬态速度场以及出口温度分布实验数据相互吻合,表明采用并行计算和动态内存分配在保证计算结果正确的前提下,大幅度降低计算机内存和大涡模拟计算时间。     

航空发动机过渡态最优控制规律设计的新方法

针对航空发动机过渡态最优控制规律的设计问题,提出了一种新的方法——动态稳定法:在发动机动态特性计算模型的基础上,通过额外提取其中所有的状态量变化率,使其共同工作方程组的偏差趋于0,从而让过渡态仿真稳定下来,此时的稳态参数值即为对应过渡态工作点的各项参数;然后分别根据指定的物理约束条件,通过简单的静态迭代优化即可直接建立相应的控制规律;最后将这些控制规律通过取大/小的方式进行合并,从而获得所需过渡态最优控制规律.以某型涡扇发动机最优加速控制规律的改进设计为例,结果表明该方法具有设计精度高、实现简单以及快速直观等优点.      

COIL化学非平衡流动的一种解耦方法研究

提出了一种计算氧碘化学激光器化学非平衡流动的解耦方法.求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分混合反应流场.空间离散格式为AUSM up格式,用四步龙格-库塔方法作显式时间推进.假设混合气体为热力学完全而热值非完全气体,化学反应模型采用有限速率反应模型.使用一种"点隐式"方法来处理化学源项的刚性问题.计算结果表明本文方法可以有效解决计算中的刚性问题,与现有隐式算法相比减少了大量的矩阵运算,公式推导简单,易于实现和改进.     

超燃冲压发动机等截面隔离段流动特性的LBM模拟

为了探寻研究超燃冲压发动机隔离段内复杂流动的新方法,采用耦合的双分布函数格子-Boltzmann方法(LBM)对不同条件下隔离段内的流动进行了研究,分析了出口反压、进口马赫数、进口附面层厚度对隔离段流动特性的影响.结果表明:随着出口反压、进口附面层厚度的增加,以及进口马赫数的减小,可使激波串起始位置前移.结果说明耦合的双分布函数LBM有望成为研究超燃冲压发动机内复杂流场的一种有效的数值模拟手段.