气液两相流混合模型代数重构方法

为提高国家数值风洞工程气液两相流软件混合模型求解器的自由界面模拟分辨率，增强多介质流体界面问题的模拟能力，基于代数重构方法，应用双曲正切函数界面捕捉格式对体积分数进行了插值重构，实现两相流混合模型与自由界面代数重构方法的耦合，抑制原混合模型下自由界面的非物理耗散现象，保持界面清晰尖锐；引入预处理技术，解决低速流动中由于特征值量级差距导致的计算刚性问题，提高收敛速度；应用双时间进行非定常时间推进。最后，采用方块对流、Zalesak盘刚性旋转和溃坝算例进行考核验证，数值结果显示，耦合了双曲正切函数界面捕捉格式的两相流混合模型在保持质量守恒的同时，获得了更高的自由界面分辨率，以及更丰富的多相流动细节。     

气液同轴离心式喷嘴喷雾流场数值模拟

采用计算流体动力学（ＣＦＤ）方法对气液同轴离心式喷嘴冷态液雾两相流场进行了数值模拟，并且与试验结果进行了比较。为了比较，计算中采用的入口边界条件和液滴质量平均直径（ＭＭＤ）及其分布规律都是由试验确定的。研究了气液喷注压降对喷雾流动过程和流强分布的影响规律。结果表明，数值模型和计算程序能较好地模拟气液同轴离心式喷嘴两相喷雾流动过程，喷雾的存在对流场流型有很大的影响。     

固体火箭发动机内气粒两相流动的SPH-FVM耦合方法数值模拟

为解决现有数值模型在求解固体火箭发动机内颗粒的燃烧流动等问题时所遇到的难题,基于拟流体模型,采用光滑粒子流体动力学方法 (SPH)求解离散颗粒相,追踪颗粒运动轨迹,采用有限体积方法 (FVM)对气体连续相进行描述,捕捉流场特性,两相间通过曳力、压力梯度、热传导、体积分数等参量进行耦合,建立了两种不同坐标系下方法间的耦合框架。耦合方法模拟了喷气推进实验室(JPL)喷管中气粒两相流动过程,得到了气相和颗粒相的参数分布规律,与相关实验结果及离散颗粒模型(DPM)数值模拟结果对比,吻合较好。针对Φ315mm实验发动机工作状态和结构特点,进行了发动机燃烧室内过载条件下的气粒两相流动数值模拟,分析了纵、横加速度载荷对发动机燃烧室内粒子场和聚集带的影响,与已有的数值模拟结果基本吻合。结果表明,新方法求解发动机气粒两相流问题准确可靠,适用于发动机内更为复杂的两相流问题数值模拟。     

基于气体射流的气液两相流动减阻特性

为了明确射流表面减阻特性,基于气液两相流动原理,建立具有射流功能的气液两相流动减阻计算模型,采用数值计算方法,通过RNG k-ε湍流模型对射流表面不同射流方向角、以及不同气体射流速度情况下的流场进行数值模拟及实验验证,分析气体射流对壁面黏性阻力、减阻率及流场的影响,研究气体射流的气液两相的减阻特性。结果表明:气体射流表面能够减小壁面黏性阻力,具有较好的减阻效果;四种射流方向角下,射流表面减阻率均在气体射流速度为5m/s时达到极大值,并且当射流方向角在xoz平面为-30°时,减阻率最大,为7.18%;数值模拟的减阻率效果略高于实验结果;气体射流的气液两相流动对壁面边界层进行了有效控制,减小了壁面的黏性剪应力和雷诺应力,导致壁面的黏性阻力减小,使得射流表面具有减阻性能。     

离心式喷嘴内流场特性的数值模拟

基于单流体模型及气液两相混合分数的概念，建立了离心式喷嘴内流场的数学模型，并进行了数值计算，计算结果给出了离心式喷嘴内部液膜与气涡共存的流场结构及气液交界面的几何形状，计算得到的壁面压力分布、出口处液膜内轴向速度分布以及雾角大小都与实验结果吻合的很好，与常用的跟踪气液界面的ALE方法相比，本文建立的数值模拟方法具有简单，计算量小和易于实现的特点。     

车载燃气轮机进气砂尘分离器设计与验证

对车载燃气轮机进气砂尘分离器的设计方法进行了研究,根据设计要求,提出了一套涡旋管粒子分离器的设计方案和设计流程,并据此设计了隔栅式涡旋管和空气滤清器总成.对涡旋管的三维气固两相流动进行数值模拟.以空气滤清器总成为试验对象,进行了气动试验和砂尘分离试验.数值模拟和试验结果表明:进气砂尘分离器的设计方法能够达到预期的设计效果,设计的分离器满足设计要求,所采用的数值模拟方法可以较准确地模拟涡旋管的三维气固两相流动,计算结果可信.      

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

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

超声速气流中横向煤油射流的数值模拟

为了对超声速气流中液体横向射流进行数值模拟,采用Eulerian Lagrangian方法描述气液两相流动现象,用随机轨道模型追踪液滴的运动历程,并考虑气流可压缩性、流场不均匀性及液滴变形对液滴运动的影响,开发了相应的计算程序。计算结果发现,液体射流与气流之间存在着强烈的相互作用,液滴在进入气流中不久就破碎成很小的子液滴,受煤油液雾的影响气流速度、温度急剧下降,同时在喷孔上游出现一道弓形激波。与液雾结构的纹影图像和液雾穿透的实验结果对比显示,数值计算结果基本合理。     

气液交叉射流动力学行为的基础研究

为了模拟燃油在高速喷嘴内的快速蒸发混合过程,基于Eulerian—Eulerian方法,用RNGk一8湍流模型建立了高速喷嘴内两相流流动的数学模型,数值模拟了高速喷嘴内气液两相交叉射流的流动规律。改变气液动量比和韦伯数等条件下的两相流场计算结果表明：动量比增大3倍时横向穿透深度略有增大;而较大的韦伯数能起到强化横向射流的破碎效果。计算结果与同条件下的试验结果相符合。结果为研究HiTAC快速速喷嘴内可燃混合气的快速形成奠定了基础。     

三维管道超声速气固两相流动的CFD/DSMC仿真

采用NND格式求解气相N-S方程，采用DSMC方法模拟固相颗粒运动，数值模拟了三维管道超声速气固两相流场，分析了固体颗粒对流场结构产生的影响。计算结果表明，将DSMC方法与CFD方法相结合是解决气固两相流动问题的一种有效途径。     

低速下水滴撞击固体表面运动模拟

以水滴撞击飞机部件造成结冰问题为背景,运用level set界面追踪方法,把不可压黏性方程组中的表面张力用level set函数表示,并利用2阶变密度投影算法求解二维不可压Navier-Stokes方程组,数值模拟了低速下的水气界面追踪问题.为了能够处理密度比为1000的水气界面问题,在界面处引入了一定宽度的光滑过渡带宽,进而消除了由于界面两边大密度比引起的数值震荡,由此较为成功地模拟了大水滴在不同表面张力时撞击平板后的形态变化及其运动轨迹.数值实验结果表明本算法可以有效地处理界面处密度比和黏性很大的不可压多相流问题.考虑界面处的表面张力,本算法可模拟不同表面张力下的界面运动变化.      

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

Solving a radiation problem with forward speed using a lifting surface method with a Green''s function

The progress of seakeeping computations requires development of computating codes for unsteady flows around a ship or its elements. In this paper, we present a method of calculation concerning waves radiated by an oscillating surface-piercing flat plate with forward speed, with a yaw angle. By use of Green's third identity, the problem is transformed into the resolution of a Fredholm integral equation of the first kind by a panel method using Green's function. The Green's diffraction-radiation function with forward speed is used. Its numerical values are calculated by an adaptative integration procedure to reduce the computation time. The present method permits determination of the pressure jump distribution across the plate, the total forces and moments. The results obtained are compared with other numerical methods in hydrodynamics and in aerodynamics, and with experimental data obtained in a water tank.     

Comparative analysis of mixing chambers for two multiphase flows of oppositely charged particles

In this paper, we formulated the criteria for evaluating the effectiveness of mixing chambers, needed in obtaining polymeric materials with reproducible properties. The results of comparative analysis of organizing the methods for mixing two multiphase flows of oppositely charged particles are presented. This analysis, carried out using CFD programs, shows that the mixing chamber construction in which the flows being mixed are directed at an angle to each other, and additional gas flow inlets are provided in the lateral wall, is the most efficient.     

多介质流动改进型Level-Set重新初始化方法

根据Level-Set函数在网格点处的值,预估出间断面的位置,并对间断面附近网格点上的LS值作相应调整,使其保持距离函数的性质,能够避免由于重新初始化造成间断面位置移动的问题。数值计算采用高精度WENO格式,对于多个可压缩流问题进行了数值模拟,表明了该方法的有效性。     

数值模拟中的边界处理方法研究

为适应于模拟超声速流场中的一些特殊物理现象,研究了不同物理现象对边界处理的要求。针对激波、旋涡和声波等流场中典型的流动现象,采用常用的外推、特征和缓冲三类边界处理方法分别进行数值模拟。计算结果表明超声速流场边界采用外推技术,只适于计算带有激波的边界。特征方法还能适于旋涡运动的计算。声波的计算只能采用缓冲层技术。     

一种适用于颗粒非规则分布的阻止SPH数值断裂的方法

 〖HT5”H〗摘〓要〖HTSS〗: 光滑粒子流体动力学(SPH)方法已广泛应用于航空航天领域中大变形的碰撞计算。数值断裂是SPH计算中一种常见问题,可严重影响计算结果。把适用于颗粒规则分布的阻止SPH数值断裂的方法推广到颗粒非规则分布的情形,得到了一种适用于颗粒非规则分布的阻止数值断裂的加颗粒方法。然后将该方法应用于存在数值断裂的悬臂梁弯曲仿真计算,验证了其有效性。最后把该方法应用于弹丸撞击飞机蒙皮的非规则颗粒模型计算,并和规则颗粒模型得到的结果及试验结果进行比较。对比结果表明：对于非规则分布的颗粒模型,该加颗粒方法可有效地阻止数值断裂,提高计算精度,并可较为普遍使用。     

水锤效应影响因素及防护结构的数值研究

燃油箱作为飞机上易损性最高的部件,当其遭受高速射弹袭击时会产生破坏严重的液压水锤效应,直接威胁飞机安全,对水锤效应问题的数值研究具有重要意义。采用CEL 方法对水锤效应问题进行数值模拟,首先同已有试验结果进行对比,验证数值计算模型的准确性;然后分析射弹入射速度、充液率对水锤效应的影响;最后对水锤效应的防护设计进行数值分析,优选防护结构设置形式。结果表明：水锤效应的破坏威力来自于冲击波,射弹动能越大、充液率越高,油箱破坏性越大;设置防护挡板能够在一定程度上降低冲击波对油箱结构的破坏程度,空气防护挡板结构的防护效果最佳。     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。