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

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

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

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

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

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

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

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

小推力液体火箭发动机燃烧室内流场数值模拟

为研究不同直径的推进剂液滴碰撞火箭发动机燃烧室壁面后对流场的影响,在FLUENT基础上通过UDF程序加入液滴撞壁模型,对一种双组元自燃推进剂发动机燃烧室内推进剂的蒸发、流动、燃烧过程进行分析计算,得出燃烧室的流场,液滴分布,液滴撞壁后轨迹的变化以及发动机的各项性能参数,符合规律的数值模拟结果,可为发动机的设计和试验提供依据。     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

液滴撞击超疏水表面的能量耗散机制

针对飞机表面易结冰部位设计超疏水表面,可以大幅度减轻对高能耗防/除冰技术的依赖程度,进而提高飞机的燃油经济性。主要通过实验研究与数值模拟的手段,分析讨论了液滴撞击分级粗糙结构超疏水表面过程中的能量耗散机制。以Ti6Al4V为基体经过喷砂处理形成微米级粗糙结构,然后在1mol/L的低浓度NaOH溶液中水热生长一层一维纳米线,构建出微/纳米复合粗糙结构并氟化修饰获得超疏水表面。通过场发射扫描电镜(FE-SEM)观察了微观形貌的变化规律,利用动态视频接触角测量仪表征试样表面液滴表观接触角与接触角滞后。基于气液两相流动界面追踪的复合Level set-VOF方法,实现了液滴撞击超疏水表面过程的数值模拟。采用高速摄像技术记录了撞击液滴在超疏水表面的运动过程,实验验证了模拟方法与铺展计算模型的正确性,并详细讨论了液滴运动过程中的能量耗散问题,分析表明液滴撞击过程中的能量耗散主要取决于超疏水表面的动态润湿特性和润湿界面模型。     

大密度比双气泡在孔板结构微通道内上升行为的格子Boltzmann方法模拟

基于格子Boltzmann两相流大密度比模型模拟了孔板结构微通道内双气泡在浮力作用下的上升过程,主要研究E?tv?s数（Eo）、气泡相对大小、气泡之间的距离以及气泡和孔板间的距离对气泡变形、合并的动力学行为以及气泡上升速度和气泡剩余质量的运动特性的影响。研究发现,随着Eo数的增大,气泡在通过孔板通道时形变越严重,表现为上部气泡和下部气泡在合并过程中所夹带的液泡数量和质量同时增加,且气泡在通过通道的过程中会发生多次接触、合并与破裂;数值结果还表明,随着Eo数的增大,气泡达到顶端的时间增加而气泡穿过孔板的质量减小。另一方面,当上方气泡的尺寸大于下方气泡的尺寸时,两气泡在合并的过程中夹带的液泡数量更少,气泡穿过孔板时更迟缓但能够穿过孔板的气泡质量增多。此外,对于不同的气泡间距离和不同的气泡与孔板之间的距离,发现上下气泡之间的距离过大或者过小时,在气泡的合并过程中都不容易夹带液泡,且气泡穿过孔板的质量随着两气泡之间距离以及上方气泡与孔板之间距离的减小而增加。     

障碍物导致火焰变形的数值模拟

工业部大量存在着障碍物。实践证明，障碍物对加速燃烧和诱导爆炸具有非常重要的影响。为了深入研究障碍物与火焰之间的关系本文基于湍流的k-ε模型和改进的EBU－Arrhenius燃烧模型，考虑了障碍物对流动的附加作用，通过修改方程的源项，建立了湍流加速火焰使火焰形状改变的理论模型。本文选用Simple格式，壁面边界层采用壁面函数法处理，模拟了火焰在含有障碍物的三维空间中的传播现象，其数值模拟得到的火焰形状与实验结果做了对比，结果表明本计算是成功的。最后，还分析了导致火焰变形的原因。     

复合材料整体化结构固化变形的数值模拟

介绍了复合材料整体化结构的概念,提出了复合材料结构件固化变形的具体分析方法,最后运用商业软件ABAQUS对树脂基碳纤维增强复合材料T700/QY8911在整体共固化过程中的变形进行了数值模拟.     

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.     

大叶尖间隙下叶栅性能数值模拟

利用NUMECA软件对某线性叶栅的三维流场进行了数值模拟,对比研究了小叶尖间隙和大叶尖间隙时叶尖泄漏涡流的形成和发展,探讨了叶尖间隙大小对叶栅流场和气动性能的影响。研究表明,随着叶尖间隙的增大,叶尖泄漏射流发展成为叶尖泄漏涡,涡流范围不断增大,涡流强度增大趋缓,涡流使得气流偏转减小。从小间隙逐渐增大到大间隙,总压损失与叶片载荷均增大,而在大叶尖间隙时,总压损失增加并不显著,叶片载荷增大趋缓。其结论为进一步揭示叶尖间隙涡流的流动机理以及工业燃气轮机的优化设计提供了参考。     

平面自由剪切层三维离散涡丝的数值模拟

本文用离散涡丝方法对平面自由剪切层进行数值模拟，给出了初始展向涡丝的随空间发展所形成的成对反转流向涡以及涡卷起和涡丝合并的图案，分析了流向涡产生的根源。得到了剪切层的速度剖面，给出了不同速度差与对流速度比之下涡丝的演变特点。     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

细长锥体有攻角绕流的非对称分岔及其涡结构的数值模拟

本文用ＮＮＤ格式计算了具有复杂涡结构的细长锥体有攻角绕流问题，着重研究了非对称分离的分岔现象以及分岔前后涡结构的演变，数值计算结果表明，大约在ａ／（２β）＝１．０时，Ｎ－Ｓ方程出现分岔解，同时还给出了分岔前后涡结构及物理最分布的特征。     

管道中超声速旋涡运动的数值模拟

本文利用ＥＮＮ格式通过求解ＮＳ方程，数值模拟了超声速旋涡在管道中沿其轴向的演变规律。结果表明，当流动处于顺压区，在涡轴附近，截面流线为不稳定的螺旋点形态；当流动处于逆压区，在涡轴附近，截面流线为稳定的螺旋点形态；当从顺压区过渡到逆压区，涡心附近出现不稳定极限环；反之，出现稳定极限环。数值模拟结果同张涵信的拓扑分析理论十分相符。     

涡流冷却推力室流场特征与性能仿真

针对2000N气氢/气氧涡流冷却推力室,采用三维全尺寸计算模型开展了仿真研究,得到了流场速度分布特点,验证了涡流冷却推力室内具有双向涡旋结构,内外涡流分界面约占涡流冷却推力室圆柱段半径的86%,燃烧区域约占涡流冷却推力室圆柱段半径的70%.分析表明:外层涡流主要受来流速度与涡流冷却推力室几何参数影响,内层涡流在黏性、燃烧等作用下室压、密度稳定.侧壁温度平均为388K,比冲效率达92%以上,仿真结果与试验对比一致.      

基于数值模拟的应力成形极限图的应用研究

根据极限应变到极限应力的转换关系,采用 Hill48塑性准则和 Hill79塑性准则建立了应力成形极限图的数学计算模型。并基于板料成形有限元软件 Dynaform,实现了用应力成形极限图作为判据对板料成形进行分析。     

跨,超声速三角翼背风区旋涡运动的数值模拟

采用杂交通量分裂的ＮＮＤ格式模拟了跨，超声速三角翼背风区的旋涡运动结构。结果表明，旋涡的截面结构依赖于涡轴上的速度与密度的乘积沿轴向的导数λ＝「（１／ρ）（δρｕ１／δｌ）」，当它大于零时，涡心附近的截面流线向里转，而当它小于零时，涡心附近的截面流线向外转，λ由大于零变到小于零，涡心附近将出现一个稳定的极限环，多一次变号，将多出现一个极限环。     

用二阶矩两相湍流模型模拟鼓泡床内气泡—液体湍流两相流动

建立了二阶矩气液两相湍流模型，模拟了不同工况下二维矩形断面气液鼓泡床中气泡－液体湍流两相流动，给出了气泡和液体的速度场、气泡体积分数和两相雷诺应力分布，基本工况的模拟结果与美国俄亥俄州立大学PIV测量结果符合较好。文中研究了不同气体表观速度对两相流场的影响。模拟结果显示了鼓泡床内液体的回流流动和气泡的上升运动、各向异性的上湍流，气泡湍流脉动比液体的强，以及气泡体积分数和两相湍流强度随着气体表观速度的增大而增大等规律。