The Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media

Recently developed transport equations for two-phase flow through porous media usually have a second term that has been included to account properly for interfacial coupling between the two flowing phases. The source and magnitude of such coupling is not well understood. In this study, a partition concept has been introduced into Kalaydjian's transport equations to construct modified transport equations that enable a better understanding of the role of interfacial coupling in two-phase flow through natural porous media. Using these equations, it is demonstrated that, in natural porous media, the physical origin of interfacial coupling is the capillarity of the porous medium, and not interfacial momentum transfer, as is usually assumed. The new equations are also used to show that, under conditions of steady-state flow, the magnitude of mobilities measured in a countercurrent flow experiment is the same as that measured in a cocurrent flow experiment, contrary to what has been reported previously. Moreover, the new equations are used to explicate the mechanism by which a saturation front steepens in an unstabilized displacement, and to show that the rate at which a wetting fluid is imbibed into a porous medium is controlled by the capillary coupling parameter, . Finally, it is argued that the capillary coupling parameter, , is dependent, at least in part, on porosity. Because a clear understanding of the role played by interfacial coupling is important to an improved understanding of two-phase flow through porous media, the new transport equations should prove to be effective tools for the study of such flow.     

Multiple time scales of fluvial processes — theory and applications

Fluvial processes comprise water flow, sediment transport and bed evolution, which normally feature distinct time scales. The time scales of sediment transport and bed deformation relative to the flow essentially measure how fast sediment transport adapts to capacity region in line with local flow scenario and the bed deforms in comparison with the flow, which literally dictates if a capacity based and/or decoupled model is justified. This paper synthesizes the recently developed multiscale theory for sediment-laden flows over erodible bed, with bed load and suspended load transport, respectively. It is unravelled that bed load transport can adapt to capacity sufficiently rapidly even under highly unsteady flows and thus a capacity model is mostly applicable, whereas a non-capacity model is critical for suspended sediment because of the lower rate of adaptation to capacity. Physically coupled modelling is critical for fluvial processes characterized by rapid bed variation. Applications are outlined on very active bed load sediment transported by flash floods and landslide dam break floods.     

Some issues related to the modeling of interfacial areas in gas–liquid flows I. The conceptual issues

Two-phase flow modeling has been under constant development for the past forty years. Actually there exists a hierarchy of models which extends from the homogeneous model valid for two-phase flows where the phases are strongly coupled to the two-fluid model valid for two-phase flows where the phases are a priori weakly coupled. However the latter model has been used extensively in computer codes because of its potential in handling many different physical situations.The two-fluid model is based on the balance equations for mass, momentum and energy, averaged in a certain sense and expressed for each phase and for the interface between the phases. The difficulty in using the two-fluid model stems from the closure relations needed to arrive at a complete set of partial differential equations describing the flow. These closure relations should supply the information lost during the averaging of the balance equations and should specify in particular the interactions of mass, momentum and energy between the phases. Another requirement for the interaction terms is that they should satisfy the interfacial balance equations. Some of these terms such as the added mass term or the lift force term do not depend on the interfacial area but some others do, such as the mass transfer term, the drag term or the heat flux term. It is then necessary to model the interfacial area in order to evaluate the corresponding fluxes. Another benefit resulting from the modeling of the interfacial area would be to replace the usual static flow pattern maps which specify the flow configuration by a dynamic follow-up of the flow pattern. All these reasons explain why so much effort has been put during the past twenty years on the modeling and measurement of the interfacial area in two-phase flows.This article contains two parts. The first one deals with the conceptual issues and has the following objectives:

1.
to give precise definitions of the interfacial area concentrations;
2.
to explain the origin of the interfacial area concentration transport equation suggested by M. Ishii in 1975;
3.
to explain some paradoxical behaviors encountered when calculating the interfacial area concentration transport velocity.

     

The prediction of stratified two-phase flow with a two-equation model of turbulence

A two-equation model is applied to a stratified two-phase flow system to predict turbulent transport mechanisms in both phases.In the present analysis, the effects of interfacial waves on the flow field are formulated in terms of boundary conditions for the gas-liquid interface. For the gas phase, the wavy interface has such flow separation effects as a rough surface in a single-phase flow. While for the liquid phase, the waves generate turbulant energy which is transported progressively toward a lower wall region. The analytical results are in good agreement with available data regarding pressure drop, holdup and velocity profiles.     

A first order relaxation model for the prediction of the local interfacial area density in two-phase flows

Many energy production and chemical processes involve vapor/liquid two-phase flows. Mass and energy are often exchanged between the vapor and the liquid phases, and the fluid mechanics of the two-phase system is strongly influenced by the exchange of momentum between each phase. Significantly, the transport of mass, energy and momentum between the phases takes place across interfaces. Therefore the interfacial area density (i.e. the interfacial area per unit volume) has to be accurately known in order to make reliable predictions of the interfacial transfers. Indeed, the interfacial area density must be known for both steady and transient two-phase flows. It is the purpose of this paper to present a first order relaxation model which is derived from the Boltzmann transport equation, and which accurately describes the evolution of interfacial area density for bubbly flows. In particular, the local, instantaneous interfacial area densities and volume fractions are predicted for vertical flow of a vapor/liquid bubbly flow involving both bubble clusters and individual bubbles.     

Diffusion effects on the interfacial tension of immiscible polymer blends

Most methods of measuring the interfacial tension between two immiscible polymers are based on the analysis of the shape that a drop of one polymer immersed in the other one exhibits under the action of flow or gravity. In such a situation, the small, yet nonzero mutual solubility between the two polymers acts toward mass transfer between the drop and the surrounding fluid. In this work, diffusion effects on the interfacial tension of the pair polyisobutylene/polydimethylsiloxane have been investigated by drop deformation under shear flow. When the drop was made of polyisobutylene, drop size decreased with time due to diffusion. Drop shrinkage was associated with a significant increase in interfacial tension, until an apparent plateau value was reached. The effect was attributed to a selective migration of molecular weights, which would act to enrich the drop with higher molar mass material. To support such an interpretation, drop viscosity was evaluated by drop shape analysis and it was actually found to increase with time. In some cases, the ratio between drop and continuous phase viscosity became higher than the critical value for drop breakup in shear flow. Upon inverting the phases (i.e., when the drop was made of polydimethylsiloxane), no significant transient effects were observed. In the light of these results, the problem of what are the correct values of interfacial tension and viscosity ratio for a polymer blend of a certain composition will also be discussed. Received: 25 January 1999 Accepted: 24 May 1999     

A Dynamic Pore Network Model for Oil Displacement by Wettability-Altering Surfactant Solution

A dynamic pore network model, capable of predicting the displacement of oil from a porous medium by a wettability-altering and interfacial tension reducing surfactant solution, is presented. The key ingredients of the model are (1) a dynamic network model for the displacement of oil by aqueous phase taking account of capillary and viscous effects, (2) a simulation of the transport of surfactant through the network by advection and diffusion taking account of adsorption on the solid surface, and (3) the coupling of these two by linking the contact angle and interfacial tension appearing in the dynamic network simulation to the local concentration of surfactant computed in the transport simulation. The coupling is two-way: The flow field used to advect the surfactant concentration is that associated with the displacement of oil by the injected aqueous phase, and the surfactant concentration influences the flow field through its effect on the capillarity parameters. We present results obtained using the model to validate that it reproduces the displacement patterns observed by other authors in two-dimensional networks as capillary number and mobility ratio are varied, and to illustrate the effects of surfactant on displacement patterns. A mechanism is demonstrated whereby in an initially mixed-wet medium, surfactant-induced wettability alteration can lead to stabilization of displacement fronts.     

Interfacial structure of air-water two-phase flow in a relatively large pipe

Miniaturized four-sensor conductivity probes are used to study flow structure development in air-water bubbly flow, cap-bubbly flow, and transition to slug flow. The measurements are performed at three different elevations in a vertical round pipe with an inner diameter of 101.6 mm. The time-averaged local void fraction, interfacial velocity, and bubble number frequency are measured by the conductivity probes. Also, the interfacial area concentration and averaged bubble Sauter mean diameter are obtained. A detailed representation of the flow structure is revealed by investigating the acquired data. Furthermore, comparisons of the data at three elevations demonstrate the development of the interfacial structure along the flow direction due to bubble interactions and hydrodynamic effects. This may provide the community with a better knowledge about two-phase flow in a relatively large pipe. In addition, these data can also serve as an experimental database for investigation of the interfacial area transport in large-pipe two-phase flow. Published online: 19 November 2002 This work was performed under the auspices of the U.S. Nuclear Regulatory Commission through the Institute of Thermal-hydraulics.     

A finite volume shock‐capturing solver of the fully coupled shallow water‐sediment equations

This paper describes a numerical solver of well‐balanced, 2D depth‐averaged shallow water‐sediment equations. The equations permit variable horizontal fluid density and are designed to model water‐sediment flow over a mobile bed. A Godunov‐type, Harten–Lax–van Leer contact (HLLC) finite volume scheme is used to solve the fully coupled system of hyperbolic conservation laws that describe flow hydrodynamics, suspended sediment transport, bedload transport and bed morphological change. Dependent variables are specially selected to handle the presence of the variable density property in the mathematical formulation. The model is verified against analytical and semi‐analytical solutions for bedload transport and suspended sediment transport, respectively. The well‐balanced property of the equations is verified for a variable‐density dam break flow over discontinuous bathymetry. Simulations of an idealised dam‐break flow over an erodible bed are in excellent agreement with previously published results, validating the ability of the model to capture the complex interaction between rapidly varying flow and an erodible bed and validating the eigenstructure of the system of variable‐density governing equations. Flow hydrodynamics and final bed topography of a laboratory‐based 2D partial dam breach over a mobile bed are satisfactorily reproduced by the numerical model. Comparison of the final bed topographies, computed for two distinct sediment transport methods, highlights the sensitivity of shallow water‐sediment models to the choice of closure relationships. Copyright © 2016 John Wiley & Sons, Ltd.     

Particle velocity field measurements in a near-wall flow using evanescent wave illumination

Evanescent waves from the total internal reflection of a 488 nm argon-ion laser beam at a glass-water interface were used to measure velocity fields in creeping rotating Couette flow within 380 nm of the stationary solid surface. Images of fluorescent 300 and 500 nm diameter polystyrene and silica particles suspended in water recorded at 30 Hz were processed using cross-correlation particle image velocimetry to determine the two in-plane velocity components with an in-plane spatial resolution of 40Ꮀ µm over a 200 µm (h)쏦 µm (v) field of view. The results are in reasonable agreement with the exact solution for the corresponding single-phase Stokesian flow. These data are, to our knowledge, the first velocity field measurements with this small out-of-plane spatial resolution (in all cases less than 380 nm), and the first such measurements in this interfacial or near-wall region. This paper describes the novel experimental diagnostic technique used to obtain these results.     

Rheology of polymer blends: linear model for viscoelastic emulsions

 Kerner's model for flow of composite elastic media is extended to an emulsion of viscoelastic phases with interfacial tension undergoing deformations of small amplitude. A privileged internal structure inside the suspended drops is discussed in terms of fluid circulation across the interface. It is shown that for usual drop radius and interfacial tension values of emulsions, the rheological behavior predicted by the model, with very simple expression for the complex shear modulus, is quantitatively similar to that predicted by Palierne's model. Predictions of the model are compared with experimental data obtained on a polystyrene/polyethylene blend sheared in a small-amplitude oscillatory mode. Received: 10 August 1998 Accepted: 18 December 1998     

Effect of gravity on axial development of bubbly flow at low liquid Reynolds number

In view of the great importance of two geometrical parameters such as void fraction and interfacial area concentration to the accurate two-phase flow analysis at microgravity conditions, axial developments of flow parameters such as void fraction, interfacial area concentration, bubble Sauter mean diameter, and bubble number density were measured in bubbly flow at microgravity and low liquid Reynolds number conditions where the gravity effect on the flow parameters were pronounced. A total of seven data sets were acquired in the flow range of the void fraction from 1.01% to 3.36% and the liquid Reynolds number from 1,400 to 4,750. The measurements were also performed in the similar flow range at normal gravity conditions. The transport mechanisms of the flow parameters are discussed in detail based on the data measured at normal and microgravity conditions, and the drift-flux model developed at microgravity conditions are compared with the measured data.An erratum to this article can be found at     

One-group interfacial area transport equation and its sink and source terms in narrow rectangular channel

The characteristics of two-phase flow in a narrow rectangular channel are expected to be different from those in other channel geometries, because of the significant restriction of the bubble shape which, consequently, may affect the heat removal by boiling under various operating conditions. The objective of this study is to develop an interfacial area transport equation with the sink and source terms being properly modeled for the gas–liquid two-phase flow in a narrow rectangular channel. By taking into account the crushed characteristics of the bubbles a new one-group interfacial area transport equation was derived for the two-phase flow in a narrow rectangular channel. The random collisions between bubbles and the impacts of turbulent eddies with bubbles were modeled for the bubble coalescence and breakup respectively in the two-phase flow in a narrow rectangular channel. The newly-developed one-group interfacial area transport equation with the derived sink and source terms was evaluated by using the area-averaged flow parameters of vertical upwardly-moving adiabatic air–water two-phase flows measured in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm. The flow conditions of the data set covered spherical bubbly, crushed pancake bubbly, crushed cap-bubbly and crushed slug flow regimes and their superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. Good agreement with the average relative deviation of 9.98% was obtained between the predicted and measured interfacial area concentrations in this study.     

Numerical Simulation of Three-Dimensional Bubble Oscillations by a Generalized Vortex Method

A numerical method is implemented for simulating the simultaneous three-dimensional volume and shape oscillations of a compressible vapor or gas bubble suspended in an inviscid ambient fluid in the presence of interfacial tension. The flow generated by the bubble expansion, contraction, and deformation is represented by an interfacial distribution of potential dipoles supplemented by a point source situated inside the bubble, accounting for changes in the bubble volume. The mathematical formulation is completed by setting the strength of the point source proportional to the integral of the density of the double-layer potential over the interface. The motion of marker points distributed over the interface is computed using a boundary-element implementation of Baker's generalized vortex method in which the normal component of the interfacial velocity is computed in terms of tangential derivatives of the vector potential associated with the dipoles, whereas the tangential component of the interfacial velocity is computed in terms of the surface gradient of the scalar harmonic potential. The density of the double-layer distribution is computed by solving an integral equation of the second kind using an iterative method, while the evolution of the interfacial distribution of the harmonic potential is computed using Bernoulli's equation for irrotational flow. The onset of interfacial irregularities due to numerical instabilities is prevented by truncating the Fourier–Legendre spectrum of the interfacial distribution of the harmonic potential. With smoothing implemented, the numerical method is capable of describing simultaneous volume and shape oscillations for an indefinite period of time. Received 7 September 2001 and accepted 30 April 2002 Published online 30 October 2002 RID="*" ID="*" This research was supported by a grant provided by NASA. Communicated by J.R. Blake     

Suspended load and bed-load transport of particle-laden gravity currents: the role of particle–bed interaction

The development of particle-enriched regions (bed-load) at the base of particle-laden gravity currents has been widely observed, yet the controls and relative partitioning of material into the bed-load is poorly understood. We examine particle-laden gravity currents whose initial mixture (particle and fluid) density is greater than the ambient fluid, but whose interstitial fluid density is less than the ambient fluid (such as occurs in pyroclastic flows produced during volcanic eruptions or when sediment-enriched river discharge enters the ocean, generating hyperpycnal turbidity currents). A multifluid numerical approach is employed to assess suspended load and bed-load transport in particle-laden gravity currents under varying boundary conditions. Particle-laden flows that traverse denser fluid (such as pyroclastic flows crossing water) have leaky boundaries that provide the conceptual framework to study suspended load in isolation from bed-load transport. We develop leaky and saltation boundary conditions to study the influence of flow substrate on the development of bed-load. Flows with saltating boundaries develop particle–enriched basal layers (bed-load) where momentum transfer is primarily a result of particle–particle collisions. The grain size distribution is more homogeneous in the bed-load and the saltation boundaries increase the run-out distance and residence time of particles in the flow by as much as 25% over leaky boundary conditions. Transport over a leaky substrate removes particles that reach the bottom boundary and only the suspended load remains. Particle transport to the boundary is proportional to the settling velocity of particles, and flow dilution results in shear and buoyancy instabilities at the upper interface of these flows. These instabilities entrain ambient fluid, and the continued dilution ultimately results in these currents becoming less dense than the ambient fluid. A unifying concept is energy dissipation due to particle–boundary interaction: leaky boundaries dissipate energy more efficiently at the boundary than their saltating counterparts and have smaller run-out distance.

---

     

细颗粒泥沙净冲刷和输移的大涡模拟研究

在传统水沙输移数值模拟研究中一般采用雷诺时均模拟技术(Reynolds-averaged simulation,RANS).与RANS相比,大涡模拟技术(large eddy simulation,LES)能够更加精确反映细部流动结构,计算机的发展使得采用LES探讨水流和泥沙运动规律成为可能.本文尝试给出净冲刷条件下悬沙计算的边界条件,采用动态亚格子模式对循环槽道和长槽道中的水流运动和泥沙输移进行了三维大涡模拟研究.利用直接数值模拟(direct numerical simulation,DNS)结果对LES模型进行了率定,计算结果符合良好,在此基础上初步探讨了泥沙浓度、湍动强度和湍动通量等的分布特征.结果表明,净冲刷条件下输沙平衡时泥沙浓度符合Rouse公式分布,单向流动中泥沙浓度沿着流向逐渐增大.泥沙浓度湍动强度和湍动通量都在近底部达到最大值,沿着垂向迅速减小.湍动黏性系数和扩散系数基本上在水深中间处达到最大.湍动Schmidt数沿着水深方向不是常数,在近底部和自由水面附近较大,在水深中间处较小.     

Some issues related to the modeling of interfacial areas in gas–liquid flows,II. Modeling the source terms for dispersed flows

The first part of this paper dealt with the conceptual issues encountered in the definition of the interfacial area and in the derivation of a transport equation. The second part has the following objectives: 1. to address the closure issues for the source terms appearing in the transport equation when dealing with a bubbly flow in a vertical pipe; 2. to provide a list of open questions to be answered before introducing a transport equation for the interfacial area concentration in a thermal–hydraulic computer code.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

Three-dimensional mathematical model and its application in the neighborhood of the Three Gorges Reservoir dam in the Yangtze River

The calculation of flow and sediment transport is one of the most important tasks in river engineering. The task is particularly difficult because a number of complex physical phenomena should be accounted for more realistically in a model with a predictive power. Three-dimensional calculations of river flow and suspended sediment transport are performed in this paper with application in the Three Gorges Reservoir in the Yangtze River. A period of 76 years after the dam is built is simulated and the results are compared with laboratory measurements obtained by Tsinghua University whereby the model is verified and calibrated. Generally speaking, the calculated results agree well with the experiments, demonstrating that the present model can be used for flow and sediment transport prediction in rivers. The project supported by the National Natural Science Foundation of China (50009004)