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

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

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

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

可控流场尺度预混湍流燃烧器及其火焰结构分析

为了研究单一湍流场参数对预混湍流火焰结构的影响,以及拓宽湍流场的强度和尺度范围,发展了一套可变结构的预混湍流燃烧器。采用恒温型热线风速仪标定流场,得到了一系列湍流参数。流场标定结果表明:该燃烧器能显著拓宽湍流强度和尺度范围,并能利用不同几何结构产生多种可控流场,实现研究单一湍流参数对湍流燃烧速度和火焰结构影响的目的。选用有代表性的15种湍流孔板组合结构,利用OH-PLIF燃烧激光诊断技术,开展了湍流燃烧实验,结果表明:湍流强度的增大（1 < u'/S_L,0 < 10）使得湍流火焰分区扩展到了薄层反应区,火焰面破碎程度明显增强,孤岛结构明显增多。高宏观雷诺数下,积分尺度的增长对湍流燃烧速度起抑制作用,可能存在临界宏观雷诺数Re_c,能够表现流体惯性力占主导地位的程度,决定积分尺度对湍流燃烧速度的影响效果。积分尺度能量大,扰动能力强,故积分尺度越大,火焰体积越大;但过高的湍流强度会使火焰面褶皱更加剧烈,小尺度叠加在大尺度上的程度增强,最终也使火焰体积显著增大,掩盖了积分尺度对火焰体积的影响,说明积分尺度（表征大尺度）不如湍流强度（表征叠加小尺度的程度）对火焰放热率影响大。     

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

本文采用直接过滤的Navier-Stokes方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似,但仅保留级数的一阶导数项。这样产生的误差相当于丢失了模型的部分耗散作用,本文用一种动力学非线性亚格子模型来补偿丢失的耗散影响。本文根据Caylay-Hamilton定理导出了一个非线性亚格子模型,模型中包含的系数由动力学模式确定。与传统的动力学Smagorinsky模式相比,这种动力学非线性模型(DNM)稳定性更好,且不需要在统计均匀方向进行统计平均来计算模型系数,因此减少了计算开销。本文用这种非线性亚格子模型对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验值、计算值及理论值进行了对比。结果表明,本文模型可以有效地模拟可压缩湍流流场。     

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

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

零质量射流对翼型增升影响的数值模拟研究（下）

针对安装了激振器的NACA0012二维翼型,采用数值模拟的方式对零质量射流技术进行研究。流场解算器采用了我院计算所自主研发的ENSMBV1．2计算平台：计算结果给出了湍流模型对零质量射流数值模拟技术的影响,并且分析了激振器参数对翼型气动力特性的影响规律。     

应用三维粒子动态分析仪测量三分支联接流场

应用三维粒子动态分析仪 (3D PDA)测得了矩形截面三分支联接的三维速度分布。在测量中 ,经过理论分析与对比试验 ,选择了蚊香烟雾作为散射粒子。实验结果表明 ,在只有支管进气的情况下 ,在总管的封闭端和支管下游总管壁面处产生了回流 ;以及气流从支管进入总管后由于截面扩张在总管横截面上所诱发的二次流 ,并且二次流现象只发生在支管下游附近的总管横截面上。     

PDPA技术在圆管射流测量中的应用

采用激光相位多普勒技术（ＰＤＰＡ）和传统热线对圆管射流进行试验研究,试验结果表明：ＰＤＰＡ技术具有精度高、动态响应快、不干扰流场等特点;自由剪切湍流的扩展率ｄｙｃ／ｄｘ＝０．０８０１,完全自模拟特性区在ｘ／Ｄ＞８。     

翼型边界层及近场尾流的湍流特性与气动噪声特性相关规律试验研究

变湍流度试验结果表明 ,湍流度对翼型边界层和尾流的湍流特性、转捩情况、气动噪声特性及其它们间的相关特性有显著影响。这种影响有一定规律。笔者得出了苦干有价值的半经验关系式。