蒸汽-水冷凝泡状流相界面浓度输运方程

通过对冷凝消失项的合理模化,并考虑连续相湍流导致的气泡聚合和破裂现象,建立适用于蒸汽水冷凝泡状流的相界面浓度输运方程。过冷条件下蒸汽气泡冷凝过程可分为:传热控制区和惯性控制区。在传热控制区气泡尺寸逐渐减小,在整个冷凝中所占时间较长,而在惯性控制区,气泡在很短的时间内因为周围液体的惯性迅速消失。利用所建立的相界面浓度输运方程对轴向相界面浓度进行了预测,预测结果与实验数据比较表明吻合很好。     

下降管中稠密两相湍流的数值模拟

本文将统一二阶矩两相湍流模型和颗粒动力学理论结合,推导并封闭了稠密两相流考虑颗粒间碰撞的统一二阶矩两相湍流模型。该模型用颗粒动力学理论模拟颗粒之间的碰撞,用各向异性的统一二阶矩模型考虑气相和颗粒相的湍流脉动,并用输运方程描述气固两相湍流之间的相互作用。最后用该模型对下降管中的气粒两相流动进行了模拟,模拟所得的沿径向颗粒浓度分布和速度分布与实验吻合较好,预报结果也反应出了颗粒雷诺应力的各向异性。     

方截面鼓泡床气液两相瞬态数值研究

在双流体模型中引入界面浓度输运方程,利用界面浓度和气泡平均Sauter直径模化各相间作用力。引入一附加 湍动能输运方程模化气泡诱导引起的液相湍流。利用该模型对方截面鼓泡床内气液两相流进行三维瞬态数值模拟。计算结 果表明该模型能较好得模拟方截面鼓泡床内气液两相流时均和瞬态流动特征。     

碰撞聚合引起液滴数目变化的生灭过程模型

液滴碰撞聚合是气体/液滴两相流动中的重要现象.为研究碰撞聚合引起的液滴数目变化,考虑液相仅包含大小两类尺寸液滴的情形,计算不同尺寸的两个液滴间的碰撞概率,对固定观察区域内液滴的碰撞聚合、迁入及迁出因素作随机分析,建立描述液滴数目变化的生灭过程模型.由无条件概率模型方程,得出模型方程的相应平稳分布.以雨滴下落过程的质量增长为例,说明所得分布在气体液滴两相流动中的应用.     

基于拉氏概率密度的两相湍流二阶矩模型

用颗粒运动的拉氏分析和PDF方法,改进了颗粒相的二阶矩模型。由拉氏两相运动的随机微分方程出发,采用随机过程分析和信号分析法得到湍流两相流动的PDF输运方程,双流体模型方程和两相脉动速度相关的基本模式的封闭式,和用其它方法导出的方程与封闭式的结果一致,对封闭式作了重要的改进,在分析颗粒轨道上的流体湍流作用时间时,全面地引入拉氏分析的轨道穿越效应、惯性效应、连续效应和湍流的各向异性。     

颗粒Reynolds正应力轨道模型及后台阶两相流动数值模拟

用随机过程理论建立气固两相耦合脉动量Lagrange方程,并建立了气固两相流随机颗粒轨道模型中颗粒Reynolds正应力的Lagrange方程.将新的模型运用于各向同性的湍流衰减的流场中,模拟颗粒的湍流扩散特性,与文[1]中的模型和实验结果作比较,并使用新模型模拟了后台阶两相流动.     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

基于小关联时间两相湍流色噪声的近似模型及数值模拟

从一般高斯型色噪声模型出发,通过泛函导数,应用小关联时间,近似计算多维色噪声,得到有效Fokker-Planck方程.将其应用到两相湍流中得到颗粒相的概率密度函数输运方程,从而得到颗粒相的二阶矩模型.将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,并将求解结果与实验结果进行对比分析.     

湍流边界层中重粒子弥散的随机模型

在重粒子轨道模型的基础上,引入了Saffman力,并考虑了粒子-固壁碰撞和粒子-粒子碰撞的影响,建立了重粒子运动方程,耦合湍流脉动的随机方程,发展了重粒子弥散的随机模型,并在湍流边界层中考察该模型.将数值计算结果与实验结果进行比较,同时考察了Saffman力和粒子碰撞对计算结果的影响.     

快速流化床内相分布现象的初步动力学分析

借助于双流体模型建立了快速流化床内充分发展状态颗粒相与气相平均湍流运动方程组.两相联合径向动量方程的形式揭示了床内径向空隙率分布对于当地湍流结构的依赖关系,并解释了两相流中相分布现象存在的一般性.两相联合动量方程的数值求解示范结果表明床内空隙率分布的变化趋势与实验基本吻合,为模型的合理性提供了一种佐证.     

垂直管内油水两相流局部相分布特性实验研究

应用双头电导探针测量系统,对垂直上升管内油水两相分散流局部相分布特性进行了系统测量。得到了油水两相分散流的局部含油率分布类型图。研究结果表明低折算水速和低折算油速条件下,局部含油率在实验段截面上呈抛物线型局部分布特征,局部最大值出现在实验段中心区域。随折算油速增大,油滴受到横向力如升力的作用,逐渐向实验段壁面区域迁移,形成局部含油率的壁面峰值分布特性。当折算水速大于0．8 m／s时,局部含油率在实验段截面上呈均匀分布。     

Droplet transport and coalescence kinetics in emulsions subjected to acoustic fields

A method to aid the separation of the oil phase from aqueous emulsions using a low-intensity, resonant ultrasonic field has recently been developed. The density and compressibility difference between the dispersed and continuous phases within the emulsion results in a net force on the oil drops that pushes them toward the pressure antinodes of the standing-wave field, where coalescence subsequently occurs. A trajectory model is developed to predict the relative motion of drops subjected to the acoustic field. Such trajectories are sensitive to the physical properties and relative size of interacting drops, the initial configuration of the drops, and acoustic field parameters. Model predictions are validated by comparing experimentally observed trajectories with those predicted by the model. The modeling approach is then extended to determine the temporal evolution of the size of the region surrounding a target drop cleared by coalescence as a function of physical and acoustic field parameters. These results form the basis of a population balance model that attempts to track the size-evolution of a drop population coalescing under the influence of an acoustic field.     

差压法测量两相流相含率的研究

本文通过差压式测量方法,对气液两相流中含气率及油水两相流中含水率的测量进行了深入的试验研究。通过水平及垂直上升段的结构组合,测量气液两相流中的含气率值,在对流型测量与分析的基础上,建立了不同流型下垂直上升段与水平段摩擦压降间的修正关系式,减小了测量误差;通过采用一缓流的垂直下降段,测量油水两相流中的含水率值,获得了好的结果。     

倾斜均质表面上等径水滴聚合过程的可视化实验

对倾斜均匀表面上等直径水滴的聚合过程及特性进行了可视化实验研究,获得了水滴直径和表面倾角等参数对液滴聚合过程中液滴液桥半径、接触角和接触线变化特性的影响,分析了水滴聚合对其运动的影响.实验结果表明:表面倾角越大,下滑的临界半径越小;液滴的直径越大,液滴聚合后越容易下滑;液滴聚合可以加快液滴的运动,使下滑临界半径减小.     

Interfacial Area Transport Equation for Bubble Coalescence and Breakup: Developments and Comparisons

Bubble coalescence and breakup play important roles in physical-chemical processes and bubbles are treated in two groups in the interfacial area transport equation (IATE). This paper presents a review of IATE for bubble coalescence and breakup to model five bubble interaction mechanisms: bubble coalescence due to random collision, bubble coalescence due to wake entrainment, bubble breakup due to turbulent impact, bubble breakup due to shearing-off, and bubble breakup due to surface instability. In bubble coalescence, bubble size, velocity and collision frequency are dominant. In bubble breakup, the influence of viscous shear, shearing-off, and surface instability are neglected, and their corresponding theory and modelling are rare in the literature. Furthermore, combining turbulent kinetic energy and inertial force together is the best choice for the bubble breakup criterion. The reviewed one-group constitutive models include the one developed by Wu et al., Ishii and Kim, Hibiki and Ishii, Yao and Morel, and Nguyen et al. To extend the IATE prediction capability beyond bubbly flow, two-group IATE is needed and its performance is strongly dependent on the channel size and geometry. Therefore, constitutive models for two-group IATE in a three-type channel (i.e., narrow confined channel, round pipe and relatively larger pipe) are summarized. Although great progress in extending the IATE beyond churn-turbulent flow to churn-annual flow was made, there are still some issues in their modelling and experiments due to the highly distorted interface measurement. Regarded as the challenges to be addressed in the further study, some limitations of IATE general applicability and the directions for future development are highlighted.     

A level-set continuum method for two-phase flows with insoluble surfactant

A level-set continuum surface force method is presented to compute two-phase flows with insoluble surfactant. Our method recasts the Navier–Stokes equations for a two-phase flow with insoluble surfactant as “one-fluid” formulation. Interfacial transport and interfacial jump conditions are treated using the level-set method and the discrete Dirac function. Based on the density-weighted projection method, a stable semi-implicit scheme is used to decouple the velocity components in solving the regularized Navier–Stokes equations. It allows numerical simulations for a wide range of viscosity ratios and density ratios.Numerical simulations on single drop deformation in a 2D shear flow are presented. Simulations on two drop interaction shows that surfactants can play a critical role in preventing drop coalescence. A fully 3D simulation demonstrating the physical interactions of multiple surfactant-laden drops is presented.     

柔性管内油水环状流流动结构数值模拟

海洋中柔性管道受到外界猛烈撞击,撞击物的不同会导致管道内流道结构发生不同的改变.为探寻流道结构变化对环状流动结构的影响规律,采用N-S方程和湍流模型描述流体运动,用VoF模型与CSF模型追踪油水界面,建立数值方程与分析模型,分析流道不同形式的结构变形、不同程度变形和不同流速下油水环状流在管内的湍动能和漩涡结构演变,及环状流流型的演变和压力场的演变.分析结果表明:凹陷变形程度越大,漩涡的强度及面积越大,V形变形影响远大于○形变形影响,环状流的流型发生改变时仍然保持着水相包裹油相的稳定状态,最大压缩区域的油水界面压强也随之增大,且增速变大.所得结论对研究油水环状流在柔性管道运输中的稳定性具有一定指导意义.      

Modeling Three-Dimensional Multiphase Flow Using a Level Contour Reconstruction Method for Front Tracking without Connectivity

Three-dimensional multiphase flow and flow with phase change are simulated using a simplified method of tracking and reconstructing the phase interface. The new level contour reconstruction technique presented here enables front tracking methods to naturally, automatically, and robustly model the merging and breakup of interfaces in three-dimensional flows. The method is designed so that the phase surface is treated as a collection of physically linked but not logically connected surface elements. Eliminating the need to bookkeep logical connections between neighboring surface elements greatly simplifies the Lagrangian tracking of interfaces, particularly for 3D flows exhibiting topology change. The motivation for this new method is the modeling of complex three-dimensional boiling flows where repeated merging and breakup are inherent features of the interface dynamics. Results of 3D film boiling simulations with multiple interacting bubbles are presented. The capabilities of the new interface reconstruction method are also tested in a variety of two-phase flows without phase change. Three-dimensional simulations of bubble merging and droplet collision, coalescence, and breakup demonstrate the new method's ability to easily handle topology change by film rupture or filamentary breakup. Validation tests are conducted for drop oscillation and bubble rise. The susceptibility of the numerical method to parasitic currents is also thoroughly assessed.     

A speculative study of 23-order fractional Laplacian modeling of turbulence: some thoughts and conjectures

This study makes the first attempt to use the 23-order fractional Laplacian modeling of Kolmogorov -53 scaling of fully developed turbulence and enhanced diffusing movements of random turbulent particles. Nonlinear inertial interactions and molecular Brownian diffusivity are considered to be the bifractal mechanism behind multifractal scaling of moderate Reynolds number turbulence. Accordingly, a stochastic equation is proposed to describe turbulence intermittency. The 23-order fractional Laplacian representation is also used to model nonlinear interactions of fluctuating velocity components, and then we conjecture a fractional Reynolds equation, underlying fractal spacetime structures of Levy 23 stable distribution and the Kolmogorov scaling at inertial scales. The new perspective of this study is that the fractional calculus is an effective approach to modeling the chaotic fractal phenomena induced by nonlinear interactions.