二维激波斜壁问题的数值研究

1.前言 在数值模拟激波构造时,通常将流体看作理想流体,即对Euler方程用有限差分法求解。有限差分法的优点是通用性强,但将其用于计算激波问题时,会在不连续波面附近产生数值衰减、过冲或欠冲,难以进行精确计算。FCT(Flux-Corrected Transport Method)法是Boris等开发的一种新的计算方法,它能充分发挥有限差分法的长处,可有效地计算激渡问题。但是,能对二维复杂形状的绕流进行数值研究还是在Thompson等     

固体火箭喷管排气中的粒子分布

1.前言 现代的固体火箭发动机为提高推进剂的能量特性及抑制不稳定燃烧，常常在推进剂中添加金属粉末，最常见的是加金属铝粉。在发动机工作时，铝粉燃烧形成凝聚相的氧化铝粒子。这些粒子在喷管中的流动过程中，温度、速度等方面的变化滞后于燃·气本身的变     

关于风洞中用尖劈和粗糙元模拟大气边界层的讨论

尖劈和粗糙元广泛地应用于风洞试验中的大气边界层模拟，该技术成功模拟了不同地貌特征的平均风速和湍流度剖面。随着风工程研究的深入，了解尖劈和粗糙元模拟过程中的作用机理有助于准确地模拟各种大气边界层湍流功率谱和尺度特性。试验表明：尖劈利用其迎风平板的分离流产生湍流涡旋，迎风板的宽度决定了涡旋的大小和湍流脉动强度，同时迎风板阻塞比沿高度递减产生近似线性的风速剖面；粗糙元用于模拟实际地面的摩擦效应，调整平均风速和湍流度的剖面分布。遗憾的是，尖劈下宽上窄的结构特点决定了该技术模拟的湍流功率谱和积分尺度的高度变化律与实际大气边界层相反。基于对模拟机理的认识，异型尖劈上部形状有助于模拟大比例模型试验要求的湍流风场。     

六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

利用林胶茶人工群落减弱台风对橡胶林破坏机理的风洞模拟试验研究

本文研究台风登陆时（中性大层结）胶茶人工群落内部及上部风速的变化，湍流结构的变化及阻力的变化，来分析判断基减弱台风破坏的机理，评定最佳防风效果的林胶茶人工群落方案。试验研究证明台风吹向胶茶人工群落，台风越过和穿过林带，由于气流与树冠，枝叶的冲击与摩作用，使台风风速大大降低，使动能大大损失，使原来动能强和湍流弱的较大尺度旋涡的台风，变成旋涡尺度较小的动能弱湍流强的湍流，湍流能量的增强反映了风速递减率     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

双钟型喷管高度补偿特性的数值分析

针对双钟型喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象，应用RNGκ-ε湍流模型封闭二维粘性可压缩N-S方程组，采用二阶迎风格式及隐式Gauss-Seidel迭代方法进行耦合求解，对双钟型喷管流场进行数值模拟研究。结果表明，双钟型喷管在低空出现壁面可控流动分离，避免了过膨胀现象，性能优于大面积比钟型喷管；在高空下为完全附着流动，性能优于小面积比钟型喷管，实现了高度补偿特性。     

空腔倾向对泄漏影响的流动分析

给出了前倾斜腔，矩形空腔和后倾斜腔三种结构的四空腔迷密封实验结果。通过数值求解非定常Ｎ－Ｓ方程，预测上述三种结构的单空腔迷宫通道的内部流动，计算结果表明：后倾斜腔结构内部是稳态流动，而对于前斜腔和矩形空腔结构，流动表现为总体结构的基本稳定和局部结构的不稳定，局部不稳定的强弱顺序依次是前倾斜腔，矩形腔和后倾斜腔的迷宫结构。通过分析局部结构不稳定对减小泄漏的作用机制，说明不稳定效应的强弱是造成密封效果     

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。