Dependence of Pore-to-Core Up-scaled Reaction Rate on Flow Rate in Porous Media

Due to inherent heterogeneities in structure, mineral placement and fluid velocity in rock, bulk reaction rates realized during reactive flow through porous media may differ significantly from that predicted by laboratory-measured rate laws. In particular, rate laws determined in batch reactor experiments do not capture any of the flow dependence that will be experienced in the porous medium. Based on network flow model simulations of anorthite and kaolinite reactions in two sandstone pore networks under acidic conditions commensurate with CO2 sequestration, we compute up-scaled reaction rates at the core scale and investigate the dependence of the observed reaction rates on flow rate. For the anorthite reaction which, under these acidic conditions is far from equilibrium and dominated by pH, we find a power law dependence of reaction rate on flow rate. For the kaolinite reaction, which is near equilibrium, a more complex dependence emerges, with the up-scaled rate tending to rapidly increasing net precipitation at low-flow rates, then reversing and tending toward net dissolution at high-flow rates.     

Drying Rate Measurements in Convection- and Diffusion-Driven Conditions on a Shaly Sandstone Using Nuclear Magnetic Resonance

Carbon Capture and Storage (CCS) is one of the solutions studied to reduce greenhouse gas accumulation in the atmosphere. Depleted oil and gas reservoirs have been studied for potential storage sites but also saline aquifers that have the advantages of much larger pore volume. In this latter case, injection of large volume of anhydrous carbon dioxide will lead to a strong water desaturation of the near wellbore region because of evaporation mechanisms. Even the capillary trapped water can be removed by thermodynamical transfer of water vapor in the CO₂ phase. The extension in time and space of the dry zone will be controlled by the drying rate induced by the gas flow. Consequences of drying may induce alteration of the injectivity by salt precipitation and/or alteration of the rock fabric itself, especially for shaly sandstones in the case of clay drying. The context of CCS has raised new interests in the understanding of drying kinetic where the water vapor is evacuated by gas convection. In this study, we investigated experimentally the drying rate evolution with time on a shaly sandstone sample in two conditions of drying: convective and diffusive. In convective conditions, air is injected at different flow rates through the porous media in conditions of drying representative of a CO₂ injection site at one million ton per year. In diffusive conditions, no flow is imposed and the water vapor escape by diffusion. Drying rates dynamics in both conditions were measured by Nuclear Magnetic Resonance (NMR) and compared. We varied the temperature and the salinity in diffusive-driven drying and the gas flow rate in convective-driven drying. The water distribution in the pore network and the water saturation profiles were monitored continuously using T₂ relaxation and 1D imaging NMR techniques. For the range of temperature and air flow rate used, we show that drying rates in the two drying conditions are similar but not identical. They both present different periods characteristic of the main mechanisms for water mass transfer. Drying rate has a power law dependence on the temperature, as predicted by thermodynamic, and drying rate was found proportional to the flow rate in convective drying. Presence of salt has a complex effect: an increase of the drying rate at early stage of drying followed by a strong decrease for the remaining time of drying.     

Heat dispersion in stationary mixed convection flow about horizontal surfaces in porous media

An analysis has been performed to study the influence of velocity dependent dispersion on transverse heat transfer in mixed convection flow above a horizontal wall of prescribed temperature in a saturated porous medium. The Boussinesq approximation and boundary layer analysis were used to numerically obtain gravity affected temperature and velocity distributions within the frames of Darcy's law and a total thermal diffusivity tensor comprising both of constant coefficient heat conduction and velocity proportional mechanical heat dispersion. Dependending on Pe_∞, the molecular Peclét number basing on the effective thermal diffusivity and the velocity of the oncoming flow, density coupling has distinct influences on heat transfer rates between the wall surface and the porous medium flow region. For small Peclét numbers, when heat conduction is the prevailing mechanism, wall heat fluxes are the higher the larger the density difference between the oncoming and the near wall fluid is. The opposite is true for larger Peclét numbers, when mechanical heat dispersion is the main cause of heat spreading. For Pe_∞ tending to infinity these wall heat fluxes approach finite maximum values in the total heat diffusivity model, they grow beyond any limit if only constant coefficient heat conduction is considered. Thus, the inclusion of mechanical heat dispersion effects yields physically more realistic predictions. Received on 18 September 1996     

Application of X-ray CT investigation of CO<Subscript>2</Subscript>–brine flow in porous media

A clear understanding of two-phase flows in porous media is important for investigating CO₂ geological storage. In this study, we conducted an experiment of CO₂/brine flow process in porous media under sequestration conditions using X-ray CT technique. The flow properties of relative permeability, porosity heterogeneity, and CO₂ saturation were observed in this experiment. The porous media was packed with glass beads having a diameter of 0.2 mm. The porosity distribution along the flow direction is heterogeneous owing to the diameter and shape of glass beads along the flow direction. There is a relationship between CO₂ saturation and porosity distribution, which changes with different flow rates and fractional flows. The heterogeneity of the porous media influences the distribution of CO₂; moreover, gravity, fractional flows, and flow rates influence CO₂ distribution and saturation. The relative permeability curve was constructed using the steady-state method. The results agreed well with the relative permeability curve simulated using pore-network model.     

Equilibrium of a crack in a porous medium with injection of the filtering liquid

In connection with the exploitation of petroleum deposits, the article discusses the equilibrium of a porous medium with a crack under conditions of plane deformation, with the steady-state filtration of a liquid injected into the porous medium through a crack. It is assumed that the crack, which has initial zero dimensions, can become wider and longer with a rise in the pressure. The displacement of the sides of the crack is determined on the basis of the theory of elasticity, taking account of the deformation properties of a saturated porous medium. The stress and the displacement are expressed in terms of two analytical Muskhelishvili functions and the complex filtration potential. A change in the volume of the porous medium leads to a discontinuity of the displacements at the feed contour, and to distortion in the filtration region. For a circular stratum, the dimensions of the crack and the mass flow rate of the liquid are determined in the first approximation. The region of values of the pressure in which there exists a stable equilibrium state of the open crack and a steady-state flow of the liquid is found.     

Extension of the Darcy�CForchheimer Law for Shear-Thinning Fluids and Validation via Pore-Scale Flow Simulations

Flow of non-Newtonian fluids through porous media at high Reynolds numbers is often encountered in chemical, pharmaceutical and food, as well as petroleum and groundwater engineering, and in many other industrial applications. Under the majority of operating conditions typically explored, the dependence of pressure drops on flow rate is non-linear and the development of models capable of describing accurately this dependence, in conjunction with non-trivial rheological behaviors, is of paramount importance. In this work, pore-scale single-phase flow simulations conducted on synthetic two-dimensional porous media are performed via computational fluid dynamics for both Newtonian and non-Newtonian fluids and the results are used for the extension and validation of the Darcy?CForchheimer law, herein proposed for shear-thinning fluid models of Cross, Ellis and Carreau. The inertial parameter ?? is demonstrated to be independent of the viscous properties of the fluids. The results of flow simulations show the superposition of two contributions to pressure drops: one, strictly related to the non-Newtonian properties of the fluid, dominates at low Reynolds numbers, while a quadratic one, arising at higher Reynolds numbers, is dependent on the porous medium properties. The use of pore-scale flow simulations on limited portions of the porous medium is here proposed for the determination of the macroscale-averaged parameters (permeability K, inertial coefficient ?? and shift factor ??), which are required for the estimation of pressure drops via the extended Darcy?CForchheimer law. The method can be applied for those fluids which would lead to critical conditions (high pressures for low permeability media and/or high flow rates) in laboratory tests.     

Temperature dependent mechanical behaviour of PVDF: Experiments and numerical modelling

The mechanical behaviour of Polyvinylidene Fluoride (PVDF) is analysed. To this end, tensile tests are performed on both smooth and notched specimens, for several values of the notch radius in order to set specific values of the stress triaxiality ratio in the net section. Tests were performed at various temperatures and at various strain rates. Experimental data together with fracture surface examinations by SEM allow the dependence of deformation and void growth processes on strain rate and temperature to be investigated. This experimental work was carried out in order to test the mechanics of porous media model. For each investigated temperature, constitutive relations take both porosity and strain rate sensitivity into account. The model is proposed for deformation leading to crazing. The material coefficients are optimised by imposing a continuous dependence on temperature.     

On the Domain of Validity of Brinkman’s Equation

An increasing number of articles are adopting Brinkman’s equation in place of Darcy’s law for describing flow in porous media. That poses the question of the respective domains of validity of both laws, as well as the question of the value of the effective viscosity μ ^e which is present in Brinkman’s equation. These two topics are addressed in this article, mainly by a priori estimates and by recalling existing analyses. Three main classes of porous media can be distinguished: “classical” porous media with a connected solid structure where the pore surface S _p is a function of the characteristic pore size l _p (such as for cylindrical pores), swarms of low concentration fixed particles where the pore surface is a function of the characteristic particle size l _s , and fiber-made porous media at low solid concentration where the pore surface is a function of the fiber diameter. If Brinkman’s 3D flow equation is valid to describe the flow of a Newtonian fluid through a swarm of fixed particles or fibrous media at low concentration under very precise conditions (Lévy 1983), then we show that it cannot apply to the flow of such a fluid through classical porous media.     

Dissociation rate constants for oxygen at temperatures up to 11000 K

Measurements of the absorption and vibrational temperature evolution of molecular oxygen behind the shock wave front in regimes with temperatures near the front varying over the range 4000–11 000 K made it possible to determine the dissociation rate constant for O₂ molecules under both thermal equilibrium and nonequilibrium conditions. The dependence of the dissociation rate constant on the ratio T _v /T, that is, on the deviation from the thermal equilibrium conditions, is demonstrated experimentally. The corresponding expression for the two-temperature dissociation rate constant is proposed.     

Pressure effects on viscosity and flow stability of polyethylene melts during extrusion

In the present work, the effects of pressure on the viscosity and flow stability of four commercial grade polyethylenes (PEs) have been studied: linear-low-density polyethylene copolymer, high-density polyethylene, metallocene polyethylenes with short-chain branches (mPE-SCB), and metallocene polyethylenes with long chain branching (mPE-LCB). The range of shear rates considered covers both stable and unstable flow regimes. “Enhanced exit-pressure” experiments have been performed attaining pressures of the order of 500×10⁵ Pa at the die exit. The necessary experimental conditions have been clearly defined so that dissipative heating can be neglected and pressure effects isolated. The results obtained show an exponential increase in both shear and entrance-flow pressure drop with mean pressure when shear rate is fixed and as long as flow is stable. These pressure effects are described by two pressure coefficients, β_S under shear and, β_E under elongation, that are calculated using time–pressure superposition and that are independent of mean pressure and flow rate. For three out of four PE, pressure coefficient values can be considered equal under shear and under elongation. However, for the mPE-LCB, the pressure coefficient under elongation is found to be about 30% lower than under shear. Flow instabilities in the form of oscillating flows or of upstream instabilities appear at lower shear rates as mean pressure increases. Nevertheless, the critical shear stress at which they are triggered remains independent of mean pressure. Moreover, it is found that the β_S values obtained for stable flows do not differ much from the values obtained during upstream instability regimes, and differ really from pressure effects observed under oscillating flow and slip conditions.     

Effects of anisotropy in permeability on the two-phase flow and heat transfer in a porous cavity

This paper reports on the results of a numerical study of convection flow and heat transfer in a rectangular porous cavity filled with a phase change material under steady state conditions. The two vertical walls of the cavity are subject respectively to temperatures below and above the melting point of the PCM while adiabatic conditions are imposed on the horizontal walls. The porous medium is characterized by an anisotropic permeability tensor with the principal axes arbitrarily oriented with respect to the gravity vector. The problem is governed by the aspect ratioA, the Rayleigh numberRa, the anisotropy ratioR and the orientation angle θ of the permeability tensor. Attention is focused on these two latter parameters in order to investigate the effects of the anisotropic permeability on the fluid flow and heat transfer of the liquid/solid phase change process. The method of solution is based on the control volume approach in conjunction with the Landau-transformation to map the irregular flow domain into a rectangular one. The results are obtained for the flow field, temperature distribution, interface position and heat transfer rate forA=2.5,Ra=40, 0≤θ≤π, 0.25≤R≤4. It was found that the equilibrium state of the solid/liquid phase change process may be strongly influenced by the anisotropy ratioR as well as by the orientation angle θ of the permeability tensor. First, for a given set of parametersA,Ra andR, there exists an optimum orientation θ_max for which the flow strength, the liquid volume and the heat transfer rate are maximum. There also exists an orientation θ_min=θ_max+π/2 for which these quantities are minimum. Second, when an anisotropic medium is oriented along the optimum direction θ_max, an increase of the permeability component along that direction will increase the flow and heat transfer rate in a same order while an increase of the other permeability component only has a negligible effect. For the parameter ranges considered in the present study, it was found that the optimum direction is lying between the gravity vector and the dominant flow direction.     

Steady and Unsteady Heat Transfer in a Channel Partially Filled with Porous Media Under Thermal Non-Equilibrium Condition

Steady and pulsatile flow and heat transfer in a channel lined with two porous layers subject to constant wall heat flux under local thermal non-equilibrium (LTNE) condition is numerically investigated. To do this, a physical boundary condition in the interface of porous media and clear region of the channel is derived. The objective of this work is, first, to assess the effects of local solid-to-fluid heat transfer (a criterion indicating on departure from local thermal equilibrium (LTE) condition), solid-to-fluid thermal conductivity ratio and porous layer thickness on convective heat transfer in steady condition inside a channel partially filled with porous media; second, to examine the impact of pulsatile flow on heat transfer in the same channel. The effects of LTNE condition and thermal conductivity ratio in pulsatile flow are also briefly discussed. It is observed that Nusselt number inside the channel increases when the problem is tending to LTE condition. Therefore, careless consideration of LTE may lead to overestimation of heat transfer. Solid-to-fluid thermal conductivity ratio is also shown to enhance heat transfer in constant porous media thickness. It is also revealed that an increase in the amplitude of pulsation may result in enhancement of Nusselt number, while Nusselt number has a minimum in a certain frequency for each value of amplitude.     

Air-water two-phase flow in a helically coiled tube

Air-water flow has been studied in a helically coiled tube. The flow pattern transition between stratified and annular flow was examined, and a series of measurements were then taken in the annular flow regime. Local values of the liquid film thickness and liquid film flowrate around the tube periphery were obtained. The variations of these values around the periphery was similar. For most of the cases studied the liquid film flow rate was greatest on the inside of tbe bend, but in some results a subsidiary peak at the outside position was also obtained. There was little net entrained flow because of the centrifugal forces tending to deposit drops very quickly. Attempts to use correlations developed in vertical annular flow at a local position on the tube periphery were not very successful.     

PIV measurements combined with the motion tracking technique to analyze flow around a moving porous structure

To gain a better understanding of the fluid–structure interaction and especially when dealing with a flow around an arbitrarily moving body, it is essential to develop measurement tools enabling the instantaneous detection of moving deformable interface during the flow measurements. A particularly useful application is the determination of unsteady turbulent flow velocity field around a moving porous fishing net structure which is of great interest for selectivity and also for the numerical code validation which needs a realistic database. To do this, a representative piece of fishing net structure is used to investigate both the Turbulent Boundary Layer (TBL) developing over the horizontal porous moving fishing net structure and the turbulent flow passing through the moving porous structure. For such an investigation, Time Resolved PIV measurements are carried out and combined with a motion tracking technique allowing the measurement of the instantaneous motion of the deformable fishing net during PIV measurements. Once the two-dimensional motion of the porous structure is accessed, PIV velocity measurements are analyzed in connection with the detected motion. Finally, the TBL is characterized and the effect of the structure motion on the volumetric flow rate passing though the moving porous structure is clearly demonstrated.     

Effect of molecular weight and shear on phase diagram of PS/PVME blend

The molecular weight dependence of the lower critical solution temperature of polystyrene (PS) and poly(vinyl methyl ether) blends was studied by laser light transmission. The temperature of phase separation was found to decrease with increasing PS molecular weight. In the steady shear flow conditions and near the critical temperature, shear was found to enhance fluctuations of concentration at relatively small shear rates, whereas it suppresses such fluctuations at high shear rates. The shift in Flory-interaction parameter Δ_χ was calculated from experimental data and its sign was used to predict shear-induced mixing or shear-induced demixing under flow field. The obtained experimental results were compared to Criado-Sancho et al. and Clarke-Mcleish models.     

Heat Transfer Study of the Hencken Burner Flame

The Hencken burner flame is often used in combustion laser diagnostics as a calibration flame because of its near adiabatic condition. For a fast burning H₂ flame, it can tolerate high flow rate and the flame is indeed near adiabatic; however, for a slow burning CH₄ flame, the flow rate is not always high enough to maintain near adiabatic conditions. The heat transfer of the H₂ and CH₄ Hencken burner flames are studied numerically and experimentally. Three heat loss mechanisms are analyzed: the burner surface radiation, the hot gas radiation, and the convection heat transfer between the main flow and the co-flow. The surface radiation produces negligible temperature drop while the gas radiation and the convection heat loss contribute significant temperature drop. Reducing the co-flow rate can decrease the convection heat loss slightly. The temperature drop caused by the heat loss is inversely proportional to the main flow rate. Increasing the burner size and running the flame premixed mode can increase the flow rate and reduce the temperature deviation from the adiabatic equilibrium value. Based on the heat loss and temperature drop analysis, suggestions are given to maintain the flame at near adiabatic conditions.     

Flow properties of oven-curing high solids white enamels

For the characterization of the rheological behaviour of a white, high solids, oven-curing enamel for household equipment with regard to its industrial application, attention was focused on the flow aspects involved in laying down process and film formation phenomena, i.e. viscosity at very high and very low shear rate, structural build-up in rest conditions, presence and magnitude of elastic components. Hence, investigations were aimed at (1) the determination of the equilibrium flow curve; the application of the Shangraw-Grim-Mattocks model led to the evaluation of parameters to describe infinite shear-rate viscosity and yield stress; (2) obtaining information about particle aggregation state in shear conditions. Quemada model gave indication that, even at very high shear, particle aggregates are stable; (3) the determination of time-dependent behaviour: elastic components were found to be almost negligible; the Trapeznikov-Fedotova procedure allowed thixotropic build-up in rest conditions to be evaluated, as concerns both amount and kinetics. Remarkable flow features found were: differences in the temperature dependence of viscosity at low shear rate and of yield value for enamels formulated with different pigment volume concentrations and the strong effect of the pigment volume concentration on the initial rate of thixotropic build-up in rest conditions.     

Experimental investigation of integral and local flow field properties on rotating porous disc

Flow field through rotating porous disc was investigated with experimental fluid dynamics and compared with computational fluid dynamics. Open cell aluminum metal foam with 88% porosity was used. On rotating porous disc, integral measurements of static pressure difference in dependence of air volume flow rate were performed. Local measurements of velocity profiles close to disc circumference were performed with hot-wire anemometry. The airflow visualization method using smoke generator and digital camera was performed. Flow structures through porous disc were visualized at three different air volume flow rates. Numerical simulation of homogenous rotating porous disc was performed. Experimental and numerical results were compared. The results showed appropriate comparison of integral and local properties. The presented experimental approach can be used for the investigation and understanding of flow field phenomena on rotating porous materials. The proposed conclusions can be applied for variable applications on rotating porous materials.     

Numerical and experimental investigation of convective drying in unsaturated porous media with bound water

The heat and mass transfer in an unsaturated wet cylindrical porous bed packed with quartz particles was investigated theoretically for relatively low convective drying rates. Local thermodynamic equilibrium was assumed in the mathematical model describing the multi-phase flow in the unsaturated porous media using the energy and mass conservation equations to describe the heat and mass transfer during the drying. The drying model included convection and capillary transport of the free water, diffusion of bound water, and convection and diffusion of the gas. The numerical results indicated that the drying process could be divided into three periods, the temperature rise period, the constant drying rate period and the decreasing drying rate period. The numerical results agreed well with the experimental data verifying that the mathematical model can evaluate the drying performance of porous media for low drying rates. The effects of drying conditions such as the ambient temperature, the relative humidity, and the velocity of the drying air, on the drying process were evaluated by numerical solution.     

Cotransport of Suspended Colloids and Nanoparticles in Porous Media

The objective of this study is to develop a model for cotransport of colloids and nanoparticles (NPs) in porous media under two particle capture mechanisms. The particle capture rate is proportional to the capture probability, which is a function of retained concentration, called the filtration function. Laboratory bench-scale experiments of individual transport of NPs and colloidal-size kaolinite clay particles through packed columns produced breakthrough curves (BTCs) that monotonically increased with time and stabilised at some value lower than the injected concentration. We discuss the filtration function that corresponds to BTCs stabilising at the concentration lower than the injected value. This so-called binary filtration function incorporates two particle capture mechanisms. The analytical transport model with a binary filtration function was capable to fit successfully BTCs obtained from individual transport experiments using kaolinite and NPs conducted by Chrysikopoulos et al. (Transp Porous Med 119(1):181–204, 2017). Assuming that the electrostatic particle–solid matrix interaction and the fraction of the solid matrix surface area occupied by a single attached particle (kaolinite or NP) are the same for individual transport of either kaolinite particles or NPs and for simultaneous cotransport of kaolinite particles and NPs, the proposed binary filtration function was extended for the cotransport case. Although the breakthrough data from cotransport experiments with kaolinite particles and NPs have six degrees of freedom, the developed cotransport model successfully matches the BTCs by tuning two constants only. This validates the developed model for cotransport of two colloidal populations with different attachments and straining rates.