可压缩钝楔边界层转捩到湍流的直接数值模拟

直接数值模拟了来流Mach数为6的钝楔边界层在特定扰动(音速点附近壁面吹吸)下的转捩到湍流的整个过程. 分析了平均速度剖面, 脉动速度均方根及剪切应力等统计量, 并与不可压理论及实验结果进行了比较. 展示了转捩过程中的涡结构并分析了压力梯度对转捩的影响.     

可压缩气体合成喷射流数值模拟

为了更深入地研究合成喷射流性能与激励器结构参数的关系,在考虑气体可压缩性质的条件下,对合成喷射流进行了理论分析,得出了合成喷射流性能与激励器结构参数的函数关系式,并采用动边界条件模拟实际压电膜的振动,对合成喷外部射流场和激励器腔体内部射流场进行了联合数值模拟.计算结果与已有的实验结果及理论分析结果吻合较好,验证了该数值模拟方法的可行性.     

高超声速圆锥边界层转捩反转数值研究

为研究高超声速边界层转捩中存在的转捩反转现象,基于雷诺平均的N-S方程,采用γ转捩模型对高超声速圆锥体绕流进行数值模拟,得到了圆锥边界层转捩位置随攻角、头部钝度的变化规律,分析了圆锥边界层转捩反转现象的流动机理.结果表明：当圆锥头部钝度增大到一定值时,圆锥边界层转捩位置随攻角增大呈现迎风面提前、背风面推迟的反转现象;圆锥边界层转捩出现攻角反转现象时,迎风面摩阻系数明显大于背风面摩阻系数且转捩前后摩阻变化更大,此时应重点关注迎风面的热防护;在零攻角条件下,随着圆锥头部钝度从很小值逐渐增大,圆锥边界层转捩位置呈现出“N”型反转现象;圆锥头部钝度增大到一定值时出现的边界层转捩反转现象的可能原因是随钝度增大头部激波后熵层高度显著增大引起的.     

气体轴承的动态特性分析及实验研究

针对气体静压主轴回转误差求解不准确的问题,提出了一种利用气膜阻尼模型预测气体轴承动态特性的计算方法,揭示了气膜阻尼影响气体轴承动态特性的内在机制。首先通过微扰动法和分离变量法对气体轴承的气膜阻尼系数建模并求解分析;然后将气膜阻尼系数引入轴承转子系统进行动力学建模,计算分析阻尼系数对气体轴承动态特性的影响;最后搭建气体轴承回转误差测量试验台,得到主轴实时回转跳动误差信号并与理论结果进行对比。实验结果表明:与未考虑气膜阻尼特性的回转跳动误差相比,考虑气膜阻尼特性计算的主轴回转跳动误差值更加接近实验测量值,考虑气膜阻尼特性使主轴跳动误差率减小了4.93%～8%。研究结果可为气体轴承动态特性的预测和精度控制提供理论参考。     

用扰动方程研究可压缩边界层中扰动的演化

本文对传统的可压缩流动中的通量分裂方法进行了修正，使其可以用于扰动方程的计算。利用这一方法，对来流马赫数Ma为4．5、雷诺数Re为10^5的可压缩平板边界层中扰动的演化进行了数值模拟。对于给定的基本流，在计算域入口加入T-S波扰动，研究扰动的空间发展演化。对于小幅值的扰动，计算得到的扰动幅值演化和扰动法向分布与线性稳定性理论的结果符合的很好。对于有限幅值的扰动，结果表明，各次谐波无论幅值还是扰动分布都与完整的Navier-Stocks方程计算的结果符合的很好，但由扰动方程得到的平均流修正比N-S方程得到的大一些，这也许是由于在N-S方程计算时，基本流与零阶谐波一起计算而产生的误差所至。     

雷诺气体润滑方程的适用范围

就前人对可压缩雷诺方程基本假设适用范围的研究结果进行了归纳,指出常规亚音速气体轴承在正常的工作气膜间隙下,雷诺方程的基本假设是成立的.但供气孔、节流气腔以及气膜中可能出现的边界层发展区、超音速区的气体流动则并不满足雷诺方程.     

可压缩效应对平板湍流边界层影响的直接数值模拟

采用高精度紧致型差分格式求解三维可压缩Navier-Stokes方程,通过在局部平板上引入周期性吹吸气小扰动的方法,直接数值模拟来流马赫数M∞为2.25的空间发展的可压缩平板湍流边界层.所得流场统计特征与相关理论和试验符合较好.在此基础上研究可压缩效应对平均流动特征以及湍能的生成和耗散特征的影响,指出在M∞=2.25的情况下,Morkovin假设基本成立,但可压缩效应使得拟序结构有明显差别.与不可压相比上抛、下扫对不同区域的作用有所不同.在近壁区,压力-膨胀项和压力-速度相关项仍是不可忽略的量,这导致声的产生和可压缩效应对湍能产生抑制作用.平板上声压的分布表明,随流动向下游的发展,边界层转捩过程激发了由物面脉动压力所产生的偶极子声源.     

可压缩气体合成射流微喷优化设计

借助Fluent软件,基于标准k-ε湍流模型和PISO算法,采用动网格方法构造动边界条件,在可压缩模型下,对二维、黏性、非定常气体合成射流微喷的频率特性及影响射流性能的3个关键因素(驱动频率、振动膜峰值速度和喷口尺寸)进行数值模拟,为合成射流微喷的优化设计提供了依据.结果表明,在可压缩模型下,合成射流微喷的频率特性仿真结果与理论计算结果较吻合,验证了其在可压缩气体模型下仿真的必要性;同时,对任一给定的合成射流微喷,驱动频率和喷口宽度存在最优值.     

推力波箔轴承的理论分析和试验研究

本文对半径比为0.3、0.5和0.7三种尺寸的推力波箔轴承进行了理论计算。并对其中的二种进行了参数优化。计算结果绘成了便于设计应用的图线,并对设计图线各参数间的关系作了讨论。文中还介绍了该类轴承的试验装置和试验方法,给出了一组试验数据。并将理论结果与试验结果作了比较。最后用实例说明了设计图线的使用方法。     

Gao-Yong湍流模型对边界层转捩的适用性研究

准确预测边界层转捩对深入理解湍流边界层发展和飞行器翼型优化工业应用等具有重要意义。为提升Gao-Yong(G-Y)湍流模型对平板边界层发展及转捩特性预测的准确性与适用性，该文提出了采用G-Y漂移位矢Reynolds数(Re_T)表征当地湍流脉动强度的方法，建立实度系数(Cs)与Re_T间的函数关系，解决G-Y湍流模型中Cs取值间断、分布离散的问题。基于开源OpenFOAM软件开发了改进G-Y湍流模型的求解器，对经典平板边界层和转捩实验开展计算流体力学(computational fluid dynamics,CFD)数值模拟并验证。研究结果表明，相比k-ω等传统雷诺平均(Reynolds averaged Naviers-Stokes,RANS)模型，改进后的G-Y湍流模型不仅能准确模拟边界层流动特征，还能合理预测边界层发展过程及转捩位置，准确性较高。结果表明，经过改进的G-Y湍流模型可进一步用于对边界层转捩特性的研究。     

Analytical solutions to a compressible boundary layer problem with heat transfer

The problem of momentum and heat transfer in a compressible boundary layer behind a thin expansion wave was solved by the application of the similarity transformation and the shooting technique. Utilizing the analytical expression of a two-point boundary value problem for momentum transfer, the energy boundary layer solution was represented as a function of the dimensionless velocity, and as the parameters of the Prandtl number, the velocity ratio, and the temperature ratio.     

双城气田扶余油层高分辨率等时地层格架建立

建立高分辨率等时地层格架对认识地下单砂体展布规律、研究储层非均质性及分析注采矛盾的岩性油气藏勘探和开发具有重要的地质意义。通过对双城气田扶余油层进行高分辨率层序地层学研究,制定了沉积单元划分原则、对比方法;优选出了分层标准井,建立了封闭主干对比剖面。首次进行了约54000井层次的整体、统一、精细的沉积单元划分与对比,把研究区扶余油层划分为3砂组17小层39个沉积单元,从而建立起了高分辨率等时地层格架。     

理想气体条件下平行圆盘止推气体轴承承载力特性研究

在高供气压力和大间隙下,平行圆盘止推气体轴承间隙中呈现出超音速主流区与边界层的复杂作用而导致高度复杂的流场。为探讨超音速主流区的特性,采用绝热理想气体和简化轴承模型分析环境背压不变时气膜间隙中流动状态与供气总压的依赖关系,推导出以供气总压和马赫数为主要变量的不同流态下轴承承载力的计算公式。结果表明,随着供气总压的变化,间隙内可能出现4种不同的流动状态,具体由供气总压与3个特征滞止压强的比值来决定;3个特征滞止压强与出口背压的比值仅取决于圆盘的内外半径比。当供气总压在环境背压和第2特征滞止压强之间时,随着供气总压的增大,轴承承载力从0降至最低负值;当供气总压大于第2特征滞止压强时,轴承承载力随着供气总压的提高而线性增大。     

大涡结构在有压梯度边界层流中的演化机理

采用直接数值模拟方法,研究大涡结构在不同压力梯度作用下的非线性演化特征.计算结果表明:在逆压梯度条件下,大涡结构幅值快速增长,形成很强的流向涡结构,雷诺应力迅速增加,这可能导致猝发现象的产生;与此相反,在顺压梯度边界层流中,大涡结构很难被激发起来,易于流体流动趋于稳定状态.这些结论与某些理论和试验结果一致.     

滑动轴承油膜压力及合金层应力分布

利用有限差分法求解滑动轴承油膜压力的分布;以油膜压力为载荷,建立滑动轴承的三维有限元分析模型,得出滑动轴承合金层应力应变的分布.研究结果表明:滑动轴承应力和应变的分布取决于油膜压力的分布和梯度变化,应力和应变的分布与油膜压力的分布相同,但剪应力的峰值位于滑动轴承中截面合金层与钢背的结合处:应变在压力梯度最大时方向将发生改变.     

最大刚度下单孔静压圆盘止推气体轴承的参数优化

根据已知承载力求最大静刚度的设计准则,采用不可压缩流简化方法,推导出显含圆盘止推轴承结构参数的等效刚度计算式.图解法分析结果表明,轴承的最大刚度点位于设计平面内稳定性约束条件和供气孔直径约束条件的交点;数值优化计算的结果与图解法一致,计算方法可行.     

高速轿车二维风沙流边界层数学模型

基于普朗特的边界层思想,以不可简化成稀相的风沙流体为研究对象,以风沙流的一般方程为基础,在基本假设的条件下,建立了高速轿车车身顶部过流表面的边界层微分方程并给出沙流在近壁表面以滑移为特征的边界条件。又考虑到固体颗粒对边界层的扰动作用,引入了扰动因子。通过对边界层厚度的积分并利用相关条件,得到了边界层动量积分方程。     

低度欠膨胀超音速圆盘止推气体轴承缝隙射流的数值模拟

采用SSTk-ω湍流模型,对高压圆盘气体轴承出流气体形成的低度欠膨胀超音速圆盘缝隙射流流场进行数值模拟。结果表明,低度欠膨胀超音速圆盘射流起始段1的流场波系,与低度欠膨胀超音速轴对称射流类似;起始段2中由于黏性作用逐渐深入核心区,使核心区宽度持续减小。进入亚音速基本段之后,核心区消失,射流对称面上速度持续下降,直至出现滞止区。与此同时,缝隙射流发生规则分离,两股气流偏离滞止区向两侧流动,在射流左右两侧的远场形成规则对称的两个大尺度环状漩涡。