建筑物对污染物扩散影响的数值与风洞模拟研究

采用k-ε(RNG)与RSM湍流模型对处于方形建筑物不同位置污染源所排放污染物的扩散规律以及建筑物对流场的影响进行模拟,并且与相应的风洞试验结果进行了比较.流场分析结果表明:数值模拟能够较好地模拟建筑物前方迎风侧停滞回流、顶部回流以及后方空腔区等.浓度场分析结果表明:建筑物前方迎风侧以及顶部回流区污染物扩散的数值模拟结果与风洞试验结果基本相等,而在建筑物后方空腔区污染物的数值模拟结果略高于风洞试验结果.综合分析并与风洞试验结果相比较,RSM模型能够较好地模拟污染物浓度场以及建筑物周围流场的变化规律.     

高速流非平行边界层稳定性研究

导出了可压缩流的抛物化稳定性方程（PSE）。针对高速流动,特别是超声速和高超声速流动的非平行边界层稳定性问题进行了研究。引入高效的边界层变换、全流场高精度的差分格式及预估校正迭代和推进求解法来求解PSE方程,使得PSE方法中至关重要的正规化条件得到了满足,确保了数值计算的稳定。采用高马赫数下对稳定性起支配作用的第二模式,研究了高速流边界层稳定性的演变和特征,分析了流动的非平行性、压缩性,以及壁面冷却等因素对流动稳定性的影响．所得结果与相关实验数据吻合较好。     

基于CBS有限元的流动-传热-变形耦合计算方法

研究了基于CBS有限元法和常规有限元法相结合的流动-传热-变形耦合计算方法。在该方法中,流体流动和传热采用CBS有限元方法计算,固体变形采用常规的有限元方法计算,实现了流体域和固体域统一的有限元网格划分,简化了变形过程中的网格生成和不同网格间的数据交换。然后,依据此方法发展了计算程序,并通过算例分析,校验了计算方法的可行性和程序的计算能力。     

超临界、大折转角扩压叶型的气动性能

采用数值模拟的方法,对高亚声进口条件下,具有60°几何折转角的2D扩压叶型的气动性能进行了研究.通过对不同几何参数条件下叶型气动性能的对比,总结了中弧线进口曲率、最大厚度位置以及最大厚度对叶型性能的影响.在此基础之上,对在进口马赫数达到0.83以上,仍能保持较低损失以及攻角特性的超临界流场特点进行了分析,总结了该种叶型损失机理以及流场结构特点.结果表明,本文的造型能够实现更高超临界马赫数下的大折转角扩压流动.叶型损失主要来源于激波与附面层的相互作用.     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

高焓电弧风洞试验热化学非平衡流场数值模拟

针对高焓电弧风洞内部流动的热化学非平衡效应及气体组分和振动能量冻结效应导致的试验数据外推困难问题,基于高焓风洞喷管/试验段/试验模型一体化数值模拟的思路,通过数值求解三维热化学非平衡Navier-Stokes方程,开展了FD-15高焓电弧风洞典型运行状态下流场的数值模拟,与典型试验状态的气动热数据进行了对比验证,研究了试验数据外推飞行条件的方法及有效性问题,分析了提高驻室总压对试验数据外推的影响。研究表明:（1）风洞试验段来流离解度高,热化学非平衡效应及其冻结现象严重;（2）热流校核试验测量数据位于一体化数值模拟的完全催化热流和非催化热流之间,分布合理,验证了计算方法和程序的正确性;（3）试验模型安放位置对模型表面压力和热流存在影响,模型与喷管出口的距离越大,模型表面压力和热流越低;（4）当驻室总压较低时,通过双尺度模拟准则（模拟飞行条件总焓和双尺度参数ρ_∞L）外推热流失效,使用部分模拟准则（模拟飞行条件总焓和驻点压力）外推热流也会出现较大差异,在非催化条件下这一现象更加明显;（5）当驻室总压较高时,使用双尺度模拟准则或部分模拟准则外推飞行条件,产生的热流差异明显减小。     

氮气辐射强度的激波管测量与验证

超高速飞行器的热防护设计必须考虑激波层内高温气体发射与吸收的辐射能量,需要有效的辐射加热评估手段。相应飞行条件下的光谱辐射强度地面实验测量是验证数值模型和方法、理解高焓流动的重要手段。基于燃烧驱动激波管,发展辐射强度标定技术,针对富氮气环境,开展高温气体光谱辐射强度的高分辨定量化测试,掌握辐射特征,为数值验证提供基础数据。实验获得了激波速度5.70和6.20km/s条件下的气体光谱辐射强度精细结构,数据表明激波波后的非平衡过程对辐射强度存在很大影响。通过求解耦合化学反应动力学模型的Navier-Stokes方程和辐射特性模型,得到对应实验条件下的流场参数和辐射强度,计算结果和实验数据符合很好,验证了数值模拟方法。     

湍流燃烧模拟中化学反应的加速算法研究进展

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法,讨论了动态自适应化学（Dynamic Adaptive Chemistry,DAC）方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中,使用DAC简化可以加速化学反应计算。然而,在并行燃烧数值模拟中,处理器核心的负载极度不平衡,加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心,更有利于整体计算效率的提高。在同等精度下,相比于隐式格式耦合DAC和MTS加速方法,Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。     

一种非介入式高超声速边界层不稳定波的测量方法

地面风洞实验是开展高超声速边界层转捩研究的主要手段之一,但是目前可用于高超声速边界层三维空间测量的实验技术仍极为缺乏,且已有测量技术的动态响应频率普遍较低。基于光的折射和干涉原理,搭建了一套非介入式聚焦激光差分干涉仪测量系统（Focused Laser Differential Interferometry,FLDI）,可有效获取三维流场空间点的密度变化。在马赫数为8的常规高超声速风洞中,使用FLDI开展了来流雷诺数107/m、7°半锥角尖锥标模边界层的不稳定波测量实验。结果显示FLDI成功捕获到频率在327 kHz的第二模态不稳定波及其谐波（645 kHz）。通过与PCB测试结果进行对比,FLDI的高信噪比、高解析频率（本文实验有效解析频率1.5 MHz）、高空间分辨率（沿流向小于1 mm）等优点得以体现。鉴于FLDI的高时空分辨率等优良特性,其可用于高超声速边界层不稳定波行为以及感受性等问题的研究,为深入认识高超声速边界层转捩机制以及感受性问题提供了有效手段。     

湍流边界层外区超大尺度相干结构相位平均波形

用IFA300恒温式热线风速仪和X形二分量热线探针精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、法向速度分量的时间序列信号。通过对不同速度分量进行子波变换,得到了近壁区和外边界区流向和法向脉动速度分量的分尺度湍动能随尺度的分布,在外边界区发现流向脉动速度和法向脉动速度的能量最大尺度随法向坐标远离壁面增大,相干结构的平均猝发周期也相应增加,说明在外区存在超大尺度的相干结构。采用条件采样和相位平均技术,提取了近壁区和外边界区能量最大尺度相干结构流向及法向速度的喷射和扫掠的条件相位平均波形,发现流向和法向脉动速度及雷诺应力喷射和扫掠事件的条件相位波形较近壁区都发生了改变。雷诺应力负的幅值显著降低,雷诺应力表现为正负交替的波动现象。