无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

超声速平板边界层旁路转捩直接数值模拟

采用具有七阶精度的有限差分格式和基于Euler方程组的特征边界条件,通过直接数值模拟研究了由固定相位周期性吹吸扰动诱导的、Ma∞=2.25、Re∞=2.5×107/m(635 000/inch)空间发展平板边界层的旁路转捩过程,并建立了相应的直接数值模拟数据库。数值结果表明:扰动相位的随机性并不是转捩发生的必要条件,固定相位周期性吹吸扰动即可诱发边界层的转捩;当采用固定相位时,随着吹吸频率的增加,转捩位置向下游移动;当扰动频率足够大时,周期性吹吸无法诱导出边界层的转捩。随后,在所得数据库中提取了速度及热力学变量的平均值,与已有的数据库和理论分析吻合良好;速度脉动均方根与公认的不可压缩平板边界层实验基本吻合;热力学变量脉动均方根和雷诺应力与已有数据库的统计结果及实验数据具有相同的分布特征。对比发现,不同转捩的诱导因素对平均转捩位置和热力学变量均方根、雷诺应力等高阶统计量的影响较大。     

抽吸腔中心隔板对弯曲壁面边界层抽吸性能影响研究

为了研究抽吸腔内隔板对边界层抽吸性能的影响,在弯曲壁面上均匀分布22个抽吸槽,并在抽吸腔内布置中心隔板,改变抽吸腔出口大小以及入射激波的位置。数值计算采用基于有限体积法的二阶迎风格式来离散二维可压N-S方程,湍流模型采用标准k-ε模型,从流场结构、抽吸流量、弯曲壁面的表面摩擦阻力系数、平均总压恢复系数以及平均马赫数等方面对流场进行了分析。结果表明在弯曲壁面抽吸腔内布置隔板会对抽吸效果产生影响,并且影响程度随抽吸腔出口大小以及激波位置的改变而改变。抽吸腔出口非节流时:布置隔板后,激波位于弯曲壁面中段、后段,抽吸效果均无明显变化,进气道性能亦无明显改变;激波位于弯曲壁面前段时,平均总压恢复系数增加8.33%,质量流量增加3.27%,抽吸效果有所改善,进气道性能也有所增强。抽吸腔出口节流时:布置隔板后,与非节流时相反,当激波位于弯曲壁面前段时,抽吸效果无明显改变,进气道性能不变;当激波位于弯曲壁面中段时,弯曲壁面分离泡的长度减小近2倍,质量流量增加66.15%,抽吸效果显著增强,进气道性能变好;当激波位于弯曲壁面后段时,分离泡长度增加了5倍,分离泡高度增加了近2倍,平均总压恢复系数降低8.33%,同时质量流量也减少1.83%,抽吸效果变差,进气道性能恶化。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

不同钝体宽度下的驻涡燃烧室排放试验

为了研究空气流量分配对驻涡燃烧室对排放特性的影响,了解对驻涡燃烧室内污染物生成的过程及其影响因素,设计了一个能够改变中心钝体宽度、仅凹腔供油的驻涡燃烧室.在常压下对该驻涡燃烧室进行了排放特性试验,进口温度保持200℃.试验中,燃烧室进口马赫数为0.15~0.3.影响排放的因素主要包括雾化质量、凹腔当量比以及与进口马赫数相关的驻留时间等.总体来说雾化质量、凹腔当量比的提高对降低CO和HC的排放是有利的,但是这会使NOx排放增加.在低凹腔当量比时,CO排放曲线变化下降比较平缓,甚至出现上升趋势,而HC排放曲线比较陡峭.这是由于HC的消耗速度比CO消耗速度快,随着凹腔当量比的增加,供油压力提高,燃油雾化粒径变小,燃油蒸发时间缩短,使HC排放快速减少,中间产物CO大量产生而来不及消耗.凹腔当量比进一步上升时,由于燃烧温度的提高,使得CO排放快速减少.在燃烧室内燃烧过程中,NOx的形成和消耗是非常复杂的过程,目前只能作一些定性的分析,而CO和HC的反应过程相对简单.通过对不同钝体槽宽下,具有相似凹腔前壁流量的工况的比较,发现CO和HC的形成主要受凹腔内工作状况影响,而NOx的形成过程更复杂,主流也对其产生着重要的影响.      

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

CARDC 激波风洞 TSP 技术研究进展

从基本原理、关键技术和验证应用三个方面总结了近两年在中国空气动力研究与发展中心激波风洞中开展的温敏涂层(TSP)技术相关研究工作。通过解决快速响应温敏发光材料研制、模型研制、数据处理等一系列关键技术,完成图像采集系统、光学系统及标定系统的配套和系统集成,建立了一套适于激波风洞试验的高速 TSP 测量及标定系统。该技术可在激波风洞试验中获取模型被测面温敏涂层的发光图像,基于该图像可以直接观察模型表面热流分布和捕捉峰值热流的准确位置。结合温敏发光材料的物性参数标定数据,能够实现对模型表面热流的定量测量。不同于传统的传感器点热流测量技术只能得到模型表面有限数量的离散点的热流值,TSP 技术能够以高空间分辨率得到较大面积区域的详细热流分布信息,可更加全面的测量模型外表面的热环境,并且可以据此进一步分析和辨别边界层流态以及确定边界层转捩位置。试验对比表明,TSP 技术的测量结果与点热流传感器的测量结果具有良好的一致性。目前该技术已趋于成熟,在Φ2 m 和Φ0.6 m 激波风洞上成功应用于边界层转捩研究、局部干扰区热环境研究和复杂外形飞行器热环境研究等领域,已成为激波风洞除点测热技术之外又一重要测热技术。     

大气延迟误差对InSAR数据处理影响的定量分析

从重复轨道InSAR测量基本原理出发, 详细给出了相位测量误差对InSAR测高、 双轨法D-InSAR形变测量、三轨法D-InSAR形变测量、四轨法D-InSAR形变测量影响的近似关系式以及大气延迟误差对相位测量影响的近似关系式; 以此为基础分别推导出了大气延迟误差对InSAR测高、双轨法D-InSAR形变测量、三轨法D-InSAR形变测量和四轨法D-InSAR形变测量影响的近似关系式, 同时以ERS-1星载系统为例进行模拟实验, 分析讨论了大气延迟误差对InSAR测高、双轨法D-InSAR形变测量、三轨法D-InSAR形变测量以及四轨法D-InSAR形变测量的影响, 从而得出了大气延迟误差对InSAR数据处理影响的结论.