纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.      

纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.     

黏性不可压缩流体流动前沿的数值模拟

提出了模拟注射成型中黏性、不可压缩流体流动前沿的新方法. 将Hele-Shaw流动应用于非 等温条件下的黏性、不可压缩流体，建立了流动分析模型，用充填因子的输运方程描述流动 前沿. 应用高阶Taylor展开式计算每一时间步长的充填因子，用Galerkin方法导出了计算 充填因子各阶导数的递推公式. 给出了时间增量的选取方法，证明了它的稳定性. 针对Han 设计的试验模具，用相同的材料及工艺条件模拟充填过程，比较了传统方法和该方法的模 拟结果与实验结果的差异. 算例分析表明，该方法可以有效地提高注射成型中流动前沿的 模拟精度和计算效率.     

虚拟质量力对酸性气体-钻井液两相流波速的影响

基于两流体模型、酸性气体和钻井液状态方程,考虑酸性气体与钻井液相间虚拟质量力、粘性剪切力、相间动量交换及狭义相间阻力等条件,建立酸性气体与钻井液两相中压力波传播速度的数学模型,依据小扰动原理,对波速模型求解,得到关于波数K的波速方程。结果表明,在一定范围内,随空隙率、频率的增大,虚拟质量力对波速的影响显著增强;在高空隙率下,压强增大,虚拟质量力对波速的影响减弱;增大流体的密度或不可压缩性,均可使两相压力波速增大;延长气液交换时间或减小波动频率使相间有足够时间进行动量交换,两相压力波波速随之减小。     

Basset力研究进展与应用分析

Basset力为两相流中颗粒与流体存在相对加速度时所产生的一种非恒定气动力 ,以往对其进行了大量的实验研究、理论分析和数值模拟。本文通过对Basset力研究文献的综述 ,分析并归纳了各种不同两相流动问题中该力的影响情况 ,得到的结论是 :对于气泡在液体内的流动问题 ,当气泡运动的脉动频率很大或很小时可以忽略Basset力对其运动的影响 ;对于固体或液体颗粒在气体中的运动问题可以忽略Basset力的影响 ;而其它两相流动问题则需要根据具体问题的特点来决定是否考虑Basset力的影响 ,其中需要考察 :流体与颗粒密度差别 ,颗粒尺寸 ,流动特征时间和颗粒运动弛豫时间 ,相对加速度 (或减速度 )等因素。本文还探讨了Basset力研究的发展方向。     

基于水平集方法的毛管上升特征研究

毛管上升现象与许多行业密切相关，系统地对此现象进行研究具有重大意义。与传统理论研究方法不同，本文使用N-S方程耦合水平集方法模拟毛管气液上升行为。通过与简化条件的解析解对比，验证了模拟方法的可靠性。此外，详细地研究了毛管振荡现象，并分析了影响毛管振荡行为的主要因素。结果表明，水平集方法能够精确地表征毛管振荡现象，与数值解相比具有更高的精度。毛管长度的增加能够减弱液柱振荡，主要归结于非湿相气体的粘滞力作用；湿相密度和湿相粘度同样对毛管振荡现象影响显著。湿相密度越大，惯性力越大，促进了毛管振荡；而湿相粘度变大，会增大粘滞力作用，因此减弱了毛管振荡现象。毛管振荡是由多种影响因素共同控制的，流体的惯性力是造成毛管振荡的主要原因，而粘滞力是减弱毛管振荡行为的主要因素，使液柱振荡逐渐衰减，并稳定至平衡高度。     

两相正激波松弛流的数值研究

对气、固两相正激波松弛流动进行了数值研究,方程中考虑了压力梯度及非匀速运动引起的附加质量力。给出了完整的计算方法,对三种不同固相容积比、不同固相颗粒直径和不同马赫数分别进行了数值计算,所得计算结果给出了气、固两相的速度、温度、压力、密度、固相容积比及熵值沿流向分布情况。文章对气、固两相参数在松弛区中的变化规律进行了讨论。     

两相正激波松弛流的数值研究

对气、固两相正激波松弛流动进行了数值研究,方程中考虑了压力梯度及非匀速运动引起的附加质量力。给出了完整的计算方法,对三种不同固相容积比、不同固相颗粒直径和不同马赫数分别进行了数值计算,所得计算结果给出了气、固两相的速度、温度、压力、密度、固相容积比及熵值沿流向分布情况。文章对气、固两相参数在松弛区中的变化规律进行了讨论。     

湍流两相流的脉动速度联合PDF输运方程

概率密度函数（PDF）的方法是构造两相湍流模型的一种重要的方法.构建气体-颗粒速度联合PDF输运方程的关键是颗粒所见气体微团速度的Langevin方程.首先由N-S方程出发,精确推导出颗粒所见气体微团脉动速度的Langevin方程,进而通过理论分析表明,对比通常采用的颗粒所见气体微团瞬时速度的Langevin方程而言,采用前者能有效地减少关联量的统计偏差.最后,给出了颗粒-气体脉动速度的联合PDF输运方程.     

基于气液相界面捕捉的统一气体动理学格式

气液相界面运动的研究无论是在科学还是工程领域都是非常重要的. 其中, 非平衡流动的计算尤其受到关注. 基于此, 我们构造了捕捉气液相界面的统一气体动理学格式. 由于统一气体动理学格式将自由输运和粒子碰撞耦合起来更新宏观物理量和微观分布函数, 故而可以求解非平衡流动. 具体思路是, 通过将范德瓦尔斯状态方程所表达的非理想气体效应引入统一气体动理学格式之中来捕捉气液相界面, 两相的分离与共存通过范德瓦尔斯状态方程描述. 由于流体在椭圆区域是不稳定的, 因此气液相界面可以通过蒸发和凝结过程自动捕捉. 如此, 一个锋锐的相界面便可以通过数值耗散和相变而得到. 利用该方法得到麦克斯韦等面积律(Maxwell construction)对应的数值解, 并与其相应的理论解相比较, 二者符合良好. 而后, 通过对范德瓦尔斯状态方程所描述的液滴表面张力进行数值计算, 验证了Laplace定理. 此外, 通过模拟两个液滴的碰撞融合过程, 进一步证明了该格式的有效性. 但是, 由于范德瓦尔斯状态方程的特性, 其所构造的格式仅适用于液/气两相密度比小于5的情况.      

A Steady-State Upscaling Approach for Immiscible Two-Phase Flow

The paper presents a model for computing rate-dependent effective capillary pressure and relative permeabilities for two-phase flow, in 2 and 3 space-dimensions. The model is based on solving the equations for immiscible two-phase flow at steady-state, accounting for viscous and capillary forces, at a given external pressure drop. The computational performance of the steady-state model and its accuracy is evaluated through comparison with a commercial simulator ECLIPSE. The properties of the rate-dependent effective relative permeabilities are studied by way of computations using the developed steady-state model. Examples presented show the dependence of the effective relative permeabilities and capillary pressures, which incorporate the effects of fine scale wettability heterogeneity, on the external pressure drop, and thereby on the dimensionless macro-scale capillary number. The effective relative permeabilities converge towards the viscous limit functions as the capillary number tends to infinity. Special cases, when the effective relative permeabilities are rate-invariant, are also studied. The applicability of the steady-state upscaling algorithm in dynamic displacement situations is validated by comparing fine-gridded simulations in heterogeneous reservoirs against their homogenized counterparts. It is concluded that the steady-state upscaling method is able to accurately predict the dynamic behavior of a heterogeneous reservoir, including small scale heterogeneities in both the absolute permeability and the wettability.     

Relative Permeability Analysis of Tube Bundle Models,Including Capillary Pressure

The analytical equations for calculating two-phase flow, including local capillary pressures, are developed for the bundle of parallel capillary tubes model. The flow equations that are derived were used to calculate dynamic immiscible displacements of oil by water under the constraint of a constant overall pressure drop across the tube bundle. Expressions for averaged fluid pressure gradients and total flow rates are developed, and relative permeabilities are calculated directly from the two-phase form of Darcy's law. The effects of pressure drop and viscosity ratio on the relative permeabilities are discussed. Capillary pressure as a function of water saturation was delineated for several cases and compared to a steady-state mercury-injection drainage type of capillary pressure profile. The bundle of serial tubes model (a model containing tubes whose diameters change randomly at periodic intervals along the direction of flow), including local Young-Laplace capillary pressures, was analyzed with respect to obtaining relative permeabilities and macroscopic capillary pressures. Relative permeabilities for the bundle of parallel tubes model were seen to be significantly affected by altering the overall pressure drop and the viscosity ratio; relative permeabilities for the bundle of serial tubes were seen to be relatively insensitive to viscosity ratio and pressure, and were consistently X-like in profile. This work also considers the standard Leverett (1941) type of capillary pressure versus saturation profile, where drainage of a wetting phase is completed in a step-wise steady fashion; it was delineated for both tube bundle models. Although the expected increase in capillary pressure at low wetting-phase saturation was produced, comparison of the primary-drainage capillary pressure curves with the pseudo-capillary pressure profiles, that are computed directly using the averaged pressures during the displacements, revealed inconsistencies between the two definitions of capillary pressure.     

A mathematical model for pressure drop of two-phase dry-plug flow in circular mini/micro channels

Dry-plug flow is a variation of two-phase plug flow that occurs in small scale channels and refers to the dry wall conditions at gas portions of the flow regime. Previous experimental studies found a significant increase in pressure drop in this flow regime from the two-phase wet-plug flow regime. In this work an analytical model for the pressure drop of this flow regime is developed and phenomena that influence pressure drop are examined. Unlike previous models, the proposed model seems to be applicable to a wide range of capillary numbers and give good estimations for all static contact angles. Contact angle hysteresis turned out to play a major role in inducing pressure drop in this flow regime. Model’s predictions were in good agreement with previous experimental data. Finally the model is applied to very low capillary number region and pressure drop predictions for this region are presented.     

Discontinuous Galerkin finite element method applied to the coupled unsteady Stokes/Cahn-Hilliard equations

Two-phase flows driven by the interfacial dynamics are studied by tracking implicitly interfaces in the framework of the Cahn-Hilliard theory. The fluid dynamics is described by the Stokes equations with an additional source term in the momentum equation taking into account the capillary forces. A discontinuous Galerkin finite element method is used to solve the coupled Stokes/Cahn-Hilliard equations. The Cahn-Hilliard equation is treated as a system of two coupled equations corresponding to the advection-diffusion equation for the phase field and a nonlinear elliptic equation for the chemical potential. First, the variational formulation of the Cahn-Hilliard equation is presented. A numerical test is achieved showing the optimal order in error bounds. Second, the variational formulation in discontinuous Galerkin finite element approach of the Stokes equations is recalled, in which the same space of approximation is used for the velocity and the pressure with an adequate stabilization technique. The rates of convergence in space and time are evaluated leading to an optimal order in error bounds in space and a second order in time with a backward differentiation formula at the second order. Numerical tests devoted to two-phase flows are provided on ellipsoidal droplet retraction, on the capillary rising of a liquid in a tube, and on the wetting drop over a horizontal solid wall.     

Dynamic Network Model of Mobility Control in Emulsion Flow Through Porous Media

Modeling the flow of emulsion in porous media is extremely challenging due to the complex nature of the associated flows and multiscale phenomena. At the pore scale, the dispersed phase size can be of the same order of magnitude of the pore length scale and therefore effective viscosity models do not apply. A physically meaningful macroscopic flow model must incorporate the transport of the dispersed phase through the porous material and the changes on flow resistance due to drop deformation as it flows through pore throats. In this work, we present a dynamic capillary network model that uses experimentally determined pore-level constitutive relationships between flow rate and pressure drop in constricted capillaries to obtain representative transient macroscopic flow behavior emerging from microscopic emulsion flow at the pore level. A parametric analysis is conducted to study the effect of dispersed phase droplet size and capillary number on the flow response to both emulsion and alternating water/emulsion flooding in porous media. The results clearly show that emulsion flooding changes the continuous-phase mobility and consequently flow paths through the porous media, and how the intensity of mobility control can be tuned by the emulsion characteristics.     

Simulation and transport phenomena of a ternary two-phase flow

A chemical flood model for a three-component (petroleum, water, injected chemical) two-phase (aqueous, oleic) system is presented. It is ruled by a system of nonlinear partial differential equations: the continuity equation for the transport of each of its components and Darcy's equation for the two-phase flow. The transport mechanisms considered are ultralow interfacial tension, capillary pressure, dispersion, adsorption, and partition of the components between the fluid phases (including solubilization and swelling).The mathematical model is numerically solved in the one-dimensional case by finite differences using an explicit and direct iterative procedure for the discretization of the conservation equations. Numerical results are compared with Yortsos and Fokas' exact solution for the linear waterflood case including capillary pressure effects and with Larson's model for surfactant flooding. The effects of the above-mentioned transport mechanisms on concentration profiles and on oil recovery are also analyzed.     

Numerical prediction of droplet dynamics in turbulent flow,using the level set method

In the present article, the droplet dynamics in turbulent flow is numerically predicted. The modelling is based on an interfacial marker-level set (IMLS) method, coupled with the Reynolds-averaged Navier–Stokes (RANS) equations to predict the dynamics of turbulent two-phase flow. The governing equations for time-dependent, two-dimensional and incompressible two-phase flow are described in both phases and solved separately using a control volume approach on structured cell-centred collocated grids. The topological changes of the interface are predicted by applying the level set approach. The kinematic and dynamic conditions on the interface separating the two phases are satisfied. The numerical method proposed is validated against a well-known computational fluid dynamics problem. Further, the deformation and breakup of a single droplet either suddenly moved in air or exposed to turbulent stream are numerically investigated. In general, the developed numerical method demonstrates remarkable capability in predicting the characteristics of complex turbulent two-phase flows.     

流场中聚合物共混体系液滴形变的理论模型

讨论了两相聚合物共混体系中，悬浮于另一种牛顿（或粘弹）液体中的牛顿（或粘弹）液滴的形变理论模型．影响液滴形变的主要因素有两相的组成、粘度比和弹性比、动态界面张力、临界界面张力系数，外流场形式及其强度．对于两相均为牛顿流体的体系，理论预测能够与实验相符；对于两相（或其中一相）为粘弹流体的体系，由于弹性的影响而使液滴形变的研究变得复杂，理论模型尚需完善．建立完整的液滴形变理论模型还需深入研究界面层、微观分子形变、液滴之间及液滴和连续相介质之间的相互作用对液滴形变的影响     

PTV investigation of phase interaction in dispersed liquid–liquid two-phase turbulent swirling flow

An investigation of dispersed liquid–liquid two-phase turbulent swirling flow in a horizontal pipe is conducted using a particle tracking velocimetry (PTV) technique and a shadow image technique (SIT). Silicone oil with a low specific gravity is used as immiscible droplets. A swirling motion is given to the main flow by an impeller installed in the pipe. Fluorescent tracer particles are applied to flow visualization. Red/green/blue components extracted from color images taken with a digital color CCD camera are used to simultaneously estimate the liquid and droplet velocity vectors. Under a relatively low swirl motion, a large number of droplets with low specific gravity tend to accumulate in the central region of the pipe. With increasing droplet volume fraction, the liquid turbulence intensity in the axial direction increases while that in the wall-normal direction decreases in the central region of the pipe. In addition, the turbulence modification in the present flow is strongly dependent on the droplet Reynolds number; however, the interaction of droplet-induced turbulences is significant due to vortex shedding, particularly at high droplet Reynolds numbers and higher droplet volume fraction.     

A probabilistic method for field solutions in two-phase flows

Solid-particle motion and related transport phenomena in two-phase flow are fluctuating processes in space and time. A deterministic method can describe only partially the intrinsic physics of these processes. In this paper, the fluctuations of the flow parameters are modelled by considering the spatial correlations, and a probabilistic computational method for two-phase flow is presented. The probabilistic governing equations have been discretized in space using a finite volume method, and then solved by applying the Neumann expansion method. This last method is time efficient, and its convergence can be guaranteed even for large fluctuations. A liquid-solid particle mixture flow in a circular pipe is taken as an example. Computational results illustrate the merit of the probabilistic approach for the prediction of two-phase flow phenomena.