超声速气流中液滴破碎的大涡模拟

为了研究超声速气流中的液滴破碎过程,采用基于CLSVOF界面追踪方法的大涡模拟方法。该方法中用可压流和不可压流求解器分别计算超声速气流流场和液体流场,并且在分别计算气相或液相时,另一相作为该相计算时的边界条件参与计算,从而实现可压流求解器和不可压求解的融合。通过去散度化液体速度外推解决界面处大液气密度比所带来的动量误差大的问题。利用该方法对来流Ma=1.358的超声速气流中的水滴破碎进行了数值模拟。模拟结果展示了超声速气流中液滴破碎的几何结构变化细节和基本破碎形态,且与实验吻合较好。数值模拟结果揭示了Rayleigh-Taylor不稳定和气动剪切力对液滴破碎的作用过程。破碎过程中,液滴在垂直来流方向上的最大无量纲宽度(Dmax/D0)约为4.5,无量纲破碎时间(tbU∞/D0)约为29.85。     

大型商用运输机机翼增升构型水滴撞击特性计算

基于欧拉-拉格朗日数值模拟框架,发展了一套适合求解三维多体模型水滴撞击特性的计算方法.采用粒子统计与面积比率计算相结合的方式,实现了水滴收集率的高效精确预测.在常规粒径条件下,该方法计算结果与国外标模试验数据吻合良好.应用上述算法,开展了某大型商用运输机增升构型的水滴撞击特性数值模拟研究,计算结果显示:水滴在缝翼背面缝...     

过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.      

背压对撞击式喷嘴雾化特性影响研究

为研究背压对撞击式喷嘴雾化特性的影响,将压力的变化等效为气体密度的变化,基于一种树形自适应加密算法,通过直接数值求解不可压Navier-Stokes方程组实现了不同背压条件下射流撞击雾化的数值模拟。首先将数值模拟结果与试验数据进行对比,验证数值模拟的有效性,在此基础上开展了高背压条件下雾化过程的数值模拟。结果表明,随着背压的提高,气动力相应增强,液膜的破碎更加剧烈,一次雾化区域的液滴数密度增大,雾场由稀疏向稠密发展;液膜在向下游运动过程中波动速度的幅值逐渐增大,并且随着背压的提高,液膜的波动由线性向非线性转变;背压增大导致破碎长度减小,经过参数修正得到了液膜破碎长度的经验公式,并与试验数据进行了对比;背压对液滴尺寸分布规律没有显著影响,但随着背压的提高,同一时刻大液滴所占的比例提高,整个雾场的Sauter平均直径有增加的趋势,当背压从0.1MPa增大到1MPa时,雾场的Sauter平均直径由155.5μm增大到166.9μm;背压增大,液滴粒径分布的均匀度指数减小,液滴尺寸分布更加不均匀。     

基于SPH方法的射流撞击仿真

为了研究液体推进剂的撞击雾化问题,采取光滑粒子流体动力学方法 (SPH)进行了射流撞击形成液膜过程的数值模拟。修正了基于连续表面力模型的表面张力算法,推导了考虑粘性与表面张力作用的SPH形式的控制方程,对撞击雾化模型进行了合理简化,采用水作为推进剂模拟液。得到了不同撞击速度、角度下的液膜图像及液膜的形成过程图像。证明了SPH方法适于解决存在自由表面及大变形的射流撞击问题。     

幂律型流体雾化SPH方法数值分析

幂律型流体的雾化过程存在复杂的界面运动,用传统网格法很难精确追踪运动界面.为研究幂律型流体的雾化特性,运用SPH方法对典型的双股幂律型流体撞击雾化问题进行三维数值分析.根据文献实验,采用SPH方法模拟获得与实验条件相应的数值结果,对比后表明,二者雾化角相当吻合,数值结果还成功捕捉到液膜向液丝的破碎过程及网状的液丝分布状态,验证了方法的有效性.分析了射流速度和撞击角度对雾化角的影响,得到雾化角随着射流速度和撞击角度的增加而增大;对雾化后的速度场进行数值分析后表明,撞击点附近,惯性力的作用使速度场变化剧烈;在流体远离撞击点的过程中,粘性耗散作用使速度场趋于稳定,但速度大小小于初始撞击速度.     

考虑遮蔽区影响的旋翼三维水滴撞击特性计算新方法

针对直升机旋翼三维黏性流场特有的复杂环境,建立了一种基于欧拉法的旋翼三维水滴撞击特性计算的新方法。首先,在旋翼桨叶嵌套网格的基础上,发展了一套用于预测旋翼绕流流场的计算流体力学(CFD)模拟方法。然后,为克服传统直升机旋翼二维水滴撞击特性计算方法的不足,充分考虑旋翼流场的三维效应,在嵌套网格中基于欧拉法求解旋翼三维水滴撞击流场。其中,为解决尾流等区域的密度脉冲现象所引起的稳定性和收敛性问题,提出并建立了遮蔽区扩散模型。该模型通过判断遮蔽区变量,在计算域中动态生成遮蔽区域,并随迭代步数逐渐扩散。最后,通过与NACA0012翼型及国外UH-1H桨叶的试验和计算结果的对比,验证了旋翼三维水滴撞击特性计算新方法的可靠性,并进行了温度和水滴当量直径(MVD)对旋翼三维水滴撞击特性的影响分析。结果表明:遮蔽区扩散模型的加入,使二维情况的计算时间减少了22%,并增加了三维情况的计算稳定性,显著提高了旋翼三维水滴撞击特性的计算效率;沿着旋翼桨叶展向位置增大的方向,旋翼桨叶剖面水滴撞击范围有所增大,最大水滴局部收集系数呈先增加后减少再增加的变化趋势,其变化幅度接近50%;旋翼桨叶表面的水滴撞击区域和水滴局部收集系数随水滴当量直径的增加而增加。     

旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。     

预处理方法在低速流场数值模拟中的应用

自20世纪70年代以来，随着高分辨率差分格式和高计算效率时间推进格式的建立和发展，可压缩流动数值模拟取得了很大进展。并广泛采用时间推进法。但是。对于低速不可压流动，采用时间推进法效果并不太好。目前对低速不可压流动的数值模拟主要有三种途径：（1）基于压力的算法，如著名的SIMPLE算法（压力偶合方程组的半隐式方法）；     

发动机旋转表面水滴撞击特性数值研究

介绍了发动机旋转机械内部过冷水滴轨迹的三维数值计算方法,阐述了旋转坐标系下气流及粒子的运动规律。采用ANSYS-CFX软件及其粒子输运模型,对某发动机风扇转子叶片外围空气及水滴流场进行了数值模拟。利用水滴的速度矢量、撞击区域等参数表征水滴的撞击特性,获得了转速、水滴直径等对风扇叶片表面水滴撞击特性的影响:水滴撞击区域集中在风扇叶片迎风面叶盆侧,且水滴撞击区域随着转速的增加而减小;水滴在叶片表面的撞击范围随着水滴直径的增大而减小。     

SPH方法对气液两相流自由界面运动的追踪模拟

气液两相流动是自然界中常见的流动现象,对其进行数值模拟时要求必须能够准确跟踪界面运动变形。本文利用光滑粒子流体动力学方法,结合微可压缩模型(SCM),引入界面控制方法XSPH速度修正以及Van Derwaals状态方程修正,对典型的气液两相流动如二维溃坝、气泡上浮等问题进行了数值模拟,分析了空气和水相互作用机理以及界面运动规律,并同实验结果和其他数值计算结果进行了比较。结果表明,该方法在模拟多相介质界面运动问题准确有效,可用于处理更为复杂的多相流动工程问题。     

燃气排导箱三维场的数值计算

对燃气排导箱内的典型三维稳态流动过程进行了粘性不可压数值计算,采用较为成熟的ＳＩＭＰＬＥ方法。所发展的计算程序可求解质量、动量和能量守恒方程,ｋ－ε双方程湍流模型方程,以及热辐射通量方程,并进行了程序的算例验证。三维燃气排导箱内的计算结果也得到了一些实验测量的数据的验证,得出一些有工程实用价值的结果和结论。     

新型等离子体激励器的数值模拟

为了对新型多电极对等离子体激励器进行可靠的数值模拟,通过求解带有体积力源项的Navier—Stokes方程,模拟单电极对等离子体激励器在二维不可压流中对周围流场的影响,计算结果表明所使用唯象模型的参数基本正确。在此基础上,改进原有模型,对新型双极性等离子体激励器诱导流动进行计算。通过与实验结果的对比,发现改进后的模型能较好地模拟新型双极性等离子体激励器的控制作用,可为今后新型等离子体激励器的应用提供可靠的数值模拟方法。     

喷雾液滴撞击液膜影响参数及流动机理分析

液滴撞击壁面是喷雾冷却中最常见的现象。使用计算流体力学软件中的Volume of Fluid （VOF）模型对液滴撞击壁面过程进行了数值模拟研究，研究了不同参数的液滴撞击壁面液膜过程的水花形态特点，针对水花高度和直径等参数进行了分析，通过运动间断理论和无量纲参数Re和We的变化分析了不同参数对液膜的流动特征和水花形态的影响及水花产生机理。结果表明：大液滴撞击液膜时能够改善液膜的流动状态；过大的液滴速度会产生大量飞溅及表面干涸；液滴撞击薄液膜时，液膜的流动性较好，当0.6<h*<1.2时，液膜的流动性不随液膜厚度改变。     

LES of film cooling for different jet fluids

The present paper investigates the impact of the velocity and density ratio on the turbulent mixing process in gas turbine blade film cooling.A cooling fluid is injected from an inclined pipe at α=30° into a turbulent boundary layer profile at a freestream Reynolds number of Re∞=400000.This jet-in-a-crossflow(JICF) problem is investigated using large-eddy simulations(LES).The governing equations comprise the Navier-Stokes equations plus additional transport equations for several species to simulate a non-reacting gas mixture.A variation of the density ratio is simulated by the heat-mass transfer analogy,i.e.,gases of different density are effused into an an air crossflow at a constant temperature.An efficient large-eddy simulation method for low subsonic flows based on an implicit dual time-stepping scheme combined with low Mach number preconditioning is applied.The numerical results and experimental velocity data measured using two-component particle-image velocimetry (PIV) are in excellent agreement.The results show the dynamics of the flow field in the vicinity of the jet hole,i.e.,the recirculation region and the inclination of the shear layers,to be mainly determined by the velocity ratio.However,evaluating the cooling efficiency downstream of the jet hole the mass flux ratio proves to be the dominant similarity parameter,i.e.,the density ratio between the fluids and the velocity ratio have to be considered.      

Direct Numerical Simulation of Secondary Breakup of Liquid Drops

This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional (2D) unsteady axisymmetric Navier-Stokes equations for both the drop and the ambient gas flow. The level set method is applied to capturing the liquid/gas interface. Through calculation are obtained four typical breakup modes-oscillation, bag breakup, sheet stripping breakup and shear breakup governed by four non-dimensional numbers which are gas Weber number (Weg), liquid Reynolds number (Rel), gas Reynolds number (Reg) and density ratio (γ) . Their effects upon each mode are analyzed. The results indicate that among the four numbers, Weg is of the highest importance with Rel, Reg and γfollowing up. By widening the range of the density ratio up to 1 000, the breakup mode is discovered to be so complicated that a new one called multimode breakup mode turns up. This mode contains the shearing breakup and piercing breakup, which successively happen. The calculation results agree well with what is observed from the experiments.     

自适应直角坐标网格方法在对流换热数值模拟中的应用

复杂几何边界是对流换热计算面临的主要难点之一。自适应直角坐标网格方法可以比较方便的处理复杂几何边界,但其应用目前主要集中于可压无粘流的Euler方程、可压粘性流的Navier-Stokes方程的外部绕流计算。本文应用自适应直角坐标网格方法计算了两个不可压内流对流换热算例,并把计算结果同商用CFD软件STAR-CD的计算结果进行了对比。自适应直角坐标网格方法在本文的两个算例中得到了同STAR-CD一致的结果。这说明采用该方法数值模拟对流换热问题是适用的。      

幂律流体液膜破裂的线性稳定性分析

为对凝胶推进剂撞击式喷嘴产生的液膜的破裂行为进行研究,使用线性稳定性分析方法并结合幂定律本构方程对幂律流体液膜的稳定性进行了研究,并对流体物性参数、流变参数等对液膜流动稳定性的影响规律进行了分析.结果发现:稠度系数、流动指数和表面张力越大,液膜表面扰动波的增长率越小,波长越大.而较高的液膜速度、气液密度比和液膜厚度都使扰动波增长率变大,但波长随液膜速度和气液密度比的增大而减小,而液膜厚度的情况却相反.理论计算与实验值的对比表明,线性稳定性分析能较准确地预测液膜的破碎长度,但预测的表面波长却大于实验值.      

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

基于N-S方程的地面效应模拟

模拟了低速情况下二维翼型的地面效应问题，其中翼面流动利用一种改进的、基于不可压N－S方程的离散涡方法来模拟，地面模拟则利用一系列连续分布的涡板块来实现。利用上述模型研究了地面对翼型表面流动和压力分布的影响，并对有地面存在条件下的尾涡发展规律进行了讨论。计算结果表明，该方法能够有效地模拟低速情况下翼型的近地面流动，具有一定的工程应用价值。