直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.      

高温高压超声速气液两相流场数值模拟

采用欧拉-拉格朗日方法对燃气-蒸汽发射动力装置内高温高压超声速燃气中横向喷雾的气液两相流进行了数值模拟研究,气相采用RK(Runge-Kutta)-AUSM+(advection upstream splitting method+)格式求解Navier-Stokes方程,液相应用颗粒轨道模型,两相之间的耦合通过在气相...     

三维管道超声速气固两相流动的CFD/DSMC仿真

采用NND格式求解气相N-S方程，采用DSMC方法模拟固相颗粒运动，数值模拟了三维管道超声速气固两相流场，分析了固体颗粒对流场结构产生的影响。计算结果表明，将DSMC方法与CFD方法相结合是解决气固两相流动问题的一种有效途径。     

航线网络中基于最小费用流的航线选择问题

针对航线网络中航线选择优化问题,分析了空中交通航线网络结构和空中交通流的特点,给出了空中交通航线网络的数学描述,考虑从求解最小费用流问题的角度寻求最优航线,建立了航线网络中求解最小费用流问题的数学模型。借鉴了连续时间网络最小费用流问题中结点到达费用的概念,并结合简单算例的求解验证了航线网络中用最小费用流方法选择航线的可行性,具有良好的计算性能。     

考虑汽化的喷嘴高温气流注水两相流动规律

针对某火箭发动机考虑汽化的注水气液两相流问题，研究了液态水的汽化机理。根据不同环境压力下水的饱和温度建立了汽化模型，编制适合于液态水的专用汽化相变求解程序，并通过添加汽化组分源项及热源项的方法将汽化相变计算嵌入到多相流场控制方程中，实现了考虑液态水汽化相变的三维多相流场求解。结合经典算例进行了对比，特征点的两相流场温度计算值和试验的误差控制在8.5%以内，验证了程序三维计算的可靠性。距火箭喷嘴不同距离横截面的两相质量转化率曲线显示，质量转化率从最低值上升到最高值，然后降低到最低值，揭示了考虑相变的气液两相相间反应转化机理。该方法可为运载火箭发射平台发射起始阶段水室的汽化降温提供参考。      

固体火箭发动机气动喉部的两相流数值模拟

本文对两相流条件下的环缝式气动喉部方案进行了数值模拟。气-粒两相流动模型采用二维轴对称欧拉-拉格朗日模型。研究了二次流不同喷射位置、角度、流率及喷嘴个数对气动喉部调节有效喉部面积大小的影响规律。模拟得到了两相流条件下气动喉部的流场特征,结果表明使二次流的喷入位置越靠近喉部、增大二次流流量或减小喷射角度都能明显增加气动喉部调节性能;由此归纳出增强气动喉部对主流扼流能力的两种主要机制,即增加挤压和流阻。     

双工况氢氧发动机燃烧与传热数值分析

应用三维湍流N-S方程以及颗粒轨道模型描述双工况氢氧发动机内部喷雾两相燃烧流动过程。两相之间的质量、能量交换由液滴蒸发模型计算，气相化学反应速率由Arrhenius公式计算。通过耦合求解气液两相的模型方程，对发动机转工况前后的三维流场进行了数值计算，并耦合计算了燃气与壁面之间的传热以确定壁面的温度和热流分布。另外还对分别采用同轴离心式喷嘴和直流式喷嘴得到的燃烧流场与燃烧效率进行了比较。计算结果表明转工况前的壁面温度与热流都比转工况后大。离心式喷嘴的雾化混合效果与燃烧效率都比直流式喷嘴好。     

大气层内模型预测静态规划拦截中制导

在临近空间机动目标拦截中,拦截器的初始动力段对中制导段乃至末端拦截性能具有重要影响。在模型预测静态规划（MPSP）理论基础上提出了一种初、中制导联合规划制导方法,旨在解决多阶段、快速、最优拦截轨迹规划和制导问题。首先,提出了一种改进的模型预测静态规划方法,该方法不仅可以满足终端约束,还可以生成最优初始状态,并在性能指标中考虑状态变量相关函数。其次,基于等效阻力模型建立了包含动力段与非动力段的两段规划模型,通过采用分段离散以及构建关机点变分关系的方法,避免了内点约束的引入,使MPSP算法可直接求解该两段规划问题。最后,结合提出的MPSP算法以及两段规划模型,实现了终端速度最优的拦截轨迹规划,结合目标预测方法,实现了对机动目标的预测拦截。仿真结果表明,本文方法可提高制导精度和终端速度,且能更好满足对机动目标的拦截需求。     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

轴承腔油气两相流数值模拟的试验验证

为探索一种适于模拟轴承腔油气两相流流动特性的数值计算模型,建立了一套简化结构的轴承腔试验装置,通过试验研究获得了不同工况对轴承腔通风口出气比例和回油口滑油温升的影响规律.应用Fluent软件,采用Mixture多相流模型和Eulerian多相流模型分别搭配标准k-ε湍流模型和RNG(renormalization group) k-ε湍流模型对选定试验工况进行数值模拟.对比结果表明:当转速为2200r/min时,使用Mixture多相流模型搭配RNG k-ε湍流模型计算的滑油温升与试验值的相对误差仅为2.72%,体现了该计算模型在高转速时的精确性,并在此基础上获得了油气两相流速度场、压力场和温度场的分布特征.      

Numerical simulation of a boiling up fluid flow in the variable cross-section channels

A physical and mathematical model of a boiling-up fluid flow in the variable crosssection channels is proposed. The numerical model is based on the solution of the Navier-Stokes nonstationary equations for compressible fluid by the implicit finite-difference method. We present the results of the numerical modeling for two-phase flow characteristics at the initial stage of the evaporation process and the steady-state quasistationary flow of a two-phase medium.     

极低压下两相流挤压油膜阻尼器油膜参数特性

针对极低供油压力工况开展实验研究,以考察气液两相流对挤压油膜阻尼器(SFD)油膜参数特性的影响。结果表明:当SFD入口气体体积分数小于0.9时,油膜阻尼随着入口气体体积分数的增加而减小,直至气体体积分数增大到0.9时,油膜阻尼是纯油状态时阻尼的60%;当气体体积分数大于0.9时,油膜阻尼大幅减小至几乎可以忽略不计。现有理论模型并不适用于极低供油压力工况。基于实验结果,找到了最符合SFD两相流动的等效黏度模型,其理论预测的油膜阻尼与不同供油压力下的实验数据吻合较好,为SFD两相流研究的模型选择提供了依据。     

低含气率影响下圆柱动静压轴承特性分析

为研究低含气率条件下,温度和气泡含量对润滑油黏度的影响,采用涡轮油和氮气制备油气两相流,在不同温度下采集含气率和动力黏度,拟合得到两相流的黏度计算模型。针对带有深浅腔的圆柱动静压轴承,基于拟合的黏度模型,结合有限元法和有限差分法联立求解油膜Reynolds方程、能量方程及两相流等效密度模型,分析了不同偏心率下承载力、摩擦力、端泄流量、刚度、阻尼系数和失稳转速随含气率变化的规律。结果表明:含气率在2.50%以内,含气率越高,黏度值越大;不同偏心率下油膜承载力、摩擦力和动力特性参数均随含气率的增加而增加,端泄流量随含气率的增加而减小;计入两相流的影响使稳定性有所提高,失稳转速最高增幅为9.14%。     

环形燃烧室两相燃烧数值研究

根据环形燃烧室设计需要,运用计算机模拟技术研究模型环形燃烧室三维有旋液雾燃烧流场,研究两相流动与混合,液雾蒸发对燃烧过程的影响。两相流动模型采用颗粒群轨道模型,气相紊流输运方程求解采用双方程k-ε模型和EBU-Arrhenius紊流燃烧模型,热辐射采用热通量法模型,本文提供计算方法可以用来研究不同性能参数对环形燃烧室两相燃烧特性影响。     

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

工作参数对航空发动机轴承腔内二相流动的影响

作为航空发动机润滑系统油气二相流的重要区域,主轴承腔的工作参数对内部二相流动的影响对于发动机润滑系统设计具有重要意义。利用DPM壁面液膜模型,采用CFD方法对某型发动机轴承腔简化模型内油气二相流进行了数值计算,计算结果与现有试验数据符合良好;给出了轴承腔在不同主轴转速及不同滑油流量下油膜厚度、空气和油膜速度的分布以及出口速度变化规律。     

N2O_HTPB固液火箭发动机喷管两相流计算

利用二维轴对称N-S方程对选用氧化亚氮/丁羟基燃料推进剂的固液混合火箭发动机的喷管两相流进行了计算。计算采用MacCormack时间推进预报校正二步格式,采用了Baldwin－Lomax代数湍流模型和两相平衡流模型。计算了三种氧燃比下四个不同喷管的喷管流场参数,并计算了喷管性能,通过比较两相流和气相流的计算结果,分析了不同氧燃比和喷管形状对喷管性能的影响,认为固液火箭发动机的性能主要受氧燃比的影响,为固液混合火箭发动机的设计提供了依据。     

航空发动机轴承腔内油气两相分布数值计算

主轴承腔是航空发动机润滑系统油气两相流的重要区域,了解滑油及密封空气在腔内的分布预测对于认知腔内两相流动具有重要的意义。利用VOF数值计算模型对某航空发动机轴承腔简化模型内油气两相流进行了数值计算,定义无量纲数β^＋用来分析腔内两相分布。计算结果与现有实验数据符合良好。给出了滑油在腔内的相界面分布,描述了腔内的油膜分布及运动,并且分析了无量纲数β＋在腔内周向分布及出口处随着转速及滑油进口流量的变化规律。     

N2O/HTPB固液火箭发动机喷管两相流计算

利用二维轴对称N-S方程对选用氧化亚氮/丁羟基燃料推进剂的固液混合火箭发动机的喷管两相流进行了计算.计算采用MacCormack时间推进预报校正二步格式,采用了Baldwin-Lomax代数湍流模型和两相平衡流模型.计算了三种氧燃比下4个不同喷管的喷管流场参数,并计算了喷管性能,通过比较两相流和气相流的计算结果,分析了不同氧燃比和喷管形状对喷管性能的影响,认为固液火箭发动机的性能主要受氧燃比的影响,为固液混合火箭发动机的设计提供了依据.