可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...     

燃烧室-涡轮耦合流动传热超大涡模拟研究

为了准确预测发动机燃烧室和涡轮耦合复杂气热环境下涡轮部件的流动和换热特性,应用基于BSL k-ω湍流模型的高精度超大涡模拟方法(VLES),以及SST k-ω雷诺平均湍流模型对雷诺数为Re=3.8×105的高压涡轮导叶在均匀进口条件以及燃烧室-涡轮耦合情况进口条件下的涡轮流动和传热进行了数值研究.耦合火焰面生成流型(F...     

导流槽构型对火箭起飞时流场和声场的影响

为了研究四喷管运载火箭起飞时火箭周围噪声环境问题,建立燃气/空气双组分的可压缩流动模型,采用2阶Roe格式、SAS(scale-adaptive simulation)湍流模型和声学类比积分法Ffowcs-Williams Hawkings(FW-H)求解三维Navier-Stokes方程。以单机火箭的噪声问题为对象开展数值模拟,并将噪声数值计算结果与试验数据对比,误差在3dB以内(相对误差小于1.6%),验证预测噪声方法的有效性,进而研究四喷管运载火箭起飞阶段采用不同导流槽构型对箭体舱段区域噪声环境的影响,结果表明:在相同噪声接收点处,单侧导流槽对应的总声压级比双侧导流槽大,两者的总声压级(OASPL)之差最大为10.7dB。另外,当采用单侧导流槽时,沿周向接收点的噪声总声压相对单侧导流槽中心截面呈对称分布,而且沿导流出口方向逐渐增大。所建立的噪声数值方法为大推力捆绑运载火箭舱段的噪声环境预测及其控制提供一定参考。     

湍流控制屏声学校准及校准精度影响因素

湍流控制屏（TCS）是航空发动机风扇噪声试验必备的试验设施。由于湍流控制屏的作用,风扇前传噪声向远场传播时会产生传递损失。通过对湍流控制屏的声学校准,可以获取湍流控制屏的声学修正量,从而实现对风扇前传噪声的修正。本文结合实际工程应用条件,分析了湍流控制屏重复安装的位置精度、校准声源重复性、稳定性、校准声源位置偏差、温湿度修正等因素对湍流控制屏声学修正量精度的影响;结合上述影响因素,对湍流控制屏160 Hz~40 kHz的声学修正量的特征进行分析,总结工程应用中湍流控制屏高精度声学校准的注意事项,提出湍流控制屏的声学校准应包含测量不确定度。     

燃油喷嘴雾化特性试验研究

为了检验某型航空发动机燃油喷嘴改进设计效果,利用相位多普勒粒子分析仪对燃油喷嘴的雾化性能参数进行试验研究。得到雾化液滴的索太尔平均直径的空间分布、轴向平均速度、脉动速度及其湍流度的分布情况。结果表明:轴向平均速度呈凹盆状分布,脉动速度呈双峰状分布;喷雾中心湍流度大,喷雾边缘湍流度小。随着供油压力增大,在相同测试截面上,喷雾的范围和中心区域粒径变大,边缘位置粒径变小。在相同供油压力下,随着与喷嘴距离的增加,喷雾范围增大,喷雾的轴向平均速度和脉动速度减小,轴向速度的湍流度波动幅度减小。     

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

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

低雷诺数仿生分离流翼型气动特性初探

为了探索适合低雷诺数微型飞行器的翼型形式,基于对自然界鸟类和昆虫滑翔飞行时翅膀形状的观察,设计出一种由前缘削尖平板和后缘圆弧翼型组合而成的仿生分离流翼型。数值研究结果表明,气流在削尖平板的前缘点强制分离,形成大范围低压分离流动,随后在后部圆弧翼上表面再附形成稳定低压涡流区,从而实现较高的气动效率和较强的抵抗大气湍流的能力。上削尖平板可以使流动分离点固定在削尖点。相对于单独平板,仿生分离流翼型的升力系数有大幅提高,迎角为4°时提高了112%。此外,仿生分离流翼型可以在较宽的迎角范围内（4°~20°）保持高升力,但是迎角增加,阻力也快速增大,因此小迎角情况下（小于4°）气动效率更优。      

马赫6柱-裙构型激波/湍流边界层干扰摩阻统计特性

为探究激波/边界层干扰表面摩阻统计特性,对马赫6柱-裙构型的激波/湍流边界层流场进行了直接数值模拟,推导了柱坐标形式的时均摩阻分解公式,将其与对流效应、流向不均匀性效应、分子黏性耗散效应、曲率效应和湍动能耗散关联起来,在此基础上对摩阻脉动和时均摩阻的统计特性进行了研究。概率密度分布结果表明干扰区内摩阻脉动概率密度明显偏离正态分布,摩阻脉动在该区域的间歇性较强,这与湍流脉动在激波干扰后存在显著内在可压缩性效应有关;谱功率密度结果表明摩阻脉动的能量在干扰前中频部分的能量占据主导,峰值对应频率约为0.14,干扰区后附近区域峰值位置转移到高频区域。采用摩阻分解公式对激波干扰前后的时均摩阻进行分解,结果表明时均摩阻的主要贡献来自分子黏性耗散效应和湍动能耗散效应;在激波干扰前分子黏性应力效应占据主导位置,干扰后湍动能耗散效应占据主导位置;激波干扰后对流效应明显增强,曲率效应的贡献在激波干扰前后基本保持不变;流向不均匀项贡献由正贡献变为负贡献,且所占比例有所上升,这主要是由压力梯度项的变化导致的。对分解后各项积分内函数沿法向分布的分析表明,分解后各项在边界层内分布呈现较大的差异,特别地,湍动能耗散项...     

湍流数据同化技术及应用

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

贴壁方柱湍流场TR-PIV实验研究

采用时变粒子图像速度场测试技术(TR-PIV),对低速循环水槽中贴壁二维方柱绕流湍流场进行了细致的测量.通过对实验得到的30000个连续瞬态速度场进行分析,得到了时均速度场和流线图谱,以及流向、法向速度分量的脉动强度场和涡量场.此外,对脉动速度时序信号进行谱分析,得到了流场中低频大尺度相干结构脱落频率.通过对瞬态速度场的分析,揭示了贴壁方柱绕流的旋涡脱落及其发展过程.