Lab-on-chip methodologies for the study of transport in porous media: energy applications

We present a lab-on-chip approach to the study of multiphase transport in porous media. The applicability of microfluidics to biological and chemical analysis has motivated much development in lab-on-chip methodologies. Several of these methodologies are also well suited to the study of transport in porous media. We demonstrate the application of rapid prototyping of microfluidic networks with approximately 5000 channels, controllable wettability, and fluorescence-based analysis to the study of multiphase transport phenomena in porous media. The method is applied to measure the influence of wettability relative to network regularity, and to differentiate initial percolation patterns from active flow paths. Transport phenomena in porous media are of critical importance to many fields and particularly in many energy-related applications including liquid water transport in fuel cells, oil recovery, and CO(2) sequestration.     

Pore-scale modeling of rarefied gas flow in fractal micro-porous media,using lattice Boltzmann method (LBM)

Due to the widespread use of rarefied gas flow in micro-porous media in industrial and engineering problems, a pore-scale modeling of rarefied gas flow through two micro-porous media with fractal geometries is presented, using lattice Boltzmann method. For this purpose, square- and circular-based Sierpinski carpets with fractal geometries are selected due to their inherent behavior for real porous media. Diffusive reflection slip model is used and developed for these porous media through this study. With this respect, the planar Poiseuille flow is selected as a benchmark and validated with the literature. The effect of Knudsen number (Kn) on the permeability is investigated and compared in each geometry. It is shown that as Knudsen number increases, the permeability will increase due to the gas slippage effect on the solid blocks. In addition, it is observed that the permeability is more sensitive to the gaseous flow behavior at the slip and beginning of transition flow regimes. At last, the permeability relationship with Knudsen number is presented with a higher coefficient of determination for both fractal geometries, showing that this relation is logarithmic.

     

Electroosmosis in homogeneously charged micro- and nanoscale random porous media

Electroosmosis in homogeneously charged micro- and nanoscale random porous media has been numerically investigated using mesoscopic simulation methods which involve a random generation-growth method for reproducing three-dimensional random microstructures of porous media and a high-efficiency lattice Poisson-Boltzmann algorithm for solving the strongly nonlinear governing equations of electroosmosis in three-dimensional porous media. The numerical modeling and predictions of EOF in micro- and nanoscale random porous media indicate that the electroosmotic permeability increases monotonically with the porosity of porous media and the increasing rate rises with the porosity as well; the electroosmotic permeability increases with the average solid particle size for a given porosity and with the bulk ionic concentration also; the proportionally linear relationship between the electroosmotic permeability and the zeta potential on solid surfaces breaks down for high zeta potentials. The present predictions agree well with the available experimental data while some results deviate from the predictions based on the macroscopic theories.     

Penetration into three-dimensional complex porous structures

The short-term uptake of a fluid by porous media is important in a number of processes, such as in coating and printing operations. We present a new model to predict short-term absorption into real pore geometries taking into account fluid properties, surface forces, and the complex pore geometry. Two assumptions are made to reduce the complexity of the situation: (1) the flow resistance between pores can be estimated from pore geometry or air permeability measurements, and (2) the volume of fluid in the constrictions between pores is small. Pores can be connected in any manner and can be in any arrangement. The absorption rates predicted by the model are compared to experimental values obtained with coating layers of plastic, kaolin, and calcium carbonate pigments. These coatings are characterized in terms of void fraction, pore size, contact angle, and permeability. The comparison is good for water and inks when the air permeability of the porous layer is used to determine the average resistance to flow in the sample. These resistance values are close to the values obtained from pore geometries estimated from particle packing simulations.     

Properties and performance of novel high-resolution/high-permeability ion-exchange media for protein chromatography

There is continued interest in the development of stationary phases for protein chromatography that can provide high resolution at elevated flow rates of the mobile phase. When using porous particles, resolution and dynamic binding capacity decline rapidly as the flow rate is increased. Monolithic columns have been developed to overcome these limitations. However, there are difficulties in manufacturing homogeneous larger scale monoliths. In this paper we investigate the morphology and performance characteristics of columns based on new ion exchangers obtained by mechanically disrupting continuous beds of acrylamido-based polymeric media. Near colloidal suspensions of loose particles obtained with this procedure can be flow-packed in ordinary chromatography columns resulting in beds of unexpectedly high hydraulic permeability. Columns up to 2.2 cm in diameter were studied with both Q and S functionalized media. The hydraulic permeability and interparticle porosity of these columns were rather high. The permeabilities of the S and Q media were 1.5 x 10(-13) and 2.4 x 10(-13) m2, respectively, while the corresponding porosities were 60 and 70%. These porosity values are similar to those of monoliths, suggesting that these particles assemble under flow to give high-porosity bridged structures. The structure of these packed beds was further characterized by embedding small packed columns in resins and obtaining sections for microscopic observation. The sections reveal the presence of small aggregates of non-porous 1-3 microm particles, surrounded by flow channels several micrometers in size. The height equivalent to a theoretical plate under isocratic and gradient elution conditions and the dynamic binding capacity were determined for several proteins and were found to be virtually independent of flow.     

Wettability control and patterning of PDMS using UV-ozone and water immersion

We demonstrate a simple method to tune and pattern the wettability of polydimethylsiloxane (PDMS) to generate microfluidic mimics of heterogeneous porous media. This technique allows one to tailor the capillary forces at different regions within the PDMS channel to mimic multi-phase flow in oil reservoirs. In this method, UV-ozone treatment is utilized to oxidize and hydrophilize the surface of PDMS. To maintain a stable surface wettability, the oxidized surfaces are immersed in water. Additionally, the use of a photomask makes it convenient to pattern the wettability in the porous media. A one-dimensional diffusive reaction model is established to understand the UV-ozone oxidation as well as hydrophobic recovery of oxidized PDMS surfaces. The modeling results show that during UV-ozone, surface oxidation dominates over diffusion of low-molecular-weight (LMW) species. However, the diffusivity of LMW species plays an important role in wettability control of PDMS surfaces.     

Modeling electroosmotic and pressure-driven flows in porous microfluidic devices: zeta potential and porosity changes near the channel walls

This work presents analytical solutions for both pressure-driven and electroosmotic flows in microchannels incorporating porous media. Solutions are based on a volume-averaged flow model using a scaling of the Navier-Stokes equations for fluid flow. The general model allows analysis of fluid flow in channels with porous regions bordering open regions and includes viscous forces, permitting consideration of porosity and zeta potential variations near channel walls. To obtain analytical solutions problems are constrained to the linearized Poisson-Boltzmann equation and a variation of Brinkman's equation [Appl. Sci. Res., Sect. A 1, 27 (1947); 1, 81 (1947)]. Cases include one continuous porous medium, two adjacent regions of different porosities, or one open channel adjacent to a porous region, and the porous material may have a different zeta potential than that of the channel walls. Solutions are described for two geometries, including flow between two parallel plates or in a cylinder. The model illustrates the relative importance of porosity and zeta potential in different regions of each channel.     

Fluctuation-induced dynamics of multiphase liquid jets with ultra-low interfacial tension

Control of fluid dynamics at the micrometer scale is essential to emulsion science and materials design, which is ubiquitous in everyday life and is frequently encountered in industrial applications. Most studies on multiphase flow focus on oil-water systems with substantial interfacial tension. Advances in microfluidics have enabled the study of multiphase flow with more complex dynamics. Here, we show that the evolution of the interface in a jet surrounded by a co-flowing continuous phase with an ultra-low interfacial tension presents new opportunities to the control of flow morphologies. The introduction of a harmonic perturbation to the dispersed phase leads to the formation of interfaces with unique shapes. The periodic structures can be tuned by controlling the fluid flow rates and the input perturbation; this demonstrates the importance of the inertial effects in flow control at ultra-low interfacial tension. Our work provides new insights into microfluidic flows at ultra-low interfacial tension and their potential applications.     

Porous media: A model system for the physics of disorder

We discuss several examples of applications of soft matter and statistical physics to various problems related to porous and fractured media. The structure of porous media often displays multiple scale features which can be analysed by neutron diffraction or electron microscopy and can be approached by fractal models. Such characteristics are also observed on the surface of natural fractures.The relation between various transport parameters such as permeability or conductivity introduce characteristic microscopic length scales which can sometimes be independently determined. In some cases, if the porous medium get clogged or if the number of flow channels or fractures is low enough so that threshold effects appear which can be analysed in terms of percolation models. Disorder physics approaches are particularly useful to analyse non miscible diphasic flows in some cases in which multiscale heterogeneities of the fluid mixture composition appear. This is for instance the case for very slow non wetting invasions and of the fast injection of a low viscosity fluid : these processes can be described respectively by the “invasion percolation” and the “diffusion limited aggregation” models. Finally tracer dispersion provides an application of random walk models to disordered systems : examples of the response of this measurement in partly saturated and double porosity media are presented.     

Effect of capillary geometry on predicting electroosmotic volumetric flowrates in porous or fibrous media

Electrokinetic-based methods are used in a variety of applications including drug delivery and separation of biomolecules, among others. Many of these applications feature a fibrous or a porous medium that can be modeled by using capillary bundle models to predict the behavior of the electroosmotic flow within the particular system. The role of geometry in predicting volumetric flowrates in porous media is investigated by modeling the electroosmotic flow in idealized capillaries of rectangular, cylindrical, and annular geometries. This is achieved by the coupling of electrostatics and continuum hydrodynamics to obtain analytical expressions that govern the electrokinetically - driven volumetric flow within these idealized capillary geometries. A previous study developed a model to compare the cylindrical and annular capillary geometries by utilizing two methods that compare the areas of the two geometries. The methods used in this previous work will also be used in the present contribution to compare the volumetric flowrates in the cylindrical and annular capillaries with a rectangular capillary. Illustrative results will be presented to aid in the understanding of the influence of the various geometrical and electrostatic parameters that arise from the analysis of these volumetric flowrates. It was found that the electroosmotic volumetric flowrates are significantly affected by the capillary geometry.     

Computational study of a latent heat thermal energy storage system enhanced by highly conductive metal foams and heat pipes

Numerical simulations are performed to analyze the thermal characteristics of a latent heat thermal energy storage system with phase change material embedded in highly conductive porous media. A network of finned heat pipes is also employed to enhance the heat transfer within the system. ANSYS-FLUENT 19.0 is used to create a transient multiphase computational model to simulate the thermal behavior of the storage unit. Copper foam is the porous medium used to enhance the heat transfer and is impregnated with the phase change material, potassium nitrate (KNO₃). The effects of the porosity of the metal foam and the quantity of heat pipes on the thermal characteristics of storage unit have been investigated. The results indicated that increasing the quantity of the embedded heat pipes leads to drastic acceleration of both charging and discharging process. Impregnating the copper foam with potassium nitrate phase change material significantly affects the total charging and discharging times of the storage unit. It was shown that the porosity of the metal foam plays a key role in the thermal behavior of the system during the charging and discharging processes.

     

The rheology of pseudoplastic fluids in porous media using network modeling

This paper considers the rheology of pseudoplastic (shear thinning) fluids in porous media. The central problem studied is the relationship between the viscometric behavior of the polymer solution and its observed behavior in the porous matrix. In the past, a number of macroscopic approaches have been applied, usually based on capillary bundle models of the porous medium. These simplified models have been used along with constitutive equations describing the fluid behavior (usually of power law type) to establish semiempirical macroscopic equations describing the flow of non-Newtonian fluids in porous media. This procedure has been reasonably successful in correlating experimental results on the flow of polymer solutions through both consolidated and unconsolidated porous materials. However, it does not allow an interpretation of polymer flow in porous media in terms of the flows on a microscopic scale; nor does it allow us to predict changes in macroscopic behavior resulting from variations at a microscopic level in the characteristics of the porous medium such as pore size distribution. In this work, we use a network approach to the modeling of non-Newtonian rheology, in order to understand some of the more detailed features of polymjer flow in porous media. This approach provides a mathematical bridge between the behavior of the non-Newtonian fluid in a single capillary and the macroscopic behavior as deduced from the pressure drop-flow rate relation across the whole network model. It demonstrates the importance of flow redistribution within the elements of the capillary network as the overall pressure gradient varies. As an example of a pseudoplastic fluid in a porous medium, we consider the flow of xanthan biopolymer. This polymer is important as a displacing fluid viscosifier in enhanced oil recovery applications and, for that reason, a considerable amount of experimental data has been published on the flow of xanthan solutions in various porous media.     

Dynamic pore network model of surface heterogeneity in brine-filled porous media for carbon sequestration

Trapping of carbon in deep underground brine-filled reservoirs is a promising approach for the reduction of atmospheric greenhouse gas emissions. However, estimation of the amount of carbon dioxide (CO(2)) that can be captured in a given reservoir and the long-term storage stability remain a challenge. One difficulty lies in the estimation of local capillary pressure effects that arise from mineral surface heterogeneity inherent in underground geological formations. As a preliminary step to address this issue, we have performed dynamic pore network modelling (PNM) simulations of two-phase immiscible flow in two-dimensional structured porous media with contact angle heterogeneity under typical reservoir conditions. We begin by characterizing the network with a single, uniform contact angle. We then present saturation patterns for networks with homogeneous and heterogeneous contact angles distributions, based on two common reservoir minerals: quartz and mica, both of which have been well-characterized experimentally for their brine-CO(2) contact angles. At lower flow rates, we found moderately higher saturations for the heterogeneous networks than for the homogeneous ones. To characterize the fingering patterns, we have introduced R as the ratio of filled throats to the total network saturation. Based on this measure, the heterogeneous networks demonstrated thicker fingering patterns than the homogeneous networks. The computed saturation patterns demonstrate the importance of considering surface heterogeneity in pore-scale modelling of deep saline aquifers.     

Study of the permeability characteristics of porous media with methane hydrate by pore network model

The permeability in the methane hydrate reservoir is one of the key parameters in estimating the gas production performance and the flow behavior of gas and water during dissociation. In this paper, a three-dimensional cubic pore-network model based on invasion percolation is developed to study the effect of hydrate particle formation and growth habit on the permeability. The variation of permeability in porous media with different hydrate saturation is studied by solving the network problem. The simulation results are well consistent with the experimental data. The proposed model predicts that the permeability will reduce exponentially with the increase of hydrate saturation, which is crucial in developing a deeper understanding of the mechanism of hydrate formation and dissociation in porous media.     

Transportation characteristics of HPAM solution in the micro-fractures

In order to improve the success rate of practical application of an organic polymer conformance control agent in a fractured low permeability reservoir, the transportation characteristics of hydrolyzed polyacrylamide (HPAM) solution in micro-fracture is systematically studied by displacement experiments using the visual model and fractured cores. The results are as follows: (1)There is no threshold pressure for HPAM solution transportation in micro-fracture, and its flow is linear at a certain range of flow rate. (2) There is an appropriate corresponding range between the viscosity of HPAM solution and fracture permeability, which can provide guidance for screening out the appropriate molecular weight or concentration to obey the permeability. (3)Flow resistance of HPAM solution in the micro-fracture is lower than that in porous media at the same water phase permeability, but the difference degree of flow resistance is reduced with the increase of permeability. (4)When the aperture of the lower permeability fractured core in the parallel combination is high, the shunt flow volume of HPAM solution is linear. When the aperture is low, the shunt flow volume is nonlinear. The experimental results can provide important guidance for optimization and application of the polymer conformance control agent in a fractured low permeability reservoir.     

Study on the flow resistance of the dispersion system of deformable preformed particle gel in porous media using LBM-DEM-IMB method

After being injected into the porous media, the dispersion system of preformed particle gel (PPG) tends to enter high permeability regions and block water channeling passages, which forces the subsequent water to turn to the low permeability regions and thus increases sweep efficiency and enhances oil recovery. However, it is still unclear about the influence factors and the mechanisms how PPG increases water flow resistance, which limits the application of PPG in more oilfields. Therefore, the paper combines the lattice Boltzmann method (LBM), the discrete element method (DEM) and the improved immersed moving boundary (IMB) method to simulate the migration of deformable PPG in porous media. On the basis, the paper quantitatively analyzes the variation law of displacement pressure across the porous media and discusses the influence factors such as the PPG diameter, elasticity modulus and the number concentration. Results indicate that, because of the friction and retention of PPG in pore-throat, the displacement pressure across the porous media during PPG flooding is much higher than that during water flooding. In other words, the existence of PPG increases the flow resistance of injected water. Besides, the displacement pressure is always fluctuant resulting from the continuous process of PPG migration, retention, deformation and remigration. Influence factor analysis shows that the incremental value and fluctuation degree of flow resistance increase with the PPG diameter, elasticity modulus and the number concentration. The study not only provides useful reference for future PPG flooding, but also benefits the development of deformable particle flow theory.