超声速湍流密度脉动预测的神经网络方法

针对气动光学效应对湍流密度脉动预测的需求,发展了超声速湍流密度脉动预测的神经网络方法,从直接数值模拟（DNS）的超声速湍流边界层的流场数据中挖掘规律,建立了包含5个隐含层的密度脉动模型。实验结果表明,所发展的神经网络方法可以很好地预测密度脉动均方值,它不仅能很好地预测训练样本,对测试样本预测的精度和稳定性也显著高于传统模型,且具有一定的泛化能力。通过特征选择和加入先验信息,确定了密度脉动模型的7个输入参数特征量,进一步提高了模型的泛化能力和实用性。     

带凹窗斜劈高速湍流气动光学效应研究

采用数值方法分析研究了带凹窗斜劈高速湍流气动光学特性受窗口深度的影响。湍流数值模拟采用耦合J-B模型的RANS/LES混合方法即DES方法,依靠双时间步LU-SGS方法开展非定常迭代求解,并且利用物理光学方法和波前重构技术等计算气动光学效应,用于分析近场波前畸变和远场光斑分布情况。模拟结果表明,在窗口较浅时,窗口前缘处产生较强激波,带来更大的波前倾斜影响,而窗口较深时,分离剪切层、旋涡等流场结构导致气动光学效应的非定常特性更为显著;高速湍流导致的气动光学效应很强,在所模拟条件下,波前倾斜、均方根光程差、峰谷值以及Strehl比分别达到69μrad、0.22λ、1.2λ和0.31,其对应的跳动值38μrad、0.04λ、0.8λ和0.43。     

三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。     

NASA TrapWing高升力标模数值模拟研究

针对NASA TrapWing高升力全展襟翼构型，采用计算流体力学（ CFD）方法进行三维复杂流场仿真模拟，考察网格尺度和湍流模型对高升力模型气动特性的影响。采用“超立方体”概念，生成绕TrapWing模型的不同网格密度高质量多块结构网格，通过求解雷诺平均Navier-Stokes方程，研究网格尺度对高升力模型气动特性的影响。在此基础上选取中等规模计算网格，考察Spalart Allmaras和Menter k-ωSST湍流模型对高升力全展构型流场模拟能力，分析湍流模型对气动特性的影响。研究结果表明：网格策略具有较好的网格收敛性；湍流模型对翼稍附近位置上翼面压力系数的预测稍有影响，SA湍流模型预测的压力系数较SST更接近实验值。确认研究工作为大型飞机增升装置数值模拟提供了一定的参考。     

钝头双锥喷流致冷流场结构及密度脉动特性

采用大涡模拟方法对钝头双锥喷流致冷流场开展了数值模拟,研究了超声速喷流混合流场结构特征及密度脉动特性。大涡模拟方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runger-Kutta方法。数值模拟清晰地捕捉到了流场波系结构,精细地预测了流动发生失稳、转捩以及发展为充分发展湍流的物理过程,直接获得了流场密度脉动特性。通过有、无喷流状态对称面流场的对比,发现超声速喷流能够有效冷却光学窗口;喷流与主流形成的混合层不稳定,很快发生失稳和转捩,形成大尺度湍流结构,进而引起强烈的密度脉动。此外,获得了钝头双锥整体模型喷流致冷流场的空间发展形态特征。      

湍流等离子体尾迹雷达散射截面的计算及其影响因素分析

针对一阶畸变波Born近似模型,深入分析了湍流内外尺度和电子数密度脉动值对雷达散射截面的影响,并且给出了关于湍流外尺度的一个经验的,能被工程上较好使用的公式.在以上分析的基础上,计算了几种高程条件下再入小钝头锥体等离子尾迹的雷达散射截面,与已有实验结果进行了对比分析.分析和计算结果表明,湍流的外尺度和局部电子数密度值对雷达散射截面值影响较大;湍流内尺度变化的影响不大.     

微观能量分析气动光学效应时变误差的方法

针对高超声速飞行器在大气层内飞行过程中采用自主天文导航时,气动光学效应导致的光学传感器成像畸变严重干扰导航精度的问题,提出一种在微观能量变化机制层面探究气动光学效应时变误差的方法。该方法基于光子传输理论分析光子与湍流分子的相互作用机制,建立光子在高速流场中的传输模型,对光子在湍流传输过程中的能量分布进行统计。通过建立微观光子体系下的气动光学效应评价函数,得到气动光学效应的时变误差描述,统一了微观能量分析与宏观几何光学之间的关系,利用光子微观方法进行数值仿真分析,并将仿真结果与风洞试验结果进行了对比,为高超声速飞行器气动光学效应的研究提供了新的研究思路。     

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景,采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型,采用高精度WENO格式进行空间离散,并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致,来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构,并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析,揭示了流场结构特征及其时空演化规律,并获得了流场密度脉动特性。      

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

高超声速尾迹雷达散射湍流场研究

为给高超声速再入尾迹亚密湍流雷达散射分析提供所需的脉动背景场参数 ,本文提出计算非平衡再入湍流尾迹脉动等离子体场的理论方法。在研究高超声速尾迹流动特征的基础上 ,推导、使用包括化学组份浓度脉动强度的k ε g湍流模型 ,用以封闭高超声速尾迹雷诺平均控制方程 ,并用全隐式有限差分求解。以M∞ =2 1 .2 6、Re∞D =1 .33× 1 0 6的小钝锥体流动为例 ,得到的结果说明 :流场的流向和径向参数分布合理 ;在转捩点后较近距离内湍流脉动影响较大 ,随着向下游流动脉动影响迅速减弱 ;本文计算尾迹湍流脉动等离子体场的方法是可行的     

射流中速度和温度的脉动关联及耗散规律

本文采用X热线-冷线组合探头及全速度偏航角标定法,对非等温圆形射流中的脉动关联量以及脉动量的耗散规律进行了实测和分析。得到了Reynolds应力、脉动热通量和Prandtl数的分布。对耗散率的测量结果证实了由满足自模条件导出的Taylor微尺度和湍动能耗散率在中心线上的变化关系是成立的,并且推导了中心线上湍流温度耗散率的表达式。本文的研究还表明,以往模拟计算中,将脉动温度场的时间尺度与脉动速度场的时间尺度之比取为常数的假定是缺乏实验依据的。     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

加糙平板湍流边界层紊动特性

本文应用偏振差动式激光测速仪对加糙平板湍流边界层紊动特性进行了研究,得到了紊动强度、三阶矩、四阶矩、概率密度分布、功率谱、自相关等特征量。结果表明,可据高阶矩的符号来对湍流边界层沿垂向进行分区,并给出了分区标准,同时发现,加糙使边界层湍流趋于均匀、各向同性,即概率密度分布更接近于高斯分布。     

边界层中湍动能和耗散能最大的尺度分量特征研究

阐述了利用离散正交小波理论分析 1 -D湍流测试信号的方法,使用热线风速仪,在 3个雷诺数下,分别测量了光滑壁面湍流边界层的速度和加速度信号。运用所述方法分析了边界层测试信号,给出了最大湍动能密度和最大耗散能密度及其对应尺度沿 y+ 的分布曲线,并且给出该尺度分量的发生频率沿 y+ 分布曲线。分析可见 :两种最大能量密度的无量纲尺度沿 y+分布具有雷诺相似性,其发生频率随着雷诺数增加而增加,但是最大湍动能密度尺度的频率沿 y+基本是一常数;无量纲的最大湍动能密度沿 y+在边界层小部分区域具有雷诺相似性,而无量纲的最大耗散能密度沿 y+ 在整个边界层内具有雷诺相似性。     

光滑通道内格栅湍流特性实验

为了更好地模拟实际涡轮叶片内冷通道的流动换热,采用格栅对光滑通道内湍流度进行控制。并通过实验的方法对不同格栅尺寸激发的湍流度进行测量。实验中,在边长为80mm×80mm的方形通道中,放置了3种不同尺寸的格栅,利用热线风速仪,得到了通道雷诺数为5000~30000范围内的格栅后下游湍流特性。研究发现:流体通过该格栅后,气流在流经格栅后较短距离内就获得了4.5%的湍流度,同时湍流度沿程呈现衰减趋势,雷诺数对湍流度的影响较小。湍流积分时间尺度与雷诺数呈负相关的关系。      

旋涡运动与湍流模型

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

湍流大涡数值模拟的现状与展望

近年来,随着计算机技术和计算流体力学(CFD)技术的发展,大涡数值模拟(LES)越来越多地应用于解决湍流流动问题.本文介绍了当前LES研究与应用的主要领域以及取得的成果,分析了LES的局限及实施的难点,并对LES未来的研究方向和重点进行了总结.     

任意转速下旋转直通道内的流场数值分析

采用 SIMPLEC算法结合 k-ε湍流模型对各种转速下矩形截面直通道内的流场进行了数值模拟。与文献中的试验值和计算值进行比较后 ,表明本文算法具有良好的收敛稳定性和准确性。计算结果显示 :流动入口条件不变提高转速时 ,在低转速范围内主流速度峰值向压力面移动 ,二次流增强 ,流动具有较强的三维特性 ;在较高转速范围内主流速度峰值向吸力面移动 ,二次流减弱 ,流动由三维特性向二维特性过渡     

低湍流度下翼型边界层及近场尾流中湍流切应力和湍流动能间的相关

简述了在西北工业大学低湍流度风洞中，用热线风速仪对翼型边界层及近场尾流中雷诺正应力及切应力等的测量结果，着重讨论了雷诺切应力及湍流动能的分布规律及相关特性。结果表明，各区段的相关规律并不相同，与常规结果有明显差异。