用双流体模型模拟湍流两相流场

首先针对颗粒浓度在两相流动中的重要性，建立了颗粒体积浓度方程，阐述了颗粒浓度在求解过程中的特点及颗粒体积浓度方程在计算方法中的贯穿。其次计算了沙层上的稀相可压气固两相流、直管内的不同压液固两相流及叶片扩压器中的密相可压力固两相流的流动，发现了颗粒浓度两种截然不同的分布，并验证了用不可压或稀相的模型来计算可压或密相的模型会带来误差，说明了在计算可压密相两相流动时，可压及密相的因素必须考虑到。     

密相、湍流、可压、气固两相流理论及在管道中的应用

就密相、湍流、可压、气固两相流的基本方程用双流体模型进行了较说尽地归纳。与单相流体相比，在两相流的基本方程中，多了一个变量（体积浓度）及两相间的耦合作用，并用CDIPSA方法对可压密相上升管道湍流流动进行了计算，结果表明：随颗粒直升的增加，使得单个颗粒脉动速度减小；但是气体及颗粒的速度差却随之增大。     

高焓流动中的可压缩颗粒求解器（第1部分）：考虑多物理效应的点力颗粒两相流理论方程

从含烧蚀颗粒的高超声速边界层到航空发动机燃烧室或火箭发动机燃烧室里的气液、气固两相燃烧等问题，都存在弥散相颗粒调制携带流体的两相耦合物理过程。尽管不可压颗粒两相流的理论、数值工具和实验都比较成熟，但可压缩颗粒两相流与不可压缩的情况是截然不同的。本文首先对可压缩颗粒两相流动力学和热力学进行理论研究，理论分析了经典的可压缩流体动力学和热力学方程，把颗粒动力学对不可压流动的调制理论修正拓展到了可压缩情况，发现了可压缩颗粒两相流特有的动量和能量调制无量纲参数；推导了考虑多物理效应的颗粒动力学方程，并将其与携带流体相的控制方程耦合。研究发现，尽管多物理效应下的颗粒动力学是复杂的，但其影响携带流体相的路径是唯一确定的，即：多物理效应的颗粒动力学仅能通过相间阻力调制携带流体相的动量，通过相间阻力做功和热对流调制携带流体相的能量。     

可压稀相两相流场的数值模拟

针对可压两相流动，在阻力系数中考虑了可压缩性及颗粒聚积所造成的影响，用TVD格式对沙层，平板及JPL喷管内的气固稀相两相流进行了数值模拟，所获结果和已有的实验数据及有关文献吻合很好，在JPL喷管中，随着质量分数的增加，由于气体和颗粒间的相互作用逐渐增强，气相的温度增大，马赫数减小，但压力变化很小，而随着颗粒直径的增加，气相轴线上的温度减小。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

管内流动的临界热流研究

本文是对管内流动的临界热流进行试验研究和理论分析.从沸腾机理分析和应用两相流的相似原理进行处理,我们推导了表达临界热流计算的准则方程.用6种不同的冷却剂,355个实验数据进行修正处理求得了常数A和指数h,m,n的数值.引入一个考虑冷却套缝隙高度变化影响的无因次修正项,上述方程成为计算管内流动临界热流的通用方程     

微分雷诺应力模型在激波分离流中的应用

由于Boussinesq涡粘性假设的湍流模型在一些流动中如激波分离流中预测失真,本文介绍一种基于完全雷诺应力方程的微分雷诺应力模型(DRSM),应用压缩性修正和改进耗散率方程构造出应用于可压流体的可实现DRSM,SSG压应力模型测试了两组参数。将其应用于ONERA M6机翼在跨声速大攻角的计算,结果表明DRSM能更好地模拟激波位置和流动分离。     

可压流双方程紊流模型及其应用

 本文提出了可压流紊流模型。所用的方法类似于K-KL双方程模型,但考虑了气流压缩性的影响。方程中的常数已由实验确定。用本模型计算了管内超音速混合流动,其理论计算曲线与实验数据相当符合,证明该模型较不可压紊流模型更为适用。     

吸气式发动机试验中空气流量的计算

为降低吸气式发动机试验中空气流量测量的计算误差，分析三种常用的测量空气流量计算方法的基础上，提出了可压流温比修正算法和无方法误差的程序计算法，分析了各种方法的适用范围和方法误差。同时，提出了空气湿度、直接加热污染和管壁受热等对物性参数和管径的修正方法。结果表明，程序计算法没有误差，在高超声速飞行的模拟试验中尤为适用；可压流温比修正法可以在工程上代替程序计算法，在广阔的压比和来流总温条件下，其方法误差可以忽略。     

亚音速扩压器附面层计算

本方法是一个比较实用的工程近似方法。其基本思想是考虑了管内附面层和核心流之间存在着强干扰,必须将附面层方程和无粘核心流方程联立求解。否则,两层之间的流动可能出现计算上的不协调,导致发散或物理解不合理。 将小曲率下的定常可压紊流普朗特附面层方程,通过沿附面层厚度方向的积分运算,得到动量积分方程和速度矩积分方程。为了利用不可压附面层的结果和参考焓概念,引用Mager变换,将此两积分方程变换成以转换物面坐标为自变量的不可压形     

可压缩流动SIMPLE算法压力修正方程边界条件处理探析

可压缩粘性流体Nivier-Stokes方程SIMPLE算法求解中压力修正方程起着关键的作用.提出从连续性方程本质出发确定压力修正方程边界条件的方法.通过对双圆弧(凸包)通道内流场的计算验证了方法的正确性.     

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

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

压缩拐角区的压力,热流工程计算方法

 <正> 引言 从压缩拐角分离流动的试验及数值计算出发,建立起一个实用范围较广精度较高的工程估算程序,方便快速地估算出压缩拐角分离区和再附区的峰值压力、峰值热流;分离区和再附点位置;分离区长度等。这是气动设计部门关注的课题。这些工作是在文献(3,4)中的经验公式的基础上,按照分离流动流谱,分区进行迭代修正计算,对二维情况给出了完整     

固体火箭发动机内轴对称两相流场的一体化数值模拟

将固体火箭发动机燃烧室和喷管作为一个整体，应用时间相关法进行了轴对称两相流场的一体化数值模拟。对气相采用ＭａｃＣｏｒｍａｃｋ显式差分格式，对颗粒相则采用特征线法进行计算。计算结果表明该方法可以对带潜入式喷管的发动机内复杂的两相流场进行有效的数值模拟。     

大涡模拟超声速两相流混合层的非线性分析

为研究超声速两相流混合机理,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构.用二维盒式滤波法处理可压缩气相N-S方程,得到亚格子形式的LES控制方程;用轨道模型仿真液滴运动,并用单液滴蒸发模型模拟液滴蒸发过程.利用伪邻近法和互信息量法确定了相空间重构的嵌入维数和时间延迟,对混合层压力时间序列进行了相空间重构.分析了仿真结果中混合层的涡结构特性和液滴对流场的反作用效果,并讨论了混合层失稳过程的非线性特征.发现混合层发展的不同阶段与其系统的非线性特征值密切相关.      

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

对颗粒相采用颗粒轨道模型,气相求解可压缩N-S方程组,计算方法采用显式Runge-Kutta时间推进法与有总变差衰减(TVD)性质的高精度MUSCL-Roe格式;自主开发了曲线坐标系下二维轴对称可压缩N-S方程组的解算器Solve2D,研究了固体火箭发动机喷管中颗粒相对流场的影响以及不同尺寸颗粒运动规律.结果表明:颗粒相对流场的影响主要表现在喷管喉部以及扩张段,和单相流场相比,沿轴线马赫数减小,且颗粒尺寸越小减少得越多;沿轴线气相温度升高,且颗粒尺寸越小温度升高越多;颗粒尺寸越小,无粒子区越小;颗粒越大与收缩段壁面碰撞越剧烈,无粒子区越大.      

水平槽道内气固两相湍流中颗粒行为的PIV实验研究

运用PIV技术对充分发展的水平槽道内稀疏湍流两相流中的颗粒行为进行了定量研究。实验中壁面雷诺数Rer=430,使用颗粒为直径60μm的聚乙烯微珠,测量了0．1％和0．5％两种质量荷载。研究发现,颗粒最大浓度出现在y＋≈10的位置上。对两种质量荷载,下壁面附近的颗粒体积分数最大值均远小于10^3量级,仍属于稀疏两相流。槽道下壁面附近,向上运动的颗粒概率大于向下运动的颗粒,两者概率的差异随着法向距离的增大而减小。在y’=10～30的范围内,颗粒相“喷射”事件（Q2）和“下扫”事件（Q4）概率分别显著地大于和小于当地流体对应值。颗粒相主要通过促进Q2且抑制Q4来实现对流体湍流的影响,颗粒的这一作用随着法向距离的增加而逐渐减弱。     

过载条件下固体发动机内流场数值模拟

应用颗粒轨道模型，连续相控制方程按二阶迎风有限体积方法进行离散，并对纵、横加速度载荷均为20g、30g和35g的固体火箭发动机燃烧室内两相流动进行了模拟。结果表明，有纵横加速度载荷的情况下，发动机燃烧室内颗粒相会形成粒子聚集流，对承载方向的装药和壁面产生严重的冲蚀，明显改变了发动机燃烧室内原有的轴对称流动形态，同时承载方向上粒子聚集流的最大密度点随横向加速度的增加而远离发动机后封头。这些结果与实验发动机试车结果有较好的一致性，为发动机绝热层异常烧蚀机理分析提供了理论依据。