DES与DDES在湍流分离中的原理与性能研究

随着工程上流动结构的日益复杂，兼具雷诺平均Navier-Stokes(RANS)方法高效率与大涡模拟（LES）高精度的分离涡模拟（DES）类混合方法成为现阶段工程中最有效的湍流模拟方法之一。围绕DES类混合方法中的DES与延迟分离涡模拟（DDES）方法开展工作，分析二者开关函数构造上的不同，研究延迟因子作用机理，并考察DES与DDES方法的求解能力。研究表明：DES与DDES方法在模拟表现上存在一定差异，DDES方法通过引入延迟因子，保护RANS求解区域，改善模化应力不足，降低了DDES方法对交界面系数CDES敏感程度；DDES方法在计算过程中容易出现过度保护，导致求解瞬时涡结构能力不如DES方法，分析与延迟因子引入比重及开关函数构造形式有关。     

气液针栓式喷注器液膜破碎过程和结构的数值研究

为了更好地认识针栓式喷注器雾化场的结构，基于网格自适应加密技术以及VOF (volume of fraction)方法追踪气液的分界面，采用realizable k-ε湍流模型模拟整个流动过程，还原了不同时刻气/液撞击的初次破碎过程，数值模拟结果与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。进一步对针栓式喷注器气/液撞击的初次破碎过程、内部流场涡结构、速度场进行分析，研究了初次破碎雾化的动力学过程和机理。研究结果表明：液桥的形成主要是由液洞的扩展和拉伸、合并而形成，而液滴主要是由中心液膜拉伸、液丝断裂以及液桥断裂而形成，液膜破碎阶段形成的涡结构是造成液膜断裂的主要原因。     

带自由表面水流的分离涡模拟

目前水动力学问题的数值求解仍然广泛采用RANS模型,特别对于大尺度的地表水流的数值模拟,该类模型计算效率较高,但其难以给出湍流场的小尺度涡结构.深入研究复杂地形条件下的水流运动,需要发展诸如LES等的高级数值模型.天然地表水流,地形、边界复杂,雷诺数通常较高,高级数值模拟受限于计算能力的限制,还难以工程应用.RANS和LES混合模型是目前具有工程应用的一种混合数值模型,其已经在CFD领域获得了长足发展.然而对于地表水流运动的数值模拟,相关工作还较少.本文将分离涡模型(DES),即一种RANS和LES的混合模型,应用于带自由表面的地表水流运动,建立了一套数值仿真模型.模型基于有限体积法,水平面内采用非结构计算网格,垂向为结构化网格,对流项离散格式采用二阶TVD格式,并行基于OpenMP语言库.算例表明DES模型有助于揭示复杂地形条件下带自由表面水流的大涡拟序结构.     

端壁组合射流对高速扩压叶栅损失特性的影响

提出了一种端壁组合射流技术以控制进口马赫数0.67的高速扩压叶栅端区流动。通过前缘射流旋涡可以增强端壁附面层与主流间的流体交换,阻碍横向二次流动,减小角区低能流体堆积;而采用角区射流注入能量能够进一步减弱吸力面侧流动分离。以上组合控制方法可较单独采用前缘或角区射流更有效减小栅内损失,提高其气动性能。当角区射流位于近吸力面侧的分离起始位置附近时,其改善栅内流动的效果最佳;远离吸力面的端壁射流则可抑制端区低能流体横向迁移及其与分离区流体间的相互作用,但其减小损失的效果弱于近吸力面侧的射流。随着射流总压比的增加,组合射流减小损失的效果先增加后减小;过大的总压比会加剧射流与来流间的掺混损失,使得叶栅气动性能恶化。当射流总压比为1.2时,损失减小最大可达12.6%,而射流流量仅相当于叶栅进口流量的0.64%。     

翅片位置对复合式减涡器减阻性能影响数值模拟

为了探索翅片-管复合式减涡器的翅片安装位置对共转盘腔径向内流压力损失的影响规律,对不同转速、翅片周向位置及安装角度下的去旋系统开展了数值研究,得到了不同工况下共转盘腔径向内流的流场结构及压力损失分布曲线。研究结果表明:减涡管能引导流体径向流入,并降低流体的旋流比;相比于管式减涡器,翅片-管复合式减涡器能明显降低盘腔内的总压损失;在不同旋转雷诺数下,翅片的周向安装位置α及安装角β均存在最佳值;在中、高旋转雷诺数下,最佳值分别为α=9°,β=30°,最佳结构下总压损失较基础模型低40%左右;改变翅片周向位置及安装角度可以明显改变气流进入减涡管的角度,在较优情况下,可以减小流体流入减涡管的阻力及在减涡管内的流动阻力,整体上减小了盘腔内总压损失。     

流场可压缩性对涡相互作用影响的数值研究

涡相互作用作为冲压发动机喷注装置的典型抽象流动现象,研究其可压缩性影响对于认识包含化学反应的真实燃料喷射场具有一定基础理论价值。基于经过数值验证的可压缩Navier-Stokes算法与压力泊松方程初始条件设置方法研究流场可压缩性对涡相互作用演化过程的影响。结果表明,以较高涡旋马赫数表征的可压缩性在改变涡对形态的同时,具有延缓涡旋相互靠近,迟滞融合进程的作用。关于涡对系统开始融合的临界条件,可压缩相互作用开始的临界展弦比与无量纲时间相比于低速涡对明显提高。为在无量纲时间意义下统一不同马赫数涡对相互作用的进程,从涡心密度随时间变化规律出发在不可压涡对特征时间的基础上,初步构建了考虑可压缩性的时间尺度修正关系。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

室火轰燃的大涡模拟

采用一套修正的适用于低马赫数浮力流的N-S方程组来描述火灾场,运用Smagorinsky亚格子模型的大涡模拟方法对室内池火的轰燃发展过程进行了数值模拟,把模拟与实验结果进行了对比和分析。结果表明,该方法能较好地预测轰燃。同时模拟结果显示火焰区向房间后部偏转稍大,近壁处温度预测的误差大一些,随着亚格子模型的不断完善和计算能力的提高,用该技术来预测轰燃现象会更加准确。     

基于轨道任务几何的"嫦娥五号"采样区选择

"嫦娥五号"(Chang'E-5,CE-5)任务飞行过程复杂、约束条件多,采样区域的选择与轨道设计密不可分.首次提出了基于轨道任务几何的采样区选择方法,以确保探测器任务时序安排紧凑、推进剂消耗少、发射窗口适应范围广、发射机会最多、并满足各个分系统要求.首先给出了工程实施代价最小的可行采样区范围,再综合考虑采样区的地形安...     

基于新型丝状电极等离子体激励器的前体涡控制

本文首次将新型丝状暴露电极DBD等离子激励器应用于大迎角下细长体非对称涡控制.丝状暴露电极的材料的选择对DBD推力以及推力效率至关重要,通过地面精细推力测量对丝状暴露电极等离子体激励器进行了优化,结果表明,本文研究材料中采用钨丝作为暴露电极,其推力效率最优;且随着电极直径从d=0.3 mm减小到d=0.08 mm,DB...