外挂物干扰流场特性数值仿真研究

飞行器研制中对于带外挂物的复杂外形气动特性分析研究一直是个难点,对外挂物与挂架不同缝隙的复杂机翼-挂架-外挂物组合外形采用结构重叠网格数值方法,快捷高质地完成了网格生成,通过求解带有k-e湍流模型的Navier-Stokes方程组,得到多个计算工况下的绕流外流场,分析了外挂物在不同工况下的气动特性,结果表明,CFD数值方法可方便快捷地对复杂干扰流场进行数值计算,可为飞行器相关设计提供技术支持。     

适用于混合网格的约束最小二乘重构方法

为强化边界条件对计算流体力学(CFD)计算中流场重构过程的影响,基于三维任意多面体混合网格,对一种约束最小二乘法进行了研究。通过对不同边界类型的透明化处理,统一了约束方程组的形式,从而简化了该方法的应用。为了充分吸收消去法和加权法各自求解约束方程组的优点,提出了混合采用两种方法的思路。针对层流平板研究了加权法中加权系数如何取值问题,数值实验表明加权系数取到5基本已经足够大。最后通过亚声速层流平板和跨声速湍流ONERA-M6算例对比了混合法、加权法与原始最小二乘法,计算结果表明:约束最小二乘法相比原始最小二乘法的流场重构更加准确,特别是对边界质量较差网格的计算优势更为突出;混合消去加权法相比单独采用加权法的计算结果也有所改善。     

叶片外换热计算程序STAN5收敛性改进

涡轮叶片外换热计算程序STAN5,采用Patankar-Spalding方法求解可压缩边界层微分方程,获得流动及温度边界层内的参数分布,进而得到叶片表面的外换热系数分布。本文分析了STAN5程序的求解方法,研究了程序在外流场存在较强逆压梯度时出现计算不收敛及错误收敛等现象的机理,并对程序进行了两点改进:修正了逆压梯度下边界层内部速度分布,建立了边界层内部沿流向的网格分布与主流压力梯度间的联系。改进后的程序有效地避免了计算发散等问题。     

高效气动优化设计方法

基于CFD方法开展气动外形优化设计通常计算量较大,采用离散共轭方法计算目标函数梯度,建立了高效的飞机气动优化设计系统。采用NURBS方法对翼剖面几何外形的扰动量进行参数化,避免了对原始外形的拟合,建立了基于NURBS方法的机翼剖面参数化方法。在复杂外形的优化设计中,动网格方法是关键技术之一,采用无限插值方法生成变形后的网格,并提出采用无限插值方法处理部件之间相贯线发生变化的情况。最后开展了翼身组合体+吊舱+挂架等复杂外形的优化设计,成功地减小了阻力,证明了方法的有效性。     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

改进的弹簧比拟非结构化网格优化方法

提出了改进的弹簧比拟非结构化网格优化方法.分别从网格的稀密和单元质量两个方面对非结构化网格进行优化.首先采用非等刚度弹簧,实现了弹簧比拟方法对网格疏密的控制;其次,分别采用修正刚度系数的顶点弹簧模型和修正初始平衡长度的棱边弹簧模型,改善了原弹簧比拟光顺算法中存在的局部网格单元过渡不够平滑的问题.该方法保留了原弹簧比拟法易于实现、数据结构简单的优点,提高了原算法对网格优化能力,算例中最短边与最长边之比从0.24提高到了0.42(越接近1越好),且对于存在交错单元的初始网格也能优化.      

气助式直接喷射喷油器动态喷雾场数值模拟及试验

为了研究小型飞机发动机气助式直接喷射(air-assisted direct injection,AADI)喷油器喷雾形态及发展规律,采用数值模拟和高速摄影技术对其动态喷雾场进行研究.以锥阀运动的物理学模型为数值计算的移动边界条件,结合动网格技术对AADI喷油器动态喷雾场的速度分布、静动压分布和不同时刻雾化颗粒的空间分布进行动态数值模拟分析.以0.2ms为步长,对3.0ms内不同锥阀开启时刻的喷雾形态、雾化颗粒和贯穿距离进行了试验和仿真对比分析.结果表明:采用CFD仿真模拟技术可对AADI喷油器动态喷雾场的速度分布、压力分布、喷雾形态等参数进行定性分析,有利于揭示喷雾的发展规律,同时通过对喷雾粒径和贯穿距离的定量对比分析,表明了仿真计算的精确性,数值模拟计算与试验结果误差小于5%,满足工程应用的要求.      

双级涡流器环形燃烧室整体流场数值模拟

在任意曲线坐标系下对包括扩压器、双级涡流器及火焰筒在内的环形燃烧室三维两相燃烧整体流场进行数值模拟。由于环形燃烧室形状复杂,采用偏微分方程法与整体分区结合法生成环形燃烧室整体网格。所用的数学模型有:k-ε紊流模型、EBU-Arrhenius紊流燃烧模型、六通量热辐射模型以及颗粒群轨道模型等。在非交错网格系下,气相采用SIMPLE算法求解,液相采用PSIC算法求解。数值分析不同燃烧室进口气流参数以及涡流器几何尺寸对燃烧室流场的影响,并将计算结果与实验数据进行了比较。结果表明本文的计算方法合理性可用于环形燃烧室的研制与优化设计。     

一种基于冗余分配的网格计算任务调度方法的研究

提出了一种基于冗余分配的网格计算任务调度方法,包含资源预测、资源调度和任务分配三部分,并将其应用于用二分法求解高次方程中,取得了运行时间短和兼顾负载平衡的应用效果,验证了该算法的有效性。     

多电飞机混合储能系统的模型研究

简要分析了适用于多电飞机电网的蓄电池—超级电容混合储能系统研究的意义和价值,建立较为完整的混合储能模型,并在模型基础上采用模糊控制等方法,设计出三级控制策略,实现稳定飞机电网电压目的.