自由射流不同速度剖面下界面稳定性的伪谱分析

建立了外围为气体,内部为液体射流的双流体模型,并用Chebyshev配点法研究了不同速度剖面自由射流的界面稳定性问题．利用非线性的坐标变换,将物理空间变换到计算空间,得到关于时间频率的广义特征值问题．通过对计算结果的分析与讨论,得出射流稳定性的基本特征及其随各种参数的变化趋势．     

热传导对气-液射流界面不稳定性的作用机理研究

研究了平面分层气-液射流在非线性温度分布条件下的界面不稳定性性质.考虑了气体的可压缩性、液体的粘性、以及气体热导率和密度随温度变化等事实.并应用正则模态方法将问题转化为四阶变系数常微分方程,用数值积分和多重打靶法对模型的空间模式进行了计算,研究了不稳定模态随各物理参量的变化趋势.计算表明模型所体现的不稳定性特征与其它模型的计算结果是一致的.同时计算还得出气体和液体的温差越小、雷诺数越大、热导率变大均将有利于液体射流有效雾化的结果.该结论与HJE.Co.Inc(Glens Falls,NY,USA)的实验数据是定性吻合的.     

比热容非常值对高超声速平板边界层稳定性的影响

在高超声速条件下,边界层中气体的温度可能很高,以致气体的比热容不再是常数而与温度有关．这时边界层中的流动稳定性如何是值得研究的问题．采用线性稳定性理论,考虑比热容与温度有关时高超声速可压缩平板边界层的稳定性,并与假定比热容为常值的情况作比较,发现对第一模态和第二模态波的中性曲线、最大增长率都有影响．因此,在高超声速情况下,比热容随温度变化是研究边界层稳定性时必须考虑的一个因素．     

对冲圆射流生成径向扩张液膜研究

采用两股互相冲击的圆射流可以形成环形的液体薄膜,液膜在径向扩展到一定的临界半径距离会破碎．数值模拟了液膜在周围气体中形成和破碎的非定常过程．考虑了液体和气体都是不可压缩Newton流体的轴对称问题．液体和气体的界面采用Level set函数来跟踪,Navier-Stokes 控制方程和物理边界条件采用有限差分格式离散求解．计算结果给出了环形液体薄膜形成并在其环形边缘处破碎,并缓慢运动的过程．液膜的厚度随着液膜在轴向的扩展会逐渐变薄,因此定义的局部Weber数会在径向逐渐减小,这里的局部Weber数定义为ρu2h/σ,其中ρ和σ分别为液体的密度和界面的张力,u和h分别为在径向某个位置的液膜的平均径向速度和半液膜厚度．数值结果表明就像实验中所观察到的那样,液膜径向扩展的过程的确会在局部Weber数趋向于1的时候终结而停止扩张．根据空间-时间线性稳定性理论,液膜的破碎最初是由正弦模式在临界局部Weber数Wec=1引起的,在临界局部Weber数小于1时会发生绝对不稳定性．在线性理论中另一个独立的模式,所谓的余弦模式,则增长比正弦模式要慢,从而会推测到正弦模式主导破碎的结论．然而,这里的数值结果却表明,余弦模式在界面波的非线性发展阶段实质的超越了正弦模式的增长,并对液膜的最终阶段的破碎起主导作用．这验证了线性理论只能够对触发时扰动波的性质进行预测,而对失稳后情况和结果的预测则不一定正确．     

高浓度悬浮固粒对射流界面的粘性稳定性影响

该文首次利用双流体模型和扰动速度势理论,推得含高浓度悬浮固粒的射流界面粘性稳定性方程和对应的固气扰动速度比值方程．通过数值计算,得到了不同雷诺数及固粒属性的射流界面粘性稳定性曲线和对应的固气扰动速度比值曲线．在分析和比较所得的粘性稳定性曲线的基础上,得到了流场雷诺数及固粒特性对射流界面粘性稳定性影响的结论．同时,通过分析所得的固气扰动速度比值曲线,得到了流场雷诺数及固粒等效斯托克斯数对固粒跟随气流的扰动 性能的影响的结论．这些结论是首次在计入气流的粘性的条件下得到的,不同于文献［８］和文献［１０］相关的囿于无粘情形的研究,对于两相射流发展的认识和工程实际中实施对两相射流场的人工控制有重要意义．     

含周期脉动的大密度比气液同轴射流的Floquet稳定性分析

基于周期脉动速度激励下气液同轴射流的数学模型,运用线性稳定性理论,采用Chebyshev配点法和Floquet理论,将含周期脉动分层流的Floquet稳定性分析扩展到大密度比的情况.研究了液铝-氮气射流的参数共振特性,分析了不同的物理参数对系统稳定性的影响,计算了实验工况并和实验结果进行了比较.     

合成射流激励器射流矢量控制的物理因素

对不同出口构型合成射流激励器进行射流矢量控制进行了数值研究,并对决定合成射流激励器射流矢量控制的物理因素进行了分析和归纳．低压区位置和面积及其压强梯度、合成射流动量分量、合成射流对主流的卷吸率是直接控制主射流矢量力和矢量角的物理因素．合成射流的3个特征参数直接影响和控制低压区的面积及其压强梯度,合成射流激励器出口台阶和出口斜喷角都对低压区位置、面积和合成射流对主流的卷吸率有影响和调节作用,合成射流激励器出口斜喷角还直接控制合成射流动量分量．基于对合成射流激励器射流矢量控制物理因素的分析,确定了控制物理因素的源变量,建立了由控制能力函数和调节功能函数组成的合成射流矢量控制初步模型,初步模型能够对源变量引起的合成射流激励器射流矢量控制效率不同作出解释,并进一步指出了进行射流矢量控制的最佳激励器是充分利用调节功能函数．     

液体射流泵理论的研究

本文运用流体力学及湍流射流理论,导出了液体射流泵基本方程、最优参数方程、汽蚀方程及装置性能方程.通过计算数学方法,利用电子计算机求出了上述方程的数值解,并经过国内外的试验资料验证.     

非线性Kelvin-Helmholtz不稳定性

本文用导数展开法对液体薄层与亚音速气流接壤时的界面稳定性作非线性分析.文中考虑了液体的表面张力与体积力,故非线性的Rayleigh-Taylor不稳定性可作为特例而导出;液体与气体均不计粘性.虽然Nayfeh[1]曾算过这一情况,但其三阶方程有遗漏(如213页的式(2.29)).同时解也不自洽(如其一阶解(2.31)并不满足他的初始条件(2.20)),此外,在截止波数附近,对行波他并未考虑.本文弥补了这些,并得出了新的结论.     

含悬浮固粒射流界面稳定性研究

利用气固两相耦合模型,理论推导出含悬浮固粒射流的稳定性方程,通过数值计算得到了两相射流稳定性特征曲线、固气扰动速度比值幅值曲线及固气相位差曲线,进而得到了关于固粒对流场中扰动增长和传播的影响及失稳过程中固粒扰动特性的结论。这些结论对于两相射流发展的认识和工程实际中实施对两相射流场的人工控制有重要意义。     

Stability of liquid films adjacent to compressible streams

An analysis is presented for the linear stability of a liquid film, adjacent to a compressible viscous gas stream. The analysis is valid for all wavelengths and liquid Reynolds numbers. The pressure and shear perturbations exerted by the gas on the liquid are calculated, using a gas model which takes into account the gas viscosity, velocity profile, and heat transfer. The results show that an inviscid, uniform stream model for the gas is inadequate unless the disturbed boundary layer is very thin. Although the present linear analysis is in fairly good agreement with the experimental observations for subsonic flow, it does not predict the observed wavelengths and wave speeds for supersonic flow.     

Multiphase transient flow model in wellbore annuli during gas kick in deepwater drilling based on oil-based mud

Transient-state gas and oil-based mud (OBM) two-phase flow in wellbore annuli will occur during gas kick. The phase behavior of influx gas and OBM will make the gas kick during OBM drilling more complicated. There are three possible cases in an annulus: only liquid flow in the entire annulus, gas and liquid two-phase flow in part of the annulus, and gas and liquid two-phase flow in the entire annulus. First, the phase behaviors of gas and OBM in wellbore annuli are studied based on the phase behavior of methane and diesel. A multiphase transient-flow model in annuli during gas kick based on OBM is then established based on gas–liquid two-phase flow theory and on flash theory in annuli. The influences of phase behavior in annuli and annular geometry are taken into account. The local flow parameters are predicted by the hydrodynamic models and the local thermodynamic parameters are predicted by the heat-transfer models in the corresponding flow pattern. The proposed model has a better performance, compared with two other models, against the published experimental data. Finally, the variation of pit gain, well-bottom hole pressure, and gas void fraction are obtained, leading to a better understanding of the occurrence and evolution mechanism of gas kick during deepwater drilling.     

Modeling and simulation of severe slugging in air-water systems including inertial effects

A mathematical model and numerical simulations corresponding to severe slugging in air-water pipeline-riser systems are presented. The mathematical model considers continuity equations for liquid and gas phases, with a simplified momentum equation for the mixture. A drift-flux model, evaluated for the local conditions in the riser, is used as a closure law. In many models appearing in the literature, propagation of pressure waves is neglected both in the pipeline and in the riser. Besides, variations of void fraction in the stratified flow in the pipeline are also neglected and the void fraction obtained from the stationary state is used in the simulations. This paper shows an improvement in a model previously published by the author, including inertial effects. In the riser, inertial terms are taken into account by using the rigid water-hammer approximation. In the pipeline, the local acceleration of the water and gas phases are included in the momentum equations for stratified flow, allowing to calculate the instantaneous values of pressure drop and void fraction. The developed model predicts the location of the liquid accumulation front in the pipeline and the liquid level in the riser, so it is possible to determine which type of severe slugging occurs in the system. A comparison is made with experimental results published in literature including a choke valve and gas injection at the bottom of the riser, showing very good results for slugging cycle and stability maps. Simulations were also made assessing the effect of different strategies to mitigate severe slugging, such as choking, gas injection and increase in separation pressure, showing correct trends.     

Well-Posedness of A Compressible Gas-Liquid Model for Deepwater Oil Well Operations

The main purpose of this paper is two-fold: (i) to generalize an existence result for a compressible gas-liquid model with a friction term recently published by Friis and Evje [SIAM J. Appl. Math., 71 (2011), pp. 2014–2047]; (ii) to derive a uniqueness result for the same model. A main ingredient in the existence part is the observation that we can consider weaker assumptions on the initial liquid and gas mass, and still obtain an existence result. Compared to the above mentioned work, we rely on a more refined application of the estimates provided by the basic energy estimate. Concerning the uniqueness result, we borrow ideas from Fang and Zhang [Nonlinear Anal. TMA, 58 (2004), pp. 719–731] and derive a stability result under appropriate constraints on parameters that determine rate of decay toward zero at the boundary for gas and liquid masses, and growth rate of masses associated with the friction term and viscous coefficient.     

On modeling evaporation

We develop a theoretical model for evaporation of a pure liquid drop on a thermally conductive solid substrate. We discuss a variety of effects regarding evaporation regime, the state of the liquid/gas interface and the content of gas phase. Then, we further consider two models: the one resulting from the one-sided non-equilibrium assumption and the other that assumes diffusion-limited regime and equilibrium at the liquid/gas interface. A single governing equation for the evolution of drop thickness is derived for both models. We show that although the model predicts qualitatively different temperature along liquid/gas and liquid/solid interface, the dynamics of the drops is almost the same.      

Analysis of a model of a cocurrent packed-bed column with periodic inlet conditions

We consider a model of a process in a cocurrent packed-bed column with the simplest kinetics of the rate of liquid holdup formation. Under study is the dependence of it on the time-periodic velocities of the liquid and gas flows in a neighborhood of the stationary values of these velocities. The principal terms of the asymptotics in time are found for the liquid holdup and the flow velocities. We detect the growth of the amplitudes of the sinusoidal oscillations of the liquid holdup and the liquid and gas velocities with respect to the column length. The constructed model is used for design of the cocurrent two-phase columns.     

CFD modeling of hydrodynamic characteristics of a gas–liquid two-phase stirred tank

A three-dimensional CFD model was developed in this work to simulate hydrodynamic characteristics of a gas–liquid two-phase stirred tank with two six-bladed turbines and four baffles, coupling of the Multiple Size Group model to determine bubble size distribution. Important hydrodynamic parameters of the multi-phase system such as volume-averaged overall and time-averaged local gas holdups and axial liquid velocities along time and transversal courses were simulated and analyzed in detail, under varied operating conditions (inlet air flow rate and impeller rotation speed). Model predictions of local transient gas holdup and liquid velocity distributions on vertical and horizontal sections of the tank were also carried out. The overall flow patterns were discussed in detail to assess the mixing. Bubble size distributions were further predicted to reveal the unique properties of gas phase. Experimental measurements of overall gas holdups and local axial liquid velocities were used to validate the developed model.     

Numerical simulation of internal flow fields of swirl coaxial injector in a hot environment

Based on the fractional volume of fluid (VOF), a pure Eulerian model for defining and capturing the gas/liquid interface is developed in this paper. This model can describe gas/liquid interface in high refinement, which is better than the original VOF methodology. To validate the proposed model and the algorithm, the computational code is employed to predict the flow performance in a cylindrical swirl injector under cold-flow condition, and the predicted results agree well with experimental measurements. Furthermore, the proposed model is used to simulate gas-liquid reacting flows inside a gas/liquid coaxial swirl injector operating in a hot environment. The turbulent combustion process is simulated with the k−ε−f−g model. The numerical simulation is carried out under actual operating condition of the coaxial injector. The injector performances, such as liquid film thickness, liquid film injection velocity, spray angle, pressure drop, are obtained based on the detailed information of the internal flow field. The predicted results also show that droplets are shed from the liquid film in the recess cup of the coaxial injector because of the large velocity gradient between the gas and liquid streams, and a burning area, which is characterized by high temperature, is present inside the injector.     

Temporal analysis of a power-law liquid sheet in the presence of acoustic oscillations

The instability of a non-Newtonian liquid sheet in the presence of acoustic oscillations is investigated theoretically. The power-law model is used to describe the viscosity of the non-Newtonian liquid. The corresponding dispersion relation is obtained by linear analysis. The effects of the mean velocity of the gas, the oscillation amplitude, the oscillation frequency, and the gas density on the instability of the power-law liquid sheet are studied. The results show that the shear-thickening liquid sheet is more unstable than Newtonian and shear-thinning liquid sheets when the effects of acoustic oscillations are considered. In particular, a second unstable region appears on the shear-thickening liquid sheet at a low oscillation frequency. Especially, for the shear-thinning liquid sheet, there is a second unstable region in the dispersion curve at a high mean gas velocity. A third unstable region appears on the shear-thinning liquid sheet at a high gas density in the presence of acoustic oscillations. The unstable range of the Newtonian liquid is always the widest among these liquids.