可压缩湍流封闭模型初探：湍流动能K方程

为解决可压缩流的湍流运动方程组封闭问题，对瞬时值的Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ运动方程进行质量加权平均处理，根据湍流封闭理论对各脉动量关联项进行模化，得到可压缩流动的湍流动能Ｋ方程，完成了建立一个适用于三维可压缩流的湍流封闭模型的第一步工作．     

牛顿流体在圆管内湍流流动剪切率的近似计算

根据剪切率的定义导出了描述湍流流动剪切率与时均流剪切率和脉动流剪切率关系的基本方程 ,该关系式与流体的流变性类型无关 ,适用于牛顿流体和非牛顿流体。根据不可压缩牛顿流体湍流剪切率与能量耗散率的关系和湍流能量耗散率的近似算法 ,得出了不可压缩牛顿流体在圆管内湍流流动剪切率及脉动流剪切率的实用算式 ,计算所得湍流剪切率及脉动流剪切率的分布特征与已知实验规律相符。     

可压缩k-ξ湍流模式理论研究

基于Boussinesq涡黏性假设,选取新的第二输运变量,推导出k-ξ湍流模式的可压缩形式。针对可压缩流动的特点以及分离流动的复杂性,在可压缩k-ξ湍流模式的基础上添加应力限制器,通过限制涡黏系数的方式来达到限制湍流剪切应力以及减小湍能增幅的效果。将添加剪切应力限制器的k-ξ可压缩湍流模式用于RAE2822翼型流场数值模拟,结果表明:压力系数及不同截面速度与实验值符合较好,且在分离与掺混的复杂流动区域,修正后的k-ξ湍流模式数值模拟结果优于k-ε模式及k-ω模式。     

牛顿流体在圆管内湍流流动剪切率的近似计算

根据剪切率的定义导出了描述湍流流动剪切率与时均流剪切率和脉动流剪切率关系的基本方程,该关系式与流体的流变性类型无关,适用于牛顿流体和非牛顿流体.根据不可压缩牛顿流体湍流剪切率与能量耗散率的关系和湍流能量耗散率的近似算法,得出了不可压缩牛顿流体在圆管内湍流流动剪切率及脉动流剪切率的实用算式,计算所得湍流剪切率及脉动流剪切率的分布特征与已知实验规律相符.     

边界层流动的建模与计算 第二版

本书是第二版，在第一版的基础上增加了可压缩边界层流动的建模与数值计算的内容，主题更加广泛，涵盖二维和三维不可压缩与可压缩层流、湍流和过渡边界层，还论述了在边界层和无粘流动之间的粘性-无粘耦合现象。     

Mach8的平板可压缩湍流边界层直接数值模拟及分析

对来流马赫数等于8,壁温等于10.03倍来流参考温度的平板可压缩湍流边界层做了直接数值模拟,计算涵盖了从层流到转捩以及最终充分发展湍流的全空间演化过程.对湍流的统计特征做了详细的分析,结果表明,在当前的计算工况下,湍流边界层核心区平均速度剖面仍然满足对数率,且卡门常数基本不变;可压缩效应明显增强,由于采用近似恢复温度的等温壁条件,使得近壁区温度较高,导致当地声速增大,使得湍流马赫数绝对值较低,造成内在压缩性效应不强,与经典强雷诺比拟相比,除在数量上产生一些偏差外,强雷诺比拟关系近似成立,且Morkovin假设依然有效;对扩展自相似性和标度率分析表明,对于平板可压缩湍流边界层而言,高马赫数流动使得其适用范围减小;压缩性效应对近壁湍动能,条带结构,涡等值面分布的影响得到分析.     

高精度中心-WENO混合格式在可压缩各向同性湍流大涡模拟中的应用

为提高可压缩湍流大涡模拟的格式精度和分辨率,通过引入当地压力脉动的感应因子和格式加权函数的取值界限,发展了一种数值耗散自适应可控的近似6阶中心-WENO(加权本质无振荡)混合格式,采用傅里叶法对离散格式数学特性进行了理论分析,并对一维激波/密度脉动干涉问题和三维可压缩各向同性湍流大涡模拟问题进行了计算。结果表明:近似6阶中心-WENO混合格式相比于5阶-WENO格式具有更小的耗散误差,且对激波和物理脉动均具有较高的分辨率;基于所发展的中心-WENO混合格式的大涡模拟计算结果与已有的直接数值模拟结果符合较好,且能够成功捕捉-5/3幂律能谱特性曲线;该格式标定了适合于可压缩流动大涡模拟的格式加权函数界限数值,为流体机械内部可压缩湍流的高精度大涡模拟研究奠定了算法基础。     

高精度气体动理学格式与湍流模拟

本文介绍了多维高阶精度气体动理学格式的构造原理及其在可压缩湍流直接数值模拟中的初步应用.通过对初始分布函数和平衡态分布的多维高阶展开,利用BGK方程的通解可以得到单元界面上随时间演化的分布函数,从而建立具有真正多维特性的高阶精度气体动理学有限体积格式.应用建立的格式对可压缩各向同性衰减湍流和槽道湍流进行了数值模拟,得到的流场结果与已有研究相吻合.研究表明,相比于已有二阶气体动理学方法,高阶动理学格式具有更高的求解精度.     

楔形体回流区的射流汇合

本文建立了两股二维不可压缩湍流射流在楔形体后汇合的数学模型,分析了汇合点位置、回流区中压力分布、射流形状和汇合点下游的流动情况,并测量了整个流场的速度分布、湍流能量分布和压力分布.     

剪切应力输运k-ω湍流模型在细长弹体可压缩分离流中的改进

为研究飞机和导弹的细长体可压缩湍流分离流问题,对剪切应力输运(SST)k-ω湍流模型进行了改进。数值模拟了马赫数2.5和0.7、攻角14°下细长旋成体的湍流分离流场,给出了原SST k-ω模型和改进的SST k-ω湍流模型的细长体背风面极限流线、分离涡的强度和位置、物面压力分布的计算结果,并与实验结果进行了对比。结果表明,对计算细长体可压缩分离流动的绕流特性,SST k-ω模型引入的Bradshaw数(雷诺切应力与湍动能之比)应由0.31修正为10/29,修正后的SST k-ω模型与原SST k-ω模型相比,所计算的分离涡的强度和位置、物面压强分布与实验结果更接近。     

高超声速流动中双尺度湍流模式的应用

论文研究了双尺度湍流模式,并对其在壁面附的近长度尺度进行了修正,选择四个基准流动－超声速和高超声速二维压缩拐角,锥柱裙组合体绕流和斜激波／平板湍流边界层干扰－进行了数据计算,数值计算和实验结果的比较表明修正双尺度湍流模式对流动分离,摩阻和热流的计算具有更好的效果。     

数值模拟入射斜激波／平板湍流边界层干扰流动

应用GAO-YONG可压缩湍流模型数值计算了入射斜激波／平板湍流边界层相互干扰现象。计算程序中的对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程。计算中包含了无分离流动、初始分离流动以及较大分离流动等多个情况,比较了平板壁面压力、法向平均速度剖面、壁面摩阻系数Cf以及壁面斯坦顿数St等的计算结果与实验值。结果发现：GAO-YONG可压缩湍流模型能够很好地预测入射斜激波／平板湍流边界层相互干扰下的无分离以及小分离流动,对高马赫数下的大分离流动也能得到较合理的结果,但再附点之后的壁面摩阻系数以及斯坦顿数的计算值不够理想。     

FAP——流动分析通用有限元程序

采用高分辨率有限元格式，建立了流动与输运过程的数值模拟系统。该系统包含范围相当广泛的流动、传热与输运问题，如不可压粘流与可压缩流，自由表面流与浅水环流，传热与相奕问题，渗流以及多相流等。由于采用非结构化网络制备了多种二维与三维单元，系统具有网络生成容易、通用性强等特点。     

二维浅水流动的涡运动学特性及动力学特征

建立了二维浅水流动的旋涡诱导速度公式，在已知旋涡场和胀量场时，由该表达式可确定流动速度场，对比分析了二维浅水流动中的涡运动与可压缩流动中涡运动的异同点。     

超临界二氧化碳BLEVE机理数值模拟

 碳捕集与贮存及提高油田采收率（CCS-EOR）技术是节能减排的有效手段。本文以超临界CO₂为研究对象,采用数值模拟的方法（使用FLUENT 软件及可压缩形式的N-S 方程,采用SIMPLEC 求解有限差分方程,选用标准k-ε湍流模型,多相流模型选用FLUENT 自带的VOF（volume of fluid）模型）,对其在注入油管内发生沸腾液体扩展蒸气云爆炸（BLEVE）事故的初期过程进行模拟计算,建立并验证合理可靠的数学模型,分析容器内流场演化过程及压力、温度的变化过程,揭示超临界CO₂ BLEVE 过程的发生机理与变化规律。     

可压密相气固两相湍流流动的数值模拟方法

在ＤＩＰＳＡＲ算法的基础上,推导了适用于可压密相气固两相湍流流动数据会计的ＣＤＩＰＳＡＲ算法,利用提出的可压密相的两个ｋ－ε低湍流雷诺数模型,用ＣＤＩＰＳＡＲ算法,对竖直上升管中可压密相气固两相湍流流动进行了数值计算,同时考虑了密度变化和体积分数变化的影响,研究结果表明,随着颗粒体积分数的增加,颗料间及颗粒与壁面间的相互作用占了主导地位,使得单个颗粒脉动速度增大,但是气体及颗粒的速度差却随之减小,     

低浓度固液两相流颗粒相本构关系的动理学分析

固液两相流的动理学理论认为固体颗粒相可用Boltzmann方程描述。从流场中的颗粒受力分析出发 ,应用动理学的 Chapman- Enskog法 ,探讨了低浓度固液两相流下的 Boltzmann方程的二阶近似解和颗粒相本构关系。结果表明 ,尽管颗粒所受相间力与两相脉动速度有关 ,但在两相脉动速度近乎无关的极限条件下 ,颗粒速度分布函数和颗粒相本构关系与快速颗粒流的相同。这样的极限条件在含粗重颗粒的两相流中可以近似实现     

MPC模块三分支联接湍流流场的数值模拟

用任意拉格朗日－欧拉（ＡＬＥ）法模拟了ＭＰＣ模块的二维定常可压缩湍流流场．以某一型号发动机的ＭＰＣ排气系统为例，对仅支管进气及支管总管都进气两种情况获得了速度场和压力场，揭示了ＭＰＣ模块三分支联接湍流流动特性，也为用ＡＬＥ法模拟ＭＰＣ模块三维流场奠定了基础．     

高温气冷堆螺旋管直流蒸汽发生器时域模型

螺旋管直流蒸汽发生器(SG)是高温气冷堆核电站的关键部件。为研究SG内氦气、水/蒸汽流动和换热的动态过程,该文建立了SG时域模型并编制了计算程序。其中,水/蒸汽两相流采用一维漂移流模型描述;氦气采用一维、可压缩流动模型描述;两侧流体与管壁的换热系数和流动阻力系数采用经验关系式计算。求解方法采用适用于动态、低速、可压缩流动的压力修正算法,以克服低Mach数造成的数值不稳定。采用此模型计算了THTR-300SG温度分布,计算结果与实验结果相比平均温度误差小于10℃。动态计算结果表明,此模型可以捕捉到规律的两相流脉动。利用此模型可以进行SG热工设计、不稳定性分析以及电站系统的工艺设计。