两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

轴对称射流矢量喷管的试验和数值模拟

对基于激波来实现推力矢量的轴对称射流矢量喷管的缩比模型进行了测力和测压试验，用推广到可计算可压缩流的SIMPLE方法其内外流场进行了数值模拟，根据试验和数值模拟结果分析了喷管主流与次流相互作用产生的复杂流场结构、二次流流量和落压比对气动矢量角的影响，在落压比3～6范围内，二次流流量和喷管主流比值增大，气动矢量角增大，两者比值相同时，落压比增大，气动矢量角减小。     

应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

扫频式射流对涡轮叶栅间隙泄漏流动影响的数值研究

为控制涡轮叶栅中叶顶间隙泄漏流动和改善涡轮气动性能,将扫频式射流器(SJA)作为一种主动流动控制方法应用在涡轮叶栅的研究中。通过非定常数值计算,分析了SJA对涡轮叶栅叶顶间隙流动的作用过程以及作用机理,并且研究了不同工况下SJA对涡轮叶顶流场改善效果以及不同频率的SJA对叶顶流场的影响。结果表明:通过在涡轮叶栅上端壁增加单个SJA装置,可以有效地延迟上端壁的流动分离,其中最佳方案射流流量仅为进口总流量的0.35%,涡轮叶栅出口截面总压损失系数减少了11.48%。存在着最佳的频率284Hz,使SJA装置对流场的作用效果最佳,有效地改善了涡轮叶栅内的间隙流动。     

应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

湍流数据同化技术及应用

近年来数据同化（DA）被引入湍流动力学研究中,通过融合实验测量和数值计算,提高了实验测量的精度和广度,改善了数值模拟的预测性能。实验观测、预测模型和同化算法是数据同化的三要素,湍流研究中的实验观测包括热线风速仪、粒子图像测速法（PIV）、压力传感器等局部测量数据,预测模型主要指流动控制方程及湍流封闭方程,而同化算法包括贝叶斯推断、集合卡尔曼滤波（EnKF）、伴随等。稳态数据同化一般结合雷诺平均Navier-Stokes （RANS）模型方程,从重新标定模型常数、修正涡黏模型方程形式误差、修正雷诺应力项等方面着手。非稳态的数据同化一般包括四维变分（4DVar）等时间连续的数据同化方式以及顺序数据同化。4DVar通过时间正向和逆向积分迭代,存储量和计算量都非常大。顺序数据同化不需要时间逆向积分,可以在若干时刻对实验观测进行间断地植入,正向求解整个系统。另外,随着人工智能的飞速发展,湍流数据同化研究也向智能化迈进。对于纯数据驱动的湍流机器学习,其缺乏物理本质的约束,而基于物理信息的机器学习在物理本质上与数据同化是一致的。     

高速可压缩来流与凹腔驻涡相互作用及其物质输运特性研究

为了深入认知凹腔驻涡与高速可压缩来流之间的相互作用过程以及来流马赫数对上述动力学过程的内在影响规律,对长深比为1.0的单凹腔基础构型在常温、常压、Ma=0.3~0.7的高速可压缩来流下进行冷态非定常数值模拟,并基于拉格朗日视角研究来流与凹腔驻涡相互作用及其物质输运特性。结果表明：凹腔上方剪切层内旋涡结构仍然具有类似自由剪切层的周期性演化过程,主流与凹腔通过三次掺混进行物质输运;随着主流马赫数的增加,旋涡运动方式及凹腔-主流物质掺混机制都发生了本质上的改变,剪切层内旋涡和凹腔主涡涡量均呈现数量级式增长,剪切层旋涡结构脱落频率近似线性增长并逐渐卷吸至凹腔中部,进入凹腔的主流流体占比降低,但扩散程度和速度加快;此外,源于凹腔的流体流出唇口线后经剪切层的夹带作用回流占比增大,降低了凹腔-主流的整体质量交换水平;驻留参数整体衰减,高马赫数时仅为低速工况的1/2,显著降低了回流区对于流体驻留作用的贡献。     

快响应PSP技术用于高超声速圆柱绕流的非定常压力测量

压敏涂料(PSP)技术是飞行器风洞实验大面积定量测压和流动显示的重要工具。为了发展非定常流动压力测量和脉冲风洞中全局压力测量能力,航天空气动力技术研究院开发了快速响应压敏涂料技术,与中科院化学所共同合作开发的压敏涂料采用PtTFPP作为发光基团,稳定性较强,持续光照下发光强度衰减为1.5%/h。采用自主研制的静态标定设备在标定腔内测量不同温度和压力条件下20mm×20mm涂料样片的表面发光强度来确定涂料的压力响应特性和Stern-Volmer公式,并设计制作了2套快响应时间动态标定设备,测得涂料典型响应时间在0.2ms。在中国航天空气动力技术研究院FD-03高超声速风洞中对平板圆柱装置进行了Ma=5的 PSP 试验,利用高速相机采集图像数据后,经过批量数据处理,采集频率为250Hz,得到了间隔时间4ms 的连续压力场数据。结合纹影图像和油流图对得到的PSP结果进行了分析,利用同时采集的测压孔数据对PSP 结果进行了比较。试验结果表明,快响应PSP技术可以更详细地显示流场结构,并能同时得到很好的定量压力数据。