基于扩展TAB模型的凝胶液滴二次雾化特性研究

为初步研究幂率型凝胶液滴的二次雾化特性,将原始的Taylor-analogy-breakup(TAB)模型扩展并应用到幂率型凝胶颗粒中,采用四阶Runge-Kutta法对凝胶液滴的振荡方程和运动方程进行了数值求解,计算了不同空气动力和物性参数条件下的初始破碎时间和临界特征。结果表明:随着相对速度和液滴直径的增加,初始破碎时间迅速降低,然后保持稳定;随液/气密度比和表面张力系数增大,初始破碎时间呈线性增长趋势;初始破碎时间随稠度系数增大而增大,而当流动指数较小时初始破碎时间变化很小,流动指数超过0.6后初始破碎时间增长迅速,二者与相关实验比较存在一定误差;随Web数增大,液滴的振荡幅度变大,达到稳定后其无量纲变形系数就越大;凝胶液滴的稠度系数越高,临界Webc数越大,液滴二次雾化能力越低。流动指数小于0.6时,临界Webc数变化较小,而后其值则迅速上升。     

横向气流中航空煤油凝胶液滴二次雾化特性实验研究

为研究航空煤油凝胶液滴在横向气流中的二次雾化特性,建立了二次雾化实验台,制备了3种胶凝剂含量下的航空煤油凝胶,在不同实验条件下(82We691,9.6×10-4Oh3.77)对液滴二次雾化过程进行了拍摄。测量了液滴变形、破碎模态、破碎时间、剥离速度及R-T(Rayleigh-Taylor)不稳定波波长等特征参数,并与航空煤油液滴进行了对比。结果表明:在实验条件,煤油和煤油凝胶液滴无量纲变形率在0.5~2.0,胶凝剂含量增加,液滴粘性效应增强,平均无量纲变形率降低;航空煤油和凝胶液滴具有两种破碎模态,同种破碎模态下,煤油和煤油凝胶的破碎过程相同;煤油和煤油凝胶液滴初始和总破碎时间随We数变化规律基本相同。不同Oh数下,3%煤油凝胶初始和总破碎时间与文献计算值相差较大,而其余3种材料液滴则吻合较好;液滴剥离速度随We数逐渐增加。We数低于200时,煤油剥离速度高于煤油凝胶,但继续增大We数时则相差不大;R-T不稳定波波长随We数增加而减小。煤油、1%和2%煤油凝胶波长实验值和理论值基本吻合,而对3%煤油凝胶,理论值则高于实验值。     

旋流雾化喷嘴雾化特性的数值模拟

为了掌握旋流雾化喷嘴的雾化特性,基于CFD软件建立了该喷嘴的几何模型进行数值模拟,在验证了网格无关性和数值模拟结果可靠性之后,通过VOF和Realizable k-ε模型对不同压差工作条件下和不同结构的喷嘴分别进行数值模拟。结果表明:旋流雾化喷嘴的雾化流量、雾化锥角和流量系数会随着进、出口压差的增大而增大;喷嘴的旋流强度会随着喷嘴进、出口压差的增大而减小;当喷嘴的出口直径和出口长度增大时,喷嘴的各项雾化参数与理论值更接近,相应的雾化效果变好。     

二维涡轮叶栅流场的直接数值模拟

采用高精度有限差分格式求解非定常N-S方程组,对低雷诺数下二维涡轮叶栅流动进行了直接数值模拟,计算了雷诺数为10000,VKI涡轮叶栅在0°,8°以及-8°攻角下的流场,对涡轮叶栅非定常流动机理做了初步的探讨。计算结果表明:在叶栅尾缘处,逆时针方向和顺时针方向的主涡交替在壁面产生,并和主流相互作用产生二次涡,而当二次涡与主流连通发生掺混时,将引起主涡被分割并从叶片表面脱落;攻角在一定范围内的变化对VKI涡轮叶片表面边界层发展影响不明显。文中还对尾迹区的统计量特性和速度亏损特性等进行了研究。      

直流互击式喷注单元雾化特性准直接数值模拟

为实现直流互击式喷注单元雾化过程的数值求解并探究结构及工作参数对雾化特性的影响规律,基于一种树形自适应加密算法与流体容积(Volume-of-Fluid,VOF)方法实现了雾化过程的准直接数值模拟。计算得到了两股射流由喷射到撞击形成液膜,液膜进一步破碎形成液丝、液滴的全过程,获得了液膜的破碎长度、液滴的Sauter平均直径、雾化频率等雾化特性参数。通过将典型算例的计算结果与试验数据进行对比验证了计算的有效性,给出了数值求解精度。探讨了撞击波的形成机理,分析了雾场液滴的尺寸分布规律,撞击夹角、孔径比、射流速度、动量比对雾化特性的影响规律。结果表明:所采用的算法可以实现多相、多尺度雾化过程的数值求解;撞击波的形成是由于两股射流撞击时惯性力的不完全对称导致的;雾场液滴的尺寸分布近似服从Rosin-Rammler分布;撞击夹角的增大与射流速度的提高导致液膜的破碎长度减小,液滴的平均粒径减小,撞击夹角增大雾化频率呈减小的趋势,射流速度提高雾化频率呈增大的趋势;动量比为1而孔径比不为1时会形成凹形液膜,雾场存在一定程度的偏斜;动量比主要影响雾场的偏斜程度。     

叶片通道内二次涡的数值模拟

采用三维Ｅｕｌｅｒ方程数值模拟叶片通道内流向二次涡,及其对叶片出口气流角的影响。算例表明：采用这种方法可以计算出流向二次涡;该涡造成的出口气流偏转角与实验值比较,在远离壁面处相差较小。此方法可用于叶栅通道内流向二次涡的初步分析     

涡轮叶尖泄漏流主动控制数值模拟研究(英文)

在涡轮设计中一直追求减小损失、提高效率,流动控制作为一种减小损失的一种有效手段,在涡轮叶尖泄漏流动方面引起广泛的关注,因此本文对在涡轮叶尖表面进行主动射流对叶尖泄漏流动的影响作了数值模拟。结合基于密度修正的采用雷诺应力湍流模型加壁面函数的三维计算流体力学程序,详细分析不同叶尖间隙高度下射流孔位置和射流流量改变对间隙泄漏流场的影响,并相应地阐述射流主动控制的机理,分析了主动射流条件下叶尖间隙流动的拓扑结构,最后计算考虑射流影响的涡轮转子效率。结果表明:在涡轮叶尖表面选择合适射流孔位置进行射流可以提高涡轮转子效率,其中大间隙下通过射流孔组合射流可以提高涡轮效率0.35%,小间隙下可以提高0.30;在射流孔区域靠近叶片表面处流场结构中鞍点数和结点数相等,若不考虑周围流场的拓扑结构,则射流孔附近的流场满足奇点分布规律。     

复合材料十字型接头力学性能实验研究与数值模拟(英文)

通过拉伸和剪切试验,研究了3种不同的树脂转型模型(RTM)成型工艺(即RTM、引入缝纫的RTM和共胶接)对复合材料十字型接头力学性能的影响。根据实验现象,分析和讨论了其破坏机理。实验结果表明,3种工艺中,RTM成型的十字型接头力学性能最优,胶接和缝纫工艺都会对接头的强度产生不同程度的影响。通过有限元分析,数值模拟了RTM十字型接头的破坏过程,与实验现象吻合良好。     

用双时间法数值模拟悬停旋翼流场(英文)

本文用改进的双时间法和重叠网格通过求解可压缩的欧拉和Ｎ－Ｓ方程计算了悬停旋翼流场。在双时间法中,我们采用牛顿松弛以加速每个真实时间步中的子迭代。高效的双时间法加上重叠网格使在微机上计算复杂的悬停旋翼流场成为可能,计算结果和实验符合得很好,计算得到的下洗流和载荷证明了我们的方法是可行的。     

民用飞机发动机舱灭火剂雾化数值模拟

针对民用飞机灭火系统的特点。使用商业软件FLUENT对民用飞机发动机舱着火后灭火剂的雾化过程进行数值仿真模拟．分析发动机舱着火状态下灭火剂喷射后舱内灭火剂浓度的分布情况,为民用飞机发动机灭火系统设计提供技术支持。     

结冰风洞中液滴相变效应数值模拟

为明晰结冰风洞中液滴相变效应,发展了基于Euler法的气液两相传质传热耦合流动计算方法,探索了3 m×2 m结冰风洞主试验段构型相变效应对液滴传热过程的影响,开展了参数影响研究,评估了试验段内液滴过冷状态。结果表明：构型内液滴经历了先蒸发后凝结两个阶段,蒸发效应促进了准一维传热阶段液滴温度的下降趋势,使液滴温度趋于湿球温度,而凝结效应则抑制了三维收缩阶段液滴温度的下降趋势,进而增大了试验段液气温差,影响液滴过冷状态;增大初始相对湿度和试验段气流速度,会导致蒸发效应减弱而凝结效应增强,进而增强了相对湿度和气流速度对液滴过冷状态的影响程度,与此相反,初始液滴尺寸的增加,则会导致蒸发效应和凝结效应均减弱,进而减弱了液滴尺寸对液滴过冷状态的影响程度;在典型工况下,小尺寸液滴（液滴直径范围为40~100 μm）在高风速时（试验段气流速度大于164 m/s）将偏离过冷状态（液气温差超过3℃）。     

液滴喷射过程流场和温度场的数值模拟

 建立了液滴喷射过程流场和温度场的计算模型,研究了在外加周期性正弦扰动作用下,石蜡液滴形成过程的流场和温度场及其变化情况,给出了不同扰动频率下所形成液滴的形态及其变化规律,得出了比较均匀的液滴流,指出扰动频率对液滴的形成及其均匀性具有重要影响。模拟结果与文献中的实验结果吻合较好。     

胶质微推进液滴加速过程的数值模拟

采用拉格朗日粒子跟踪法计算液滴的运动轨迹,对胶质推进器中带电液滴在加速区的运动过程进行了数值模拟,获得液滴速度和密度的空间分布。使用两种溶剂(TBP3,甲酰胺)作为推进剂,分析了电压、流量变化对喷雾扩散角、液滴轴向速度及对推进器性能的影响。模拟结果表明,随着电压的增加,液滴轴向速度增大,喷雾扩散角减小;随着流量的增大,液滴轴向速度和喷雾扩散角均减小;使用甲酰胺推进剂比使用TP3推进剂得到更高的推力和比冲;模拟得到的推进器性能能良好符合实验结果。     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

盐指型双扩散湍流二维直接数值模拟研究

针对海洋相关盐指型双扩散湍流开展了二维直接数值模拟研究,特别关注了四种不同物性参数组合情况下,传输特性和流动结构随瑞利数和密度比的变化趋势。结果表明：温度努塞尔特数、盐度努塞尔特数和流动雷诺数随瑞利数的整体变化规律受密度比的影响较小,但其具体数值随密度比的增加而减小;密度比的影响幅度随着组分扩散率比值,即刘易斯数的增大而减小;相关变化规律可由幂律模型描述,且基于数值模拟结果本文确定了四种物性参数组合所对应的幂律指数。对于大多数算例,流动主导结构为盐指。对具有小刘易斯数和大瑞利数的数个算例,主导结构变为类对流涡结构。对所有物性参数组合,盐指水平宽度随瑞利数变化趋势与线性稳定性分析给出的结果类似。这些结果对理解海洋盐指双扩散湍流、以及实验和模拟结果应用于真实海洋流动等方面有重要的意义。     

超声速气流中液体燃料雾化数值模拟

为了模拟超声速气流中液体射流的雾化,发展了一种液滴破碎混合模型。该模型将Kelvin-Helmhotz(K-H)模型和Rayleigh-Taylor(R-T)波动模型耦合在一起,首先计算R-T波增长模型,当R-T模型不能导致液滴破碎时,然后计算K-H波增长模型。利用该模型对来流Ma=1.94的超声速气流中的水射流的雾化进行了数值模拟。模拟结果再现了超声速气流中射流的雾化结构,计算得到的射流穿透深度、颗粒直径分布和速度分布与实验结果吻合较好。     

塞式喷管底部二次流的数值模拟研究

从 N-S方程出发 ,采用二阶精度的 NND格式对塞式喷管的流场进行了数值模拟 ,重点研究了一线性塞式喷管底部二次流在低中高空的特性。研究表明 ,塞式喷管底部回流区的消失不仅与二次流流量大小有关 ,而且与底部压强有关 ;加入二次流会使塞式喷管的推力上升 ,但在大多数情况下其比冲会下降 ,只有特定压强比条件下 ,加入少量的二次流比冲才会上升     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

翼型结冰过程中过冷液滴运动的数值模拟研究

本文采用拉各朗日法、标准四步龙格一库塔格式和Euler方程流场数值模拟等方法实现过冷液滴运动轨迹与碰撞的数值模拟，分析了液滴直径、推进时间间隔、投放位置、流场特征等因素的影响。     

气体介质中液滴破碎的LES/VOF数值模拟

为进一步揭示气体介质中液滴破碎的机理与影响因素,采用LES与VOF相结合的方法对液滴破碎过程进行了数值模拟,计算结果与实验结果符合较好,并且优于前人RANS模拟的计算结果。在此基础上,对多工况液滴破碎进行了LES/VOF模拟,分析计算结果发现,液滴破碎初期,液滴的尾涡效应影响液滴的初始变形位置,在变形过程中,液滴表面压力分布直接影响液滴的破碎状态;另外,在一定初始速度范围内,液滴破碎过程保持一致,但初始速度对液滴变形破碎时间及强度影响明显。