带挡板圆柱可压缩绕流的大涡模拟

采用大涡模拟方法数值研究了来流马赫数为0.75、雷诺数为2×105的带挡板圆柱可压缩绕流.由于带挡板圆柱绕流的被动控制,与圆柱相比,带挡板圆柱的阻力减小了33%,升力系数的振荡几乎被完全抑制.在此类流动中,由剪切效应引起的减阻大于可压缩效应的影响.通过分析挡板的被动控制机理,我们发现抑制剪切层的失稳可以提高柱体的背压,这与柱体减阻和升力系数振荡抑制密切相关.通过分析沿剪切层的对流马赫数和自激振荡现象,发现在可压缩效应和高频力效应共同作用下,带挡板圆柱的剪切层失稳得到了抑制.     

受损伴流湍流氢气射流火焰的直接数值模拟

采用直接数值模拟DNS的方法对受损伴流湍流氢气射流火焰进行了数值模拟, 采用16步的氢气燃烧详细化学反应机理, 冷的高速H₂/N₂燃料射流喷入热的低速伴随流, 伴随流由贫燃氢气预混火焰燃烧产生, 温度1045 K, 氧量较低. 化学反应源项由主程序在每一时间步长内动态调用CHEMKIN库函数获得. 计算采用消息传递MPI的并行计算方法, 采用12颗CPU在并行计算机上完成. 作为与实验对比的Faver平均结果由DNS瞬态结果做长时间的统计平均后获得. 火焰中涡结构的卷起以及发展过程均能很好地被捕捉, 可以观察到同旋向涡结构之间的相互吸引和反旋向涡结构之间的相互排斥过程, 伴随射流两侧涡结构彼此复杂的吸引、合并、挤压和撕裂过程, 湍流拟序结构由最初的轴对称模式开始向非对称模式演化. 流场中5.67 ms时刻瞬态的H, OH和H₂O分布, 表征了燃料射流自点燃过程中的详细火焰结构. 计算中获取的火焰抬升高度为9d ~ 11d, 与实验结果相吻合. 计算发现由OH和H粒子表征的火焰锋面中, 在火焰锋面转角位置, 燃烧过程得到强化, 可能与火焰面的拉伸以及较长的停留时间有关. 从湍流强度的分布曲线来看, 火焰的传播应该是从两侧向中心发展的. 这里的DNS结果可以作为今后发展更准确通用湍流燃烧模型的参考.     

生物油燃烧与污染物排放特性的数值模拟

应用CFD(computational fluid dynamics)软件FLUENT对实验室自砌小型窑炉装置内生物油的燃烧进行三维数值模拟.获得了炉膛流场、温度场、组分分布、火焰形状及污染物排放的详细信息,揭示了炉膛内部流动、燃烧及传热传质过程的特点.对不同过量空气系数下生物油的燃烧进行数值模拟,结果表明,随着过量空气系数的增加,炉膛内最高温度减小,燃烧区域增大,尾气中CO的浓度减小,但NO的浓度增加.模拟结果与试验数据吻合良好,验证了模型的可靠性,为燃烧工况的优化设计和污染物的控制提供了一定的依据.     

中性大气边界层中风力机的湍流演化及叶根载荷分析

针对一台外场试验风电机组,采用大涡模拟(large eddy simulation, LES)耦合致动线模型的方法,构建了中性大气边界层和风力机风轮的气动耦合求解模型,模拟风力机在中性大气边界层中的流场.通过连续小波分析、频谱分析和相关性分析,研究了中性大气边界层中风力机前、后的湍流演化过程及其与叶根载荷的相关性.研究发现,自然来流从风轮前1D(D为风轮直径)处运动到后1D处时,大气中的湍流强度逐渐增大;在风轮平面处出现了较强的小尺度湍流结构,这些小尺度的湍流结构在向下游运动过程中不断耗散,并在风轮后1D处能量基本耗散殆尽;叶尖位置处的高频湍流出现频率约为1.82 Hz,此频率正好与叶片通过频率相对应.风力机的叶根挥舞载荷对大气中的湍流结构响应明显,低频湍流结构对叶根挥舞载荷的低频段影响显著,高频湍流结构对叶根挥舞载荷的高频段影响明显;叶尖高频湍流结构相对于叶根高频湍流结构,频率更高,能量更大,其对叶根挥舞载荷高频段的影响更为明显;同时,叶尖高频湍流与叶根挥舞载荷的高频部分表现出了一致的周期性变化规律.     

基于时驱硬球算法与格子玻尔兹曼方法的颗粒流体系统直接数值模拟

实现了一种直接数值模拟颗粒流体系统的耦合算法, 颗粒间相互作用由时驱硬球算法描述, 而流体的控制方程采用格子玻尔兹曼方法求解, 流固耦合用浸入运动边界法实现.该方法使用欧拉网格求解流场, 拉格朗日网格跟踪颗粒, 避免了非结构化贴体网格方法需要重新划分网格的问题. 通过模拟两个圆形颗粒在黏性流体中的沉降过程, 成功地复现了经典的Drafting-Kissing-Tumbling(DKT)过程, 验证了耦合算法的有效性.