An Analysis of the Movement of Wetting and Nonwetting Fluids in Homogeneous Porous Media

The movement of wetting and nonwetting fluid flow in columns packed with glass beads is used to understand the more complicated flows in homogeneous porous media. The motion of two immiscible liquids (oil and water) is observed with different surfactants. Through dimensional analyses, fluid velocity is well correlated with interfacial tension and less dependent on particle size. In water–oil (W/O) experiments, finger pattern flows are observed if water is the displacing fluid that flows in an oil-filled porous media, whereas oil ganglia tend to form if oil is the displacing fluid in the water-wetted porous media. The results are well described by a simple model based on an earlier theory of flow in a tube.     

Micromechanics of residual saturation formation during immiscible displacement in a porous medium

The process of displacement of a nonwetting fluid has been studied experimentally on a transparent model of a porous medium for various percolation velocities in the stable front regime, when the viscosity of the displacing fluid is greater than that of the fluid displaced. The flow structures in the final displacement regime, when the nonwetting phase is distributed in the pore space in the form of individual drops or ganglia, have been visually investigated. Imbibition is numerically modeled on a two-dimensional network model with allowance for various microdisplacement mechanisms. The effect of the initial displacing phase saturation on the magnitude and structure of the residual displaced fluid saturation is demonstrated. The fractal dimensionality of the displacement boundary is measured.Novosibirsk. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.3, pp. 116–121, May–June, 1994.     

Finite-Difference Approximation for Fluid-Flow Simulation and Calculation of Permeability in Porous Media

We introduce a finite-difference method to simulate pore scale steady-state creeping fluid flow in porous media. First, a geometrical approximation is invoked to describe the interstitial space of grid-based images of porous media. Subsequently, a generalized Laplace equation is derived and solved to calculate fluid pressure and velocity distributions in the interstitial space domain. We use a previously validated lattice-Boltzmann method (LBM) as ground truth for modeling comparison purposes. Our method requires on average 17 % of the CPU time used by LBM to calculate permeability in the same pore-scale distributions. After grid refinement, calculations of permeability performed from velocity distributions converge with both methods, and our modeling results differ within 6 % from those yielded by LBM. However, without grid refinement, permeability calculations differ within 20 % from those yielded by LBM for the case of high-porosity rocks and by as much as 100 % in low-porosity and highly tortuous porous media. We confirm that grid refinement is essential to secure reliable results when modeling fluid flow in porous media. Without grid refinement, permeability results obtained with our modeling method are closer to converged results than those yielded by LBM in low-porosity and highly tortuous media. However, the accuracy of the presented model decreases in pores with elongated cross sections.     

Wetting and nonwetting fluid displacements in porous media

The understanding of simple laminar flow in tubes has often been used to interpret the more complicated flow in porous media. A study of the motion of two immiscible liquids in closed tubes with relatively large diameter (> 0.3 cm i.d), was conducted in order to examine the influence of wetting and nonwetting liquids on the flow behavior. The results indicate that the wetting properties of the fluids with regard to the tube wall have a major efffect on the formation and motion of long bubbles. A physically based model was used to predict the velocity and the conditions for no motion of bubbles and drops in tubes. These results were used to interpret the nature of oil and water flow in porous media. Experiments in which the wetting liquid was displaced by the nonwetting, or vice versa, were conducted by injecting the displacing liquid at a constant flux at the center of a two-dimensional chamber saturated with the displaced liquid. The influence of wetting-nonwetting characteristics on the quantity of liquid displaced, the shape of the interface between the two liquids, and the interpretation of the no motion radius in a closed tube to the case of a porous medium are discussed. It would appear that the no motion radius gives a good indication of the minimum width of a nonwetting penetrating finger and the maximum width of nonwetting ganglia left by drainage.     

Effects of an impermeable wall in dissipative dynamics of saturated porous media

A phase transition model for porous media in consolidation is studied. The model is able to describe the phenomenon of fluid-segregation during the consolidation process, i.e., the coexistence of two phases differing on fluid content inside the porous medium under static load. Considering pure Darcy dissipation, the dynamics is described by a Cahn–Hilliard-like system of partial differential equations (PDE). The goal is to study the dynamics of the formation of stationary fluid-rich bubbles. The evolution of the strain and fluid density profiles of the porous medium is analysed in two physical situations: fluid free to flow through the boundaries of the medium and fluid flow prevented at one of the two boundaries. Moreover, an analytic result on the position of the interface between the two phases is provided.     

Magnetohydrodynamics convective flow in a vertical slot through a porous medium

An analysis is presented with magnetohydrodynamics natural convective flow of a viscous Newtonian fluid saturated porous medium in a vertical slot. The flow in the porous media has been modeled using the Brinkman model. The fully-developed two-dimensional flow from capped to open ends is considered for which a continuum of solutions is obtained. The influence of pertinent parameters on the flow is delineated and appropriate conclusions are drawn. The asymptotic behaviour and the volume flux are analyzed and incorporated graphically for the three-parameter family of solution.     

Motion equations of the dynamic theory of consolidation from the point of view of the theory of transient pipe flows

An elastic fluid-saturated porous medium is modeled as a bundle of parallel cylindrical tubes aligned in a direction parallel to the fluid movement. The pore space is filled with viscous compressible liquid. A cell model and the theory of transient pipe flow are used to derive one-dimensional governing equations in such media. All macroscopic constants in these equations are defined by the individual material constants of the fluid and solid. The interaction force includes an additional term not found in Biot's theory.     

Experimental Conditions Favoring the Formation of Fluid Banks during Counter-Current Flow in Porous Media

The purpose of this study is to investigate factors that affect the formation of fluid banks during gravity-driven counter-current flow in porous media. To our knowledge, development of a fluid bank has been observed in only one previous counter-current flow experiment, although there are some hints of fluid banks in other experiments. We have undertaken experimental and simulation studies to confirm the presence of such banks and to delineate factors which enhance or inhibit their formation. Experiments were performed using glass bead packs and X-ray Computed Tomography to monitor saturation distribution as a function of time. The simulation approach considers saturation history at every point in the sample, defining conditions at each time point from hysteresis in capillary pressure and relative permeability. The model proved to reproduce experimental observations accurately. The experiments and associated model show that a minimal vertical sample height is needed for the development of a fluid bank. In addition, round sample boundaries and higher average nonwetting phase saturation tend to prevent the formation of a bank. The validated model can improve our ability to predict and optimize counter-current flow processes, both in the laboratory and in the field (e.g. exploration and hydrocarbon extraction).     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.     

Natural convection in a porous medium coupled across an impermeable vertical wall with film condensation

This paper presents a theoretical study of thermofluid interaction between natural convection in fluid-saturated porous medium and film condensation, coupled through an impermeable vertical wall. The two heat transfer modes are analyzed separately. The solutions are matched on the wall. The complexion of this two-fluid problem is governed by a dimensionless interaction parameter which relates the heat transfer effectiveness of the two heat transfer mechanisms. The effect of this parameter on the flow and heat transfer is documented. Results regarding the overall heat transfer coefficient are obtained for a wide range of the independent parameters. Received on 19 January 1998     

Axi-symmetric simulation of a two phase vertical thermosyphon using Eulerian two-fluid methodology

Numerical simulation of steady state operation of a vertical two phase closed thermosyphon is performed using the two-fluid methodology within Eulerian multiphase domain. A full scale axi-symmetric model is developed for computational fluid dynamics simulation of thermosyphon using ANSYS/FLUENT 13.0. The effects of evaporation, condensation and interfacial heat and mass transfer are taken into account within the whole domain. Cooling water jacket is also modelled along with the wall of thermosyphon to simulate the effect of conjugate heat transfer between the wall and fluid phase. The results obtained are presented and compared with available experimental investigations for a similar thermosyphon. It is established that two-fluid methodology can be used effectively for the purpose of simulation of two phase system like a typical thermosyphon.     

Mechanistic Model of Steady-State Two-Phase Flow in Porous Media Based on Ganglion Dynamics

Recent experimental work has shown that the pore-scale flow mechanism during steady-state two-phase flow in porous media is ganglion dynamics (GD) over a broad and practically significant range of the system parameters. This observation suggests that our conception and theoretical treatment of fractional flow in porous media need careful reconsideration. Here is proposed a mechanistic model of steady-state two-phase flow in those cases where the dominant flow regime is ganglion dynamics. The approach is based on the ganglion population balance equations in combination with a microflow network simulator. The fundamental information on the cooperative flow behavior of the two fluids at the scale of a few hundred pores is expressed through the system factors, which are functions of the system parameters and are calculated using the simulator. These system factors are utilized by the population balance equations to predict the macroscopic behavior of the process. The dependence of the conventional relative permeability coefficients not only on the wetting fluid saturation S_wbut also on the capillary number, Ca, the viscosity ratio the wettability (⁰ a, ⁰ r), the coalescence factor, Co, as well as the porous medium geometry and topology is explained and predicted on a mechanistic basis. Sample calculations have been performed for steady-state fully developed (SSFD) and steady-state nonfully developed (SSnonFD) flow conditions. The number distributions of the moving and the stranded ganglia, the mean ganglion size, the fraction of the nonwetting fluid in the form of mobile ganglia, the ratio of the conventional relative permeability coefficients and the fractional flows are studied as functions of the system parameters and are correlated with the flow phenomena at pore level and the system factors.     

A linear dynamic model for a saturated porous medium

A linear isothermal dynamic model for a porous medium saturated by a Newtonian fluid is developed in the paper. In contrast to the mixture theory, the assumption of phase separation is avoided by introducing a single constitutive energy function for the porous medium. An important advantage of the proposed model is it can account for the couplings between the solid skeleton and the pore fluid. The mass and momentum balance equations are obtained according to the generalized mixture theory. Constitutive relations for the stress, the pore pressure are derived from the total free energy accounting for inter-phase interaction. In order to describe the momentum interaction between the fluid and the solid, a frequency independent Biot-type drag force model is introduced. A temporal variable porosity model with relaxation accounting for additional attenuation is introduced for the first time. The details of parameter estimation are discussed in the paper. It is demonstrated that all the material parameters in our model can be estimated from directly measurable phenomenological parameters. In terms of the equations of motion in the frequency domain, the wave velocities and the attenuations for the two P waves and one S wave are calculated. The influences of the porosity relaxation coefficient on the velocities and attenuation coefficients of the three waves of the porous medium are discussed in a numerical example.     

On the Hydrodynamics and Heat Convection of an Impinging External Flow Upon a Cylinder with Transpiration and Embedded in a Porous Medium

This paper extends the existing studies of heat convection by an external flow impinging upon a flat porous insert to that on a circular cylinder inside a porous medium. The surface of the cylinder is subject to constant temperature and can include uniform or non-uniform transpiration. These cylindrical configurations are introduced in the analyses of stagnation-point flows in porous media for the first time. The equations governing steady transport of momentum and thermal energy in porous media are reduced to simpler nonlinear differential equations and subsequently solved numerically. This reveals the dimensionless velocity and temperature fields of the stagnation-point flow, as well as the Nusselt number and shear stress on the surface of the cylinder. The results show that transpiration on the surface of the cylinder and Reynolds number of the external flow dominate the fluid dynamics and heat transfer problems. In particular, non-uniform transpiration is shown to significantly affect the thermal and hydrodynamic responses of the system in the circumferential direction. However, the permeability and porosity of the porous medium are found to have relatively smaller influences.     

Numerical Simulation of Darcy and Forchheimer Force Distribution in a HBC Fuse

We present the breaking of a short-circuit current in a HBC fuse simulation based on an isentropic non-stationary model in a porous medium for a one dimensional geometry. The fluid flow is affected by the nature of the gas and by the morphology of the silica sand. To model the gas–silica sand interaction, we introduce two classical laws: the Darcy's law due to the viscous interaction and the Forchheimer's law due to the inertial force. Numerical simulations with realistic physical parameters have been performed using a finite volume scheme with a fractional step technique. We show the evolution of Darcy and Forchheimer forces during time and according to the position in the fuse. We place in prominent position the fact that either force is predominant in the fuse according to the time and the position which justifies a numerical treatment to cover all the situations.     

Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

LBM Simulation of Viscous Fingering Phenomenon in Immiscible Displacement of Two Fluids in Porous Media

This article presents the lattice Boltzmann simulation of viscous fingering phenomenon in immiscible displacement of two fluids in porous media. Such phenomenon generally takes place when a less viscous fluid is used to displace a more viscous fluid, and it can be found in many industrial fields. Dimensionless quantities, such as capillary number, Bond number and viscosity ratio between displaced fluid and displacing fluid are introduced to illustrate the effects of capillary force, viscous force, and gravity on the fluid behaviour. The surface wettability, which has an impact on the finger pattern, is also considered in the simulation. The numerical procedure is validated against the experiment about viscous fingering in a Hele-Shaw cell. The displacement efficiency is investigated using the parameter, areal sweep efficiency. The present simulation shows an additional evidence to demonstrate that the lattice Boltzmann method is a useful method for simulating some multiphase flow problems in porous media.     

非牛顿流体两相渗流问题的摄动解

本文用奇异摄动法结合正则摄动法求解了考虑毛管力因素时多孔介质中弱非牛顿流体的两相驱替问题,得到了分流函数和湿相饱和度的渐近解析解。所得结果同数值解和经典的牛顿流体两相渗流结果进行了比较,并着重讨论了非牛顿因素的影响。     

Coupled Processes in Charged Porous Media: From Theory to Applications

Charged porous media are pervasive, and modeling such systems is mathematically and computationally challenging due to the highly coupled hydrodynamic and electrochemical interactions caused by the presence of charged solid surfaces, ions in the fluid, and chemical reactions between the ions in the fluid and the solid surface. In addition to the microscopic physics, applied external potentials, such as hydrodynamic, electrical, and chemical potential gradients, control the macroscopic dynamics of the system. This paper aims to give fresh overview of modeling pore-scale and Darcy-scale coupled processes for different applications. At the microscale, fundamental microscopic concepts and corresponding mass and momentum balance equations for charged porous media are presented. Given the highly coupled nonlinear physiochemical processes in charged porous media as well as the huge discrepancy in length scales of these physiochemical phenomena versus the application, numerical simulation of these processes at the Darcy scale is even more challenging than the direct pore-scale simulation of multiphase flow in porous media. Thus, upscaling the microscopic processes up to the Darcy scale is essential and highly required for large-scale applications. Hence, we provide and discuss Darcy-scale porous medium theories obtained using the hybrid mixture theory and homogenization along with their corresponding assumptions. Then, application of these theoretical developments in clays, batteries, enhanced oil recovery, and biological systems is discussed.