Topological phase transition in 2D porous media flows
NASA Astrophysics Data System (ADS)
Waisbord, Nicolas; Stoop, Norbert; Kantsler, Vasily; Guasto, Jeffrey S.; Dunkel, Jorn; Guasto Team; Dunkel Team; Kantsler Team
2015-11-01
Since the establishment of Darcy's law, analysis of porous-media flows has focused primarily on linking macroscopic transport properties, such as mean flow rate and dispersion, to the pore statistics of the material matrix. Despite intense efforts to understand the fluid velocity statistics from the porous-media structure, a qualitative and quantitative connection remains elusive. Here, we combine precisely controlled experiments with theory to quantify how geometric disorder in the matrix affects the flow statistics and transport in a quasi-2D microfluidic channel. Experimentally measured velocity fields for a range of different microstructure configurations are found to be in excellent agreement with large-scale numerical simulations. By successively increasing the matrix disorder, we study the transition from periodic flow structures to transport networks consisting of extended high-velocity channels. Morse-Smale complex analysis of the flow patterns reveals a topological phase transition that is linked to a qualitative change in the physical transport properties. This work demonstrates that topological flow analysis provides a mathematically well-defined, broadly applicable framework for understanding and quantifying fluid transport in complex geometries.
Chaotic advection in 2D anisotropic porous media
NASA Astrophysics Data System (ADS)
Varghese, Stephen; Speetjens, Michel; Trieling, Ruben; Toschi, Federico
2015-11-01
Traditional methods for heat recovery from underground geothermal reservoirs employ a static system of injector-producer wells. Recent studies in literature have shown that using a well-devised pumping scheme, through actuation of multiple injector-producer wells, can dramatically enhance production rates due to the increased scalar / heat transport by means of chaotic advection. However the effect of reservoir anisotropy on kinematic mixing and heat transport is unknown and has to be incorporated and studied for practical deployment in the field. As a first step, we numerically investigate the effect of anisotropy (both magnitude and direction) on (chaotic) advection of passive tracers in a time-periodic Darcy flow within a 2D circular domain driven by periodically reoriented diametrically opposite source-sink pairs. Preliminary results indicate that anisotropy has a significant impact on the location, shape and size of coherent structures in the Poincare sections. This implies that the optimal operating parameters (well spacing, time period of well actuation) may vary strongly and must be carefully chosen so as to enhance subsurface transport. This work is part of the research program of the Foundation for Fundamental Research on Matter (FOM), which is part of Netherlands Organisation for Scientific Research (NWO). This research program is co-financed by Shell Global Solutions International B.V.
Gint2D-T2 correlation NMR of porous media
NASA Astrophysics Data System (ADS)
Zhang, Yan; Blümich, Bernhard
2015-03-01
The internal magnetic field gradient induced in porous media by magnetic susceptibility differences at material interfaces impacts diffusion measurements in particular at high magnetic field and can be used to probe the pore structure. Insight about the relationship between pore space and internal gradient Gint can be obtained from 2D Laplace NMR experiments. When measuring distributions of transverse relaxation times T2 in fluid filled porous media, relaxation and diffusion in internal gradients arise simultaneously and data are often interpreted with the assumption that one or the other parameter be constant throughout the sample. To examine this assumption we measure correlations of the distributions of Gint2D and T2 by 2D Laplace NMR for three different kinds of samples, glass beads with different bead diameters saturated with water, glass beads filled with oil and water, and a wet mortar sample. For the first two samples the cases where either the internal gradient or diffusion dominates were examined separately in order to better understand the relationship between Gint and D. These results are useful for assessing the impact of internal gradients and diffusion in unknown samples, such as the mortar sample. The experiments were performed at different magnetic field strengths corresponding to 300 MHz and 700 MHz 1H Larmor frequency to identify the impact of the magnetic field on the internal gradient. Subsequently, spatially resolved Gint2D-T2 maps were obtained to study the sample heterogeneity.
Gint2D-T2 correlation NMR of porous media.
Zhang, Yan; Blümich, Bernhard
2015-03-01
The internal magnetic field gradient induced in porous media by magnetic susceptibility differences at material interfaces impacts diffusion measurements in particular at high magnetic field and can be used to probe the pore structure. Insight about the relationship between pore space and internal gradient G(int) can be obtained from 2D Laplace NMR experiments. When measuring distributions of transverse relaxation times T(2) in fluid filled porous media, relaxation and diffusion in internal gradients arise simultaneously and data are often interpreted with the assumption that one or the other parameter be constant throughout the sample. To examine this assumption we measure correlations of the distributions of G(int)(2)D and T(2) by 2D Laplace NMR for three different kinds of samples, glass beads with different bead diameters saturated with water, glass beads filled with oil and water, and a wet mortar sample. For the first two samples the cases where either the internal gradient or diffusion dominates were examined separately in order to better understand the relationship between G(int) and D. These results are useful for assessing the impact of internal gradients and diffusion in unknown samples, such as the mortar sample. The experiments were performed at different magnetic field strengths corresponding to 300 MHz and 700 MHz (1)H Larmor frequency to identify the impact of the magnetic field on the internal gradient. Subsequently, spatially resolved Gint(2)D-T(2) maps were obtained to study the sample heterogeneity. PMID:25723135
Comparison of unstable water infiltration in porous media in 2D and 3D experiments
NASA Astrophysics Data System (ADS)
Schütz, C.; Neuweiler, I.; Lehmann, P.; Papafotiou, A.; Vontobel, P.; Hartmann, S.
2010-05-01
Water infiltration into unsaturated soil is an important process for groundwater recharge and thus for water balance of natural hydrosystems. The characteristics of infiltration patterns depend on porous media properties and initial moisture content. Infiltration fronts into soil can be unstable in layered media with fine over dry coarse material. To predict arrival times of infiltration fronts and average water content in upscaled models, it is necessary to understand occurrence of instabilities. The unstable flow behavior is not captured by standard models and finger characteristics have mostly been investigated experimentally. Most experiments in the past were carried out in 2D setups and it is not clear how the results of such studies relate to real 3D systems. The aim of this study is to compare development and finger characteristics of unstable infiltration in 2D and 3D setups. We carried out laboratory experiments on fast infiltration in 2D and 3D setups and measured water content in porous media with neutron transmission technology at the NEUTRA beam line at the Paul Scherrer Institute, Switzerland. The 2D experiments were carried out in a glass sandbox (260 mm high, 75 mm wide and 11 mm deep). For the 3D experiments aluminum cylindrical column (150 mm in height and 100 mm in diameter) were used. Both columns were filled homogeneously with coarse quartz sand (grain size 0.7 - 1.2 mm) below fine sand layer (0.1 - 0.3 mm) of 20 - 30 mm thickness. Two dimensional projection images of water content with spatial resolution of 125 microns were deduced from neutron images every 2 second. For the 3D setup water content distribution was reconstructed in 3D to monitor water content inside the fingers over time. Water content and finger-width (15 - 23 mm) were similar for 2D and 3D setups. In both cases water content was maximum when the front passes and was decreasing afterwards (indicating "overshoot" behavior). Also the water content difference between values after
2D experiments for characterizing solute dispersion in unsaturated heterogeneous porous media
NASA Astrophysics Data System (ADS)
Jimenez-Martinez, J.; De Anna, P.; Turuban, R.; Tabuteau, H.; Le Borgne, T.; Meheust, Y.
2012-12-01
The unsaturated zone plays a key role in the transfer of chemical elements from the surface to the subsurface. Yet, predicting the transport of chemical species through unsaturated porous media is still an open issue. The distribution of water and air clusters creates flow paths that are controlled by the water saturation, with the feature of a large velocity distribution. As saturation decreases, very low velocity zones in regions of trapped fluid coexist with connected fluid clusters with relatively high velocities. As a consequence the dispersion of solute elements strongly depends on the saturation degree. Numerical simulations of unsaturated flows at the pore scale are feasible, but to our knowledge no simulation of solute transport in the water phase during two-phase flow has been achieved yet. Due to technical difficulties, there also exists relatively few laboratory experiments that allow for visualization and quantification of unsaturated flow and transport at the pore scale. We have developed a two-dimensional (2D) horizontal set up, built by lithographic technique and in which a joint injection of the two phases (wetting and non-wetting) provides a controlled homogeneous saturation in the medium. The simultaneous precise measurement of the flow field, the spatial distribution of water and air, and the 2D tracer concentration field, as well as breakthrough curves at different locations, are used to investigate the relationship between the flow field complexity (velocity distribution and its correlation properties) and dispersion properties. Experimental results show non-Fickian transport behaviors, characterized by heavy tailed breakthrough curves, whose characteristics depend on the average saturation.
Foam flows in 2D porous media: intermittency and bubble fragmentation
NASA Astrophysics Data System (ADS)
Meheust, Y.; Géraud, B.; Jones, S. A.; Dollet, B.; Cantat, I.
2014-12-01
Flowing foams are used in many engineering and technical applications. A well-known application is oil recovery. Another one is the remediation of polluted soils: the foam is injected into the ground in order to mobilize chemical species present in the medium. Apart from potential interesting physico-chemical and biochemical properties, foams have peculiar flow properties that applications might benefit of. In particular, viscous dissipation arises mostly from the contact zones between the soap films and the walls, which results in peculiar friction laws allowing the foam to invade narrow pores more efficiently than Newtonian fluids would. In most experimental studies no local information of the foam structure is possible, and only global quantities such as the effective viscosity can be measured. Using two-dimensional transparent flow cells, we have previously shown that foam structural [1] and elastic [2] effects significantly impact the flow of foams in porous media. We now present an investigation of foam flow through a two-dimensional (2D) porous medium consisting of circular obstacles positioned randomly in a Hele-Shaw cell (see figure). The foam structure is recorded in time by a video camera and subsequently analyzed by image processing, which provides us with the velocity field and spatial distribution of bubble sizes. The flow exhibits a rich phenomenology, including flow irreversibility, preferential flow paths, local flow intermittency/non-stationarity (despite the imposed permanent global flow rate). Moreover, the medium impacts the nature of the flowing fluid by selecting the bubble size through bubble fragmentation. We investigate how preferential flow paths and intermittency depend on the imposed global flow rate and foam quality (the water content), and show that the spatial distribution of bubble sizes is to some extent correlated with the velocity field. We furthermore measure the evolution, along the flow direction, of the probability density
Propagator-resolved 2D exchange in porous media in the inhomogeneous magnetic field.
Burcaw, Lauren M; Hunter, Mark W; Callaghan, Paul T
2010-08-01
We present a propagator-resolved 2D exchange spectroscopy technique for observing fluid motion in a porous medium. The susceptibility difference between the matrix and the fluid is exploited to produce an inhomogeneous internal magnetic field, causing the Larmor frequency to change as molecules migrate. We test our method using a randomly packed monodisperse 100 microm diameter glass bead matrix saturated with distilled water. Building upon previous 2D exchange spectroscopy work we add a displacement dimension which allows us to obtain 2D exchange spectra that are defined by both mixing time and spatial displacement rather than by mixing time alone. We also simulate our system using a Monte Carlo process in a random nonpenetrating monodisperse bead pack, finding good agreement with experiment. A simple analytic model is used to interpret the NMR data in terms of a characteristic length scale over which molecules must diffuse to sample the inhomogeneous field distribution. PMID:20554230
2D Electrical Resistivity Tomography surveys optimisation of the solutes transports in porous media.
NASA Astrophysics Data System (ADS)
Lekmine, G.; Pessel, M.; Auradou, H.
2009-04-01
Electrical resistivity tomography applied in borehole or cross-borehole is a method often used to follow the invasion process of pollutant [Daily, 1991]. The aim of this work is to test experimentally the electrode arrays and inversion process used to obtain a spatial representation of tracer propagation in porous media. Experiments were conducted in a plexiglas container with glass beads of 166 microns in diameter. The height of the container is 275 mm, its width 85 mm and its thickness 10 mm. 21 electrodes, equally spaced, are placed along each of the lateral sides of the porous medium : these electrodes are used to perform the electrical measurements. The porous medium is lightened from behind and a video camera records the propagation of the fluids. The fluid containing the tracer (i.e the pollutant) is a water solution containing a small amount of dye together with NaCl (0.5g/l up to 2.0g/l). The medium is first saturated by a water solution containing a slight concentration of NaCl so that its density is smaller than the injected fluid. An upward flow is first established, then the denser fluid is injected at the bottom and over the full width of the medium. In this way, the flow is stabilized by gravity avoiding the development of unstable fingers. Still, the fluids are miscible and a mixing front develops during the flow: in the present study, both the determination by optical and electrical imaging of the mean position of the front and its width are of interest. The comparison of the two techniques allows to study the ability of the inversion process to quantify the solute transport. The electrical measurements are acquired by a standard multi electrode system (IRIS Instruments) and the data are inverted with the Res2Dinv software which models the 2D distribution of conductivity contrasts. The obtained bulk conductivity can be related through Archie's law to fluid conductivity by the porosity and the cementation factor which have been experimentally
Jones, Edward H; Smith, Colin C
2005-12-01
This paper describes an investigation into non-equilibrium partitioning tracer transport and interaction with non-aqueous-phase liquid (NAPL) contaminated water-saturated porous media using a two-dimensional (2-D) physical modelling methodology. A fluorescent partitioning tracer is employed within a transparent porous model which when imaged by a CCD digital camera can provide full spatial tracer concentrations and tracer breakthrough curves. Quasi one-dimensional (1-D) benchmarking tests in models packed with various combinations of clean quartz sand and NAPL are described. These modelled residual NAPL saturations, S(n), of 0-15%. Results demonstrated that the fluorescent partitioning tracer was able to detect and quantify the presence of NAPL at low flow rates. At larger flow rates and/or higher NAPL saturations, the tracer increasingly underpredicted the NAPL volume as expected and this is attributed primarily to non-equilibrium partitioning. Despite little change in permeability, change in NAPL saturations from 4% to 8% resulted in significant NAPL saturation underestimates at the same flow rates implying coalescence of NAPL into wider separated but larger ganglia. A 2-D investigation of an idealised heterogeneous residual NAPL contaminated flow field indicated little permeability change in the NAPL contaminated zone and thus little flow bypassing, leading to reduced underpredictions of NAPL saturations than for equivalent quasi 1-D cases. This was attributed to increased 'sampling' of the NAPL by the tracer. The process is clearly visually identifiable from the experimental images. This rapid and relatively inexpensive experimental method is of value in laboratory studies of partitioning tracer behaviour in porous media; in particular, the ability to observe full field concentrations makes it valuable for the study of complex heterogeneous systems. PMID:16298415
Phenrat, Tanapon; Cihan, Abdullah; Kim, Hye-Jin; Mital, Menka; Illangasekare, Tissa; Lowry, Gregory V
2010-12-01
Concentrated suspensions of polymer-modified Fe(0) nanoparticles (NZVI) are injected into heterogeneous porous media for groundwater remediation. This study evaluated the effect of porous media heterogeneity and the dispersion properties including particle concentration, Fe(0) content, and adsorbed polymer mass and layer thickness which are expected to affect the delivery and emplacement of NZVI in heterogeneous porous media in a two-dimensional (2-D) cell. Heterogeneity in hydraulic conductivity had a significant impact on the deposition of NZVI. Polymer modified NZVI followed preferential flow paths and deposited in the regions where fluid shear is insufficient to prevent NZVI agglomeration and deposition. NZVI transported in heterogeneous porous media better at low particle concentration (0.3 g/L) than at high particle concentrations (3 and 6 g/L) due to greater particle agglomeration at high concentration. High Fe(0) content decreased transport during injection due to agglomeration promoted by magnetic attraction. NZVI with a flat adsorbed polymeric layer (thickness ∼30 nm) could not be transported effectively due to pore clogging and deposition near the inlet, while NZVI with a more extended adsorbed layer thickness (i.e., ∼70 nm) were mobile in porous media. This study indicates the importance of characterizing porous media heterogeneity and NZVI dispersion properties as part of the design of a robust delivery strategy for NZVI in the subsurface. PMID:21058703
A 2-D Pore-Network Model of the Drying of Single-Component Liquids in Porous Media
Yortsos, Yanic C.; Yiotis, A.G.; Stubos, A.K.; Boundovis, A.G.
2000-01-20
The drying of liquid-saturated porous media is typically approaching using macroscopic continuum models involving phenomenological coefficients. Insight on these coefficients can be obtained by a more fundamental study at the pore- and pore-network levels. In this report, a model based on pore-network representation of porous media that accounts for various process at the pore-scale is presented. These include mass transfer by advection and diffusion in the gas phase, viscous flow in liquid and gas phases and capillary effects at the gas-liquid menisci in the pore throats.
SAGE 2D and 3D Simulations of the Explosive Venting of Supercritical Fluids Through Porous Media
NASA Astrophysics Data System (ADS)
Weaver, R.; Gisler, G.; Svensen, H.; Mazzini, A.
2008-12-01
Magmatic intrusive events in large igneous provinces heat sedimentary country rock leading to the eventual release of volatiles. This has been proposed as a contributor to climate change and other environmental impacts. By means of numerical simulations, we examine ways in which these volatiles can be released explosively from depth. Gases and fluids cooked out of country rock by metamorphic heating may be confined for a time by impermeable clays or other barriers, developing high pressures and supercritical fluids. If confinement is suddenly breached (by an earthquake for example) in such a way that the fluid has access to porous sediments, a violent eruption of a non-magmatic mixture of fluid and sediment may result. Surface manifestations of these events could be hydrothermal vent complexes, kimberlite pipes, pockmarks, or mud volcanoes. These are widespread on Earth, especially in large igneous provinces, as in the Karoo Basin of South Africa, the North Sea off the Norwegian margin, and the Siberian Traps. We have performed 2D and 3D simulations with the Sage hydrocode (from Los Alamos and Science Applications International) of supercritical venting in a variety of geometries and configurations. The simulations show several different patterns of propagation and fracturing in porous or otherwise weakened overburden, dependent on depth, source conditions (fluid availability, temperature, and pressure), and manner of confinement breach. Results will be given for a variety of 2D and 3D simulations of these events exploring the release of volatiles into the atmosphere.
NASA Astrophysics Data System (ADS)
Geistlinger, Helmut; Ataei-Dadavi, Iman; Mohammadian, Sadjad; Vogel, Hans-Jörg
2015-11-01
We study the impact of pore structure and surface roughness on capillary trapping of nonwetting gas phase during imbibition with water for capillary numbers between 10-7 and 5 × 10-5, within glass beads, natural sands, glass beads monolayers, and 2-D micromodels. The materials exhibit different roughness of the pore-solid interface. We found that glass beads and natural sands, which exhibit nearly the same grain size distribution, pore size distribution, and connectivity, showed a significant difference of the trapped gas phase of about 15%. This difference can be explained by the microstructure of the pore-solid interface. Based on the visualization of the trapping dynamics within glass beads monolayers and 2-D micromodels, we could show that bypass trapping controls the trapping process in glass beads monolayers, while snap-off trapping controls the trapping process in 2-D micromodels. We conclude that these different trapping processes are the reason for the different trapping efficiency, when comparing glass beads packs with natural sand packs. Moreover, for small capillary numbers of 10-6, we found that the cluster size distribution of trapped gas clusters of all 2-D and 3-D porous media can be described by a universal power law behavior predicted from percolation theory. This cannot be expected a priori for 2-D porous media, because bicontinuity of the two bulk phases is violated. Obviously, bicontinuity holds for the thin-film water phase and the bulk gas phase. The snap-off trapping process leads to ordinary bond percolation in front of the advancing bulk water phase and is the reason for the observed universal power law behavior in 2-D micromodels with rough surfaces.
Lappala, E.G.; Healy, R.W.; Weeks, E.P.
1987-01-01
This report documents FORTRAN computer code for solving problems involving variably saturated single-phase flow in porous media. The flow equation is written with total hydraulic potential as the dependent variable, which allows straightforward treatment of both saturated and unsaturated conditions. The spatial derivatives in the flow equation are approximated by central differences, and time derivatives are approximated either by a fully implicit backward or by a centered-difference scheme. Nonlinear conductance and storage terms may be linearized using either an explicit method or an implicit Newton-Raphson method. Relative hydraulic conductivity is evaluated at cell boundaries by using either full upstream weighting, the arithmetic mean, or the geometric mean of values from adjacent cells. Nonlinear boundary conditions treated by the code include infiltration, evaporation, and seepage faces. Extraction by plant roots that is caused by atmospheric demand is included as a nonlinear sink term. These nonlinear boundary and sink terms are linearized implicitly. The code has been verified for several one-dimensional linear problems for which analytical solutions exist and against two nonlinear problems that have been simulated with other numerical models. A complete listing of data-entry requirements and data entry and results for three example problems are provided. (USGS)
NASA Astrophysics Data System (ADS)
Peters, C. A.; Crandell, L. E.; Um, W.; Jones, K. W.; Lindquist, W. B.
2011-12-01
Geochemical reactions in the subsurface can alter the porosity and permeability of a porous medium through mineral precipitation and dissolution. While effects on porosity are relatively well understood, changes in permeability are more difficult to estimate. In this work, pore-network modeling is used to estimate the permeability of a porous medium using pore and throat size distributions. These distributions can be determined from 2D Scanning Electron Microscopy (SEM) images of thin sections or from 3D X-ray Computed Tomography (CT) images of small cores. Each method has unique advantages as well as unique sources of error. 3D CT imaging has the advantage of reconstructing a 3D pore network without the inherent geometry-based biases of 2D images but is limited by resolutions around 1 μm. 2D SEM imaging has the advantage of higher resolution, and the ability to examine sub-grain scale variations in porosity and mineralogy, but is limited by the small size of the sample of pores that are quantified. A pore network model was created to estimate flow permeability in a sand-packed experimental column investigating reaction of sediments with caustic radioactive tank wastes in the context of the Hanford, WA site. Before, periodically during, and after reaction, 3D images of the porous medium in the column were produced using the X2B beam line facility at the National Synchrotron Light Source (NSLS) at Brookhaven National Lab. These images were interpreted using 3DMA-Rock to characterize the pore and throat size distributions. After completion of the experiment, the column was sectioned and imaged using 2D SEM in backscattered electron mode. The 2D images were interpreted using erosion-dilation to estimate the pore and throat size distributions. A bias correction was determined by comparison with the 3D image data. A special image processing method was developed to infer the pore space before reaction by digitally removing the precipitate. The different sets of pore
Marsden, S.S.
1986-07-01
In 1978 a literature search on selective blocking of fluid flow in porous media was done by Professor S.S. Marsden and two of his graduate students, Tom Elson and Kern Huppy. This was presented as SUPRI Report No. TR-3 entitled ''Literature Preview of the Selected Blockage of Fluids in Thermal Recovery Projects.'' Since then a lot of research on foam in porous media has been done on the SUPRI project and a great deal of new information has appeared in the literature. Therefore we believed that a new, up-to-date search should be done on foam alone, one which would be helpful to our students and perhaps of interest to others. This is a chronological survey showing the development of foam flow, blockage and use in porous media, starting with laboratory studies and eventually getting into field tests and demonstrations. It is arbitrarily divided into five-year time periods. 81 refs.
FLUID TRANSPORT THROUGH POROUS MEDIA
Fluid transport through porous media is a relevant topic to many scientific and engineering fields. Soil scientists, civil engineers, hydrologists and hydrogeologists are concerned with the transport of water, gases and nonaqueous phase liquid contaminants through porous earth m...
Shepodd, Timothy J.
2002-01-01
Highly crosslinked monolithic porous polymer materials for chromatographic applications. By using solvent compositions that provide not only for polymerization of acrylate monomers in such a fashion that a porous polymer network is formed prior to phase separation but also for exchanging the polymerization solvent for a running buffer using electroosmotic flow, the need for high pressure purging is eliminated. The polymer materials have been shown to be an effective capillary electrochromatographic separations medium at lower field strengths than conventional polymer media. Further, because of their highly crosslinked nature these polymer materials are structurally stable in a wide range of organic and aqueous solvents and over a pH range of 2-12.
Eyler, L.L.; Budden, M.J.
1985-03-01
The objective of this work is to assess prediction capabilities and features of the MAGNUM-2D computer code in relation to its intended use in the Basalt Waste Isolation Project (BWIP). This objective is accomplished through a code verification and benchmarking task. Results are documented which support correctness of prediction capabilities in areas of intended model application. 10 references, 43 figures, 11 tables.
Dillon, J.
1999-09-01
A 2.8-liter tube-shaped combustion vessel was constructed to study flame propagation and quenching in porous media. For this experiment, hydrogen-air flames propagating horizontally into abed of 6 mm diameter glass beads were studied. Measurements of pressure and temperature along the length of the tube were used to observe flame propagation of quenching. The critical hydrogen concentration for Hz-air mixtures was found to be 11.5%, corresponding to a critical Peclet number of Pe* = 37. This value is substantially less than the value of Pe* = 65 quoted in the literature, for example Babkin et al. (1991). It is hypothesized that buoyancy and a dependence of Pe on the Lewis number account for the discrepancy between these two results.
Regeneratively Cooled Porous Media Jacket
NASA Technical Reports Server (NTRS)
Mungas, Greg (Inventor); Fisher, David J. (Inventor); London, Adam Pollok (Inventor); Fryer, Jack Merrill (Inventor)
2013-01-01
The fluid and heat transfer theory for regenerative cooling of a rocket combustion chamber with a porous media coolant jacket is presented. This model is used to design a regeneratively cooled rocket or other high temperature engine cooling jacket. Cooling jackets comprising impermeable inner and outer walls, and porous media channels are disclosed. Also disclosed are porous media coolant jackets with additional structures designed to transfer heat directly from the inner wall to the outer wall, and structures designed to direct movement of the coolant fluid from the inner wall to the outer wall. Methods of making such jackets are also disclosed.
Characterization of Porous Medium Properties Using 2D NMR
NASA Astrophysics Data System (ADS)
Sun, Boqin; Dunn, Keh-Jim
2003-03-01
We have successfully applied the concept of 2D NMR to the characterization of properties of fluid-saturated porous medium. Using a two-windowed modified CPMG pulse sequence, we were able to explore the magnetic internal filed gradient distribution within the pore space of a fluid-saturated porous medium due to magnetic susceptibility contrast between the solid matrix and pore fluid. Similar scheme is used to identify and quantify different types of pore fluids, such as oil, water, and gas, based on the contrast in their diffusion coefficients. The magic angle spinning technique (MAS) can also be applied in the 2D NMR framework for delineating the chemical shift spectra of the pore fluids in a porous medium at different T1 or T2 relaxation times. The results can be displayed in a two-dimensional plot, with one axis being the T1 or T2 relaxation times, the other axis being the internal field gradient, diffusion coefficient, or chemical shift, and the third axis being the proton population. Our preliminary laboratory work indicates that the 2D NMR approach can be a powerful tool for the characterization of properties of fluid-saturated porous medium, such as fluid typing, oil viscosity determination, surface wettability, etc.
Natural convection in porous media
Prasad, V.; Hussain, N.A.
1986-01-01
This book presents the papers given at a conference on free convection in porous materials. Topics considered at the conference included heat transfer, nonlinear temperature profiles and magnetic fields, boundary conditions, concentrated heat sources in stratified porous media, free convective flow in a cavity, heat flux, laminar mixed convection flow, and the onset of convection in a porous medium with internal heat generation and downward flow.
Steady flow through porous media
Greenkorn, R.A.
1981-07-01
The movement of materials through porous media is of interest in many disciplines: in chemical engineering - adsorption, chromatography, filtration, flow in packed columns, ion exchange, reactor-engineering; in petroleum engineering - displacement of oil with gas, water and miscible solvents including surface-active agent solutions and description of reservoirs; in hydrology - movement of trace pollutants in water systems, recovery of water for drinking and irrigation, saltwater encroachment into freshwater reservoirs; in soil physics - movement of water, nutrients, and pollutants into plants; and in biophysics. This work reviews the fundamentals of steady flow through porous media. It discusses the pseudotransport coefficients permeability, capillary pressure, and dispersion and relates these coefficients to the geometry of porous media. It discusses single-fluid flow, multifluid immiscible flow, and multifluid miscible flow including the effects of heterogeneity, nonuniformity, and anisotropy of media. 104 references.
Neeper, Donald A.
1994-01-01
Methods for distributing gases throughout the interstices of porous materials and removing volatile substances from the interstices of porous materials. Continuous oscillation of pressures and flows results in increased penetration of the interstices by flowing gases and increased transport of gaseous components out of the interstices. The invention is particularly useful in soil vapor extraction.
Neeper, D.A.
1994-02-22
Methods are presented for distributing gases throughout the interstices of porous materials and removing volatile substances from the interstices of porous materials. Continuous oscillation of pressures and flows results in increased penetration of the interstices by flowing gases and increased transport of gaseous components out of the interstices. The invention is particularly useful in soil vapor extraction. 10 figures.
Basser, Peter J.
2008-12-05
The fields of MR in Porous Media and Neuroradiology have largely developed separately during the past two decades with little appreciation of the problems, challenges and methodologies of the other. However, this trend is clearly changing and possibilities for significant cross-fertilization and synergies are now being realized.
Parametric study of boiling heat transfer in porous media
Shi, B.; Jones, B.G.; Pan, C.
1996-04-01
Detailed numerical modeling and parametric variation studies were conducted on boiling heat transfer processes in porous deposits with emphasis on applications associated with light water nuclear power reactor systems. The processes of boiling heat transfer in the porous corrosion deposits typically involve phase changes in finite volumetric regions in the porous media. The study examined such processes in two porous media configurations, without chimneys (homogeneous porous structures) and with chimneys (heterogeneous porous structures). A 1-D model and a 2-D model were developed to simulate two-phase flows with phase changes, without dry-out, inside the porous media for both structural configurations. For closure of the governing equations, an empirical correlation of the evaporation rate for phase changes inside the porous media was introduced. In addition, numerical algorithms were developed to solve the coupled nonlinear equations of mass, momentum, energy, capillary pressure, and evaporation rate. The distributions of temperature, thermodynamic saturation, liquid pressure, vapor pressure, liquid velocity, and vapor velocity were predicted. Furthermore, the effects of heat flux, system pressure, porosity, particle diameter, chimney population density, chimney radius, and crud thickness on the all superheat, critical heat flux, and minimum saturation were examined. The predictions were found to be in good agreement with the available experimental results.
Nonlinear flow in porous media
NASA Astrophysics Data System (ADS)
Rojas, Sergio Jesus
1998-07-01
Numerical solutions of the Navier-Stokes equations in two-dimensional quasi-periodic and quasi-isotropic random media were obtained to analyze the local and large scale aspects of finite Reynolds number flow. For Reynolds number less than one, the results show a first correction to Darcy's law which is cubic in the Darcy (averaged) velocity, while for Reynolds number greater than one, the results are in agreement with Forchheimer equation. That is, the correction to Darcy's law is quadratic in the average (Darcy) velocity. The cubic correction to Darcy's law support Mei and Auriault's (1991) theoretical study, based on homogenization theory. In addition, the results show support to a unifying empirical equation describing fluid flow in porous media of similar structure, first proposed by Beavers and Sparrow (1969). Also, the results show agreement, except by a multiplicative constant, with Sangani and Acrivos (1982) equation for the drag on dilute array of cylinders.
Microscale transport in porous media
Rashidi, M.; Rinker, R.
1996-04-01
In-pore transport processes in homogeneous and heterogeneous porous media have been investigated using novel 3D imaging techniques. The experimental system consists of a clear column packed with clear particles and a refractive index-matched fluid seeded with fluorescent tracers and an organic solute dye. By illuminating the porous regions within the column with a planar sheet of laser beam, flow and transport processes through the porous medium can be observed microscopically, and qualitative and quantitative in-pore transport information can be obtained at a good resolution and a high accuracy. Fluorescent images are captured and recorded at every vertical plane location while sweeping back and forth across the test section. These digitized transport images are then analyzed and accumulated over a 3D volume within the column. This paper reports on pore-scale observations of velocity, chemical concentration, and fluxes. Tests were undertaken with two separate columns. One is a rectangular column for chemical transport and bioremediation studies in aqueous heterogeneous systems and the other is a cylindrical column for flow and transport investigations in nonaqueous homogeneous systems.
Upscaling diffusion waves in porous media
NASA Astrophysics Data System (ADS)
Valdés-Parada, Francisco J.; Álvarez Ramírez, José; Ochoa-Tapia, J. Alberto
2016-04-01
The aim of this work is to derive the effective-medium equations and to estimate the related effective diffusivities for diffusion waves in porous media. Effective diffusivities are estimated within the framework of the volume averaging method, where they are obtained from the solution of the associated closure problems in 2D and 3D periodic unit cells. The results showed that the transport of diffusion waves are governed by the diffusion and co-diffusion mechanisms of harmonic waves. In addition, numerical results showed that the effective diffusivities increase with frequency, while the effective co-diffusivities display a resonance-like behavior. Our results also indicate that geometry plays a more significant effect over the predictions of the co-diffusion coefficient at moderate frequencies and it mainly influences the predictions of the direct diffusivity at low frequencies (i . e .,ω∗ ≪ 1).
Effective permeabilities for model heterogeneous porous media
Otevo, C.; Rusinek, I. ); Saez, A.E. )
1990-01-01
This paper presents a technique to evaluate effective absolute permeabilities for heterogeneous porous media. The technique is based on a perturbation analysis of the equations of motion of a slightly compressible fluid in a homogeneous porous medium at low Reynolds numbers. The effective permeabilities can be calculated once the local geometry of the heterogeneous medium is specified. The technique is used to evaluate two- and three-dimensional effective vertical permeabilities in porous media with shale intercalations, including the case in which the porous matrix is anisotropic.
Multiphase flow in fractured porous media
Firoozabadi, A.
1995-02-01
The major goal of this research project was to improve the understanding of the gas-oil two-phase flow in fractured porous media. In addition, miscible displacement was studied to evaluate its promise for enhanced recovery.
Optimal implicit 2-D finite differences to model wave propagation in poroelastic media
NASA Astrophysics Data System (ADS)
Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.
2016-08-01
Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite differences (FD) to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10 kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P and slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High-order explicit FD can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here, we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.
Optimal implicit 2-D finite differences to model wave propagation in poroelastic media
NASA Astrophysics Data System (ADS)
Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.
2016-05-01
Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite-differences to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P, slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High order explicit finite-differences (FD) can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.
Upscaling flow and transport properties in synthetic porous media
NASA Astrophysics Data System (ADS)
Jasinski, Lukasz; Dabrowski, Marcin
2015-04-01
Flow and transport through the porous media has instances in nature and industry: contaminant migration in geological formations, gas/oil extraction from proppant filled hydraulic fractures and surrounding porous matrix, underground carbon dioxide sequestration and many others. We would like to understand the behavior of propagating solute front in such medium, mainly flow preferential pathways and the solute dispersion due to the porous medium geometry. The motivation of our investigation is to find connection between the effective flow and transport properties and porous media geometry in 2D and 3D for large system sizes. The challenge is to discover a good way of upscaling flow and transport processes to obtain results comparable to these calculated on pore-scale in much faster way. We study synthetic porous media made of densely packed poly-disperse disk-or spherical-shaped grains in 2D and 3D, respectively. We use various protocols such as the random sequential addition (RSA) algorithm to generate densely packed grains. Imposed macroscopic pressure gradient invokes fluid flow through the pore space of generated porous medium samples. As the flow is considered in the low Reynolds number regime, a stationary velocity field is obtained by solving the Stokes equations by means of finite element method. Void space between the grains is accurately discretized by using body-fitting triangular or tetrahedral mesh. Finally, pure advection of a front carried by the velocity field is studied. Periodicity in all directions is applied to microstructure, flow and transport processes. Effective permeability of the media can be calculated by integrating the velocity field on cross sections, whereas effective dispersion coefficient is deduced by application of centered moment methods on the concentration field of transported solute in time. The effective parameters are investigated as a function of geometrical parameters of the media, such as porosity, specific surface area
NASA Astrophysics Data System (ADS)
Li, Y.; Kazemifar, F.; Blois, G.; Christensen, K. T.
2015-12-01
Multiphase flow of water and supercritical carbon dioxide (CO2) in porous media is central to geological sequestration of CO2 into saline aquifers. However, our fundamental understanding of the coupled flow dynamics of CO2 and water in complex geologic media still remains limited, especially at the pore scale. Recently, studies have been carried out in 2D homogeneous models with the micro-PIV technique, yielding very interesting observations of pore-scale flow transport. The primary aim of this work is to leverage this experimental protocol to quantify the pore-scale flow of water and liquid/supercritical CO2 in 2D heterogeneous porous micromodels under reservoir-relevant conditions. The goal is to capture the dynamics of this multi-phase flow in a porous matrix that mimics the heterogeneity of natural rock. Fluorescent microscopy and the micro-PIV technique are employed to simultaneously measure the spatially-resolved instantaneous velocity field in the water and quantify the instantaneous spatial configuration of both phases. The results for heterogeneous micromodels will be presented and compared with those for homogeneous micromodels, yielding valuable insight into flow processes at the pore scale in natural rock.
Numerical modeling of seismogram envelopes in 2-D random media
NASA Astrophysics Data System (ADS)
Fehler, Michael
2002-11-01
Several portions of seismograms recorded from regional earthquakes cannot be easily explained as resulting from waves propagating along deterministic paths within the Earth. For example, seismic coda, which is the tail portion of the seismogram of an earthquake recorded at distances of less than 100 km, is considered as resulting from waves that are multiply scattered from random heterogeneities in the Earth's lithosphere. At greater distances, observations that the duration of the initial arriving wave packet is much longer than the source-time duration is explained as being due to multiple forward scattering along the path between the source and the receiver. To investigate these phenomena, we use a finite difference method to numerically simulate 2-D scalar-waves that propagate through random media characterized by a von Karman autocorrelation function. Such media are considered to be appropriate models for the random component of the structure of the Earth's lithosphere. We investigate the characteristics of the resulting wavefields and compare them with those of observed seismograms.
Fluid dynamics in porous media with Sailfish
NASA Astrophysics Data System (ADS)
Coelho, Rodrigo C. V.; Neumann, Rodrigo F.
2016-09-01
In this work we show the application of Sailfish to the study of fluid dynamics in porous media. Sailfish is an open-source software based on the lattice-Boltzmann method. This application of computational fluid dynamics is of particular interest to the oil and gas industry and the subject could be a starting point for an undergraduate or graduate student in physics or engineering. We built artificial samples of porous media with different porosities and used Sailfish to simulate the fluid flow through them in order to calculate their permeability and tortuosity. We also present a simple way to obtain the specific superficial area of porous media using Python libraries. To contextualise these concepts, we analyse the applicability of the Kozeny–Carman equation, which is a well-known permeability–porosity relation, to our artificial samples.
Explosion propagation in inert porous media.
Ciccarelli, G
2012-02-13
Porous media are often used in flame arresters because of the high surface area to volume ratio that is required for flame quenching. However, if the flame is not quenched, the flow obstruction within the porous media can promote explosion escalation, which is a well-known phenomenon in obstacle-laden channels. There are many parallels between explosion propagation through porous media and obstacle-laden channels. In both cases, the obstructions play a duel role. On the one hand, the obstruction enhances explosion propagation through an early shear-driven turbulence production mechanism and then later by shock-flame interactions that occur from lead shock reflections. On the other hand, the presence of an obstruction can suppress explosion propagation through momentum and heat losses, which both impede the unburned gas flow and extract energy from the expanding combustion products. In obstacle-laden channels, there are well-defined propagation regimes that are easily distinguished by abrupt changes in velocity. In porous media, the propagation regimes are not as distinguishable. In porous media the entire flamefront is affected, and the effects of heat loss, turbulence and compressibility are smoothly blended over most of the propagation velocity range. At low subsonic propagation speeds, heat loss to the porous media dominates, whereas at higher supersonic speeds turbulence and compressibility are important. This blending of the important phenomena results in no clear transition in propagation mechanism that is characterized by an abrupt change in propagation velocity. This is especially true for propagation velocities above the speed of sound where many experiments performed with fuel-air mixtures show a smooth increase in the propagation velocity with mixture reactivity up to the theoretical detonation wave velocity. PMID:22213663
Propagation of polymer slugs through porous media
Lecourtier, J.; Chauveteau, G.
1984-09-01
This paper describes an experimental and theoretical study of the mechanisms governing polymer slug propagation through porous media. An analytical model taking into account the macromolecule exclusion from pore walls is proposed to predict rodlike polymer velocity in porous media and thus the spreading out of polydispersed polymer slugs. Under conditions where this wall exclusion is maximum, i.e. at low shear rates and polymer concentrations, the experiments show that xanthan propagation is effectively predicted by this model. At higher flow rates and polymer concentrations, the effects of hydrodynamic dispersion and viscous fingering are analyzed. A new fractionation method for determining molecular weight distribution of polymers used in EOR is proposed.
2-D solution for drying with internal vaporization of anisotropic media
Perre, P.; Passard, J.
1999-01-01
A set of physical assumptions for a straightforward and efficient simulation of the drying of a porous medium submitted to a convective and/or radiative heating is proposed. All of the parameters used are physical (liquid migration coefficient, permeability values along the thickness and the length, thermal conductivity, and external changes), but the relevant geometrical properties (length and thickness of the slab) were also included in the analysis. The 2-D pressure field generated within the medium during the drying process is obtained by using an analytical expression. The treatment of the pressure field, especially for a strongly anisotropic medium, is an important feature of the model, which allows an analytical model for such a complicated porous medium as wood to be used across a wide range of drying conditions. The computer code developed from the proposed formulation permits a complete simulation of the drying process within a few seconds on a personal computer. Different configurations have been tested for both anisotropic (wood) and isotropic (light concrete) porous media. Agreement with the experimental results is reasonable in terms of the observed physical phenomena. For instance, the model highlights dependence of the duration of the first drying rate on both material properties and drying conditions. This new model can be used for a global physical characterization of products by curve fitting and the collated information for the design of dryers.
Propagation of dissolution/precipitation waves in porous media
Novak, C.F.
1992-01-01
The transport of a chemically reactive fluid through a permeable medium is governed by many classes of chemical interactions. Dissolution/precipitation (D/P) reactions are among the interactions of primary importance because of their significant influence on the mobility of aqueous ions. In general, D/P reactions lead to the propagation of coherent waves. This paper provides an overview of the types of wave phenomena observed in one-dimensional (1D) and two-dimensional (2D) porous media for systems in which mineral D/P is the dominant type of chemical reaction. It is demonstrated that minerals dissolve in sharp waves in 1D advection-dominated transport, and that these waves separate zones of constant chemical compositions in the aqueous and mineral phases. Analytical solutions based on coherence methods are presented for solving 1D advection-dominated transport problems with constant and variable boundary conditions. Numerical solutions of diffusion-dominated transport in porous media show that sharp D/P fronts occur in this system as well. A final example presents a simple dual-porosity system with advection in an idealized fracture and solute diffusion into an adjacent porous matrix. The example illustrates the delay of contaminant release from the 2D domain due to a combination of physical retardation and chemical retardation.
Full-waveform inversion in 2D VTI media
NASA Astrophysics Data System (ADS)
Kamath, Nishant
Full-waveform inversion (FWI) is a technique designed to produce a high-resolution model of the subsurface by using information contained in entire seismic waveforms. This thesis presents a methodology for FWI in elastic VTI (transversely isotropic with a vertical axis of symmetry) media and discusses synthetic results for heterogeneous VTI models. First, I develop FWI for multicomponent data from a horizontally layered VTI model. The reflectivity method, which permits computation of only PP reflections or a combination of PP and PSV events, is employed to model the data. The Gauss-Newton technique is used to invert for the interval Thomsen parameters, while keeping the densities fixed at the correct values. Eigenvalue/eigenvector decompostion of the Hessian matrix helps analyze the sensitivity of the objective function to the model parameters. Whereas PP data alone are generally sufficient to constrain all four Thomsen parameters even for conventional spreads, including PS reflections provides better constraints, especially for the deeper part of the model. Next, I derive the gradients of the FWI objective function with respect to the stiffness coefficients of arbitrarily anisotropic media by employing the adjoint-state method. From these expressions, it is straightforward to compute the gradients for parameters of 2D heterogeneous VTI media. FWI is implemented in the time domain with the steepest-descent method used to iteratively update the model. The algorithm is tested on transmitted multicomponent data generated for Gaussian anomalies in Thomsen parameters embedded in homogeneous VTI media. To test the sensitivity of the objective function to different model parameters, I derive an an- alytic expression for the Frechet kernel of FWI for arbitrary anisotropic symmetry by using the Born approximation and asymptotic Green's functions. The amplitude of the kernel, which represents the radiation pattern of a secondary source (that source describes a perturbation
Numerical study of porous media thermoelectric converter
Kosaka, Kenichirou; Yamada, Akira
1996-12-31
Thermoelectric conversion is direct conversion technology that has characteristics of being maintenance free. However, the efficiency of the conventional bulk semiconductor thermoelectric device is about 20% for ideal theoretical calculation, and less than 5% for an actual application. The efficiency is very low because the heat conduction in the device and the Joule loss are too large compared with the Peltier heat which is changed into the electric power. The thermoelectric device made by porous media is heated by the radiation and maintains a large temperature difference by the gas which passes in the porous device. Therefore, the influence of the heat conduction in the thermoelectric device is small and the improvement of the conversion efficiency can be attempted. In this paper, the authors report the calculated results and the performance of thermoelectric converter made with porous media.
Computation of streaming potential in porous media: Modified permeability tensor
NASA Astrophysics Data System (ADS)
Bandopadhyay, Aditya; DasGupta, Debabrata; Mitra, Sushanta K.; Chakraborty, Suman
2015-11-01
We quantify the pressure-driven electrokinetic transport of electrolytes in porous media through a matched asymptotic expansion based method to obtain a homogenized description of the upscaled transport. The pressure driven flow of aqueous electrolytes over charged surfaces leads to the generation of an induced electric potential, commonly termed as the streaming potential. We derive an expression for the modified permeability tensor, K↔eff, which is analogous to the Darcy permeability tensor with due accounting for the induced streaming potential. The porous media herein are modeled as spatially periodic. The modified permeability tensor is obtained for both topographically simple and complex domains by enforcing a zero net global current. Towards resolving the complicated details of the porous medium in a computationally efficient framework, the domain identification and reconstruction of the geometries are performed using adaptive quadtree (in 2D) and octree (in 3D) algorithms, which allows one to resolve the solid-liquid interface as per the desired level of resolution. We discuss the influence of the induced streaming potential on the modification of the Darcy law in connection to transport processes through porous plugs, clays and soils by considering a case-study on Berea sandstone.
Statistical fusion of two-scale images of porous media
NASA Astrophysics Data System (ADS)
Mohebi, Azadeh; Fieguth, Paul; Ioannidis, Marios A.
2009-11-01
The reconstruction of the architecture of void space in porous media is a challenging task, since porous media contain pore structures at multiple scales. Whereas past methods have been limited to producing samples with matching statistical behavior, the patterns of grey-level values in a measured sample actually say something about the unresolved details, thus we propose a statistical fusion framework for reconstructing high-resolution porous media images from low-resolution measurements. The proposed framework is based on a posterior sampling approach in which information obtained by low-resolution (MRI or X-ray) measurements is combined with prior models inferred from high-resolution microscopic data, typically 2D. In this paper, we focus on two-scale reconstruction tasks in which the measurements resolve only the large scale structures, leaving the small-scale to be inferred. The evaluation of the results generated by the proposed method shows the strong ability of the proposed method in reconstructing fine-scale structures positively correlated with the underlying ground truth. Comparing our method with the recent method of Okabe and Blunt [12], in which the measurements are also used in the reconstruction, we conclude that our method is more robust to the resolution of the measurement, and more closely matches the underlying fine-scale field.
Finite volume hydromechanical simulation in porous media
NASA Astrophysics Data System (ADS)
Nordbotten, Jan Martin
2014-05-01
Cell-centered finite volume methods are prevailing in numerical simulation of flow in porous media. However, due to the lack of cell-centered finite volume methods for mechanics, coupled flow and deformation is usually treated either by coupled finite-volume-finite element discretizations, or within a finite element setting. The former approach is unfavorable as it introduces two separate grid structures, while the latter approach loses the advantages of finite volume methods for the flow equation. Recently, we proposed a cell-centered finite volume method for elasticity. Herein, we explore the applicability of this novel method to provide a compatible finite volume discretization for coupled hydromechanic flows in porous media. We detail in particular the issue of coupling terms, and show how this is naturally handled. Furthermore, we observe how the cell-centered finite volume framework naturally allows for modeling fractured and fracturing porous media through internal boundary conditions. We support the discussion with a set of numerical examples: the convergence properties of the coupled scheme are first investigated; second, we illustrate the practical applicability of the method both for fractured and heterogeneous media.
Finite volume hydromechanical simulation in porous media
Nordbotten, Jan Martin
2014-01-01
Cell-centered finite volume methods are prevailing in numerical simulation of flow in porous media. However, due to the lack of cell-centered finite volume methods for mechanics, coupled flow and deformation is usually treated either by coupled finite-volume-finite element discretizations, or within a finite element setting. The former approach is unfavorable as it introduces two separate grid structures, while the latter approach loses the advantages of finite volume methods for the flow equation. Recently, we proposed a cell-centered finite volume method for elasticity. Herein, we explore the applicability of this novel method to provide a compatible finite volume discretization for coupled hydromechanic flows in porous media. We detail in particular the issue of coupling terms, and show how this is naturally handled. Furthermore, we observe how the cell-centered finite volume framework naturally allows for modeling fractured and fracturing porous media through internal boundary conditions. We support the discussion with a set of numerical examples: the convergence properties of the coupled scheme are first investigated; second, we illustrate the practical applicability of the method both for fractured and heterogeneous media. PMID:25574061
Fluid Flow Within Fractured Porous Media
Crandall, D.M.; Ahmadi, G.; Smith, D.H.; Bromhal, G.S.
2006-10-01
Fractures provide preferential flow paths to subterranean fluid flows. In reservoir scale modeling of geologic flows fractures must be approximated by fairly simple formulations. Often this is accomplished by assuming fractures are parallel plates subjected to an applied pressure gradient. This is known as the cubic law. An induced fracture in Berea sandstone has been digitized to perform numerical flow simulations. A commercially available computational fluid dynamics software package has been used to solve the flow through this model. Single phase flows have been compared to experimental works in the literature to evaluate the accuracy with which this model can be applied. Common methods of fracture geometry classification are also calculated and compared to experimentally obtained values. Flow through regions of the fracture where the upper and lower fracture walls meet (zero aperture) are shown to induce a strong channeling effect on the flow. This model is expanded to include a domain of surrounding porous media through which the flow can travel. The inclusion of a realistic permeability in this media shows that the regions of small and zero apertures contribute to the greatest pressure losses over the fracture length and flow through the porous media is most prevalent in these regions. The flow through the fracture is shown to be the largest contributor to the net flow through the media. From this work, a novel flow relationship is proposed for flow through fractured media.
Chemo-hydrodynamic patterns in porous media.
De Wit, A
2016-10-13
Chemical reactions can interplay with hydrodynamic flows to generate chemo-hydrodynamic instabilities affecting the spatio-temporal evolution of the concentration of the chemicals. We review here such instabilities for porous media flows. We describe the influence of chemical reactions on viscous fingering, buoyancy-driven fingering in miscible systems, convective dissolution as well as precipitation patterns. Implications for environmental systems are discussed.This article is part of the themed issue 'Energy and the subsurface'. PMID:27597788
OPM: The Open Porous Media Initiative
NASA Astrophysics Data System (ADS)
Flemisch, B.; Flornes, K. M.; Lie, K.; Rasmussen, A.
2011-12-01
The principal objective of the Open Porous Media (OPM) initiative is to develop a simulation suite that is capable of modeling industrially and scientifically relevant flow and transport processes in porous media and bridge the gap between the different application areas of porous media modeling, including reservoir mechanics, CO2 sequestration, biological systems, and product development of engineered media. The OPM initiative will provide a long-lasting, efficient, and well-maintained open-source software for flow and transport in porous media built on modern software principles. The suite is released under the GNU General Public License (GPL). Our motivation is to provide a means to unite industry and public research on simulation of flow and transport in porous media. For academic users, we seek to provide a software infrastructure that facilitates testing of new ideas on models with industry-standard complexity, while at the same time giving the researcher control over discretization and solvers. Similarly, we aim to accelerate the technology transfer from academic institutions to professional companies by making new research results available as free software of professional standard. The OPM initiative is currently supported by six research groups in Norway and Germany and funded by existing grants from public research agencies as well as from Statoil Petroleum and Total E&P Norge. However, a full-scale development of the OPM initiative requires substantially more funding and involvement of more research groups and potential end users. In this talk, we will provide an overview of the current activities in the OPM initiative. Special emphasis will be given to the demonstration of the synergies achieved by combining the strengths of individual open-source software components. In particular, a new fully implicit solver developed within the DUNE-based simulator DuMux could be enhanced by the ability to read industry-standard Eclipse input files and to run on
Design of 2D Porous Coordination Polymers Based on Metallacrown Units.
Atzeri, Corrado; Marchiò, Luciano; Chow, Chun Y; Kampf, Jeff W; Pecoraro, Vincent L; Tegoni, Matteo
2016-05-01
A 12-metallacrown-4 (MC) complex was designed and employed as the building block in the synthesis of coordination polymers, one of which is the first permanently porous MC architecture. The connection of the four-fold symmetric MC subunits by Cu(II) nodes led to the formation of 2D layers of metallacrowns. Channels are present in the crystalline architecture, which exhibits permanent porosity manifested in N2 and CO2 uptake capacity. PMID:26951956
Experimental seismic investigations for porous media characterization
NASA Astrophysics Data System (ADS)
Senechal, P.; Bordes, C.; Barrière, J.; Garambois, S.; Brito, D.
2009-12-01
During decades, seismic investigations were based on the hypothesis of elastic media which is generally appropriate for geological imaging. Nevertheless, the necessity to improve non destructive geophysical methods for reservoir characterization implies to take into account the role of pore fluids for the seismic propagation. On the one hand, the water saturation modifies porous media bulk moduli and density and implies strong phase velocity variations. On the other hand, mechanical couplings between solid and fluid phases which can have inertial or viscous origin, generates attenuation and dispersion phenomena. In order to study these phenomena, we develop laboratory experiments for the measurement of seismic propagation in controlled porous media. We choose to work in a broad band [100-10 000 Hz] which includes Biot's frequency by using different seismic source (launched ball, pendulum, piezoelectric spherical source). These investigations are performed on various media from suspension (fresh concrete) to continuous framework (sand). Seismic propagation is recorded by accelerometers, discussed in term of frequency content, velocity analyses, attenuation and dispersion and compared to theoretical models. Moreover, seismoelectric measurements are used to study fluid displacements (filtration) associated to the seismic propagation.
Scaling theory of drying in porous media
Tsimpanogiannis, I.N.; Yortsos, Y.C.; Poulou, S.; Kanellopoulos, N.; Stubos, A.K.
1999-04-01
Concepts of immiscible displacements in porous media driven by mass transfer are utilized to model drying of porous media. Visualization experiments of drying in two-dimensional glass micromodels are conducted to identify pore-scale mechanisms. Then, a pore network approach is used to analyze the advancing drying front. It is shown that in a porous medium, capillarity induces a flow that effectively limits the extent of the front, which would otherwise be of the percolation type, to a finite width. In conjuction with the predictions of a macroscale stable front, obtained from a linear stability analysis, the process is shown to be equivalent to invasion percolation in a stabilizing gradient. A power-law scaling relation of the front width with a diffusion-based capillary number is also obtained. This capillary number reflects the fact that drying is controlled by diffusion in contrast to external drainage. The scaling exponent predicted is compatible with the experimental results of Shaw [Phys Rev. Lett. {bold 59}, 1671 (1987)]. A framework for a continuum description of the upstream drying regimes is also developed. {copyright} {ital 1999} {ital The American Physical Society}
Uncertainty quantification for porous media flows
NASA Astrophysics Data System (ADS)
Christie, Mike; Demyanov, Vasily; Erbas, Demet
2006-09-01
Uncertainty quantification is an increasingly important aspect of many areas of computational science, where the challenge is to make reliable predictions about the performance of complex physical systems in the absence of complete or reliable data. Predicting flows of oil and water through oil reservoirs is an example of a complex system where accuracy in prediction is needed primarily for financial reasons. Simulation of fluid flow in oil reservoirs is usually carried out using large commercially written finite difference simulators solving conservation equations describing the multi-phase flow through the porous reservoir rocks. This paper examines a Bayesian Framework for uncertainty quantification in porous media flows that uses a stochastic sampling algorithm to generate models that match observed data. Machine learning algorithms are used to speed up the identification of regions in parameter space where good matches to observed data can be found.
Foam Generation in Homogeneous Porous Media
Gauglitz, Phillip A.; Friedman, F.; Kam, S. I.; Rossen, W. R.
2002-10-01
In steady gas-liquid flow in homogeneous porous media with surfactant present, there is often observed a critical injection velocity or pressure gradient ?grad p min? at which ?weak? or ?coarse? foam is abruptly converted into ?strong foam,? with reduction of one to two orders of magnitude in total mobility: i.e., ?foam generation.? Earlier research on foam generation is extended here with extensive data for a variety of porous media, permeabilities, gases (N2 and C02), surfactants, and temperatures. For bead and sandpacks, ?grad p min? scales like (1/k), where k is permeability, over 2 1/2 orders of magnitude in k; for consolidated media the relation is more complex. For dense C02 foam, ?grad p min? exists but can be less than 1 psi/ft. If pressure drop, rather than flow rates, is fixed, one observes and unstable regime between stable ?strong? and ?coarse? foam regimes; in the unstable regime ?grad p? is nonuniform in space or variable in time.
Probing porous media with gas diffusion NMR.
Mair, R W; Wong, G P; Hoffmann, D; Hurlimann, M D; Patz, S; Schwartz, L M; Walsworth, R L
1999-10-18
We show that gas diffusion nuclear magnetic resonance (GD-NMR) provides a powerful technique for probing the structure of porous media. In random packs of glass beads, using both laser-polarized and thermally polarized xenon gas, we find that GD-NMR can accurately measure the pore space surface-area-to-volume ratio, S/V rho, and the tortuosity, alpha (the latter quantity being directly related to the system's transport properties). We also show that GD-NMR provides a good measure of the tortuosity of sandstone and complex carbonate rocks. PMID:11543587
Reservoir performance in viscoelastic porous media
Rago, F.M.; Ohkuma, H.; Sepehrnoori, K.; Thompson, T.W.
1982-01-01
The mass balance equations for a two-phase two-component fluid system are written for viscoelastic porous media. The resulting equations are approximated by finite differences and the resulting numerical simulator is used to conduct a sensitivity study on the effects of uniaxial viscoelastic deformation in geopressured aquifers. Results of this study indicate that viscoelastic deformation may have considerable influence on the pressure maintenance of these aquifers. A numerical model of the geopressured aquifer in Brazoria County, Texas, is constructed and the numerical simulator is used to predict the ultimate recovery of solution gas from this viscoelastic geopressured aquifer.
Probing porous media with gas diffusion NMR
NASA Technical Reports Server (NTRS)
Mair, R. W.; Wong, G. P.; Hoffmann, D.; Hurlimann, M. D.; Patz, S.; Schwartz, L. M.; Walsworth, R. L.
1999-01-01
We show that gas diffusion nuclear magnetic resonance (GD-NMR) provides a powerful technique for probing the structure of porous media. In random packs of glass beads, using both laser-polarized and thermally polarized xenon gas, we find that GD-NMR can accurately measure the pore space surface-area-to-volume ratio, S/V rho, and the tortuosity, alpha (the latter quantity being directly related to the system's transport properties). We also show that GD-NMR provides a good measure of the tortuosity of sandstone and complex carbonate rocks.
Experiments with metallic and ceramic porous media
NASA Technical Reports Server (NTRS)
Frederking, T. H. K.; Abbassi, P.; Khandhar, P. K.; Luna, Jack
1988-01-01
Work in the area of mechano-caloric phenomena was initiated during 1988 with startup in the Summer 1988 period. The ideal system utilizing He-II super-phenomena is modeled readily, within the frame of thermodynamics energetics, using the concept of an ideal superleak. The real system however uses porous media of non-ideal pore-grain ingredients. The early phase of experimental and related modeling studies is outlined for the time period from Summer 1988 to the end of 1988.
Mathematical modeling heat and mass transfer processes in porous media
NASA Astrophysics Data System (ADS)
Akhmed-Zaki, Darkhan
2013-11-01
On late development stages of oil-fields appears a complex problem of oil-recovery reduction. One of solution approaches is injecting of surfactant together with water in the form of active impurities into the productive layer - for decreasing oil viscosity and capillary forces between ``oil-water'' phases system. In fluids flow the surfactant can be in three states: dissolved in water, dissolved in oil and adsorbed on pore channels' walls. The surfactant's invasion into the reservoir is tracked by its diffusion with reservoir liquid and mass-exchange with two phase (liquid and solid) components of porous structure. Additionally, in this case heat exchange between fluids (injected, residual) and framework of porous medium has practical importance for evaluating of temperature influences on enhancing oil recovery. Now, the problem of designing an adequate mathematical model for describing a simultaneous flowing heat and mass transfer processes in anisotropic heterogeneous porous medium -surfactant injection during at various temperature regimes has not been fully researched. In this work is presents a 2D mathematical model of surfactant injections into the oil reservoir. Description of heat- and mass transfer processes in a porous media is done through differential and kinetic equations. For designing a computational algorithm is used modify version of IMPES method. The sequential and parallel computational algorithms are developed using an adaptive curvilinear meshes which into account heterogeneous porous structures. In this case we can evaluate the boundaries of our process flows - fronts (``invasion'', ``heat'' and ``mass'' transfers), according to the pressure, temperature, and concentration gradient changes.
A novel 2D porous indium coordination polymer with tunable luminescent property
NASA Astrophysics Data System (ADS)
Li, Xuejiao; Wang, Fangfang; Yang, He; Xu, Bo; Li, Cuncheng
2016-08-01
A new Indium coordination polymer [In(pda)1.5(phen)]n1 based on 1,4-phenylenediacetic acid (H2pda) and phen = 1,10-phenanthroline was obtained under hydrothermal condition and further characterized by single crystal X-ray analysis and other physicochemical studies such as infrared spectrum (IR), elemental analysis, thermogravimetric analysis (TGA) and powder X-ray diffraction (PXRD). Structure analysis reveals that complex 1 exhibits 2D porous (6,3) connected layer structure. Luminescent property of 1 was investigated both in the solid state and in different solvents and the results indicated that complex 1 demonstrates distinct solvent dependent luminescent property.
Electrochemical fabrication of 2D and 3D nickel nanowires using porous anodic alumina templates
NASA Astrophysics Data System (ADS)
Mebed, A. M.; Abd-Elnaiem, Alaa M.; Al-Hosiny, Najm M.
2016-06-01
Mechanically stable nickel (Ni) nanowires array and nanowires network were synthesized by pulse electrochemical deposition using 2D and 3D porous anodic alumina (PAA) templates. The structures and morphologies of as-prepared films were characterized by X-ray diffraction and scanning electron microscopy, respectively. The grown Ni nanowire using 3D PAA revealed more strength and larger surface area than has grown Ni use 2D PAA template. The prepared nanowires have a face-centered cubic crystal structure with average grain size 15 nm, and the preferred orientation of the nucleation of the nanowires is (111). The diameter of the nanowires is about 50-70 nm with length 3 µm. The resulting 3D Ni nanowire lattice, which provides enhanced mechanical stability and an increased surface area, benefits energy storage and many other applications which utilize the large surface area.
Viscoelastic fluid flow in inhomogeneous porous media
Siginer, D.A.; Bakhtiyarov, S.I.
1996-09-01
The flow of inelastic and viscoelastic fluids in two porous media of different permeabilities and same priority arranged in series has been investigated both theoretically and experimentally. The fluids are an oil field spacer fluid and aqueous solutions of polyacrylamide. The porous medium is represented by a cylindrical tube randomly packed with glass spheres. Expressions for the friction factor and the resistance coefficient as a function of the Reynolds number have been developed both for shear thinning and viscoelastic fluids based on the linear fluidity and eight constant Oldroyd models, respectively. The authors show that the energy loss is higher if the viscoelastic fluid flows first through the porous medium with the smaller permeability rather than through the section of the cylinder with the larger permeability. This effect is not observed for Newtonian and shear thinning fluids flowing through the same configuration. Energy requirements for the same volume flow rate are much higher than a Newtonian fluid of the same zero shear viscosity as the polymeric solution. Energy loss increases with increasing Reynolds number at a fixed concentration. At a fixed Reynolds number, the loss is a strong function of the concentration and increases with increasing concentration. The behavior of all fluids is predicted qualitatively except the difference in energy requirements.
Dynamics of clogging in drying porous media
NASA Astrophysics Data System (ADS)
Kaplan, C. Nadir; Mahadevan, L.
2014-11-01
Drying in porous media pervades a range of phenomena from brine evaporation arrested in porous bricks, causing efflorescence, i.e. salt aggregation on the surface where vapor leaves the medium, to clogging of reservoir rocks via salt precipitation when carbon dioxide is injected for geological storage. During the process of drying, the permeability and porosity of the medium may change due to the solute accumulation as a function of the particle concentration, in turn affecting the evaporation rate and the dynamics of the fluid flow imposed by it. To examine the dynamics of these coupled quantities, we develop a multiphase model of the particulate flow of a saline suspension in a porous medium, induced by evaporation. We further provide dimensional arguments as to how the salt concentration and the resulting change in permeability determine the transition between efflorescence and salt precipitation in the bulk. This research was supported by the Air Force Office of Scientific Research (AFOSR) under Award FA9550-09-1-0669-DOD35CAP and the Kavli Institute for Bionano Science and Technology at Harvard University.
Modern hardware architectures accelerate porous media flow computations
NASA Astrophysics Data System (ADS)
Kulczewski, Michal; Kurowski, Krzysztof; Kierzynka, Michal; Dohnalik, Marek; Kaczmarczyk, Jan; Borujeni, Ali Takbiri
2012-05-01
Investigation of rock properties, porosity and permeability particularly, which determines transport media characteristic, is crucial to reservoir engineering. Nowadays, micro-tomography (micro-CT) methods allow to obtain vast of petro-physical properties. The micro-CT method facilitates visualization of pores structures and acquisition of total porosity factor, determined by sticking together 2D slices of scanned rock and applying proper absorption cut-off point. Proper segmentation of pores representation in 3D is important to solve the permeability of porous media. This factor is recently determined by the means of Computational Fluid Dynamics (CFD), a popular method to analyze problems related to fluid flows, taking advantage of numerical methods and constantly growing computing powers. The recent advent of novel multi-, many-core and graphics processing unit (GPU) hardware architectures allows scientists to benefit even more from parallel processing and built-in new features. The high level of parallel scalability offers both, the time-to-solution decrease and greater accuracy - top factors in reservoir engineering. This paper aims to present research results related to fluid flow simulations, particularly solving the total porosity and permeability of porous media, taking advantage of modern hardware architectures. In our approach total porosity is calculated by the means of general-purpose computing on multiple GPUs. This application sticks together 2D slices of scanned rock and by the means of a marching tetrahedra algorithm, creates a 3D representation of pores and calculates the total porosity. Experimental results are compared with data obtained via other popular methods, including Nuclear Magnetic Resonance (NMR), helium porosity and nitrogen permeability tests. Then CFD simulations are performed on a large-scale high performance hardware architecture to solve the flow and permeability of porous media. In our experiments we used Lattice Boltzmann
MULTIGRID HOMOGENIZATION OF HETEROGENEOUS POROUS MEDIA
Dendy, J.E.; Moulton, J.D.
2000-10-01
This is the final report of a three-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL); this report, however, reports on only two years research, since this project was terminated at the end of two years in response to the reduction in funding for the LDRD Program at LANL. The numerical simulation of flow through heterogeneous porous media has become a vital tool in forecasting reservoir performance, analyzing groundwater supply and predicting the subsurface flow of contaminants. Consequently, the computational efficiency and accuracy of these simulations is paramount. However, the parameters of the underlying mathematical models (e.g., permeability, conductivity) typically exhibit severe variations over a range of significantly different length scales. Thus the numerical treatment of these problems relies on a homogenization or upscaling procedure to define an approximate coarse-scale problem that adequately captures the influence of the fine-scale structure, with a resultant compromise between the competing objectives of computational efficiency and numerical accuracy. For homogenization in models of flow through heterogeneous porous media, We have developed new, efficient, numerical, multilevel methods, that offer a significant improvement in the compromise between accuracy and efficiency. We recently combined this approach with the work of Dvorak to compute bounded estimates of the homogenized permeability for such flows and demonstrated the effectiveness of this new algorithm with numerical examples.
The significance of biofilms in porous media
Rittmann, B.E. )
1993-07-01
The recent literature contains conflicting claims about the characteristics of attached bacteria in subsurface porous media and how these characteristics affect permeability reduction. Some claim that the bacteria form continuous biofilms that restrict the pore size, while others claim that bacteria are attached in patchy aggregates that accumulate in pore throats. This contribution applies a recently developed tool from biofilm kinetics, the normalized surface loading, to interpret a wide range of experimental data from porous media experiments and biological filtration. The normalized surface loading is the actual substrate flux (i.e., rate of removal per unit surface area) divided by the minimum flux capable of supporting a deep biofilm. The analyses show that biofilms are continuous for normalized surface loadings greater than 1.0, but appear to become discontinuous for values less than about 0.25. For the low-load situation, distinguishing between continuous and discontinuous biofilms is not important when the modeling goal is prediction of substrate removal. However, the distinction is more critical when the modeling goal is to describe the spatial contribution of biomass and permeability loss. 62 refs., 1 fig., 1 tab.
Critical condensate saturation in porous media
Wang, X.; Mohanty, K.K.
1999-06-15
The understanding of gas and condensate flow in porous media is critical to the optimum exploitation of gas-condensate reservoirs. Critical condensate saturation and relative permeabilities are the key parameters for the evaluation of possible recovery strategies. This work is aimed at developing a mechanistic network model for the critical condensate saturation in which phase trapping and connectivity in the pore corners are critically examined. Porous media are modeled by networks of pore bodies interconnected by pore throats. Bodies and throats are characterized by their connectivity, shapes, and radii distributions. Pore-level laws are identified from micromodel experiments with near-critical fluids. A nonzero critical condensate saturation can be obtained in the absence of contact angle hysteresis due to the converging-diverging nature of the throats. The critical saturation at which the condensate flows is found to be a function of pore geometry, water saturation, and interfacial tension (or the Bond number). The modified sphere-pack model underpredicts the critical condensate saturation of typical sandstones. The cubic model adequately predicts the critical saturation and its experimentally observed trends.
Fundamentals of foam transport in porous media
Kovscek, A.R.; Radke, C.J.
1993-10-01
Foam in porous media is a fascinating fluid both because of its unique microstructure and because its dramatic influence on the flow of gas and liquid. A wealth of information is now compiled in the literature describing foam generation, destruction, and transport mechanisms. Yet there are conflicting views of these mechanisms and on the macroscopic results they produce. By critically reviewing how surfactant formulation and porous media topology conspire to control foam texture and flow resistance, we attempt to unify the disparate viewpoints. Evolution of texture during foam displacement is quantified by a population balance on bubble concentration, which is designed specifically for convenient incorporation into a standard reservoir simulator. Theories for the dominant bubble generation and coalescence mechanisms provide physically based rate expressions for the proposed population balance. Stone-type relative permeability functions along with the texture-sensitive and shear-thinning nature of confined foam complete the model. Quite good agreement is found between theory and new experiments for transient foam displacement in linear cores.
2D seismic reflection tomography in strongly anisotropic media
NASA Astrophysics Data System (ADS)
Huang, Guangnan; Zhou, Bing; Li, Hongxi; Zhang, Hua; Li, Zelin
2014-12-01
Seismic traveltime tomography is an effective method to reconstruct underground anisotropic parameters. Currently, most anisotropic tomographic methods were developed under the assumption of weak anisotropy. The tomographic method proposed here can be implemented for imaging subsurface targets in strongly anisotropic media with a known tilted symmetry axis, since the adopted ray tracing method is suitable for anisotropic media with arbitrary degree. There are three kinds of reflection waves (qP, qSV and qSH waves) that were separately used to invert the blocky abnormal body model. The reflection traveltime tomographiy is developed here because a surface observation system is the most economical and practical way compared with crosswell and VSP. The numerical examples show that the traveltimes of qP reflection wave have inverted parameters {{c}11},{{c}13},{{c}33} \\text{and} {{c}44} successfully. Traveltimes of qSV reflection wave have inverted parameters {{c}11},{{c}33} \\text{and} {{c}44} successfully, with the exception of the {{c}13}, since it is less sensitive than other parameters. Traveltimes of qSH reflection wave also have inverted parameters {{c}44} \\text{and} {{c}66} successfully. In addition, we find that the velocity sensitivity functions (derivatives of phase velocity with respect to elastic moduli parameters) and raypath illuminating angles have a great influence on the qualities of tomograms according to the inversion of theoretical models. Finally, the numerical examples confirm that the reflection traveltime tomography can be applied to invert strongly anisotropic models.
Proceedings of heat transfer and flow in porous media
Somerton, C.W.
1990-01-01
The topic of heat transfer and flow in porous media continues to be the focus of considerable research efforts. Certainly, this is partly due to the wide application of porous materials in engineering systems as well as the novel application of a porous media model to a variety of engineering problems. The work presented in this volume deals with such applications as papermaking, insulation materials, heat pipes, buried heating systems, tumor treatment, and cooling of microelectronics. This volume contains a nice mixture of experimental and computational approaches to problems and should provide the reader with a sense of the current state-of-the-art in porous media research.
Stability of Miscible Displacements Across Stratified Porous Media
Shariati, Maryam; Yortsos, Yanis C.
2000-09-11
This report studied macro-scale heterogeneity effects. Reflecting on their importance, current simulation practices of flow and displacement in porous media were invariably based on heterogeneous permeability fields. Here, it was focused on a specific aspect of such problems, namely the stability of miscible displacements in stratified porous media, where the displacement is perpendicular to the direction of stratification.
Simulations on porous media with nanofluids-initial study
NASA Astrophysics Data System (ADS)
Zeidan, D.; Alnaief, M.; Saghir, M. Ziad; Touma, R.
2016-06-01
Numerical simulations of natural convection heat and mass transfer in a square cavity using Comsol Multiphysics 5.0 software are presented. The effective thermal conductivity of nanofluid in porous media is computed using glass beads as porous media. It is observed that the heat and mass transfer rate increases with the increase of temperature variation as well as nanoparticle volume concentration.
Resistivity inversion in 2-D anisotropic media: numerical experiments
NASA Astrophysics Data System (ADS)
Wiese, Timothy; Greenhalgh, Stewart; Zhou, Bing; Greenhalgh, Mark; Marescot, Laurent
2015-04-01
Many rocks and layered/fractured sequences have a clearly expressed electrical anisotropy although it is rare in practice to incorporate anisotropy into resistivity inversion. In this contribution, we present a series of 2.5-D synthetic inversion experiments for various electrode configurations and 2-D anisotropic models. We examine and compare the image reconstructions obtained using the correct anisotropic inversion code with those obtained using the false but widely used isotropic assumption. Superior reconstruction in terms of reduced data misfit, true anomaly shape and position, and anisotropic background parameters were obtained when the correct anisotropic assumption was employed for medium to high coefficients of anisotropy. However, for low coefficient values the isotropic assumption produced better-quality results. When an erroneous isotropic inversion is performed on medium to high level anisotropic data, the images are dominated by patterns of banded artefacts and high data misfits. Various pole-pole, pole-dipole and dipole-dipole data sets were investigated and evaluated for the accuracy of the inversion result. The eigenvalue spectra of the pseudo-Hessian matrix and the formal resolution matrix were also computed to determine the information content and goodness of the results. We also present a data selection strategy based on high sensitivity measurements which drastically reduces the number of data to be inverted but still produces comparable results to that of the comprehensive data set. Inversion was carried out using transversely isotropic model parameters described in two different co-ordinate frames for the conductivity tensor, namely Cartesian versus natural or eigenframe. The Cartesian frame provided a more stable inversion product. This can be simply explained from inspection of the eigenspectra of the pseudo-Hessian matrix for the two model descriptions.
Review of enhanced vapor diffusion in porous media
Webb, S.W.; Ho, C.K.
1998-08-01
Vapor diffusion in porous media in the presence of its own liquid has often been treated similar to gas diffusion. The gas diffusion rate in porous media is much lower than in free space due to the presence of the porous medium and any liquid present. However, enhanced vapor diffusion has also been postulated such that the diffusion rate may approach free-space values. Existing data and models for enhanced vapor diffusion, including those in TOUGH2, are reviewed in this paper.
Utilization of Porous Media for Condensing Heat Exchangers
NASA Technical Reports Server (NTRS)
Tuan, George C.
2006-01-01
The use of porous media as a mean of separating liquid condensate from the air stream in condensing heat exchangers has been explored in the past inside small plant growth chambers and in the Apollo Command Module. Both applications used a cooled porous media made of sintered stainless steel to cool and separate condensation from the air stream. However, the main issues with the utilization of porous media in the past have been the deterioration of the porous media over long duration, such as clogging and changes in surface wetting characteristics. In addition, for long duration usage, biofilm growth from microorganisms on the porous medial would also be an issue. In developing Porous Media Condensing Heat Exchangers (PMCHX) for future space applications, different porous materials and microbial growth control methods will need to be explored. This paper explores the work performed at JSC and GRC to evaluate different porous materials and microbial control methods to support the development of a Porous Media Condensing Heat Exchanger. It outlines the basic principles for designing a PMCHX and issues that were encountered and ways to resolve those issues. The PMCHX has potential of mass, volume, and power savings over current CHX and water separator technology and would be beneficial for long duration space missions.
Mechanics of fluids in porous media
NASA Astrophysics Data System (ADS)
Bear, Jacob; Corapcioglu, M. Yavuz
Transport of quantities such as mass component of a phase and/or heat occurs in fields as diversified as petroleum reservoir engineering, groundwater hydraulics, soil mechanics, industrial filtration, water purification, wastewater treatment, soil drainage and irrigation, and geothermal energy production. In all these areas, scientists, engineers, and planners make use of mathematical models; these models describe the relevant transport processes that occur within controlled porous medium domains and enable forecasting of the future behavior of these domains in response to planned activities. The mathematical models, in turn, are based on the understanding of phenomena, often within the void space, and on theories that relate these phenomena to measurable quantities.Because of the pressing needs in areas of practical interest such as the development of groundwater energy storage and geothermal energy production, a vast amount of research in all these fields has contributed, especially in the last two decades, to our understanding and ability to describe transport phenomena in porous media. In recent years these research efforts have been significantly accelerated, attracting scientists from many disciplines. The practical needs of solving boundary value problems in heterogeneous domains, irregular boundaries, coupled phenomena and multiple dependent variables led to the development of a variety of powerful numerical techniques. The realization that fields are highly heterogeneous and that the degree of heterogeneity depends on the scale of the problem led to the introduction of stochastic concepts as an additional tool for the description of phenomena.
Reconstruction of three-dimensional porous media using a single thin section.
Tahmasebi, Pejman; Sahimi, Muhammad
2012-06-01
The purpose of any reconstruction method is to generate realizations of two- or multiphase disordered media that honor limited data for them, with the hope that the realizations provide accurate predictions for those properties of the media for which there are no data available, or their measurement is difficult. An important example of such stochastic systems is porous media for which the reconstruction technique must accurately represent their morphology--the connectivity and geometry--as well as their flow and transport properties. Many of the current reconstruction methods are based on low-order statistical descriptors that fail to provide accurate information on the properties of heterogeneous porous media. On the other hand, due to the availability of high resolution two-dimensional (2D) images of thin sections of a porous medium, and at the same time, the high cost, computational difficulties, and even unavailability of complete 3D images, the problem of reconstructing porous media from 2D thin sections remains an outstanding unsolved problem. We present a method based on multiple-point statistics in which a single 2D thin section of a porous medium, represented by a digitized image, is used to reconstruct the 3D porous medium to which the thin section belongs. The method utilizes a 1D raster path for inspecting the digitized image, and combines it with a cross-correlation function, a grid splitting technique for deciding the resolution of the computational grid used in the reconstruction, and the Shannon entropy as a measure of the heterogeneity of the porous sample, in order to reconstruct the 3D medium. It also utilizes an adaptive technique for identifying the locations and optimal number of hard (quantitative) data points that one can use in the reconstruction process. The method is tested on high resolution images for Berea sandstone and a carbonate rock sample, and the results are compared with the data. To make the comparison quantitative, two sets
Principles of single-phase flow through porous media
Liu, Shijie; Masliyah, J.H.
1996-12-31
Porous media are both permeable and dispersive to a traversing fluid. Flow of a single-phasefluid in porous media is not only of practical interest but also of fundamental significance in characterizing the porous media. In this chapter, the characteristics of porous media are introduced from both fundamental and application points of view. A continuum approach is used. The volume-averaged equations are used to describe the flow, where the momentum dispersion has been neglected. The relations between Darcy`s law-Brinkinans equation and the volume-averaged Navier-Stokes equation are described. The Forchheimer hypothesis, Ergun equation, and Liu-Afacan-Alashyah equation are briefly described in terms of coupling of the viscous and inertial effects oil the single-phaseflow in porous media. Discussions are provided on the concept and modeling of areal porosity, tortuosity, permeability, and shear factor. A curved passage model is discussed in terms of the shear factor and pressure-drop modeling for flow through porous media. Bounding wall effects are discussed through a simple approach. Examples of flow simulations in porous media (i.e., slightly compressible flow in oil reservoirs and incompressible flow infixed beds) are provided. 110 refs., 22 figs., 4 tabs.
Macroscopic properties of fractured porous media
NASA Astrophysics Data System (ADS)
Jean Francois, T.; Adler, P.; Bogdanov, I.; Mourzenko, V.
2006-12-01
There are two basic problems to be addressed. The first one is to solve precisely the partial differential equations which govern the phenomena which occur in these media and which are of interest in a large number of applications. The second one is to define a methodology in order to be able to estimate the macroscopic properties of real media by using quantities which are easily measurable on the field. Two major steps are needed for the numerical solution (1). First, an unstructured tetrahedral mesh of the fractures and of the porous matrix located in between is constructed; second, the governing partial equations are discretized and solved, in a finite volume formulation. A brief overview of the various problems which have been addressed so far, will be given: single and two phase flows, unsteady flows around a well, dispersion of an active and a passive solute, mechanical properties. This set of codes enabled us to cope with the second basic problem. Our approach is based on the systematic use of the excluded volume of fractures (which will be defined). The number of fractures per unit volume can be replaced by the number ρ ' of fractures per excluded volume. When numerical results such as the percolation threshold, the macroscopic permeability are plotted as functions of ρ ' they become independent of the shapes of the fractures which is a decisive simplification. Then, we show how ρ ' can be estimated from measurements performed along lines (2), planes, and galleries. It is interesting to notice that many stereological relations are actually independent of the size and shapes of fractures provided that they are convex (3); such a property adds a lot of generality to our methodology. Some tentative applications of the methodology are given and they show that the estimations are always in good agreement with the data. References (1) I.I. Bogdanov, V.V. Mourzenko, J.-F. Thovert, P.M. Adler, Effective permeability of fractured porous media in steady
Confinement properties of 2D porous molecular networks on metal surfaces
NASA Astrophysics Data System (ADS)
Müller, Kathrin; Enache, Mihaela; Stöhr, Meike
2016-04-01
Quantum effects that arise from confinement of electronic states have been extensively studied for the surface states of noble metals. Utilizing small artificial structures for confinement allows tailoring of the surface properties and offers unique opportunities for applications. So far, examples of surface state confinement include thin films, artificial nanoscale structures, vacancy and adatom islands, self-assembled 1D chains, vicinal surfaces, quantum dots and quantum corrals. In this review we summarize recent achievements in changing the electronic structure of surfaces by adsorption of nanoporous networks whose design principles are based on the concepts of supramolecular chemistry. Already in 1993, it was shown that quantum corrals made from Fe atoms on a Cu(1 1 1) surface using single atom manipulation with a scanning tunnelling microscope confine the Shockley surface state. However, since the atom manipulation technique for the construction of corral structures is a relatively time consuming process, the fabrication of periodic two-dimensional (2D) corral structures is practically impossible. On the other side, by using molecular self-assembly extended 2D porous structures can be achieved in a parallel process, i.e. all pores are formed at the same time. The molecular building blocks are usually held together by non-covalent interactions like hydrogen bonding, metal coordination or dipolar coupling. Due to the reversibility of the bond formation defect-free and long-range ordered networks can be achieved. However, recently also examples of porous networks formed by covalent coupling on the surface have been reported. By the choice of the molecular building blocks, the dimensions of the network (pore size and pore to pore distance) can be controlled. In this way, the confinement properties of the individual pores can be tuned. In addition, the effect of the confined state on the hosting properties of the pores will be discussed in this review article.
Confinement properties of 2D porous molecular networks on metal surfaces.
Müller, Kathrin; Enache, Mihaela; Stöhr, Meike
2016-04-20
Quantum effects that arise from confinement of electronic states have been extensively studied for the surface states of noble metals. Utilizing small artificial structures for confinement allows tailoring of the surface properties and offers unique opportunities for applications. So far, examples of surface state confinement include thin films, artificial nanoscale structures, vacancy and adatom islands, self-assembled 1D chains, vicinal surfaces, quantum dots and quantum corrals. In this review we summarize recent achievements in changing the electronic structure of surfaces by adsorption of nanoporous networks whose design principles are based on the concepts of supramolecular chemistry. Already in 1993, it was shown that quantum corrals made from Fe atoms on a Cu(1 1 1) surface using single atom manipulation with a scanning tunnelling microscope confine the Shockley surface state. However, since the atom manipulation technique for the construction of corral structures is a relatively time consuming process, the fabrication of periodic two-dimensional (2D) corral structures is practically impossible. On the other side, by using molecular self-assembly extended 2D porous structures can be achieved in a parallel process, i.e. all pores are formed at the same time. The molecular building blocks are usually held together by non-covalent interactions like hydrogen bonding, metal coordination or dipolar coupling. Due to the reversibility of the bond formation defect-free and long-range ordered networks can be achieved. However, recently also examples of porous networks formed by covalent coupling on the surface have been reported. By the choice of the molecular building blocks, the dimensions of the network (pore size and pore to pore distance) can be controlled. In this way, the confinement properties of the individual pores can be tuned. In addition, the effect of the confined state on the hosting properties of the pores will be discussed in this review article
Kan, Kan; Wang, Lei; Yu, Peng; Jiang, Baojiang; Shi, Keying; Fu, Honggang
2016-05-21
Two-dimensional (2D) quasi-ordered nitrogen-enriched porous carbon (QNPC) nanohybrids, with the characteristics of an ultrathin graphite nanosheet framework and thick quasi-ordered nitrogen-doped carbon cladding with a porous texture, have been synthesized via an in situ polymerization assembly method. In the synthesis, the expandable graphite (EG) is enlarged by an intermittent microwave method, and then aniline monomers are intercalated into the interlayers of the expanded EG with the assistance of a vacuum. Subsequently, the intercalated aniline monomers could assemble on the interlayer surface of the expanded EG, accompanied by the in situ polymerization from aniline monomers to polyaniline. Meanwhile, the expanded EG could be exfoliated to graphite nanosheets. By subsequent pyrolysis and activation processes, the QNPC nanohybrids could be prepared. As supercapacitor electrodes, a typical QNPC12-700 sample derived from the precursor containing an EG content of 12%, with a high level of nitrogen doping of 5.22 at%, offers a high specific capacitance of 305.7 F g(-1) (1 A g(-1)), excellent rate-capability and long-term stability. Notably, an extremely high energy density of 95.7 Wh kg(-1) at a power density of 449.7 W kg(-1) in an ionic liquid electrolyte can be achieved. The unique structural features and moderate heteroatom doping of the QNPC nanohybrids combines electrochemical double layer and faradaic capacitance contributions, which make these nanohybrids ideal candidates as electrode materials for high-performance energy storage devices. PMID:27122446
Fractal Analysis of Stress Sensitivity of Permeability in Porous Media
NASA Astrophysics Data System (ADS)
Tan, Xiao-Hua; Li, Xiao-Ping; Liu, Jian-Yi; Zhang, Lie-Hui; Cai, Jianchao
2015-12-01
A permeability model for porous media considering the stress sensitivity is derived based on mechanics of materials and the fractal characteristics of solid cluster size distribution. The permeability of porous media considering the stress sensitivity is related to solid cluster fractal dimension, solid cluster fractal tortuosity dimension, solid cluster minimum diameter and solid cluster maximum diameter, Young's modulus, Poisson's ratio, as well as power index. Every parameter has clear physical meaning without the use of empirical constants. The model predictions of permeability show good agreement with those obtained by the available experimental expression. The proposed model may be conducible to a better understanding of the mechanism for flow in elastic porous media.
Experimental tracking of the evolution of foam in porous media
Cohen, D; Patzek, T.W.; Radke, C.J.
1996-07-01
The authors discuss the experiments that have been done to track the effects of diffusion of gas in foams trapped in porous media. They describe several types of experiments and discuss the difficulties that prevent quantitative results from being obtained in most cases. However, the experiments do help them understand the physics and diffusion-driven coarsening of foams trapped in porous media. This understanding is necessary to simulate the behavior of these foams and predict the mobilization characteristics of foam in porous media. At the end of this paper, they compare the trends and predictions resulting from the experimental work to the predictions of the models which are presented elsewhere.
Molecular diffusion in porous media by PGSE ESR.
Talmon, Yael; Shtirberg, Lazar; Harneit, Wolfgang; Rogozhnikova, Olga Yu; Tormyshev, Victor; Blank, Aharon
2010-06-21
Diffusion in porous media is a general subject that involves many fields of research, such as chemistry (e.g. porous catalytic pallets), biology (e.g. porous cellular organelles), and materials science (e.g. porous polymer matrixes for controlled-release and gas-storage materials). Pulsed-gradient spin-echo nuclear magnetic resonance (PGSE NMR) is a powerful technique that is often employed to characterize complex diffusion patterns inside porous media. Typically it measures the motion of at least approximately 10(15) molecules occurring in the milliseconds-to-seconds time scale, which can be used to characterize diffusion in porous media with features of approximately 2-3 mum and above (in common aqueous environments). Electron Spin Resonance (ESR), which operates in the nanoseconds-to-microseconds time scale with much better spin sensitivity, can in principle be employed to measure complex diffusion patterns in porous media with much finer features (down to approximately 10 nm). However, up to now, severe technical constraints precluded the adaptation of PGSE ESR to porous media research. In this work we demonstrate for the first time the use of PGSE ESR in the characterization of molecular restricted diffusion in common liquid solutions embedded in a model system for porous media made of sub-micron glass spheres. A unique ESR resonator, efficient gradient coils and fast gradient current drivers enable these measurements. This work can be further extended in the future to many applications that involve dynamical processes occurring in porous media with features in the deep sub-micron range down to true nanometric length scales. PMID:20372729
Physical Properties of Fractured Porous Media
NASA Astrophysics Data System (ADS)
Mohammed, T. E.; Schmitt, D. R.
2015-12-01
The effect of fractures on the physical properties of porous media is of considerable interest to oil and gas exploration as well as enhanced geothermal systems and carbon capture and storage. This work represents an attempt to study the effect fractures have on multiple physical properties of rocks. An experimental technique to make simultaneous electric and ultrasonic measurements on cylindrical core plugs is developed. Aluminum end caps are mounted with ultrasonic transducers to transmit pules along the axis of the cylinder while non-polarizing electrodes are mounted on the sides of the core to make complex conductivity measurements perpendicular to the cylinder axis. Electrical measurements are made by applying a sinusoidal voltage across the measurement circuit that consist of a resister and the sample in series. The magnitude and phase of the signal across the sample is recorded relative to the input signal across a range of frequencies. Synthetic rock analogs are constructed using sintered glass beads with fractures imbedded in them. The fracture location, size and orientation are controlled and each fractured specimen has an unfractured counterpart. Porosity, Permeability, electrical conductivity and ultrasonic velocity measurements are conducted on each sample with the complex electrical conductivities recorded at frequencies from 10hz to 1 Mhz. These measurements allow us to examine the changes induced by these mesoscale fractures on the embedding porous medium. Of particular interest is the effect of fracture orientation on electrical conductivity of the rock. Seismic anisotropy caused by fractures is a well understood phenomenon with many rock physics models dedicated to its understanding. The effect of fractures on electrical conductivity is less well understood with electrical anisotropy scarcely investigated in the literature. None the less, using electrical conductivity to characterize fractures can add an extra constraint to characterization based
Foam Flows in Analog Porous Media
NASA Astrophysics Data System (ADS)
Meheust, Y.; Géraud, B.; Jones, S. A.; Cantat, I.; Dollet, B.
2015-12-01
Foams have been used for decades as displacing fluids for EOR and aquifer remediation, and more recently as carriers of chemical amendments for the remediation of the vadose zone. Apart from various interesting physico-chemical and biochemical properties, foams are better injection fluids due to their low sensitivity to gravity and their peculiar rheology: for foams with bubbles on the order of at least the typical pore size, viscous dissipation arises mostly from the contact zones between the soap films and the walls. In most experimental studies no local information of the foam structure is possible, and only global quantities such as the effective viscosity can be measured. We investigate foam flow through a two-dimensional porous medium consisting of circular obstacles positioned randomly in a horizontal transparent Hele-Shaw cell. The local foam structure is recorded in situ, which provides a measure of the spatial distribution of bubble velocities and sizes at regular time intervals. The flow exhibits a rich phenomenology including preferential flow paths and local flow intermittency/non-stationarity despite the imposed permanent global flow rate. Moreover, the medium selects the bubble size distribution through lamella division-triggered bubble fragmentation. Varying the mean bubble size of the injected foam, its water content, and mean velocity, we characterize those processes systematically and show that the distributions of bubble sizes and velocities are to some extent correlated. We furthermore measure the evolution, along the flow direction, of the distribution of bubble sizes, and measure the efficiency of bubble fragmentation as a function of the control parameters. The bubble fragmentation can be modeled numerically and to some extent analytically, based on statistical measures inferred from the experimental data. This study sheds new light on the local rheology of foams in porous media and opens the way towards quantitative characterization of the
Foam Transport in Porous Media - A Review
Zhang, Z. F.; Freedman, Vicky L.; Zhong, Lirong
2009-11-11
Amendment solutions with or without surfactants have been used to remove contaminants from soil. However, it has drawbacks such that the amendment solution often mobilizes the plume, and its movement is controlled by gravity and preferential flow paths. Foam is an emulsion-like, two-phase system in which gas cells are dispersed in a liquid and separated by thin liquid films called lamellae. Potential advantages of using foams in sub-surface remediation include providing better control on the volume of fluids injected, uniformity of contact, and the ability to contain the migration of contaminant laden liquids. It is expected that foam can serve as a carrier of amendments for vadose zone remediation, e.g., at the Hanford Site. As part of the U.S. Department of Energy’s EM-20 program, a numerical simulation capability will be added to the Subsurface Transport Over Multiple Phases (STOMP) flow simulator. The primary purpose of this document is to review the modeling approaches of foam transport in porous media. However, as an aid to understanding the simulation approaches, some experiments under unsaturated conditions and the processes of foam transport are also reviewed. Foam may be formed when the surfactant concentration is above the critical micelle concentration. There are two main types of foams – the ball foam (microfoam) and the polyhedral foam. The characteristics of bulk foam are described by the properties such as foam quality, texture, stability, density, surface tension, disjoining pressure, etc. Foam has been used to flush contaminants such as metals, organics, and nonaqueous phase liquids from unsaturated soil. Ball foam, or colloidal gas aphrons, reportedly have been used for soil flushing in contaminated site remediation and was found to be more efficient than surfactant solutions on the basis of weight of contaminant removed per gram of surfactant. Experiments also indicate that the polyhedral foam can be used to enhance soil remediation. The
Structure of drying fronts in three-dimensional porous media.
Shokri, Nima; Sahimi, Muhammad
2012-06-01
Evaporation in a three-dimensional (3D) porous medium, a sand column saturated by water, was studied using synchrotron x-ray tomography. Three-dimensional images of the medium with a resolution of 7 μm were obtained during the evaporation. The entire column was scanned seven times, resulting in nearly 10(4) 2D cross sections and illustrating the spatial distribution of air, liquid, and solid phases at the pore scale. The results were analyzed in order to gain new insights and better understanding of the characteristics of the drying front that was formed when the liquid-filled pores were invaded by air, as well as the structure of the liquid phase as it was dried. The analysis indicates that the liquid phase has a self-similar fractal structure, with its fractal dimension D(f) in all the cross sections being a function of the water content or saturation. In addition, D(f) for the 3D liquid structure, as well as its density correlation function, were computed using the 3D images. A crossover length scale ξ was identified that separates the fractal regime from the compact geometry. For length scales r>ξ, the density correlation function approaches asymptotically the water content of the porous medium. The drying front is shown to be rough and multi-affine, rather than self-affine. Its properties were also computed using the 3D images. The roughness characteristics agree with those for imbibition in porous media, but not with those of fracture surfaces and crack lines. PMID:23005211
Gas-phase diffusion in porous media: Comparison of models
Webb, S.W.
1998-09-01
Two models are commonly used to analyze gas-phase diffusion in porous media in the presence of advection, the Advective-Dispersive Model (ADM) and the Dusty-gas Model (DGM). The ADM, which is used in TOUGH2, is based on a simple linear addition of advection calculated by Darcy`s law and ordinary diffusion using Fick`s law with a porosity-tortuosity-gas saturation multiplier to account for the porous medium. Another approach for gas-phase transport in porous media is the Dusty-Gas Model. This model applies the kinetic theory of gases to the gaseous components and the porous media (or dust) to combine transport due to diffusion and advection that includes porous medium effects. The two approaches are compared in this paper.
Quantitative Permeability Prediction for Anisotropic Porous Media
NASA Astrophysics Data System (ADS)
Sheng, Q.; Thompson, K. E.
2012-12-01
Pore-scale modeling as a predictive tool has become an integral to both research and commercial simulation in recent years. Permeability is one of the most important of the many properties that can be simulated. Traditionally, permeability is determined using Darcy's law, based on the assumption that the pressure gradient is aligned with the principal flow direction. However, a wide variety of porous media exhibit anisotropic permeability due to particle orientation or laminated structure. In these types of materials, the direction of fluid flow is not aligned with the pressure gradient (except along the principal directions). Thus, it is desirable to predict the full permeability tensor for anisotropic materials using a first-principles pore-scale approach. In this work, we present a fast method to determine the full permeability tensor and the principal directions using a novel network modeling algorithm. We also test the ability of network modeling (which is an approximate method) to detect anisotropy in various structures. Both computational fluid dynamics (CFD) methods and network modeling have emerged as effective techniques to predict rock properties. CFD models are more rigorous but computationally expensive. Network modeling involves significant approximations but can be orders-of-magnitude more efficient computationally, which is important for both speed and the ability to model larger scales. This work uses network modeling, with simulations performed on two types of anisotropic materials: laminated packings (with layers of different sized particles) and oriented packings (containing particles with preferential orientation). Pore network models are created from the porous media data, and a novel method is used to determine the permeability tensor and principal flow direction using pore network modeling. The method is verified by comparing the calculated principal directions with the known anisotropy and also by comparing permeability with values from CFD
Spatially resolved D-T(2) correlation NMR of porous media.
Zhang, Yan; Blümich, Bernhard
2014-05-01
Within the past decade, 2D Laplace nuclear magnetic resonance (NMR) has been developed to analyze pore geometry and diffusion of fluids in porous media on the micrometer scale. Many objects like rocks and concrete are heterogeneous on the macroscopic scale, and an integral analysis of microscopic properties provides volume-averaged information. Magnetic resonance imaging (MRI) resolves this spatial average on the contrast scale set by the particular MRI technique. Desirable contrast parameters for studies of fluid transport in porous media derive from the pore-size distribution and the pore connectivity. These microscopic parameters are accessed by 1D and 2D Laplace NMR techniques. It is therefore desirable to combine MRI and 2D Laplace NMR to image functional information on fluid transport in porous media. Because 2D Laplace resolved MRI demands excessive measuring time, this study investigates the possibility to restrict the 2D Laplace analysis to the sum signals from low-resolution pixels, which correspond to pixels of similar amplitude in high-resolution images. In this exploratory study spatially resolved D-T2 correlation maps from glass beads and mortar are analyzed. Regions of similar contrast are first identified in high-resolution images to locate corresponding pixels in low-resolution images generated with D-T2 resolved MRI for subsequent pixel summation to improve the signal-to-noise ratio of contrast-specific D-T2 maps. This method is expected to contribute valuable information on correlated sample heterogeneity from the macroscopic and the microscopic scales in various types of porous materials including building materials and rock. PMID:24607821
Perspectives on Porous Media MR in Clinical MRI
NASA Astrophysics Data System (ADS)
Sigmund, E. E.
2011-03-01
Many goals and challenges of research in natural or synthetic porous media are mirrored in quantitative medical MRI. This review will describe examples where MR techniques used in porous media (particularly diffusion-weighted imaging (DWI)) are applied to physiological pathologies. Tissue microstructure is one area with great overlap with porous media science. Diffusion-weighting (esp. in neurological tissue) has motivated models with explicit physical dimensions, statistical parameters, empirical descriptors, or hybrids thereof. Another clinically relevant microscopic process is active flow. Renal (kidney) tissue possesses significant active vascular / tubular transport that manifests as "pseudodiffusion." Cancerous lesions involve anomalies in both structure and flow. The tools of magnetic resonance and their interpretation in porous media has had great impact on clinical MRI, and continued cross-fertilization of ideas can only enhance the progress of both fields.
Development of correlations to predict biopolymer mobility in porous media
Hejri, S.; Willhite, G.P.; Green, D.W. )
1991-02-01
This paper describes the flow and rheological behavior of biopolymer solutions in sandpacks over a wide range of permeability and frontal advance rates. Empirical correlations were developed to estimate polymer mobility in porous media. The correlations are based on porous medium properties, polymer concentration, and rheological parameters for the polymer derived from steady-shear measurements.
Characterization of an impinging jet into porous media
NASA Astrophysics Data System (ADS)
Wang, Cong; Alhani, Salwan; Gharib, Morteza
2015-11-01
In this work, characteristic behavior of a liquid jet into porous hydrophobic / hydrophilic particle media is investigated. In porous media, the capillary effect becomes significant, especially when the jet Reynolds Number is low. To analyze the cavity creation phenomena, the effect of jet's diameter, speed and acceleration as well as particles' size are carefully studied. Such knowledge of fluid behavior will provide guidance for medicine injection process. This work is supported by Caltech GALCIT STEM program.
Unsteady vortical structures in porous media flows
NASA Astrophysics Data System (ADS)
Finn, Justin; Apte, Sourabh; Wood, Brian
2011-11-01
The pore scale character of moderate Reynolds number, inertial flow through mono-disperse packed beds of spheres is examined using numerical experiments. Direct numerical simulations are performed for flow through (i) a periodic, 3 × 3 × 6 simple cubic arrangement at Rep = 529 , and (ii) a realistic randomly packed tube containing 326 spheres with dtube /dsp = 5 . 96 at Rep = 600 . At these Reynolds numbers, unsteady vortical regions are dominant features at the pore scale, and can have a profound effect on permeability and dispersion properties at the macro-scale. Despite similar Reynolds numbers and mean void fractions, the vortical structures observed in these two flows are remarkably different. The flow through the arranged packing is characterized by spatially and temporally periodic vortex-ring like structures, while the flow through the random packing contains many elongated helical vortices and a wider spectrum of space and time scales. The sensitive dependence of flow length and time scales and the local pore geometry is investigated using the DNS data. Funding: NSF project #0933857, Inertial Effects in Flow Through Porous Media.
Wave simulation in anisotropic, saturated porous media
Carcione, J.M.
1995-12-31
Porous media are anisotropic due to bedding, compaction and the presence of aligned microcracks and fractures. Here, I assume that the skeleton (and not the solid itself) is anisotropic. The rheological model also includes anisotropic tortuosity and permeability. The poroelastic equations are based on a transversely isotropic extension of Biot`s theory, and the problem is of plane strain type, i.e., two dimensional, and describes qP - qS propagation. In the high-frequency case, the (two) viscodynamic operators are approximated by Zener relaxation functions, that allow a close differential formulation of Biot`s equation of motion. The propagation is solved numerically, with a direct grid method and a time splitting integration algorithm, allowing the solution of the stiff part of the differential equations in closed analytical form. Snapshots in sandstone show that three waves propagate when the fluid is ideal (zero viscosity): the fast compressional and shear waves and the slow compressional wave. Anisotropic tortuosity has not a major influence on the faster modes, but significantly affects the slow wavefront. On the other hand, when the fluid is viscous, the slow wave becomes diffusive and appears as a static mode at the source location.
Survival of MEOR systems in porous media
Bryant, R.S.; Douglas, J.
1986-03-01
It was determined that many potential environmental hazards from microbial enhancement of oil recovery (MEOR) technology are extensions of EOR environmental problems. One unique concern of MEOR is that interactions between injected microorganisms and molasses, and those microorganisms already present in the reservoir might pose novel problems. It appeared necessary to obtain information regarding the indigenous microbial flora and the survival characteristics of a MEOR system after it encounters simulated reservoir microbial conditions. We found from laboratory studies that adventitious microbial species present in the molasses overgrew several of the indigenous microorganisms in Berea sandstone cores and also overgrew the injected MEOR microorganisms. From the results of these studies, we conclude that: (1) The introduction of nutrients into a petroleum reservoir could stimulate the growth of indigenous microorganisms. (2) Microorganisms present in injected non-sterile nutrients can overgrow both injected microbes and indigenous microorganisms. (3) Spore-forming bacteria cannot survive time periods in porous media of 4-20 weeks if there are other bacterial types present. 7 refs., 7 figs., 6 tabs.
Energy dissipation during sublimation from porous media
NASA Astrophysics Data System (ADS)
Keller, H. U.; Skorov, Yu. V.
2002-09-01
Several physical processes during the sublimation from and inside porous media have been investigated in detail in a series our papers (Skorov et al., 1999, Icarus 140, 173, Skorov et al., 2001, Icarus 153, 180, Davidsson and Skorov, 2002, Icarus 156, 223, Davidsson and Skorov, 2002, Icarus, in press) in order to analyse the gas production of cometary nuclei . New features are the absorption of the irradiation within the uppermost layers of the pores (rather than on the surface), taking into account the gas pressure of the coma, and temperature dependent condensation and sublimation coefficients. Detailed kinetic calculations revealed deviations from the canonical gasdynamic models. We will summarize the impact of these new calculations on the physics of sublimation from a cometary nucleus. The absorption of the irradiation below the surface leads to a decrease of sublimation flux near the subsolar point but to an increase near the evening terminator and nightside of a rotating nucleus. More absorbed energy is available to be transferred into the interior of the nucleus. This effect and consequences for the development of cometary nuclei will be discussed.
2D quasi-ordered nitrogen-enriched porous carbon nanohybrids for high energy density supercapacitors
NASA Astrophysics Data System (ADS)
Kan, Kan; Wang, Lei; Yu, Peng; Jiang, Baojiang; Shi, Keying; Fu, Honggang
2016-05-01
Two-dimensional (2D) quasi-ordered nitrogen-enriched porous carbon (QNPC) nanohybrids, with the characteristics of an ultrathin graphite nanosheet framework and thick quasi-ordered nitrogen-doped carbon cladding with a porous texture, have been synthesized via an in situ polymerization assembly method. In the synthesis, the expandable graphite (EG) is enlarged by an intermittent microwave method, and then aniline monomers are intercalated into the interlayers of the expanded EG with the assistance of a vacuum. Subsequently, the intercalated aniline monomers could assemble on the interlayer surface of the expanded EG, accompanied by the in situ polymerization from aniline monomers to polyaniline. Meanwhile, the expanded EG could be exfoliated to graphite nanosheets. By subsequent pyrolysis and activation processes, the QNPC nanohybrids could be prepared. As supercapacitor electrodes, a typical QNPC12-700 sample derived from the precursor containing an EG content of 12%, with a high level of nitrogen doping of 5.22 at%, offers a high specific capacitance of 305.7 F g-1 (1 A g-1), excellent rate-capability and long-term stability. Notably, an extremely high energy density of 95.7 Wh kg-1 at a power density of 449.7 W kg-1 in an ionic liquid electrolyte can be achieved. The unique structural features and moderate heteroatom doping of the QNPC nanohybrids combines electrochemical double layer and faradaic capacitance contributions, which make these nanohybrids ideal candidates as electrode materials for high-performance energy storage devices.Two-dimensional (2D) quasi-ordered nitrogen-enriched porous carbon (QNPC) nanohybrids, with the characteristics of an ultrathin graphite nanosheet framework and thick quasi-ordered nitrogen-doped carbon cladding with a porous texture, have been synthesized via an in situ polymerization assembly method. In the synthesis, the expandable graphite (EG) is enlarged by an intermittent microwave method, and then aniline monomers are
Effects of Cognitive Styles on 2D Drafting and Design Performance in Digital Media
ERIC Educational Resources Information Center
Pektas, Sule Tasli
2010-01-01
This paper investigates the interactions between design students' cognitive styles, as measured by Riding's Cognitive Styles Analysis, and performance in 2D drafting and design tasks in digital media. An empirical research revealed that Imager students outperformed Verbalisers in both drafting and creativity scores. Wholist-Analytic cognitive…
Enhancing nZVI mobility in porous media using humate
NASA Astrophysics Data System (ADS)
Schmid, Doris; Micic Batka, Vesna; Gondikas, Andreas; Velimirovic, Milica; von der Kammer, Frank; Hofmann, Thilo
2016-04-01
The limited transport of nanoscale zero-valent iron (nZVI) particles in porous media is a major drawback for its use in groundwater remediation. Among other factors, transport of nZVI particles might be negatively affected by mineralogical and physical heterogeneities of the aquifer matrix. Carbonate minerals and iron oxides, for instance, provide positively charged patches which would further increase particle attachment to the sand grains. This study does assess the potential of sodium humate, a salt of humic acids, to enhance the mobility of nZVI particles. Humate is a non-toxic, inexpensive material extracted from natural oxidized lignite and obtained in commercial grade, which makes it advantageous for field applications. Humate is expected to shield the positively charged patches of the sand grains and consequently enhance nZVI mobility in porous media. In this study the humate was injected into an aquifer prior to injection of the nZVI particles. The potential of humate for enhancing the mobility of nZVI particles was tested in an array of columns packed with heterogeneous natural porous media of different mineralogical composition and sediment texture. The results demonstrated that without pre-injection of humates only limited mobility of nZVI particles can be obtained in all tested porous media. After the pre-injection of low concentration of humate (10 mg/L) the mobility of nZVI particles (1 g/L) was enhanced in all tested porous media. The magnitude of this enhancement was depended on the properties of the porous media. The largest improvement of nZVI mobility was observed for homogeneous quartz. This material had also the highest porosity (~ 40%), good sorting, and therefore a higher permeability compared to the other porous media tested. It is assumed that the higher permeability of this porous medium allowed an optimal distribution of humate, resulting in an approximately 6-fold enhancement of nZVI mobility. In carbonate-rich porous medium with a
Capillary thermomechanics in serially porous media, with implications for randomly porous media
NASA Astrophysics Data System (ADS)
Miller, Robert D.
Visions of a single mobile substance present as two rival phases more or less cleanly segregated by capillarity between rival strata of a serially porous medium (a ``discontinuum'') imply explicit testable equations for externally measured capillary thermo-osmotic pressures and capillary thermo-osmosis, with implications for thermomechanical consolidation, dilation, and cracking. Underlying equations assume fluid phases governed by the laws of surface tension and viscous flow, moderated by an appropriate form of the Clapeyron equation. Derived phenomenological coefficients in macroscopic equations for steady coupled transports of mass and heat include only path-length-weighted fluid and heat conductances for rival domains and the heat of phase transformation. Expressions emphasize the phase-specific nature of Onsager's reciprocity principle and apply to serial media held within permeameters set up for measuring either ``isothermal'' or ``adiabatic'' mass transport or held within sealed containers intended for measurements of ``thermal conductivities.'' Results clarify unmet challenges facing modelers of similar processes and attributes in randomly porous media.
Using TRINET for simulating flow and transport in porous media
Najita, J.; Doughty, C.
1998-08-01
The finite element model TRINET calculates transient or steady-state fluid flow and solute transport on a lattice composed of one-dimensional finite elements (i.e., pipes) of porous medium. TRINET incorporates an adaptive gridding algorithm to minimize numerical dispersion for transport calculations. Although TRINET was originally developed to study fracture networks, the primary interest here is in applying TRINET more generally to simulate transport in porous media (or a fractured medium being treated as an effective continuum). This requires developing expressions to relate TRINET inputs to equivalent parameters used to describe flow and transport in homogeneous porous media. In this report, the authors briefly describe the basic TRINET formulation for flow and transport, present TRINET equivalences for porous medium parameters, and compare TRINET to analytical solutions using the proposed porous medium equivalents.
Homogenization of two fluid flow in porous media
Daly, K. R.; Roose, T.
2015-01-01
The macroscopic behaviour of air and water in porous media is often approximated using Richards' equation for the fluid saturation and pressure. This equation is parametrized by the hydraulic conductivity and water release curve. In this paper, we use homogenization to derive a general model for saturation and pressure in porous media based on an underlying periodic porous structure. Under an appropriate set of assumptions, i.e. constant gas pressure, this model is shown to reduce to the simpler form of Richards' equation. The starting point for this derivation is the Cahn–Hilliard phase field equation coupled with Stokes equations for fluid flow. This approach allows us, for the first time, to rigorously derive the water release curve and hydraulic conductivities through a series of cell problems. The method captures the hysteresis in the water release curve and ties the macroscopic properties of the porous media with the underlying geometrical and material properties.
Porous media heat transfer for injection molding
Beer, Neil Reginald
2016-05-31
The cooling of injection molded plastic is targeted. Coolant flows into a porous medium disposed within an injection molding component via a porous medium inlet. The porous medium is thermally coupled to a mold cavity configured to receive injected liquid plastic. The porous medium beneficially allows for an increased rate of heat transfer from the injected liquid plastic to the coolant and provides additional structural support over a hollow cooling well. When the temperature of the injected liquid plastic falls below a solidifying temperature threshold, the molded component is ejected and collected.
Transport of human adenoviruses in porous media
NASA Astrophysics Data System (ADS)
Kokkinos, Petros; Syngouna, Vasiliki I.; Tselepi, Maria A.; Bellou, Maria; Chrysikopoulos, Constantinos V.; Vantarakis, Apostolos
2015-04-01
Groundwater may be contaminated with infective human enteric viruses from various wastewater discharges, sanitary landfills, septic tanks, agricultural practices, and artificial groundwater recharge. Coliphages have been widely used as surrogates of enteric viruses, because they share many fundamental properties and features. Although a large number of studies focusing on various factors (i.e. pore water solution chemistry, fluid velocity, moisture content, temperature, and grain size) that affect biocolloid (bacteria, viruses) transport have been published over the past two decades, little attention has been given toward human adenoviruses (hAdVs). The main objective of this study was to evaluate the effect of pore water velocity on hAdV transport in water saturated laboratory-scale columns packed with glass beads. The effects of pore water velocity on virus transport and retention in porous media was examined at three pore water velocities (0.39, 0.75, and 1.22 cm/min). The results indicated that all estimated average mass recovery values for hAdV were lower than those of coliphages, which were previously reported in the literature by others for experiments conducted under similar experimental conditions. However, no obvious relationship between hAdV mass recovery and water velocity could be established from the experimental results. The collision efficiencies were quantified using the classical colloid filtration theory. Average collision efficiency, α, values decreased with decreasing flow rate, Q, and pore water velocity, U, but no significant effect of U on α was observed. Furthermore, the surface properties of viruses and glass beads were used to construct classical DLVO potential energy profiles. The results revealed that the experimental conditions of this study were unfavorable to deposition and that no aggregation between virus particles is expected to occur. A thorough understanding of the key processes governing virus transport is pivotal for public
Postbreakthrough behavior in flow through porous media
NASA Astrophysics Data System (ADS)
López, Eduardo; Buldyrev, Sergey V.; Dokholyan, Nikolay V.; Goldmakher, Leo; Havlin, Shlomo; King, Peter R.; Stanley, H. Eugene
2003-05-01
We numerically simulate the traveling time of a tracer in convective flow between two points (injection and extraction) separated by a distance r in a model of porous media, d=2 percolation. We calculate and analyze the traveling time probability density function for two values of the fraction of connecting bonds p: the homogeneous case p=1 and the inhomogeneous critical threshold case p=pc. We analyze both constant current and constant pressure conditions at p=pc. The homogeneous p=1 case serves as a comparison base for the more complicated p=pc situation. We find several regions in the probability density of the traveling times for the homogeneous case (p=1) and also for the critical case (p=pc) for both constant pressure and constant current conditions. For constant pressure, the first region IP corresponds to the short times before the flow breakthrough occurs, when the probability distribution is strictly zero. The second region IIP corresponds to numerous fast flow lines reaching the extraction point, with the probability distribution reaching its maximum. The third region IIIP corresponds to intermediate times and is characterized by a power-law decay. The fourth region IVP corresponds to very long traveling times, and is characterized by a different power-law decaying tail. The power-law characterizing region IVP is related to the multifractal properties of flow in percolation, and an expression for its dependence on the system size L is presented. The constant current behavior is different from the constant pressure behavior, and can be related analytically to the constant pressure case. We present theoretical arguments for the values of the exponents characterizing each region and crossover times. Our results are summarized in two scaling assumptions for the traveling time probability density; one for constant pressure and one for constant current. We also present the production curve associated with the probability of traveling times, which is of
Dispersion of solutes in porous media
NASA Astrophysics Data System (ADS)
Hunt, A. G.; Skinner, T. E.; Ewing, R. P.; Ghanbarian-Alavijeh, B.
2011-04-01
A recently introduced theory of solute transport in porous media is tested by comparison with experiment. The solute transport is predicted using an adaptation of the cluster statistics of percolation theory to critical path analysis together with knowledge of how the structure of such percolation clusters affects the time of transport across them. Only the effects of a single scale of medium heterogeneity are incorporated, and a minimal amount of information regarding the structure of the medium is required. This framework is used to find effectively the distributions of solute velocities and travel distances and thus generate arrival time distributions. The comparison with experiment focuses on the dispersivity (the ratio of the second to the first moment of the spatial solute distribution). The predictions of the theory in the absence of diffusion are verified by comparing with over 2200 experiments over length scales from a few microns to 100 km. At larger length scales (centimeters on up) about 95% of the data lie within our predicted bounds. At smaller length scales approximately 99.8% of the data lie where we predict. These comparisons are not trivial as the typical values of the dispersivity increase by ten orders of magnitude over ten orders of magnitude of length scale. Noteworthy is that the classical advection-dispersion (ADE) equation predicts that the dispersivity should be independent of length scale! This agreement with experiment requires rethinking of the relevance of diffusion and multi-scale heterogeneity and would also appear to signal the complete inappropriateness of using the classical ADE or any of its derivatives to model solute transport.
Liquid flow and distribution in unsaturated porous media
NASA Technical Reports Server (NTRS)
Alexander, J. Iwan
2004-01-01
Flow and transport in permeable or porous media and microchannels occurs in a variety of situations in micro- and reduced-gravity environments, many of them associated with environmental control and life support systems. While the role of gravity is limited, due to the typically small size scales associated permeable media, gravity, at the very least, affects the overall disposition of fluid in a macroscopic system. This presentation will discuss examples where the absence of gravity affects flow and phase distribution in selected examples of unsaturated flow and transport of heat and mass in porous media and microchannels that are pertinent to spacecraft systems.
Experimental Evidence of Helical Flow in Porous Media.
Ye, Yu; Chiogna, Gabriele; Cirpka, Olaf A; Grathwohl, Peter; Rolle, Massimo
2015-11-01
Helical flow leads to deformation of solute plumes and enhances transverse mixing in porous media. We present experiments in which macroscopic helical flow is created by arranging different materials to obtain an anisotropic macroscopic permeability tensor with spatially variable orientation. The resulting helical flow entails twisting streamlines which cause a significant increase in lateral mass exchange and thus a large enhancement of plume dilution (up to 235%) compared to transport in homogenous media. The setup may be used to effectively mix solutes in parallel streams similarly to static mixers, but in porous media. PMID:26588388
Impact of saturation on dispersion and mixing in porous media
NASA Astrophysics Data System (ADS)
Jimenez-Martinez, J.; De Anna, P.; Turuban, R.; Tabuteau, H.; Le Borgne, T.; Meheust, Y.
2014-12-01
In partially saturated porous media, the spatial distribution of wetting (e.g., water) and non-wetting (e.g., air) phases causes the liquid flow to be focused onto narrow and complex flow paths, leaving large volumes of wetting fluid trapped in between non-wetting phase clusters. The impact of the resulting highly heterogeneous wetting fluid velocity distributions on the dispersion and mixing of a solute in this wetting phase is critical for predicting reactive transport processes that take place in partially saturated porous media. We study the dependence of dispersion and mixing on the saturation degree using a 2D experimental setup consisting of cylindrical grains built using soft lithography. The joint injection of the two phases (wetting and non-wetting) provides a controlled homogeneous saturation in the medium. The simultaneous measurement of the flow velocity field, the spatial distribution of the wetting and non-wetting phases, and the tracer concentration field are used to investigate the relationship between the flow field complexity induced by desaturation and dispersion/mixing properties. We analyze the temporal behavior of the mean concentration gradient and the scalar dissipation rate, which quantify the temporal variation of the concentration variability and the potential for mixing-controlled chemical reactivity. The formation of preferential flowpaths in unsaturated flows is found to have an important impact on the mixing behavior. While the mean concentration gradient decays in time for saturated flow following the classical diffusive smoothing of concentration gradients, the creation of highly channelized finger structures in unsaturated flows induces persistently large concentration gradients which decay slowly in time. The highly resolved concentration field images show that this effect is due to i) a drastic increase of the surface available for creating concentration gradients across the finger boundaries, ii) the existence of dead-ends with
Chau, Jessica Furrer; Or, Dani
2006-11-01
The effect of drainage front morphology on gaseous diffusion through partially saturated porous media is analyzed using the lattice Boltzmann method (LBM). Flow regimes for immiscible displacement in porous media have been characterized as stable displacement, capillary fingering, and viscous fingering. The dominance of a flow regime is associated with the relative magnitudes of gravity, viscous, and capillary forces, quantifiable via the Bond number Bo, capillary number Ca, and their difference, Bo-Ca . Forced drainage from an initially saturated two-dimensional (2D) porous medium was simulated and the resulting flow patterns were analyzed and compared with theoretical predictions and experimental results. The LBM simulations reproduced expected flow morphologies for a range of drainage velocities and gravitational forces (i.e., a range of capillary and Bond numbers). Furthermore, measures of drainage front width as a function of the dimensionless difference Bo-Ca correspond well with scaling laws derived from percolation theory. Effects of flow morphology on residual fluid entrapment and gaseous diffusion were assessed by running LBM diffusion simulations through the partially saturated domain for a range of water contents. The effective diffusion coefficient as a function of water content was estimated for three regimes: stable drainage front, capillary fingering, and viscous fingering. Significant reductions in gaseous diffusion coefficient were found for viscous fingering relative to stable displacement, and to a lesser extent for capillary fingering, indicating that wetting phase distribution with a high degree of fingering in the 2D domain severely restricts connectivity of gas diffusion pathways through the medium. The study lends support for the use of LBM in design and management of fluids in porous media under variable gravity, and enhances the understanding of the role of dynamic fluid behavior on macroscopic transport properties of partially saturated
A pore scale study on turbulent combustion in porous media
NASA Astrophysics Data System (ADS)
Jouybari, N. F.; Maerefat, M.; Nimvari, M. E.
2016-02-01
This paper presents pore scale simulation of turbulent combustion of air/methane mixture in porous media to investigate the effects of multidimensionality and turbulence on the flame within the pores of porous media. In order to investigate combustion in the pores of porous medium, a simple but often used porous medium consisting of a staggered arrangement of square cylinders is considered in the present study. Results of turbulent kinetic energy, turbulent viscosity ratio, temperature, flame speed, convective heat transfer and thermal conductivity are presented and compared for laminar and turbulent simulations. It is shown that the turbulent kinetic energy increases from the inlet of burner, because of turbulence created by the solid matrix with a sudden jump or reduction at the flame front due to increase in temperature and velocity. Also, the pore scale simulation revealed that the laminarization of flow occurs after flame front in the combustion zone and turbulence effects are important mainly in the preheat zone. It is shown that turbulence enhances the diffusion processes in the preheat zone, but it is not enough to affect the maximum flame speed, temperature distribution and convective heat transfer in the porous burner. The dimensionless parameters associated with the Borghi-Peters diagram of turbulent combustion have been analyzed for the case of combustion in porous media and it is found that the combustion in the porous burner considered in the present study concerns the range of well stirred reactor very close to the laminar flame region.
Diffusion, Dispersion, and Uncertainty in Anisotropic Fractal Porous Media
NASA Astrophysics Data System (ADS)
Monnig, N. D.; Benson, D. A.
2007-12-01
Motivated by field measurements of aquifer hydraulic conductivity (K), recent techniques were developed to construct anisotropic fractal random fields, in which the scaling, or self-similarity parameter, varies with direction and is defined by a matrix. Ensemble numerical results are analyzed for solute transport through these 2-D "operator-scaling" fractional Brownian motion (fBm) ln(K) fields. Contrary to some analytic stochastic theories for monofractal K fields, the plume growth rates never exceed Mercado's (1967) purely stratified aquifer growth rate of plume apparent dispersivity proportional to mean distance. Apparent super-stratified growth must be the result of other demonstrable factors, such as initial plume size. The addition of large local dispersion and diffusion does not significantly change the effective longitudinal dispersivity of the plumes. In the presence of significant local dispersion or diffusion, the concentration coefficient of variation CV={σc}/{\\langle c \\rangle} remains large at the leading edge of the plumes. This indicates that even with considerable mixing due to dispersion or diffusion, there is still substantial uncertainty in the leading edge of a plume moving in fractal porous media.
Surface effects and pressure threshold of wet porous media
Lei, S.Y.; Jia, L.Q.; Wang, B.X.; Xia, C.M.
1997-07-01
The surface effects, a hysteresis switch effect, of the wet saturated porous media resulting from gas flow into the free surface were investigated experimentally. The experiments showed that there is a pressure threshold occurring in a very thin layer of wet porous media when gas is forced into the media. The gas enters and keeps flowing through the media only when the gas pressure reaches the threshold. When the pressure drops down somewhere, the gas flow will be stopped and the resistance to gas being forced into the media will recover. For screened sand naturally deposited in water whose particle size ranges from 0.10mm to 0.45mm and size rate is 1:1.25, the pressure threshold is p{sub th} = 4.94d{sup {minus}0.772} Pa. After defining threshold size of the porous media as d{sub th} = 4{sigma}/p{sub th}, the relative size can be described as d{sub th}/d = 0.0549d{sup {minus}0.237}, which is getting larger as grain getting finer. The variation of d{sub th}/d with d agrees well with that of the dimensionless permeability {radical}k/d for the screened sands. The threshold size of the porous media can unify all the experimental results of the permeability-porosity, including that of full screened soil sand.
SPH numerical simulation of fluid flow through a porous media
NASA Astrophysics Data System (ADS)
Klapp-Escribano, Jaime; Mayoral-Villa, Estela; Rodriguez-Meza, Mario Alberto; de La Cruz-Sanchez, Eduardo; di G Sigalotti, Leonardo; Inin-Abacus Collaboration; Ivic Collaboration
2013-11-01
We have tested an improved a method for 3D SPH simulations of fluid flow through a porous media using an implementation of this method with the Dual-Physics code. This improvement makes it possible to simulate many particles (of the order of several million) in reasonable computer times because its execution on GPUs processors makes it possible to reduce considerably the simulation cost for large systems. Modifications in the initial configuration have been implemented in order to simulate different arrays and geometries for the porous media. The basic tests were reproduced and the performance was analyzed. Our 3D simulations of fluid flow through a saturated homogeneous porous media shows a discharge velocity proportional to the hydraulic gradient reproducing Darcy's law at small body forces. The results are comparable with values obtained in previous work and published in the literature for simulations of flow through periodic porous media. Our simulations for a non saturated porous media produce adequate qualitative results showing that a non steady state is generated. The relaxation time for these systems were obtained. Work partially supported by Cinvestav-ABACUS, CONACyT grant EDOMEX-2011-C01-165873.
a Fractal Network Model for Fractured Porous Media
NASA Astrophysics Data System (ADS)
Xu, Peng; Li, Cuihong; Qiu, Shuxia; Sasmito, Agus Pulung
2016-04-01
The transport properties and mechanisms of fractured porous media are very important for oil and gas reservoir engineering, hydraulics, environmental science, chemical engineering, etc. In this paper, a fractal dual-porosity model is developed to estimate the equivalent hydraulic properties of fractured porous media, where a fractal tree-like network model is used to characterize the fracture system according to its fractal scaling laws and topological structures. The analytical expressions for the effective permeability of fracture system and fractured porous media, tortuosity, fracture density and fraction are derived. The proposed fractal model has been validated by comparisons with available experimental data and numerical simulation. It has been shown that fractal dimensions for fracture length and aperture have significant effect on the equivalent hydraulic properties of fractured porous media. The effective permeability of fracture system can be increased with the increase of fractal dimensions for fracture length and aperture, while it can be remarkably lowered by introducing tortuosity at large branching angle. Also, a scaling law between the fracture density and fractal dimension for fracture length has been found, where the scaling exponent depends on the fracture number. The present fractal dual-porosity model may shed light on the transport physics of fractured porous media and provide theoretical basis for oil and gas exploitation, underground water, nuclear waste disposal and geothermal energy extraction as well as chemical engineering, etc.
Modeling isothermal and non-isothermal flows in porous media
NASA Astrophysics Data System (ADS)
Mohseni Languri, Ehsan
2011-12-01
A complete understanding of the physics of flow and heat transfer phenomena in porous media is vital for accurate simulation of flow processes in industrial applications. In one such application pertaining to liquid composite molding (LCM) for manufacturing polymer composites, the fiber preforms used in LCM as reinforcements are limited not only to the single-scale porous media in the form of random fiber-mats, but also include dual-scale porous media in the form of woven or stitched fiber-mats. The conventional flow physics is not able to model the resin filling process in LCM involving the dual-scale porous media. In this study, the flow in dual-scale porous media is studied in order to predict the permeability of these fiber mats. The effect of aspect ratio of the fiber preform on the accuracy and flow during permeability estimation in single- and dual-scale porous media is analyzed experimentally and numerically. Flow of liquid in a free channel bounded on one side by porous medium is studied next, and two well-known boundary conditions of stress continuity and stress jump at the interface of the two regions are evaluated numerically. A point-wise solution for Stokes flow through periodic and non periodic porous media (made of cylindrical particles) adjacent to the free channel is presented using the Imite element based CFD software COMSOL. The efficacy of the two interfacial conditions is evaluated after volume averaging the point-wise velocity using a long averaging volume, also called the representative elementary volume or REV, and then comparing such a volume-averaged velocity profile with the available analytical solution. The investigation is carried out for five different porosities at three different Reynolds numbers to cover a wide range of applications. The presence of randomly-placed cylinders during the creation of non-periodic porous media damps out spatial fluctuations in the averaged velocity observed in periodic porous media. The analytical
Discrete Particle Model for Porous Media Flow using OpenFOAM at Intel Xeon Phi Coprocessors
NASA Astrophysics Data System (ADS)
Shang, Zhi; Nandakumar, Krishnaswamy; Liu, Honggao; Tyagi, Mayank; Lupo, James A.; Thompson, Karten
2015-11-01
The discrete particle model (DPM) in OpenFOAM was used to study the turbulent solid particle suspension flows through the porous media of a natural dual-permeability rock. The 2D and 3D pore geometries of the porous media were generated by sphere packing with the radius ratio of 3. The porosity is about 38% same as the natural dual-permeability rock. In the 2D case, the mesh cells reach 5 million with 1 million solid particles and in the 3D case, the mesh cells are above 10 million with 5 million solid particles. The solid particles are distributed by Gaussian distribution from 20 μm to 180 μm with expectation as 100 μm. Through the numerical simulations, not only was the HPC studied using Intel Xeon Phi Coprocessors but also the flow behaviors of large scale solid suspension flows in porous media were studied. The authors would like to thank the support by IPCC@LSU-Intel Parallel Computing Center (LSU # Y1SY1-1) and the HPC resources at Louisiana State University (http://www.hpc.lsu.edu).
Acoustically driven filtration of particulate suspensions in porous media
NASA Astrophysics Data System (ADS)
Gupta, Sanjay
1997-12-01
A novel method of filtration of liquid suspensions containing micron to millimeter size particles has been developed. A resonant ultrasonic field, applied across a highly porous medium, has been used to trap fine particles inside the large pores (relative to the particle size) of the medium. Three types of porous media, unconsolidated bed of 3 mm glass beads, consolidated open pore aluminum mesh, and reticulated polyester polyurethane foam were investigated as the test media. Reasonable filtration efficiencies were achieved for model aqueous suspensions of 325 mesh polystyrene particles in all three porous media. The expected trends of filtration performance with respect to suspension flow rate, its concentration, and the acoustic field intensity were confirmed. The Filtration phenomena was found to be limited by non-physical saturation of porous media. At saturation, the particles collected inside the media were found to exhibit macroscopic vibrations which allows them to escape with the carrier fluid. The highly porous POLY foam (95% porosity) was found to be the best media for suspension studied in terms of the duration of particle retention and percentage filtration efficiencies. The aluminum mesh performed slightly poorer. The unconsolidated porous media collected the least amount of solids. A simple theoretical development based on particle trajectory around an infinitely long cylindrical fiber, in the presence of acoustic field, has been initiated. In principle, the new filtration method is similar to high gravity magnetic separation but the acoustic method has a wider scope due to inherent acoustic contrast present in most suspensions. The low pressure drop, ease of operation, amenability to large scale operation and reasonable filtration efficiency make the new method highly attractive and suitable for practical applications.
Modeling microbial processes in porous media
NASA Astrophysics Data System (ADS)
Murphy, Ellyn M.; Ginn, Timothy R.
The incorporation of microbial processes into reactive transport models has generally proceeded along two separate lines of investigation: (1) transport of bacteria as inert colloids in porous media, and (2) the biodegradation of dissolved contaminants by a stationary phase of bacteria. Research over the last decade has indicated that these processes are closely linked. This linkage may occur when a change in metabolic activity alters the attachment/detachment rates of bacteria to surfaces, either promoting or retarding bacterial transport in a groundwater-contaminant plume. Changes in metabolic activity, in turn, are controlled by the time of exposure of the microbes to electron acceptors/donor and other components affecting activity. Similarly, metabolic activity can affect the reversibility of attachment, depending on the residence time of active microbes. Thus, improvements in quantitative analysis of active subsurface biota necessitate direct linkages between substrate availability, metabolic activity, growth, and attachment/detachment rates. This linkage requires both a detailed understanding of the biological processes and robust quantitative representations of these processes that can be tested experimentally. This paper presents an overview of current approaches used to represent physicochemical and biological processes in porous media, along with new conceptual approaches that link metabolic activity with partitioning of the microorganism between the aqueous and solid phases. Résumé L'introduction des processus microbiologiques dans des modèles de transport réactif a généralement suivi deux voies différentes de recherches: (1) le transport de bactéries sous forme de colloïdes inertes en milieu poreux, et (2) la biodégradation de polluants dissous par une phase stationnaire de bactéries. Les recherches conduites au cours des dix dernières années indiquent que ces processus sont intimement liés. Cette liaison peut intervenir lorsqu
Modeling microbial processes in porous media
NASA Astrophysics Data System (ADS)
Murphy, Ellyn M.; Ginn, Timothy R.
The incorporation of microbial processes into reactive transport models has generally proceeded along two separate lines of investigation: (1) transport of bacteria as inert colloids in porous media, and (2) the biodegradation of dissolved contaminants by a stationary phase of bacteria. Research over the last decade has indicated that these processes are closely linked. This linkage may occur when a change in metabolic activity alters the attachment/detachment rates of bacteria to surfaces, either promoting or retarding bacterial transport in a groundwater-contaminant plume. Changes in metabolic activity, in turn, are controlled by the time of exposure of the microbes to electron acceptors/donor and other components affecting activity. Similarly, metabolic activity can affect the reversibility of attachment, depending on the residence time of active microbes. Thus, improvements in quantitative analysis of active subsurface biota necessitate direct linkages between substrate availability, metabolic activity, growth, and attachment/detachment rates. This linkage requires both a detailed understanding of the biological processes and robust quantitative representations of these processes that can be tested experimentally. This paper presents an overview of current approaches used to represent physicochemical and biological processes in porous media, along with new conceptual approaches that link metabolic activity with partitioning of the microorganism between the aqueous and solid phases. Résumé L'introduction des processus microbiologiques dans des modèles de transport réactif a généralement suivi deux voies différentes de recherches: (1) le transport de bactéries sous forme de colloïdes inertes en milieu poreux, et (2) la biodégradation de polluants dissous par une phase stationnaire de bactéries. Les recherches conduites au cours des dix dernières années indiquent que ces processus sont intimement liés. Cette liaison peut intervenir lorsqu
Transport of Polycyclic Aromatic Hydrocarbons in Unsaturated Porous Media
NASA Astrophysics Data System (ADS)
Chahal, Maninder; Flury, Markus
2016-04-01
Polycyclic aromatic hydrocarbons (PAHs) are complex organic molecules containing 2 or more fused benzene rings. Being hydrophobic and non-polar, PAHs tend to partition to the organic matter in the soil from bulk aqueous phase. Though transport of these contaminants has been well studied in saturated environment, interactive mechanisms of these fluorescent compounds in unsaturated (identified by presence of air-water interface) porous media is still not well understood. We studied is the transport of fluoranthene in unsaturated porous media as facilitated by moving air-water interfaces. Confocal microscopy was used to visualize the interactions of fluoranthene particles in a glass channel packed with quartz glass beads. The packed glass channel was used to mimic a porous media and effects of an advancing and receding capillary fringe on the detachment of fluoranthene.
Permeability of mono- and bi-dispersed porous media
NASA Astrophysics Data System (ADS)
Byon, C.; Kim, S. J.
2013-04-01
In this study, the permeability of mono- and bi-dispersed porous media is considered. The effects of the particle size distribution and the packing structure of particles on the permeability are investigated experimentally and analytically. Both experimental and analytic results suggest that the particlesize distribution is close to the log-normal distribution, and the permeability of the mono-dispersed porous media quasi-linearly decreases as the range of the particle size distribution increases. On the other hand, the effect of packing structure of particles on the permeability is shown to be negligible.The permeability of the bidispersed porous media quasi-linearly decreases as the range of cluster size increases, and nearly independent of the particle size distribution. The present model is valid over the range of parameters typically found in heat transfer applications.
Fractal analysis of permeability near the wall in porous media
NASA Astrophysics Data System (ADS)
Liang, Mingchao; Yu, Boming; Li, Li; Yang, Shanshan; Zou, Mingqing
2014-01-01
In this paper, a fractal model for permeability of porous media is proposed based on Tamayol and Bahrami's method and the fractal theory for porous media. The proposed model is expressed as a function of the mean particle diameter, the length along the macroscopic pressure drop in the medium, porosity, fractal dimensions for pore space and tortuous capillaries, and the ratio of the minimum pore size to the maximum pore size. The relationship between the permeability near the wall and the dimensionless distance from the wall under different conditions is discussed in detail. The predictions by the present fractal model are in good agreement with available experimental data. The present results indicate that the present model may have the potential in comprehensively understanding the mechanisms of flow near the wall in porous media.
Experimental investigation of the permeability for unconsolidated porous media
Lei, S.Y.; Jia, L.Q.; Xia, C.M.; Zheng, G.Y.
1997-07-01
A device was constructed to investigate the permeability of unconsolidated media at low flow rate and small pressure drop. The stability and reliability of the device have been verified through repeated experiments on a given porous medium. The experimental investigation on the porous media demonstrated that the permeability-porosity relation is unique for a given medium. Experiments with the narrow screened sands show that conventional hydrodynamics theory and dimension analysis can not be applied satisfactorily in the study of the capillary porous media. For screened sand whose particle size ranges from 0.10mm to 0.45mm and size rate is 1:1.25, the permeability can be estimated from formula k = 4.89 x 10{sup {minus}4} d{sup 1.465} {phi}{sup 4.69} where k and d are limited in m{sup 2} and m, respectively.
Examining Asphaltene Solubility on Deposition in Model Porous Media.
Lin, Yu-Jiun; He, Peng; Tavakkoli, Mohammad; Mathew, Nevin Thunduvila; Fatt, Yap Yit; Chai, John C; Goharzadeh, Afshin; Vargas, Francisco M; Biswal, Sibani Lisa
2016-08-30
Asphaltenes are known to cause severe flow assurance problems in the near-wellbore region of oil reservoirs. Understanding the mechanism of asphaltene deposition in porous media is of great significance for the development of accurate numerical simulators and effective chemical remediation treatments. Here, we present a study of the dynamics of asphaltene deposition in porous media using microfluidic devices. A model oil containing 5 wt % dissolved asphaltenes was mixed with n-heptane, a known asphaltene precipitant, and flowed through a representative porous media microfluidic chip. Asphaltene deposition was recorded and analyzed as a function of solubility, which was directly correlated to particle size and Péclet number. In particular, pore-scale visualization and velocity profiles, as well as three stages of deposition, were identified and examined to determine the important convection-diffusion effects on deposition. PMID:27532331
A volume-balance model for flow on porous media
NASA Astrophysics Data System (ADS)
Malaga, Carlos; Mandujano, Francisco; Becerra, Julian
2015-11-01
Volume-balance models are used by petroleum engineers for simulating multiphase and multicomponent flow phenomena in porous media and the extraction process in oil reservoirs. In these models, mass conservation equations and Darcy's law are supplemented by a balance condition for the pore and fluid volumes. This provides a pressure equation suitable for simulating a compressible flow within a compressible solid matrix. Here we present an alternative interpretation of the volume-balance condition that includes the advective transport within a consolidated porous media. We obtain a modified equation for the time evolution of the pressure field. Numerical tests for phase separation under gravity are presented for multiphase three dimensional flow in heterogeneous porous media. The authors acknowledge funding from Fondo Sectorial CONACYT-SENER grant number 42536 (DGAJ-SPI-34-170412-217).
Linking Colloid Deposit Morphology and Clogging in Porous Media
NASA Astrophysics Data System (ADS)
Roth, E. J.; Mont-eton, M. E.; Mays, D. C.
2012-12-01
Innovations in the field of groundwater remediation have been hampered by delivery limitations in the porous media. For example, colloid deposits (comprising clays or silts) can cause a detrimental reduction in permeability, or clogging, which is problematic for groundwater remediation as well as granular media filtration and aquifer storage and recovery. During remediation, clogging creates preferential pathways in the media, leading to localized rather than spatially extensive contaminant treatment. Consequentially, remediation efforts become more expensive, less effective, and take a very long time. This presentation describes ongoing research investigating the link between colloid deposit morphology and clogging in porous media. As described by Darcy's Law, the velocity of fluid flow through porous media is proportional to permeability, which depends, in part, on porosity. However, changes in permeability are not in accord with changes in porosity as predicted by the Kozeny-Carman equation. It is hypothesized that unmeasured aspects of colloid deposit morphology could be the cause of this anomaly. Colloidal phenomena have important and dynamic effects on the permeability of natural porous media, and several lines of evidence suggest a correlation between clogging in porous media and the fractal dimension of colloid deposits. Here, a custom-built static light scattering apparatus is used to measure the fractal dimension of colloid deposits in refractive index matched porous media within a flow column. The media in our flow column is Nafion, which becomes essentially invisible when saturated by a solution of isopropanol and water. Polystyrene microspheres are then added to the influent through the column as a surrogate for natural colloids. Light from a laser is passed through the column, scattering from the deposited colloids, but not from the index matched Nafion. The resulting intensity of scattered light is measured as a function of scattering angle, and then
Permeability mapping in porous media by magnetization prepared centric-scan SPRITE
NASA Astrophysics Data System (ADS)
Romanenko, Konstantin V.; Balcom, Bruce J.
2011-02-01
The ability of porous media to transmit fluids is commonly referred to as permeability. The concept of permeability is central for hydrocarbon recovery from petroleum reservoirs and for studies of groundwater flow in aquifers. Spatially resolved measurements of permeability are of great significance for fluid dynamics studies. A convenient concept of local Darcy's law is suggested for parallel flow systems. The product of porosity and mean velocity images in the plane across the average flow direction is directly proportional to permeability. Single Point Ramped Imaging with T 1 Enhancement (SPRITE) permits reliable quantification of local fluid content and flow in porous media. It is particularly advantageous for reservoir rocks characterized by fast magnetic relaxation of a saturating fluid. Velocity encoding using the Cotts pulsed field gradient scheme improves the accuracy of measured flow parameters. The method is illustrated through measurements of 2D permeability maps in a capillary bundle, glass bead packs and composite sandstone samples.
Combustion and heat transfer in porous media
Sathe, S.B.; Peck, R.E.; Tong, T.W.
1990-06-01
The objective of the present study is to generate fundamental knowledge about heat transfer and combustion in porous radiant burners (PRBs) in order to improve their performance. A theoretical heat transfer and combustion model is developed to study the characteristics of PRBs. The model accounts for non-local thermal equilibrium between the solid and gas phases. The solid is assumed to absorb, emit and scatter radiant energy. Combustion is modeled as a one-step global reaction. It is revealed that the flame speed inside the porous medium is enhanced compared to the adiabatic flame speeds due to the higher conductivity of the solid compared to the gas as well as due to radiative preheating of the reactants. The effects of the properties of the porous material on the flame speeds, radiative outputs and efficiencies were investigated. To improve the radiative output from the burner, it is desirable that the porous layer has an optical thickness of about ten. The radiative output and the efficiency is higher for lower scattering albedo. The heat transfer coupling between the solid and gas phases should be high enough to ensure local thermal equilibrium, by choosing a fine porous matrix. Higher solid phase conduction enhances the flame speed and the radiative output. Experiments are performed on a ceramic foam to verify the theoretical findings. The existence of the two stability regions was verified experimentally.
Experimental investigation of contaminant transport in porous media. Research report
Wang, J.C.; Booker, J.R.; Carter, J.P.
1998-10-01
When numerical methods are applied to simulate a real contaminant transport problem, the values of a number of key parameters such as porosity, hydrodynamic dispersion coefficient or dispersivity and Darcy velocity or seepage velocity or seepage velocity are needed. In this paper, two different experimental programs, involving two types of column test and a well simulation test, were designed to demonstrate that the theory developed to explain contaminant transport in porous media is capable of representing the actual phenomenon of contaminant migration in soil. It is demonstrated that experiments can also be carried out to determine the properties necessary to model a real case of contaminant migration in porous media.
Percolation models for boiling and bubble growth in porous media
Yortsos, Y.C.
1991-05-01
We analyze the liquid-to-vapor phase change in single-component fluids in porous media at low superheats. Conditions typical to steam injection in porous media are taken. We examine nucleation, phase equilibria and their stability, and the growth of vapor bubbles. Effects of pore structure are emphasized. It is shown that at low supersaturations, bubble growth can be described as a percolation process. In the absence of spatial gradients, macroscopic flow properties are calculated in terms of nucleation parameters. A modification of gradient percolation is also proposed in the case of spatial temperature gradients, when solid conduction predominates. 22 refs., 10 figs., 1 tab.
An analytical model for permeability of isotropic porous media
NASA Astrophysics Data System (ADS)
Yang, Xiaohu; Lu, Tian Jian; Kim, Tongbeum
2014-06-01
We demonstrate that permeability of isotropic porous media e.g., open-cell foams can be analytically presented as a function of two morphological parameters: porosity and pore size. Adopting a cubic unit cell model, an existing tortuosity model from the branching algorithm method is incorporated into a generalized permeability model. The present model shows that dimensionless permeability significantly increases as the porosity of isotropic porous media and unifies the previously reported data in a wide range of porosity (ɛ=0.55-0.98) and pore size (Dp=0.254 mm-5.08 mm).
Heat Transfer through Porous Media in Static Superfluid Helium
Baudouy, B.; Juster, F.-P.; Allain, H.; Maekawa, R.
2006-04-27
Heat transfer through porous media in static saturated superfluid helium is investigated for porous media with different thickness, porosity and pore size. For large pore diameter, data are analyzed with the tortuosity concept in the pure Gorter-Mellink regime. It is shown that the tortuosity is constant over the temperature range investigated. For smaller pore diameter, the analysis reveals that the permeability is temperature-dependent in the Landau regime. In the intermediate regime, a model, including Landau and Gorter-Mellink regime, predicts a constant tortuosity within 10% but falls short predicting correctly the experimental data over the entire range of temperature.
Evaluation of bacterial detachment rates in porous media
Peyton, B.M.; Hooker, B.S.; Skeen, R.S.; Cunningham, A.B.; Lundman, R.W.
1994-05-01
The ability of published biomass detachment rate expressions to describe experimental data obtained from porous media reactors using Pseudomonas aeruginosa grown aerobically on glucose was evaluated. A first-order rate expression on attached biomass concentration best reflected effluent substrate concentration for combined data sets. Detachment rate coefficient k{sub d1} was dependent on initial substrate concentration. Simulation of porous media reactor experiments indicated that responses using higher influent substrate concentrations possessed greater sensitivity to variations in k{sub d1}. Simulations of field bioremediation systems suggest the use of accurate biofilm development kinetics is important in the prediction of well bore biofouling.
A hierarchically porous anatase TiO2 coated-WO3 2D IO bilayer film and its photochromic properties.
Li, Hua; Wu, Huazhong; Xiao, Jiajia; Su, Yanli; Robichaud, Jacques; Brüning, Ralf; Djaoued, Yahia
2016-01-18
A hierarchically porous anatase TiO2 coated-WO3 2D inverse opal (IO) bilayer film was fabricated on ITO glass using a layer by layer route with a hierarchically porous TiO2 top layer and an ordered super-macroporous WO3 2D IO bottom layer. This novel TiO2 coated-WO3 2D IO bilayer film was evaluated for photochromic applications. PMID:26576930
Bounds on Transport Coefficients of Porous Media
Berryman, J G
2005-03-21
An analytical formulation of conductivity bounds by Bergman and Milton is used in a different way to obtain rigorous bounds on the real transport coefficients (electrical conductivity, thermal conductivity, and/or fluid permeability) of a fluid-saturated porous medium. These bounds do not depend explicitly on the porosity, but rather on two formation factors--one associated with the pore space and the other with the solid frame. Hashin-Shtrikman bounds for transport in random polycrystals of porous-material laminates will also be discussed.
Exit flows from highly porous media
NASA Astrophysics Data System (ADS)
Hall, M. J.; Hiatt, J. P.
1994-02-01
This paper presents laser velocimetry measurements of the streamwise component of mean velocities and turbulence intensities measured downstream from the exit plane of porous ceramic foams through which air is flowed. The recent development and commercial availability of porous ceramic foams has lead to their application in many fields. Their uses have extended to combustion, high-temperature fluid filtering, biotechnology, and as support matrix for catalysts. These applications have created an interest in their pore scale fluid mechanics, both within the porous matrix and along surfaces open to flow. One emerging application is porous ceramic burners which combust liquid or gaseous fuels within the pore matrix or along the surface of the ceramic. The ceramic foams have pore sizes ranging from 4 to 12 pores per cm (ppcm) and porosities of 85%. Mean velocities between 0.3 and 1.5 m/s were examined. Radial distributions of mean velocities show a jet-like structure through the pores, with local mean velocities reaching maximum values over two times the area mean velocity. Negative mean velocities were often observed between pores, suggesting that recirculation zones are present above the web-like struts surrounding the pores. Levels of turbulence intensities normalized by the area mean velocity ranged from 0.05 to 0.6 for the various flow rates and pore sizes. Turbulence intensities were found to increase with increasing pore size for a given flow rate.
Pressure diffusion waves in porous media
Silin, Dmitry; Korneev, Valeri; Goloshubin, Gennady
2003-04-08
Pressure diffusion wave in porous rocks are under consideration. The pressure diffusion mechanism can provide an explanation of the high attenuation of low-frequency signals in fluid-saturated rocks. Both single and dual porosity models are considered. In either case, the attenuation coefficient is a function of the frequency.
NASA Astrophysics Data System (ADS)
Rupert, Y. K.
2015-12-01
The remediation and monitoring of soils and groundwater contaminated with organic compounds is an important goal of many environmental restoration efforts. This laboratory research focuses on combining two innovative geophysical methods: nuclear magnetic resonance (NMR) and spectral induced polarization (SIP) to assess their suitability to characterize and quantify organic contaminants in porous media. Toluene, a light non-aqueous phase liquid (LNAPL), and ethoxy-nonafluorobutane, an engineered dense non-aqueous phase liquid (DNAPL), have been selected as representative organic contaminants. Low-field NMR relaxation time (T2) measurements and diffusion-relaxation (D-T2) correlation measurements, as well as low frequency SIP measurements (<10 kHz) are performed to quantify the amount of these two organic compounds in the presence of water in three types of porous media (sands, clay, and various sand-clay mixtures). The T2, D-T2, and SIP measurements are made on water, toluene, and the synthetic DNAPL in each porous media to understand the effect of different porous media on the NMR and SIP responses in each fluid. We then plan to make measurements on water-organic mixtures with varied concentrations of organic compounds in each porous medium to resolve the NMR and SIP response of the organic contaminants from that of water and to quantify the amount of organic contaminants. Building a relationship between SIP and NMR signatures from organic contaminants not only provides a fundamental yet important petrophysical relationship, but also builds a framework for continued investigation into how these two methods synergize. This will also provide spatially dense information about organic contaminated natural sediments at scales that will improve the quantitative characterization and remediation of contaminated sites.The remediation and monitoring of soils and groundwater contaminated with organic compounds is an important goal of many environmental restoration efforts
Effects of microstructure on flow properties of fibrous porous media at moderate Reynolds number
NASA Astrophysics Data System (ADS)
Tamayol, A.; Wong, K. W.; Bahrami, M.
2012-02-01
In this study, effects of microstructure on the viscous permeability and Forchheimer coefficient of monodispersed fibers are investigated. The porous material is represented by a unit cell which is assumed to be repeated throughout the medium. Based on the orientation of the fibers in the space, fibrous media are divided into three categories: one-, two-, and three-directional (1D, 2D, and 3D) structures. Parallel and transverse flow through square arrangements of 1D fibers, simple 2D mats, and 3D simple cubic structures are solved numerically over a wide range of porosity (0.35 < ɛ < 0.95) and Reynolds number (0.01 < Re < 200). The results are used to calculate the permeability and the inertial coefficient of the considered geometries. An experimental study is performed; the flow coefficients of three different ordered tube banks in the moderate range of Reynolds number (0.001 < Re < 15) are determined. The numerical results are successfully compared with the present and the existing experimental data in the literature. The results suggest that the permeability and Forchheimer coefficient are functions of porosity and fiber orientation. A comparison of the experimental and numerical results with the Ergun equation reveals that this equation is not suitable for highly porous materials. As such, accurate correlations are proposed for determining the Forchheimer coefficient in fibrous porous media.
Rupture dynamics and ground motions from earthquakes in 2-D heterogeneous media
NASA Astrophysics Data System (ADS)
Bydlon, Samuel A.; Dunham, Eric M.
2015-03-01
We perform 2-D simulations of earthquakes on rough faults in media with random heterogeneities (with von Karman distribution) to study the effects of geometric and material heterogeneity on the rupture process and resulting high-frequency ground motions in the near-fault region (out to ˜20 km). Variations in slip and rupture velocity can arise from material heterogeneity alone but are dominantly controlled by fault roughness. Scattering effects become appreciable beyond ˜3 km from the fault. Near-fault scattering extends the duration of incoherent, high-frequency ground motions and, at least in our 2-D simulations, elevates root-mean-square accelerations (i.e., Arias intensity) with negligible reduction in peak velocities. We also demonstrate that near-fault scattering typically occurs in the power law tail of the power spectral density function, quantified by the Hurst exponent and another parameter combining standard deviation and correlation length.
NASA Astrophysics Data System (ADS)
Jiménez-Martínez, Joaquín.; Anna, Pietro de; Tabuteau, Hervé; Turuban, Régis; Borgne, Tanguy Le; Méheust, Yves
2015-07-01
Porous media in which different fluid phases coexist are common in nature (e.g., vadose zone and gas-oil reservoirs). In partially saturated porous media, the intricate spatial distributions of the wetting and nonwetting phases causes their flow to be focused onto preferential paths. Using a novel 2-D experimental setup allowing pore-scale measurement of concentration fields in a controlled unsaturated flow, we highlight mechanisms by which mixing of an invading fluid with the resident fluid is significantly enhanced when decreasing saturation. The mean scalar dissipation rate is observed to decrease slowly in time, while under saturated conditions it decays rapidly. This slow decrease is due to sustained longitudinal solute fingering, which causes concentration gradients to remain predominantly transverse to the average flow. Consequently, the effective reactivity is found to be much larger than under saturated conditions. These results provide new insights into the role that multiphase flows play on mixing/reaction in porous media.
Towards aeroacoustic sound generation by flow through porous media.
Hasert, Manuel; Bernsdorf, Joerg; Roller, Sabine
2011-06-28
In this work, we present single-step aeroacoustic calculations using the Lattice Boltzmann method (LBM). Our application case consists of the prediction of an acoustic field radiating from the outlet of a porous media silencer. It has been proved that the LBM is able to simulate acoustic wave generation and propagation. Our particular aim is to validate the LBM for aeroacoustics in porous media. As a validation case, we consider a spinning vortex pair emitting sound waves as the vortices rotate around a common centre. Non-reflective boundary conditions based on characteristics have been adopted from Navier-Stokes methods and are validated using the time evolution of a Gaussian pulse. We show preliminary results of the flow through the porous medium. PMID:21576161
Overlimiting current and deionization shocks in porous media
NASA Astrophysics Data System (ADS)
Bazant, Martin Z.; Deng, Daosheng S.; Mani, Ali; Dydek, E. Victoria
2011-11-01
Salt transport in bulk electrolytes occurs by diffusion and convection, but in microfluidic devices and porous media, surface conduction and electro-osmotic flow also contribute to ionic fluxes. The classical theory of electrokinetic phenomena in porous media assumes linear response to a small voltage, where the electrolyte concentration is only weakly perturbed. When a large voltage or concentration gradient is imposed, some surprising nonlinear electrokinetic phenomena result from the competition between bulk and interfacial transport in a microstructure. At constant voltage, the microstructure can sustain an over-limiting current (exceeding diffusion limitation) without any hydrodynamic or chemical instability. At constant current, a ``deionization shock'' can propagate through the microstructure, leaving behind a macroscopic region depleted of ions and particles. This talk will present experimental evidence for surface-driven overlimiting current and deioniziation shocks in porous glass frits, interpreted with the help of mathematical models, and applications to water deionization by ``shock electrodialysis.'' Supported by the MIT Energy Initiative.
Moisture Content and Migration Dynamics in Unsaturated Porous Media
NASA Technical Reports Server (NTRS)
Homsy, G. M.
1993-01-01
Fundamental studies of fluid mechanics and transport in partially saturated soils are presented. Solution of transient diffusion problems in support of the development of probes for the in-situ measurement of moisture content is given. Numerical and analytical methods are used to study the fundamental problem of meniscus and saturation front propagation in geometric models of porous media.
Acoustic Wave Monitoring of Biofilm Development in Porous Media
Biofilm development in porous media can result in significant changes to the hydrogeological properties of subsurface systems with implications for fluid flow and contaminant transport. As such, a number of numerical models and simulations have been developed in an attempt to qua...
Fractal and Multifractal Models Applied to Porous Media - Editorial
Technology Transfer Automated Retrieval System (TEKTRAN)
Given the current high level of interest in the use of fractal geometry to characterize natural porous media, a special issue of the Vadose Zone Journal was organized in order to expose established fractal analysis techniques and cutting-edge new developments to a wider Earth science audience. The ...
Complete dissolution of trichloroethylene in saturated porous media
Imhoff, P.T.; Arthur, M.H.; Miller, C.T.
1998-08-15
Porous media containing trichloroethylene (TCE) trapped at residual saturation in otherwise water-saturated porous media were flushed with water to asses the dissolution rate of TCE as TCE volumetric fractions approached zero. Careful attention to column design and experimental methods limited the effect of column materials on effluent concentrations. Effluent concentration measurements during TCE dissolution are presented for a glass bead porous medium, a mixed sand, and a treated soil. Effluent concentrations were measured as they decreased below 5 {micro}g/L, the maximum allowable contaminant level, in the glass bead and mixed sand media. Effluent concentrations from columns packed with treated soil were measured down to 20 {micro}g/L. Solvent extraction of the treated soil after the dissolution experiments revealed that extremely small quantities of TCE were retained in this medium. Results from parallel experiments on columns exposed to only aqueous TCE suggest that TCE remaining in the treated soil columns was sorbed to the porous medium. Existing power-law models were capable of describing TCE dissolution in these media, if the exponent on the TCE volume fraction was modified appropriately.
Colloid adhesive parameters for chemical heterogeneous porous media
Technology Transfer Automated Retrieval System (TEKTRAN)
A simple modeling approach was developed to calculate colloid adhesive parameters for chemically heterogeneous porous media. The area of the zone of electrostatic influence between a colloid and solid-water interface (Az) was discretized into a number of equally sized grid cells to capture chemical...
Microscopic interfacial phenomena during flow in porous media
Miksis, M.J.; Ida, M.P.
1996-12-31
A fundamental process during any multiphase flow in porous media is the breaking apart of one of the phases into smaller components. Here the authors investigate this breaking process as applied to a thin liquid film. They study the breaking of both a two dimensional planar film and a cylindrical thread of liquid using both analytical and numerical methods.
Autocatalytic Reaction and Flow in Porous Media
NASA Astrophysics Data System (ADS)
Atis, Severine; Auradou, Harold; Talon, Laurent; Salin, Dominique
2011-11-01
Universités Pierre et Marie Curie, Paris Sud and CNRS. Laboratoire FAST, Bâtiment 502, UPS, 91405 ORSAY Cedex France. Coupling between autocatalytic reaction front and simple hydrodynamic flows leads to front patterns revealing the underlying flow field. Flow of a passive tracer, i. e. a dye, through the complex flow field of a porous medium leads to the so-called hydrodynamic dispersion which accounts for the mixing process inside the medium. We have performed experiments and numerical simulations of the propagation of reaction front in a porous medium. We have analyzed the dependences of the shape and velocity of the stationary fronts with the flow rate for a flow either in the same direction than the chemical front propagation or opposite to it. We determine the structure and characteristics of the front as well as the velocity distribution measured along the front. As a result this active, chemical, tracer allows to access to the characteristic of the complex flow field of the porous medium.
Porous media for catalytic renewable energy conversion
NASA Astrophysics Data System (ADS)
Hotz, Nico
2012-05-01
A novel flow-based method is presented to place catalytic nanoparticles into a reactor by sol-gelation of a porous ceramic consisting of copper-based nanoparticles, silica sand, ceramic binder, and a gelation agent. This method allows for the placement of a liquid precursor containing the catalyst into the final reactor geometry without the need of impregnating or coating of a substrate with the catalytic material. The so generated foam-like porous ceramic shows properties highly appropriate for use as catalytic reactor material, e.g., reasonable pressure drop due to its porosity, high thermal and catalytic stability, and excellent catalytic behavior. The catalytic activity of micro-reactors containing this foam-like ceramic is tested in terms of their ability to convert alcoholic biofuel (e.g. methanol) to a hydrogen-rich gas mixture with low concentrations of carbon monoxide (up to 75% hydrogen content and less than 0.2% CO, for the case of methanol). This gas mixture is subsequently used in a low-temperature fuel cell, converting the hydrogen directly to electricity. A low concentration of CO is crucial to avoid poisoning of the fuel cell catalyst. Since conventional Polymer Electrolyte Membrane (PEM) fuel cells require CO concentrations far below 100 ppm and since most methods to reduce the mole fraction of CO (such as Preferential Oxidation or PROX) have CO conversions of up to 99%, the alcohol fuel reformer has to achieve initial CO mole fractions significantly below 1%. The catalyst and the porous ceramic reactor of the present study can successfully fulfill this requirement.
Diffusion Driven Combustion Waves in Porous Media
NASA Technical Reports Server (NTRS)
Aldushin, A. P.; Matkowsky, B. J.
2000-01-01
Filtration of gas containing oxidizer, to the reaction zone in a porous medium, due, e.g., to a buoyancy force or to an external pressure gradient, leads to the propagation of Filtration combustion (FC) waves. The exothermic reaction occurs between the fuel component of the solid matrix and the oxidizer. In this paper, we analyze the ability of a reaction wave to propagate in a porous medium without the aid of filtration. We find that one possible mechanism of propagation is that the wave is driven by diffusion of oxidizer from the environment. The solution of the combustion problem describing diffusion driven waves is similar to the solution of the Stefan problem describing the propagation of phase transition waves, in that the temperature on the interface between the burned and unburned regions is constant, the combustion wave is described by a similarity solution which is a function of the similarity variable x/square root of(t) and the wave velocity decays as 1/square root of(t). The difference between the two problems is that in the combustion problem the temperature is not prescribed, but rather, is determined as part of the solution. We will show that the length of samples in which such self-sustained combustion waves can occur, must exceed a critical value which strongly depends on the combustion temperature T(sub b). Smaller values of T(sub b) require longer sample lengths for diffusion driven combustion waves to exist. Because of their relatively small velocity, diffusion driven waves are considered to be relevant for the case of low heat losses, which occur for large diameter samples or in microgravity conditions, Another possible mechanism of porous medium combustion describes waves which propagate by consuming the oxidizer initially stored in the pores of the sample. This occurs for abnormally high pressure and gas density. In this case, uniformly propagating planar waves, which are kinetically controlled, can propagate, Diffusion of oxidizer decreases
Ginn, T.R.; Cushman, J.H.; Murphy, E.M.; Fletcher, M.
1998-06-01
'Dynamic microbial attachment/detachment occurs in subsurface systems in response to changing environmental conditions caused by contaminant movement and degradation. Understanding the environmental conditions and mechanisms by which anaerobic bacteria partition between aqueous and solid phases is a critical requirement for designing and evaluating in-situ bioremediation efforts. This interdisciplinary research project will provide fundamental information on the attachment/detachment dynamics of anaerobic bacteria in heterogeneous porous media under growth and growth-limiting conditions. Experiments will provide information on passive and active attachment/detachment mechanisms used by growing anaerobes capable of reductive dechlorination. Theoretical representations of these attachment/detachment mechanisms will be incorporated into existing flow and transport models that incorporate heterogeneity effects and can be used to predict behavior at field scales. These mechanistic-based models will be tested against experimental data provided through controlled laboratory experiments in heterogeneous porous media in large (meter-scale) 2-D flow cells. In addition to a mechanistic-based predictive model, this research will lead to new theories for the transient spatial distribution of microbial populations and contaminant plumes in heterogeneous porous media, improving the capability for designing staged remediation strategies for dealing with mixed contaminants.'
Ginn, T.R.; Boone, D.R.; Fletcher, M.M.; Friedrich, D.M.; Murphy, E.M.
1997-06-01
'Dynamic microbial attachment/detachment occurs in subsurface systems in response to changing environmental conditions caused by contaminant movement and degradation. Understanding the environmental conditions and mechanisms by which anaerobic bacteria partition between aqueous and solid phases is a critical requirement for designing and evaluating in situ bioremediation efforts. This interdisciplinary research project will provide fundamental information on the attachment/detachment dynamics of anaerobic bacteria in heterogeneous porous media under growth and growth-limiting conditions. Experiments will provide information on passive and active attachment/detachment mechanisms used by growing anaerobes capable of reductive dechlorination. Theoretical representations of these attachment/detachment mechanisms will be incorporated into existing groundwater flow and contaminant transport models that incorporate heterogeneity effects and can be used to predict behavior at field scales. These mechanistic-based models will be tested against experimental data provided through controlled laboratory experiments in heterogeneous porous media in large (meter-scale) 2-D flow cells. In addition to a mechanistic-based predictive model, this research will lead to new theories for the transient spatial distribution of microbial populations and contaminant plumes in heterogeneous porous media, improving the capability for designing staged remediation strategies for dealing with mixed contaminants.'
NASA Astrophysics Data System (ADS)
Gao, Mingliang; Teng, Qizhi; He, Xiaohai; Zuo, Chen; Li, ZhengJi
2016-01-01
Three-dimensional (3D) structures are useful for studying the spatial structures and physical properties of porous media. A 3D structure can be reconstructed from a single two-dimensional (2D) training image (TI) by using mathematical modeling methods. Among many reconstruction algorithms, an optimal-based algorithm was developed and has strong stability. However, this type of algorithm generally uses an autocorrelation function (which is unable to accurately describe the morphological features of porous media) as its objective function. This has negatively affected further research on porous media. To accurately reconstruct 3D porous media, a pattern density function is proposed in this paper, which is based on a random variable employed to characterize image patterns. In addition, the paper proposes an original optimal-based algorithm called the pattern density function simulation; this algorithm uses a pattern density function as its objective function, and adopts a multiple-grid system. Meanwhile, to address the key point of algorithm reconstruction speed, we propose the use of neighborhood statistics, the adjacent grid and reversed phase method, and a simplified temperature-controlled mechanism. The pattern density function is a high-order statistical function; thus, when all grids in the reconstruction results converge in the objective functions, the morphological features and statistical properties of the reconstruction results will be consistent with those of the TI. The experiments include 2D reconstruction using one artificial structure, and 3D reconstruction using battery materials and cores. Hierarchical simulated annealing and single normal equation simulation are employed as the comparison algorithms. The autocorrelation function, linear path function, and pore network model are used as the quantitative measures. Comprehensive tests show that 3D porous media can be reconstructed accurately from a single 2D training image by using the method proposed
Heat and mass transfer in unsaturated porous media. Final report
Childs, S.W.; Malstaff, G.
1982-02-01
A preliminary study of heat and water transport in unsaturated porous media is reported. The project provides background information regarding the feasibility of seasonal thermal energy storage in unconfined aquifers. A parametric analysis of the factors of importance, and an annotated bibliography of research findings pertinent to unconfined aquifer thermal energy storage (ATES) are presented. This analysis shows that heat and mass transfer of water vapor assume dominant importance in unsaturated porous media at elevated temperature. Although water vapor fluxes are seldom as large as saturated medium liquid water fluxes, they are important under unsaturated conditions. The major heat transport mechanism for unsaturated porous media at temperatures from 50 to 90/sup 0/C is latent heat flux. The mechanism is nonexistent under saturated conditions but may well control design of unconfined aquifer storage systems. The parametric analysis treats detailed physical phenomena which occur in the flow systems study and demonstrates the temperature and moisture dependence of the transport coefficients of importance. The question of design of an unconfined ATES site is also addressed by considering the effects of aquifer temperature, depth to water table, porous medium flow properties, and surface boundary conditions. Recommendations are made for continuation of this project in its second phase. Both scientific and engineering goals are considered and alternatives are presented.
Particle dispersion and deposition in porous media: a computational perspective
NASA Astrophysics Data System (ADS)
Boccardo, Gianluca; Crevacore, Eleonora; Sethi, Rajandrea; Marchisio, Daniele
2015-11-01
This work investigates particle dispersion in porous media, which is of central relevance in a number of applications ranging from groundwater remediation tochemical engineering. The challenge lies in studying the complex fluid dynamics behavior arising at the microscale (very difficult to observe experimentally) and obtaining transport models to be employed at the macroscopic scale of interest. While a wealth of studies have approached this problem, the case of particle transport with a concurrent heterogeneous chemical reaction (e.g.: particle deposition) still lacks a satisfactory description, especially when considering a polydisperse population of solid particles. Moreover, the oft-used simplified descriptions of the porous medium (via array of spheres or similar strategies) fail to fully take into account the effect of the packing structure. Our novel approach relies on an ``in-silico'' procedure where many 3-D realistic porous media models are constructed via rigid-body simulations and fluid flowand particle transport are then investigated through computational fluid dynamics. The results evidence the need for a deeper look, afforded by these methodology, into the influence of the features of realistic porous media on particle transport and deposition.
Modeling transport phenomena in porous media
Bear, J.
1996-12-31
The paper reviews the continuum approach to modelling the transport of mass, momentum and energy, of phases and of their components in a porous medium domain. The review begins with the definition of a porous medium, making use of the concept of a Representative Elementary Volume (REV) as a tool for overcoming the effect of the microscopic heterogeneity resulting from the presence of a solid matrix and a void space. The microscopic and macroscopic levels of description are defined. By averaging the description of a transport phenomenon at the microscopic level over an REV, using certain {open_quote}averaging rules{close_quote}, the macroscopic or continuum description of the same phenomenon is obtained. This methodology is first introduced in general terms for any extensive quantity, and then demonstrated for the transport of mass, momentum and energy. In the process of deriving the macroscopic models, expressions are presented also for the advective, dispersive and diffusive fluxes of extensive quantities that appear in them, in terms of averaged, measurable values of state variables.
Statistical mechanics of unsaturated porous media.
Xu, Jin; Louge, Michel Y
2015-12-01
We explore a mean-field theory of fluid imbibition and drainage through permeable porous solids. In the limit of vanishing inertial and viscous forces, the theory predicts the hysteretic "retention curves" relating the capillary pressure applied across a connected domain to its degree of saturation in wetting fluid in terms of known surface energies and void space geometry. To avoid complicated calculations, we adopt the simplest statistical mechanics, in which a pore interacts with its neighbors through narrow openings called "necks," while being either full or empty of wetting fluid. We show how the main retention curves can be calculated from the statistical distribution of two dimensionless parameters λ and α measuring the specific areas of, respectively, neck cross section and wettable pore surface relative to pore volume. The theory attributes hysteresis of these curves to collective first-order phase transitions. We illustrate predictions with a porous domain consisting of a random packing of spheres, show that hysteresis strength grows with λ and weakens as the distribution of α broadens, and reproduce the behavior of Haines jumps observed in recent experiments on an ordered pore network. PMID:26764701
Statistical mechanics of unsaturated porous media
NASA Astrophysics Data System (ADS)
Xu, Jin; Louge, Michel Y.
2015-12-01
We explore a mean-field theory of fluid imbibition and drainage through permeable porous solids. In the limit of vanishing inertial and viscous forces, the theory predicts the hysteretic "retention curves" relating the capillary pressure applied across a connected domain to its degree of saturation in wetting fluid in terms of known surface energies and void space geometry. To avoid complicated calculations, we adopt the simplest statistical mechanics, in which a pore interacts with its neighbors through narrow openings called "necks," while being either full or empty of wetting fluid. We show how the main retention curves can be calculated from the statistical distribution of two dimensionless parameters λ and α measuring the specific areas of, respectively, neck cross section and wettable pore surface relative to pore volume. The theory attributes hysteresis of these curves to collective first-order phase transitions. We illustrate predictions with a porous domain consisting of a random packing of spheres, show that hysteresis strength grows with λ and weakens as the distribution of α broadens, and reproduce the behavior of Haines jumps observed in recent experiments on an ordered pore network.
Curating Media Learning: Towards a Porous Expertise
ERIC Educational Resources Information Center
McDougall, Julian; Potter, John
2015-01-01
This article combines research results from a range of projects with two consistent themes. Firstly, we explore the potential for curation to offer a productive metaphor for the convergence of digital media learning across and between home/lifeworld and formal educational/system-world spaces--or between the public and private spheres. Secondly, we…
MULTIPHASE FLOW AND TRANSPORT IN POROUS MEDIA
Multiphase flow and transport of compositionally complex fluids in geologic media is of importance in a number of applied problems which have major social and economic effects. n petroleum reservoir engineering efficient recovery of energy reserves is the principal goal. nfortuna...
Slip effects associated with Knudsen transport phenomena in porous media
NASA Technical Reports Server (NTRS)
Frederking, T. H. K.; Hepler, W. A.; Khandhar, P. K.
1988-01-01
Porous media used in phase separators and thermomechanical pumps have been the subject of characterization efforts based on the Darcy permeability of laminar continuum flow. The latter is not always observed at low speed, in particular at permeabilities below 10 to the -9th/squared cm. The present experimental and theoretical studies address questions of slip effects associated with long mean free paths of gas flow at room temperature. Data obtained are in good agreement, within data uncertainty, with a simplified asymptotic Knudsen equation proposed for porous plugs on the basis of Knudsen's classical flow equation for long mean free paths.
NASA Astrophysics Data System (ADS)
Jourabian, Mahmoud; Farhadi, Mousa; Rabienataj Darzi, Ahmad Ali
2016-07-01
In this study, the melting process of ice as a phase-change material (PCM) saturated with a nickel-steel porous matrix inside a horizontal elliptical tube is investigated. Due to the low thermal conductivity of the PCM, it is motivated to augment the heat transfer performance of the system simultaneously by finding an optimum value of the aspect ratio and impregnating a metallic porous matrix into the base PCM. The lattice Boltzmann method with a double distribution function formulated based on the enthalpy method, is applied at the representative elementary volume scale under the local thermal equilibrium assumption between the PCM and porous matrix in the composite. While reducing or increasing the aspect ratio of the circular tubes leads to the expedited melting, the 90° inclination of each elliptical tube in the case of the pure PCM melting does not affect the melting rate. With the reduction in the porosity, the effective thermal conductivity and melting rate in all tubes promoted. Although the natural convection is fully suppressed due to the significant flow blockage in the porous structure, the melting rates are generally increased in all cases.
Theoretical relationships between reactivity and permeability for monomineralic porous media
Smith, R.W.; Schafer, A.L.; Tompson, A.F.B.
1996-08-01
The release of contaminants to the subsurface has lead to potential or actual contamination of groundwater resources. The remediation of these plumes requires improved capability for the prediction of contaminant fate and transport. Key to improving predictive capabilities is an increased understanding of natural chemical (or reactive) heterogeneity in subsurface media and how chemical and physical heterogeneity are related. Although the effects of physical heterogeneity on transport has received significant attention, similar investigations of chemical heterogeneity, and the correlation of chemical and physical heterogeneity are in their infancy. Theoretical considerations for unconsolidated porous media composed of sand-sized grains of a single reactive phase suggest that chemical reactivity and permeability are inversely related. This correlation is consistent with measured permeability, porosity, and surface area data for unconsolidated sands. The correlation arises because both reactivity and permeability are functions of the surface area of the porous media. This relationship suggests that for a heterogeneous porous media, significant contaminant transport will occur in zones of minimal reactivity.
Propagation of nitromethane detonations in porous media
NASA Astrophysics Data System (ADS)
Lee, J. J.; Frost, D. L.; Lee, J. H. S.; Dremin, A.
1995-06-01
The characteristics of the propagation of a detonation in chemically sensitized nitromethane in a dense porous medium are investigated. By introducing liquid NM+15% (by weight) DETA into densely packed beds of solid spherical glass beads 66μm to 2.4 mm in diameter, a highly heterogeneous explosive mixture is obtained. The critical (i.e., failure) charge diameter of this mixture is systematically measured in unconfined charges over a wide range of bead sizes. Velocity measurements are also made for the various charges. It is found that there exists a critical bead size above which the critical diameter decreases with increasing bead size and below which it decreases with decreasing bead size. This result indicates an abrupt change in the mechanism of propagation at the critical bead size. Velocity measurements further support this by emphasizing the different behavior above and below the critical point.
Attenuation of Shocks through Porous Media
NASA Astrophysics Data System (ADS)
Lind, Charles A.; Cybyk, Bohdan Z.; Boris, Jay P.
1998-11-01
Structures designed to mitigate the effects of blast and shock waves are important for both accidental and controlled explosions. The net effect of these mitigating structures is to reduce the strength of the transmitted shock thereby reducing the dynamic pressure loading on nearby objects. In the present study, the attenuation of planar blast and shock waves by passage through structured media is numerically studied with the FAST3D model. The FAST3D model is a state-of-the-art, portable, three-dimensional computational fluid dynamics model based on Flux-Corrected Transport and uses the Virtual Cell Embedding algorithm for simulating complex geometries. The effects of media placement, spacing, orientation, and area blockage are parametrically studied to enhance the understanding of the complex processes involved and to determine ways to minimize the adverse effects of these blast waves.
Horizontal flow and capillarity-driven redistribution in porous media.
Doster, F; Hönig, O; Hilfer, R
2012-07-01
A recent macroscopic mixture theory for two-phase immiscible displacement in porous media has introduced percolating and nonpercolating phases. Quasi-analytic solutions are computed and compared to the traditional theory. The solutions illustrate physical insights and effects due to spatiotemporal changes of nonpercolating phases, and they highlight the differences from traditional theory. Two initial and boundary value problems are solved in one spatial dimension. In the first problem a fluid is displaced by another fluid in a horizontal homogeneous porous medium. The displacing fluid is injected with a flow rate that keeps the saturation constant at the injection point. In the second problem a horizontal homogeneous porous medium is considered which is divided into two subdomains with different but constant initial saturations. Capillary forces lead to a redistribution of the fluids. Errors in the literature are reported and corrected. PMID:23005535
Pore-scale simulation of laminar flow through porous media
NASA Astrophysics Data System (ADS)
Piller, M.; Casagrande, D.; Schena, G.; Santini, M.
2014-04-01
The experimental investigation of flow through porous media is inherently difficult due to the lack of optical access. The recent developments in the fields of X-ray micro-tomography (micro-CT hereafter), digital sample reconstruction by image-processing techniques and fluid-dynamics simulation, together with the increasing power of super-computers, allow to carry out pore-scale simulations through digitally-reconstructed porous samples. The scientific relevance of pore-scale simulations lies in the possibility of upscaling the pore-level data, yielding volume-averaged quantities useful for practical purposes. One of the best-known examples of upscaling is the calculation of absolute and relative permeability of reservoir rocks. This contribution presents a complete work-flow for setting up pore-scale simulations, starting from the micro-CT of a (in general small) porous sample. Relevant applications are discussed in order to reveal the potential of the proposed methodology.
Effects of capillarity on microscopic flow in porous media
Miksis, M.J.
1993-01-01
Central theme of this proposal is to study effects of capillarity on motion of a fluid interface and to apply these results to flow in porous media. Here we report on several problems considered this year, the second year of the grant. In particular we have developed a numerical code to study the dynamics of a gas bubble in a pore in order to examine the fundamental mechanism for the generation of a foam in a porous material, we have started an investigation of the stability of a foam lamella in order to understand the stability of foam flow in a porous material and we have derived systematically a slip coefficient for flow over a rough surface, e.g., as in a pore. In addition we report on work on several other problems.
Oil ganglion dynamics in flow through porous media
NASA Astrophysics Data System (ADS)
Ng, K. M.
1980-12-01
A model is formulated to study the transient behavior of oil ganglion populations during immiscible displacement in oil recovery processes. The model is composed of three components: a suitable model for granular porous media; a stochastic simulation method capable of predicting the fate of solitary ganglia and two coupled population balance equations for studying the dynamics of oil ganglion populations. The porous medium model proposed is a network of interconnected unit cells of the constricted tube type. The permeability of this model is determined by both the statistical analysis the network analysis. It is demonstrated that fluids in different portions of a porous medium interact with one another. Two methods are used in the simulation of the motion of solitary ganglion. It is found that the velocity of an oil blob decreases as its viscosity increases.
Competition among flow, dissolution, and precipitation in porous media
Rege, S.D.; Fogler, H.S. . Dept. of Chemical Engineering)
1989-07-01
A theoretical and experimental study has been carried out on flow, dissolution, and precipitation in porous media. Flow experiments were performed on linear carbonate cores using acidic ferric chloride solutions. Dissolution of the carbonate by the acid causes an increase in the solution pH, thereby precipitating ferric hydroxide. This precipitate plugs up the pore throats in the medium and increases the resistance to fluid flow. Fluctuations in the permeability ratio were observed during core flood experiments, confirming the competition between channel formation due to dissolution and pore plugging due to precipitation. The evolution of the pore structure was characterized. A network model has also been developed to describe flow and reaction in porous media. The model was used to simulate the ferric chloride system, and pressure oscillations predicted by the model show identical trends to those observed experimentally. Additionally, the evolution of pores in the network were graphically represented.
On growth and flow: bacterial biofilms in porous media
NASA Astrophysics Data System (ADS)
Durham, William; Leombruni, Alberto; Tranzer, Olivier; Stocker, Roman
2011-11-01
Bacterial biofilms often occur in porous media, where they play pivotal roles in medicine, industry and the environment. Though flow is ubiquitous in porous media, its effects on biofilm growth have been largely ignored. Using patterned microfluidic devices that simulate unconsolidated soil, we find that the structure of Escherichia coli biofilms undergoes a self-organization mediated by the interaction of growth and flow. Intriguingly, we find that biofilm productivity peaks at intermediate flow rates, when the biofilm is irrigated by a minimum number of preferential flow channels. At larger and smaller flow rates, fluid flows more uniformly through the matrix, but productivity drops due to removal by shear and reduced nutrient transport, respectively. These dynamics are correctly predicted by a simple network model. The observed tradeoff between growth and flow may have important consequences on biofilm-mediated processes such as biochemical cycling, antibiotic resistance and water filtration.
Dissipative particle dynamics model for colloid transport in porous media
Pan, W.; Tartakovsky, A. M.
2013-08-01
We present that the transport of colloidal particles in porous media can be effectively modeled with a new formulation of dissipative particle dynamics, which augments standard DPD with non-central dissipative shear forces between particles while preserving angular momentum. Our previous studies have demonstrated that the new formulation is able to capture accurately the drag forces as well as the drag torques on colloidal particles that result from the hydrodynamic retardation effect. In the present work, we use the new formulation to study the contact efficiency in colloid filtration in saturated porous media. Note that the present model include all transport mechanisms simultaneously, including gravitational sedimentation, interception and Brownian diffusion. Our results of contact efficiency show a good agreement with the predictions of the correlation equation proposed by Tufenkji and EliMelech, which also incorporate all transport mechanisms simultaneously without the additivity assumption.
A numerical approach for groundwater flow in unsaturated porous media
NASA Astrophysics Data System (ADS)
Quintana, F.; Guarracino, L.; Saliba, R.
2006-07-01
In this article, a computational tool to simulate groundwater flow in variably saturated non-deformable fractured porous media is presented, which includes a conceptual model to obtain analytical expressions of water retention and hydraulic conductivity curves for fractured hard rocks and a numerical algorithm to solve the Richards equation. To calculate effective saturation and relative hydraulic conductivity curves we adopt the Brooks-Corey model assuming fractal laws for both aperture and number of fractures. A standard Galerkin formulation was employed to solve the Richards' equation together with a Crank-Nicholson scheme with Richardson extrapolation for the time discretization.The main contribution of this paper is to group an analytical model of the authors with a robust numerical algorithm designed to solve adequately the highly non-linear Richards' equation generating a tool for porous media engineering.
Plume dynamics in Hele-Shaw porous media convection.
Ecke, Robert E; Backhaus, Scott
2016-10-13
Mass transport in multi-species porous media is through molecular diffusion and plume dynamics. Predicting the rate of mass transport has application in determining the efficiency of the storage and sequestration of carbon dioxide. We study a water and propylene-glycol system enclosed in a Hele-Shaw cell with variable permeability that represents a laboratory analogue of the general properties of porous media convection. The interface between the fluids, tracked using an optical shadowgraph technique, is used to determine the mass transport rate, the spatial separation of solutal plumes, and the velocity and width characteristics of those plumes. One finds that the plume dynamics are closely related to the mass transport rate.This article is part of the themed issue 'Energy and the subsurface'. PMID:27597786
Resistance absorption of some groundwater tracers in porous media
NASA Astrophysics Data System (ADS)
Jafari, Fateme
2010-05-01
Absorption of tracer to the aquifer material is among the most important factors which should be considered when a tracing program is considered. In this study, the absorption of the tracer into the porous media is analyzed experimentally for some of the most important and applied tracers as uranine, rhodamine B, eosin, potassium permanganate, sodium chloride and potassium chloride. For each tracer, effect of initial tracer concentration and percentage of fine grain sediments on tracer absorption in porous media is analyzed. According to the final results, rhodamine B and potassium permanganate have the less resistance against absorption to aquifer material, whilst eosin and uranine are the most resistant tracers among the examined ones. Key Words: Tracer, Absorption, Aquifer, Column Method
Electrokinetic coupling in unsaturated porous media
Revil, A.; Linde, N.; Cerepi, A.; Jougnot, D.; Matthai, S.; Finsterle, S.
2007-02-27
We consider a charged porous material that is saturated bytwo fluid phases that are immiscible and continuous on the scale of arepresentative elementary volume. The wetting phase for the grains iswater and the nonwetting phase is assumed to be an electricallyinsulating viscous fluid. We use a volume-averaging approach to derivethe linear constitutive equations for the electrical current density aswell as the seepage velocities of the wetting and nonwetting phases onthe scale of a representative elementary volume. These macroscopicconstitutive equations are obtained by volume-averaging Ampere's lawtogether with the Nernst Planck equation and the Stokes equations. Thematerial properties entering the macroscopic constitutive equations areexplicitly described as functions of the saturation of the water phase,the electrical formation factor, and parameters that describe thecapillary pressure function, the relative permeability function, and thevariation of electrical conductivity with saturation. New equations arederived for the streaming potential and electro-osmosis couplingcoefficients. A primary drainage and imbibition experiment is simulatednumerically to demonstrate that the relative streaming potential couplingcoefficient depends not only on the water saturation, but also on thematerial properties of the sample, as well as the saturation history. Wealso compare the predicted streaming potential coupling coefficients withexperimental data from four dolomite core samples. Measurements on thesesamples include electrical conductivity, capillary pressure, thestreaming potential coupling coefficient at various level of saturation,and the permeability at saturation of the rock samples. We found verygood agreement between these experimental data and the modelpredictions.
Denitrification and chemotaxis of Pseudomonas stutzeri KC in porous media.
Roush, Caroline J; Lastoskie, Christian M; Worden, R Mark
2006-01-01
Chemotaxis is an important mechanism by which microorganisms are dispersed in porous media. A vigorous chemotactic response to concentration gradients formed by microbial consumption of chemoattractants can accelerate transport of bacteria to highly contaminated regions of soils and sediments, enhancing the efficiency of in situ bioremediation operations. Although chemotaxis plays a key role in establishment of biodegradation zones in the subsurface, the effects of physical heterogeneity on bacterial motility are poorly understood. To investigate the influence of porous media heterogeneity on microbial chemotaxis, swarm plate migration experiments were conducted using Pseudomonas stutzeri strain KC, a denitrifying bacterium used for in situ biodegradation of carbon tetrachloride in groundwater. Swarm plate measurements indicate that strain KC is strongly chemotactic toward both acetate and nitrate. A three-component mathematical model was developed to describe the migration of strain KC. Estimates of chemotactic sensitivity were obtained in the homogeneous (agar) phase and in a heterogeneous medium of aquifer solids extracted from the Schoolcraft bioremediation field site in western Michigan. Interestingly, the motility of strain KC is significantly larger in the porous medium than in the aqueous phase. We hypothesize that chemotactic response is enhanced within the heterogeneous medium because chemoattractant gradients formed by nitrate consumption are larger in the confined spaces of the porous medium than in unconfined agar solution. PMID:16760079
Normalized inverse characterization of sound absorbing rigid porous media.
Zieliński, Tomasz G
2015-06-01
This paper presents a methodology for the inverse characterization of sound absorbing rigid porous media, based on standard measurements of the surface acoustic impedance of a porous sample. The model parameters need to be normalized to have a robust identification procedure which fits the model-predicted impedance curves with the measured ones. Such a normalization provides a substitute set of dimensionless (normalized) parameters unambiguously related to the original model parameters. Moreover, two scaling frequencies are introduced, however, they are not additional parameters and for different, yet reasonable, assumptions of their values, the identification procedure should eventually lead to the same solution. The proposed identification technique uses measured and computed impedance curves for a porous sample not only in the standard configuration, that is, set to the rigid termination piston in an impedance tube, but also with air gaps of known thicknesses between the sample and the piston. Therefore, all necessary analytical formulas for sound propagation in double-layered media are provided. The methodology is illustrated by one numerical test and by two examples based on the experimental measurements of the acoustic impedance and absorption of porous ceramic samples of different thicknesses and a sample of polyurethane foam. PMID:26093413
Numerical method of crack analysis in 2D finite magnetoelectroelastic media
NASA Astrophysics Data System (ADS)
Zhao, Minghao; Xu, Guangtao; Fan, Cuiying
2010-04-01
The present paper extends the hybrid extended displacement discontinuity fundamental solution method (HEDD-FSM) (Eng Anal Bound Elem 33:592-600, 2009) to analysis of cracks in 2D finite magnetoelectroelastic media. The solution of the crack is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy the prescribed boundary conditions on the boundary of the domain and on the crack face. The Crouch fundamental solution for a parabolic element at the crack tip is derived to model the square root variations of near tip fields. The extended stress intensity factors are calculated under different electric and magnetic boundary conditions.
Reconstruction of Porous Media with Multiple Solid Phases
Losic; Thovert; Adler
1997-02-15
A process is proposed to generate three-dimensional multiphase porous media with fixed phase probabilities and an overall correlation function. By varying the parameters, a specific phase can be located either at the interface between two phases or within a single phase. When the interfacial phase has a relatively small probability, its shape can be chosen as granular or lamellar. The influence of a third phase on the macroscopic conductivity of a medium is illustrated. PMID:9056372
Unsteady-State Flow of Flexible Polymers in Porous Media.
Zitha, Pacelli L. J.; Chauveteau, Guy; Léger, Liliane
2001-02-15
In this paper we report an investigation of the unsteady-state flow of polymer solutions through granular porous media. The experiments were performed using high-molecular-weight nonionic and anionic polyacrylamides dissolved in water containing NaCl and model porous media obtained by packing silicon carbide (SiC) grains having a narrow grain size distribution. Before injection in porous media, the polymer solutions were carefully filtered according to a method that was proved to be efficient in removing any possibly remaining microgels. The SiC grain surface was passively oxidized by a controlled thermal treatment in order to obtain a surface partially covered by a thin silica layer having adsorption properties similar to those of quartzitic sand. By packing SiC grains of different sizes, porous media having identical adsorption properties and well-known pore throats sizes can be obtained with a good reproducibility. Parameters investigated include pore size, velocity gradient, polymer concentration, and adsorption energy. A striking unsteady-state flow behavior (pressure build-up at constant flow rate) is observed when three conditions are fulfilled: (a) the velocity gradient is larger than that known to be able to induce a coil-stretch transition, (b) the polymer adsorbs on the pore surfaces, and (c) the length of stretched macromolecules is larger than the effective pore throat diameter. When one of these conditions is not satisfied the flow remains steady. These observations are interpreted by a mechanism involving the adsorption and bridging across pore restrictions of elongated chains. We propose to refer to this peculiar mode of polymer adsorption as bridging adsorption. Copyright 2001 Academic Press. PMID:11161513
Dendrite Suppression by Shock Electrodeposition in Charged Porous Media
Han, Ji-Hyung; Wang, Miao; Bai, Peng; Brushett, Fikile R.; Bazant, Martin Z.
2016-01-01
It is shown that surface conduction can stabilize electrodeposition in random, charged porous media at high rates, above the diffusion-limited current. After linear sweep voltammetry and impedance spectroscopy, copper electrodeposits are visualized by scanning electron microscopy and energy dispersive spectroscopy in two different porous separators (cellulose nitrate, polyethylene), whose surfaces are modified by layer-by-layer deposition of positive or negative charged polyelectrolytes. Above the limiting current, surface conduction inhibits growth in the positive separators and produces irregular dendrites, while it enhances growth and suppresses dendrites behind a deionization shock in the negative separators, also leading to improved cycle life. The discovery of stable uniform growth in the random media differs from the non-uniform growth observed in parallel nanopores and cannot be explained by classic quasi-steady “leaky membrane” models, which always predict instability and dendritic growth. Instead, the experimental results suggest that transient electro-diffusion in random porous media imparts the stability of a deionization shock to the growing metal interface behind it. Shock electrodeposition could be exploited to enhance the cycle life and recharging rate of metal batteries or to accelerate the fabrication of metal matrix composite coatings. PMID:27307136
Upscaling of Constitutive Relations I Unsaturated Heterogeneous Porous Media
H. H. Liu; G. S. Bodvarsson
2001-05-31
When numerical model are used for modeling field scale flow and transport processes in the subsurface, the problem of ''upscaling'' arises. Typical scales, corresponding to spatial resolutions of subsurface heterogeneity in numerical models, are generally much larger than the measurement scale of the parameters and physical processes involved. The upscaling problems is, then, one of assigning parameters to gridblock scale based on parameter values measured on small scales. The focus of this study is to develop an approach to determine large-scale (upscaled) constitutive relations (relationships among relative permeability, capillary pressure and saturation) from small-scale measurements for porous media for a range of air entry values that are typical for the tuff matrix in the unsaturated zone of Yucca Mountain. For porous media with large air entry values, capillary forces play a key role in determining spatial water distribution at large-scales. Therefore, a relatively uniform capillary pressure approximately exists even for a large gridblock scale under steady state flow conditions. Based on these reasoning, we developed formulations that relate upscaled constitutive relations to ones measured at core-scale. Numerical experiments with stochastically generated heterogeneous porous media were used to evaluate the upscaling formulations.
Intermittent filtration of bacteria and colloids in porous media
NASA Astrophysics Data System (ADS)
Auset, Maria; Keller, Arturo A.; Brissaud, FrançOis; Lazarova, Valentina
2005-09-01
Intermittent filtration through porous media used for water and wastewater treatment can achieve high pathogen and colloid removal efficiencies. To predict the removal of bacteria, the effects of cyclic infiltration and draining events (transient unsaturated flow) were investigated. Using physical micromodels, we visualized the intermittent transport of bacteria and other colloids in unsaturated porous media. Column experiments provided quantitative measurements of the phenomena observed at the pore scale. Tagged Escherichia coli and a conservative tracer (NaI) were introduced in an initial pulse into a 1.5 m sand column. Subsequent hydraulic flushes without tagged bacteria or tracer were repeated every 4 hours for the next 4 days, during which outflow concentrations were monitored. Breakthrough behavior between colloids and dissolved tracer differed significantly, reflecting the differences in transport processes. Advancement of the wetting front remobilized bacteria which were held in thin water films, attached to the air-water interface (AWI), or entrapped in stagnant pore water between gas bubbles. In contrast, the tracer was only remobilized by diffusion from immobile to mobile water. Remobilization led to successive concentration peaks of bacteria and tracer in the effluent but with significant temporal differences. Observations at the pore-scale indicated that the colloids were essentially irreversibly attached to the solid-water interface, which explained to some extent the high removal efficiency of microbes in the porous media. Straining, cluster filtration, cell lysis, protozoa grazing, and bacteriophage parasitism could also contribute to the removal efficiency of bacteria.
Dendrite Suppression by Shock Electrodeposition in Charged Porous Media
NASA Astrophysics Data System (ADS)
Han, Ji-Hyung; Wang, Miao; Bai, Peng; Brushett, Fikile R.; Bazant, Martin Z.
2016-06-01
It is shown that surface conduction can stabilize electrodeposition in random, charged porous media at high rates, above the diffusion-limited current. After linear sweep voltammetry and impedance spectroscopy, copper electrodeposits are visualized by scanning electron microscopy and energy dispersive spectroscopy in two different porous separators (cellulose nitrate, polyethylene), whose surfaces are modified by layer-by-layer deposition of positive or negative charged polyelectrolytes. Above the limiting current, surface conduction inhibits growth in the positive separators and produces irregular dendrites, while it enhances growth and suppresses dendrites behind a deionization shock in the negative separators, also leading to improved cycle life. The discovery of stable uniform growth in the random media differs from the non-uniform growth observed in parallel nanopores and cannot be explained by classic quasi-steady “leaky membrane” models, which always predict instability and dendritic growth. Instead, the experimental results suggest that transient electro-diffusion in random porous media imparts the stability of a deionization shock to the growing metal interface behind it. Shock electrodeposition could be exploited to enhance the cycle life and recharging rate of metal batteries or to accelerate the fabrication of metal matrix composite coatings.
Dendrite Suppression by Shock Electrodeposition in Charged Porous Media.
Han, Ji-Hyung; Wang, Miao; Bai, Peng; Brushett, Fikile R; Bazant, Martin Z
2016-01-01
It is shown that surface conduction can stabilize electrodeposition in random, charged porous media at high rates, above the diffusion-limited current. After linear sweep voltammetry and impedance spectroscopy, copper electrodeposits are visualized by scanning electron microscopy and energy dispersive spectroscopy in two different porous separators (cellulose nitrate, polyethylene), whose surfaces are modified by layer-by-layer deposition of positive or negative charged polyelectrolytes. Above the limiting current, surface conduction inhibits growth in the positive separators and produces irregular dendrites, while it enhances growth and suppresses dendrites behind a deionization shock in the negative separators, also leading to improved cycle life. The discovery of stable uniform growth in the random media differs from the non-uniform growth observed in parallel nanopores and cannot be explained by classic quasi-steady "leaky membrane" models, which always predict instability and dendritic growth. Instead, the experimental results suggest that transient electro-diffusion in random porous media imparts the stability of a deionization shock to the growing metal interface behind it. Shock electrodeposition could be exploited to enhance the cycle life and recharging rate of metal batteries or to accelerate the fabrication of metal matrix composite coatings. PMID:27307136
Superfluid helium in fully saturated porous media
Huang, K. ); Meng, H. )
1993-09-01
The flow of superfluid [sup 4]He through spongelike media at full saturation is modeled by the flow of current through an Ohmic network with random resistors. Solving Kirchhoff's equations leads to the conclusion that the superfluid critical point is a percolation threshold, with critical exponent 1.7. The fractal dimension of the percolating cluster is 2.6. These lead to a specific-heat exponent [alpha]=[minus]5.4, by the Josephson hyperscaling relation. Existing experiments apparently do not cover the critical region. Instead, they measure mean-field'' exponents, whose values for Vycor, aerogel, and xerogel can all be reproduced by choosing appropriate distribution functions for the resistors.
Effects of surface active agents on DNAPL migration and distribution in saturated porous media.
Cheng, Zhou; Gao, Bin; Xu, Hongxia; Sun, Yuanyuan; Shi, Xiaoqing; Wu, Jichun
2016-11-15
Dissolved surface active agents such as surfactant and natural organic matter can affect the distribution and fate of dense nonaqueous liquids (DNAPLs) in soil and groundwater systems. This work investigated how two common groundwater surface active agents, humic acid (HA) and Tween 80, affected tetrachloroethylene (PCE) migration and source zone architecture in saturated porous media under environmentally relevant conditions. Batch experiments were first conducted to measure the contact angles and interfacial tensions (IFT) between PCE and quartz surface in water containing different amount of surface active agents. Results showed that the contact angle increased and IFT decreased with concentration of surface active agent increasing, and Tween 80 was much more effective than HA. Five 2-D flow cell experiments were then conducted. Correspondingly, Tween 80 showed strong effects on the migration and distribution of PCE in the porous media due to its ability to change the medium wettability from water-wet into intermediate/NAPL-wet. The downward migration velocities of the PCE in three Tween 80 cells were slower than those in the other two cells. In addition, the final saturation of the PCE in the cells containing surface active agents was higher than that in the water-only cell. Results from this work indicate that the presence of surface active agents in groundwater may strongly affect the fate and distribution of DNAPL through altering porous medium wettability. PMID:27450259
Numerical investigation of nanoparticles transport in anisotropic porous media.
Salama, Amgad; Negara, Ardiansyah; El Amin, Mohamed; Sun, Shuyu
2015-10-01
In this work the problem related to the transport of nanoparticles in anisotropic porous media is investigated numerically using the multipoint flux approximation. Anisotropy of porous media properties is an essential feature that exists almost everywhere in subsurface formations. In anisotropic media, the flux and the pressure gradient vectors are no longer collinear and therefore interesting patterns emerge. The transport of nanoparticles in subsurface formations is affected by several complex processes including surface charges, heterogeneity of nanoparticles and soil grain collectors, interfacial dynamics of double-layer and many others. We use the framework of the theory of filtration in this investigation. Processes like particles deposition, entrapment, as well as detachment are accounted for. From the numerical methods point of view, traditional two-point flux finite difference approximation cannot handle anisotropy of media properties. Therefore, in this work we use the multipoint flux approximation (MPFA). In this technique, the flux components are affected by more neighboring points as opposed to the mere two points that are usually used in traditional finite volume methods. We also use the experimenting pressure field approach which automatically constructs the global system of equations by solving multitude of local problems. This approach facilitates to a large extent the construction of the global system. A set of numerical examples is considered involving two-dimensional rectangular domain. A source of nanoparticles is inserted in the middle of the anisotropic layer. We investigate the effects of both anisotropy angle and anisotropy ratio on the transport of nanoparticles in saturated porous media. It is found that the concentration plume and porosity contours follow closely the principal direction of anisotropy of permeability of the central domain. PMID:26212784
Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media
NASA Technical Reports Server (NTRS)
Alfano, Robert R. (Inventor); Cai, Wei (Inventor); Gayen, Swapan K. (Inventor)
2000-01-01
A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise,
Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media
NASA Technical Reports Server (NTRS)
Alfano, Robert R. (Inventor); Cai, Wei (Inventor); Liu, Feng (Inventor); Lax, Melvin (Inventor); Das, Bidyut B. (Inventor)
1999-01-01
A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: ##EQU1## wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise,
Influence of biofilm accumulation on flow and reactive transport through porous media
NASA Astrophysics Data System (ADS)
Sharp, R.; Adgie, M.; Cunningham, A.
2003-04-01
A series of continuous flow, porous media reactor studies were performed to characterize the development of thick biofilms in porous media and the subsequent effects on flow and reactive transport. The bioluminescent organism Vibrio fischeri was used to produce various degrees of biofilm growth within the porous media system. V. fischeri biofilm growth, distribution, and activity in the porous media reactors was evaluated using bioluminescent imaging. Bulk fluid flow and flow channel dynamics in the porous media were monitored by imaging pulses of nigrosine bulk fluid dye. Hydrodynamics of the porous media/biofilm systems were analyzed using fluorescein break through curves and head loss across the system. Bioluminescent and bulk-fluid dye imaging, along with fluoroscein break through curve analysis, provided quantitative information on the transient distribution of biofilm within the reactor and the dynamic relationship between biofilm development and porous media hydraulics. Bioluminescent and bulk-fluid dye images showed continuous creation and closure of flow channels in the biofilm/porous media matrix. Flow channel size and distribution changed with increasing degrees of biofilm growth. Bioluminescence showed not only the density and distribution of biofilm growth in the porous media, but also the rate of oxygen uptake across the flow field. Results from the breakthrough curves suggest that thin biofilms transform the well-defined plug flow regime of clean porous media to a flow with severe axial and longitudinal dispersion. As the biofilm thickens, the flow regime is transformed to dispersed plug flow with an accelerated residence time. These studies provide a better quantitative understanding of the fundamental relationship between biofilm development and bulk fluid hydrodynamics in porous media. Results demonstrate the simultaneous visualization of biofilm growth and bulk fluid flow in porous media at the meso-scale. The studies also establish a novel
NASA Astrophysics Data System (ADS)
Turuban, R.; Jimenez-Martinez, J.; De Anna, P.; Tabuteau, H.; Meheust, Y.; Le Borgne, T.
2014-12-01
As dissolved chemical elements are transported in the subsurface, their mixing with other compounds and potential reactivity depends on the creation of local scale chemical gradients, which ultimately drive diffusive mass transfer and reaction. The distribution of concentration gradients is in turn shaped by the spatial gradients of flow velocity arising from the random distribution of solid grains. We present an experimental investigation of the relationship between the microscale flow stretching properties and the effective large scale mixing dynamics in porous media. We use a flow cell that models a horizontal quasi two-dimensional (2D) porous medium, the grains of which are cylinders randomly positioned between two glass plates [de Anna et al. 2013]. In this setup, we perform both non diffusive and diffusive transport tests, by injecting respectively microsphere solid tracers and a fluorescent dye. While the dye front propagates through the medium, it undergoes in time a kinematic stretching that is controlled by the flow heterogeneity, as it encounters stagnation zones and high velocity channels between the grains. The spatial distribution of the dye can then be described as a set of stretched lamellae whose rate of diffusive smoothing is locally enhanced by kinematic stretching [Le Borgne et al., 2013]. We show that this representation allows predicting the temporal evolution of the mixing rate and the probability distribution of concentration gradients for a range of Peclet numbers. This upscaling framework hence provides a quantification of the dynamics of effective mixing from the microscale Lagrangian velocity statistics. References:[1] P. de Anna, J. Jimenez-Martinez, H. Tabuteau, R. Turuban, T. Le Borgne, M. Derrien,and Yves Méheust, Mixing and reaction kinetics in porous media : an experimental pore scale quantification, Environ. Sci. Technol. 48, 508-516, 2014. [2] Le Borgne, T., M. Dentz, E. Villermaux, Stretching, coalescence and mixing in porous
Buckling in 2D periodic, soft and porous structures: effect of pore shape and lattice pattern
NASA Astrophysics Data System (ADS)
Shan, Sicong; Bertoldi, Katia; Shim, Jongmin; Overvelde, Johannes T. B.; Kang, Sung Hoon
2013-03-01
Adaptive structures allowing dramatic shape changes offer unique opportunities for the design of responsive and reconfigurable devices. Traditional morphing and foldable structures with stiff structural members and mechanical joints remains a challenge in manufacturing at small length scales. Soft structures where the folding mechanisms are induced by a mechanical instability represent a new class of novel adaptive materials which can be easily manufactured over a wide range of length scales. More specifically, soft porous structures with deliberately designed patterns can significantly change their architecture in response to diverse stimuli, opening avenues for reconfigurable devices that change their shapes to respond to their environment. While so far only two-dimensional periodic porous structures with circular holes arranged on a square or triangular lattice have been investigated, here we investigate both numerically and experimentally the effects of pore shape and lattice pattern on the macroscopic properties of the structures. Our results show that both the pore shape and lattice pattern can be used to effectively design desired materials and pave the way for the development of a new class of soft, active and reconfigurable devices over a wide range of length scales.
Lattice Boltzmann simulation of chemical dissolution in porous media.
Kang, Qinjun; Zhang, Dongxiao; Chen, Shiyi; He, Xiaoyi
2002-03-01
In this paper, we develop a lattice Boltzmann model for simulating the transport and reaction of fluids in porous media. To simulate such a system, we account for the interaction of forced convection, molecular diffusion, and surface reaction. The problem is complicated by the evolution of the porous media geometry due to chemical reactions, which may significantly and continuously modify the hydrologic properties of the media. The particular application that motivates the present study is acid stimulation, a common technique used to increase production from petroleum reservoirs. This technique involves the injection of acid (e.g., hydrochloric acid, HCl, acetic acid, HAc) into the formation to dissolve minerals comprising the rock. As acid is injected, highly conductive channels or "wormholes" may be formed. The dissolution of carbonate rocks in 0.5M HCl and 0.5M HAc is simulated with the lattice Boltzmann model developed in this study. The dependence of dissolution process and the geometry of the final wormhole pattern on the acid type and the injection rate is studied. The results agree qualitatively with the experimental and theoretical analyses of others and substantiate the previous finding that there exists an optimal injection rate at which the wormhole is formed as well as the number of pore volumes of the injected fluid to break through is minimized. This study also confirms the experimentally observed phenomenon that the optimal injection rate decreases and the corresponding minimized number of pore volumes to break through increases as the acid is changed from HCl to HAc. Simulations suggest that the proposed lattice Boltzmann model may serve as an alternative reliable quantitative approach to study chemical dissolution in porous media. PMID:11909255
Effect of Bacterial Motility on Contaminant Mixing in Porous Media
NASA Astrophysics Data System (ADS)
Singh, R.; Olson, M. S.; Bioremediation At Drexel
2010-12-01
Groundwater flow is typically characterized by laminar flow and therefore contaminant mixing limited conditions prevail in subsurface environments. The presence of porous media introduces tortuosity to groundwater flow paths, thereby enhancing contaminant mixing. In addition, bacterial motility is reported to induce movement in their surrounding liquid, which may enhance contaminant mixing. Enhancement of chemical diffusion coefficients in bulk fluid due to bacterial random motility and chemotaxis has been already reported in literature. The aim of this study is to investigate the effect of bacterial motility on contaminant mixing in the presence of porous media. A microfluidic device was designed and fabricated using standard photolithography and soft-lithography techniques to simulate a contaminant plume in subsurface porous media due to leakage of an underground storage tank. A non-reactive conservative tracer, Dextran solution labeled with FITC (fluorescein isothiocyanate), was used as surrogate for the contaminant and the motile bacterial strain Escherichia coli HCB33 (wild type) was used for the experiments to enhance contaminant mixing. Images were obtained at various cross-sections along the device and fluorescence intensity profile distributions were analyzed to determine the transverse dispersion of the contaminant. Enhancement in contaminant mixing was assessed by comparing the contaminant transverse dispersion coefficients (Dyi) in porous media in presence of motile bacteria, immobilized bacteria, and with no bacteria. In order to quantify the contaminant dispersion coefficients under the various test conditions, experimental data obtained were fitted to concentration profiles predicted by the contaminant advection-dispersion equation for the given experimental conditions (Figure 1). The transverse dispersion coefficient values obtained in the presence of motile bacteria (Dymb)and with no bacteria (Dynb) were 2.49 x 10-4 cm2/s and 1.39 x 10-4 cm2/s
Drying patterns of porous media containing wettability contrasts.
Shokri, N; Or, D
2013-02-01
Porous media containing sharp wettability discontinuities may occur in natural systems due to depositional processes, accumulation of organic layers or modification of soil wettability following intense forest fires all of which are known to significantly modify water flow and transport processes. We studied evaporation from sand columns containing hydraulically-interacting domains with sharp wettability contrasts. We employed neutron transmission technique to map liquid phase dynamics during evaporation, and conducted laboratory experiments to evaluate evaporative fluxes affected by interactions across wettability contrast. We explained the preferential drying front displacement in the hydrophobic domain and the spatial extent of capillary flow supporting the vaporization plane using a physically-based model. The model provides description of observed liquid phase patterns and dynamics observed in neutron radiography measurements and evaporative fluxes from laboratory experiments. Our results provide new insights into evaporation induced capillary exchange and preferential liquid phase distribution during evaporation from hydraulically interacting vertical porous domains with differing wettability properties and offer opportunities for design of selectively drying of porous media in natural and engineered systems. PMID:23123032
Measurement of Interfacial Area Production and Permeability within Porous Media
Crandall, Dustin; Ahmadi, Goodarz; Smith, Duane H.
2010-01-01
An understanding of the pore-level interactions that affect multi-phase flow in porous media is important in many subsurface engineering applications, including enhanced oil recovery, remediation of dense non-aqueous liquid contaminated sites, and geologic CO2 sequestration. Standard models of two-phase flow in porous media have been shown to have several shortcomings, which might partially be overcome using a recently developed model based on thermodynamic principles that includes interfacial area as an additional parameter. A few static experimental studies have been previously performed, which allowed the determination of static parameters of the model, but no information exists concerning the interfacial area dynamic parameters. A new experimental porous flow cell that was constructed using stereolithography for two-phase gas-liquid flow studies was used in conjunction with an in-house analysis code to provide information on dynamic evolution of both fluid phases and gas-liquid interfaces. In this paper, we give a brief introduction to the new generalized model of two-phase flow model and describe how the stereolithography flow cell experimental setup was used to obtain the dynamic parameters for the interfacial area numerical model. In particular, the methods used to determine the interfacial area permeability and production terms are shown.
NASA Astrophysics Data System (ADS)
Meheust, Y.; Turuban, R.; Jimenez-Martinez, J.; De Anna, P.; Tabuteau, H.; Le Borgne, T.
2014-12-01
Pore scale characterization of flow velocities and concentration spatial distributions is a key to understanding non-Fickian transport and mixing in porous media. We present a millifluidic setup aimed at investigating those processes in transparent porous media, at the pore scale. The porous media are quasi-2D, consisting of a Hele-Shaw cell containing cylindrical grains. They are made by soft lithography from a numerical model and provide full control on the geometry (medium porosity, permeability and heterogeneity). The setup allows for the study of primary drainage/imbibition, or the joint continuous injection of two fluids (e. g. water and air). A camera records the distributions of fluid phases, the position of solid tracers, and spatially-resolved images of light emissions inside the flow cell. The pore scale velocity field is thus measured from particle tracking, while pore scale concentration fields are measured accurately in passive transport experiments, using fluorescein; both continuous injection and finite volume solute injections can be achieved. Using two chemo-luminescent liquids, the reaction of which produces photons in addition to the reaction product, we are also able to study the local production rate of the reaction product as the reactive liquids flow through the system [1]. Pressure drops across the medium are also measured. This complete characterization (phase distributions, velocity and concentration fields, pressure drops) of the system allows to explain non-Fickian behaviors and test models that upscale transport and mixing properties from pore scale data. As examples, we shall discuss the upscaling of transport from the knowledge of Lagrangian velocities and the relationships between conservative and reactive transport under mixing-limited conditions (very large Damkhöler number). Other applications include the prediction of the mixing rate from the sole knowledge of the flow stretching [2], and the characterization of mixing by
A cellular automaton simulation of contaminant transport in porous media
Freed, D.M.; Simonson, S.A.
1995-12-01
A simulation tool to investigate radionuclide transport in porous groundwater flow is described. The flow systems of interest are those important in determining the fate of radionuclides emplaced in an underground repository, such as saturated matrix flow, matrix and fracture flow in the unsaturated zone, and viscous fingering in porous fractures. The work discussed here is confined to consideration of saturated flow in porous media carrying a dilute, sorptive species. The simulation technique is based on a special class of cellular automata known as lattice gas automata (LGA) which are capable of predicting hydrodynamic behavior. The original two-dimensional scheme (that of Frisch et. al. known as the FHP model) used particles of unit mass traveling on a triangular lattice with unit velocity and undergoing simple collisions which conserve mass and momentum at each node. These microscopic rules go over to the incompressible Navier-Stokes equations in the macroscopic limit. One of the strengths of this technique is the natural way that heterogeneities, such as boundaries, are accommodated. Complex geometries such as those associated with porous microstructures can be modeled effectively. Several constructions based on the FHP model have been devised, including techniques to eliminate statistical noise, extension to three dimensions, and the addition of surface tension which leads to multiphase flow.
P-adic model of transport in porous disordered media
NASA Astrophysics Data System (ADS)
Khrennikov, Adrei Yu.; Oleschko, Klaudia
2014-05-01
The soil porosity and permeability are the most important quantitative indicators of soil dynamics under the land-use change. The main problema in the modeling of this dynamic is still poor correlation between the real measuring data and the mathematical and computer simulation models. In order to overpassed this deep divorce we have designed a new technique, able to compare the data arised from the multiscale image analices and time series of the basic physical properties dynamics in porous media studied in time and space. We present a model of the diffusion reaction type describing transport in disordered porous media, e.g., water or oil flow in a complex network of pores. Our model is based on p-adic representation of such networks. This is a kind of fractal representation. We explore advantages of p- adic representation, namely, the possibility to endow p-adic trees with an algebraic structure and ultrametric topology and, hence, to apply analysis which have (at least some) similarities with ordinary real analysis on the straight line. We present the system of two diffusion reaction equations describing propagation of particles in networks of pores in disordered media. As an application, one can consider water transport through the soil pore Networks, or oil flow through capillaries nets. Under some restrictions on potentials and rate coefficients we found the stationary regime corresponding to water content or concentration of oil in a cluster of capillaries. Usage of p-adic analysis (in particular, p-adic wavelets) gives a possibility to find the stationary solution in the analytic form which makes possible to present a clear pedological or geological picture of the process. The mathematical model elaborated in this paper (Khrennikov, 2013) can be applied to variety of problems from water concentration in aquifers to the problem of formation of oil reservoirs in disordered media with porous structures. Another possible application may have real practical
Fixation of radioactive ions in porous media with ion exchange gels
Mercer, Jr., Basil W.; Godfrey, Wesley L.
1979-01-01
A method is provided for fixing radioactive ions in porous media by injecting into the porous media water-soluble organic monomers which are polymerizable to gel structures with ion exchange sites and polymerizing the monomers to form ion exchange gels. The ions and the particles of the porous media are thereby physically fixed in place by the gel structure and, in addition, the ions are chemically fixed by the ion exchange properties of the resulting gel.
An electrical conductivity model for fractal porous media
NASA Astrophysics Data System (ADS)
Wei, Wei; Cai, Jianchao; Hu, Xiangyun; Han, Qi
2015-06-01
Archie's equation is an empirical electrical conductivity-porosity model that has been used to predict the formation factor of porous rock for more than 70 years. However, the physical interpretation of its parameters, e.g., the cementation exponent m, remains questionable. In this study, a theoretical electrical conductivity equation is derived based on the fractal characteristics of porous media. The proposed model is expressed in terms of the tortuosity fractal dimension (DT), the pore fractal dimension (Df), the electrical conductivity of the pore liquid, and the porosity. The empirical parameter m is then determined from physically based parameters, such as DT and Df. Furthermore, a distinct interrelationship between DT and Df is obtained. We find a reasonably good match between the predicted formation factor by our model and experimental data.
Modeling Multi-process Transport of Pathogens in Porous Media
NASA Astrophysics Data System (ADS)
Cheng, L.; Brusseau, M. L.
2004-12-01
The transport behavior of microorganisms in porous media is of interest with regard to the fate of pathogens associated with wastewater recharge, riverbank filtration, and land application of biosolids. This interest has fomented research on the transport of pathogens in the subsurface environment. The factors influencing pathogen transport within the subsurface environment include advection, dispersion, filtration, and inactivation. The filtration process, which mediates the magnitude and rate of pathogen retention, comprises several mechanisms such as attachment to porous-medium surfaces, straining, and sedimentation. We present a mathematical model wherein individual filtration mechanisms are explicitly incorporated along with advection, dispersion, and inactivation. The performance of the model is evaluated by applying it to several data sets obtained from miscible-displacement experiments conducted using various pathogens. Input parameters are obtained to the extent possible from independent means.
High-precision synthetic computed tomography of reconstructed porous media
NASA Astrophysics Data System (ADS)
Hilfer, R.; Zauner, Th.
2011-12-01
Multiscale simulation of transport in disordered and porous media requires microstructures covering several decades in length scale. X-ray and synchrotron computed tomography are presently unable to resolve more than one decade of geometric detail. Recent advances in pore scale modeling [Biswal, Held, Khanna, Wang, and Hilfer, Phys. Rev. E PLEEE81539-375510.1103/PhysRevE.80.041301 80, 041301 (2009)] provide strongly correlated microstructures with several decades in microstructural detail. A carefully calibrated microstructure model for Fontainebleau sandstone has been discretized into a suite of three-dimensional microstructures with resolutions from roughly 128 μm down to roughly 500 nm. At the highest resolution the three-dimensional image consists of 327683=35184372088832 discrete cubic volume elements with gray values between 0 and 216. To the best of our knowledge, this synthetic image is the largest computed tomogram of a porous medium available at present.
1D and 2D simulations of seismic wave propagation in fractured media
NASA Astrophysics Data System (ADS)
Möller, Thomas; Friederich, Wolfgang
2016-04-01
Fractures and cracks have a significant influence on the propagation of seismic waves. Their presence causes reflections and scattering and makes the medium effectively anisotropic. We present a numerical approach to simulation of seismic waves in fractured media that does not require direct modelling of the fracture itself, but uses the concept of linear slip interfaces developed by Schoenberg (1980). This condition states that at an interface between two imperfectly bonded elastic media, stress is continuous across the interface while displacement is discontinuous. It is assumed that the jump of displacement is proportional to stress which implies a jump in particle velocity at the interface. We use this condition as a boundary condition to the elastic wave equation and solve this equation in the framework of a Nodal Discontinuous Galerkin scheme using a velocity-stress formulation. We use meshes with tetrahedral elements to discretise the medium. Each individual element face may be declared as a slip interface. Numerical fluxes have been derived by solving the 1D Riemann problem for slip interfaces with elastic and viscoelastic rheology. Viscoelasticity is realised either by a Kelvin-Voigt body or a Standard Linear Solid. These fluxes are not limited to 1D and can - with little modification - be used for simulations in higher dimensions as well. The Nodal Discontinuous Galerkin code "neXd" developed by Lambrecht (2013) is used as a basis for the numerical implementation of this concept. We present examples of simulations in 1D and 2D that illustrate the influence of fractures on the seismic wavefield. We demonstrate the accuracy of the simulation through comparison to an analytical solution in 1D.
Seismic wave propagation in multiphasic complex porous media : upscaling and downscaling approaches
NASA Astrophysics Data System (ADS)
Dupuy, B.; Garambois, S.; Virieux, J.
2011-12-01
Seismic wave propagation is often used for subsurface investigation, related either to reservoir issues (oil, gas or CO2 storage) or geotechnical problems (slope stability, water resources, territory management). Indeed, near surface media are rather heterogeneous, complex and partially fluid-filled. These characteristics are more or less sensitive to seismic waves. In order to interpret efficiently seismic attributes in these media, wave propagation may require complex poroelastic theories in multiphasic media as double porosity or patchy saturated cases. Thanks to upscaling techniques, we determine homogeneized parameters leading to a two-phases media described by complex and frequency dependent parameters. An effective generalized Biot-Gassmann theory allows wave propagation in 2D heterogeneous media through a Discontinuous Galerkin finite-element approach. Taking into account the rather complex frequency-dependent rheology of these porous media, wave propagation simulation is simpler in a frequency formulation than in a time formulation, at least for 2D geometries. We illustrate two features essential for an accurate characterization of the medium. On one hand, strong waveforms differences between the complex upscaled media (double porosity, unsaturated, squirt flow models) and the equivalent saturated Biot-Gassmann media (determined by arithmetic and harmonic averages) are underlined on simple examples. This may require the use of these theories for the characterization of these media. In the other hand, after the reconstruction of macro-scale parameters such as velocities and attenuations through seismic attributes (times, amplitudes and so on) using standard visco-elastic interpretation (first step of the downscaling procedure), we show that we recover micro-scale parameters (skeleton parameters such as porosity or dry moduli, fluid saturation...) using global search techniques (neighborhood algorithm) with some a priori information (second step of the
Solute transport in dual-permeability porous media
NASA Astrophysics Data System (ADS)
Leij, Feike J.; Toride, Nobuo; Field, Malcolm S.; Sciortino, Antonella
2012-04-01
A dual-advection dispersion equation (DADE) is presented and solved to describe solute transport in structured or layered porous media with different nonzero flow rates in two distinct pore domains with linear solute transfer between them. This dual-permeability model constitutes a generalized version of the advection-dispersion equation (ADE) for transport in uniform porous media and the mobile-immobile model (MIM) for transport in media with a mobile and an immobile pore domain. Analytical tools for the DADE have mostly been lacking. An analytical solution has therefore been derived using Laplace transformation with time and modal decomposition based on matrix diagonalization, assuming the same dispersivity for both domains. Temporal moments are derived for the DADE and contrasted with those for the ADE and the MIM. The effective dispersion coefficient for the DADE approaches that of the ADE for a similar velocity in both pore domains and large values for the first-order transfer parameter, and approaches that of the MIM for the opposite conditions. The solution of the DADE is used to illustrate how differences in pore water velocity between the domains and low transfer rates will lead to double peaks in the volume- or flux-averaged concentration profiles versus time or position. The DADE is applied to optimize experimental breakthrough curves for an Andisol with a distinct intra- and interaggregate porosity. The DADE improved the description of the breakthrough data compared to the ADE and the MIM.
Calculation of the effective permeability of saturated random porous media
NASA Astrophysics Data System (ADS)
Ostvar, S.; Wood, B. D.; Apte, S.; Liburdy, J.
2014-12-01
Estimation of the effective permeability tensor is an essential part of Darcy-scale representations of flow in porous media. The permeability tensor itself is a property of the medium, and depends only on the microscale geometry of the system. Determining the functional relationships between effective permeability (or conductivity in the general sense) and the structure of the medium is an old problem, with the earliest results for ordered porous media dating the 1920's. In this presentation, we report on the results of (1) detailed theory development, and (2) computations for the effective permeability tensor in fully-saturated random sphere packs, with a focus on the computational results. The theory is developed by volume averaging the Stokes equations, and using developing appropriate closures via potential theory, and has been reported on previously. For the computations, we have adopted an immersed boundary method to fully resolve the pore-scale velocity field. From our results, we compute the hydraulic permeability for both ordered and random media, and we compare these results with existing analytical solutions for the hydraulic conductivity in periodic arrays.
Evaluation of hydrodynamic scaling in porous media using finger dimensions
NASA Astrophysics Data System (ADS)
Selker, John S.; Schroth, Martin H.
1998-08-01
The use of dimensionless scaling is ubiquitous to hydrodynamic analysis, providing a powerful method of extending limited experimetnal results and generalizing theories. Miller and Miller [1956] contributed a scaling framework for immiscible fluid flow through porous media that relied on consistency of the contact angle between systems to be compared. It is common to assume that the effective contact angle will be zero in clean sand material where water is the wetting liquid. The well-documented unstable wetting process of fingered flow is used here as a diagnostic tool for the scaling relationships for infiltration into sandy media. Through comparison of finger cross sections produced using three liquids as well as various concentrations of anionic surfactant, it is shown that the zero contact angle assumption is very poor even for laboratory cleaned silica sand: Experimental results demonstrate effective contact angles approaching 60°. Scaling was effective for a given liquid between sands of differing particle size. These results suggest that caution should be exercised when applying scaling theory to initial wetting of porous media by liquids of differing gas-liquid interfacial tensions.
Hydraulic properties of adsorbed water films in unsaturated porous media
Tokunaga, Tetsu K.
2009-03-01
Adsorbed water films strongly influence residual water saturations and hydraulic conductivities in porous media at low saturations. Hydraulic properties of adsorbed water films in unsaturated porous media were investigated through combining Langmuir's film model with scaling analysis, without use of any adjustable parameters. Diffuse double layer influences are predicted to be important through the strong dependence of adsorbed water film thickness (f) on matric potential ({Psi}) and ion charge (z). Film thickness, film velocity, and unsaturated hydraulic conductivity are predicted to vary with z{sup -1}, z{sup -2}, and z{sup -3}, respectively. In monodisperse granular media, the characteristic grain size ({lambda}) controls film hydraulics through {lambda}{sup -1} scaling of (1) the perimeter length per unit cross sectional area over which films occur, (2) the critical matric potential ({Psi}{sub c}) below which films control flow, and (3) the magnitude of the unsaturated hydraulic conductivity when {Psi} < {Psi}{sub c}. While it is recognized that finer textured sediments have higher unsaturated hydraulic conductivities than coarser sands at intermediate {Psi}, the {lambda}{sup -1} scaling of hydraulic conductivity predicted here extends this understanding to very low saturations where all pores are drained. Extremely low unsaturated hydraulic conductivities are predicted under adsorbed film-controlled conditions (generally < 0.1 mm y{sup -1}). On flat surfaces, the film hydraulic diffusivity is shown to be constant (invariant with respect to {Psi}).
Dynamic permeability of porous media by the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Adler, P.; Pazdniakou, A.
2012-04-01
The main objective of our work is to determine the dynamic permeability of three dimensional porous media by means of the Lattice Boltzmann method (LBM). The Navier-Stokes equation can be numerically solved by LBM which is widely used to address various fluid dynamics problems. Space is discretized by a three-dimensional cubic lattice and time is discretized as well. The generally accepted notation for lattice Boltzmann models is DdQq where D stands for space dimension and Q for the number of discrete velocities. The present model is denoted by D3Q19. Moreover, the Two Relaxation Times variant of the Multi Relaxation Times model is implemented. Bounce back boundary conditions are used on the solid-fluid interfaces. The porous medium is spatially periodic. Reconstructed media were used; they are obtained by imposing a porosity and a correlation function characterized by a correlation length. Real samples can be obtained by MicroCT. In contrast with other previous contributions, the dynamic permeability K(omega) which is a complex number, is derived by imposing an oscillating body force of pulsation omega on the unit cell and by deriving the amplitude and the phase shift of the resulting time dependent seepage velocity. The influence of two limiting parameters, namely the Knudsen number Kn and the discretization for high frequencies, on K(omega) is carefully studied for the first time. Kn is proportional to nu/(cs H) where nu is the kinematic viscosity, cs the speed of sound in the fluid and H a characteristic length scale of the porous medium. Several porous media such as the classical plane Poiseuille flow and the reconstructed media are used to show that it is only for small enough values of Kn that reliable results are obtained. Otherwise, the data depend on Kn and may even be totally unphysical. However, it should be noticed that the limiting value of Kn could not be derived in general since it depends very much on the structure of the medium. Problems occur at
Dynamic permeability of porous media by the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Pazdniakou, A.; Adler, P. M.
2011-12-01
The main objective of our work is to determine the dynamic permeability of three dimensional porous media by means of the Lattice Boltzmann method (LBM). The Navier-Stokes equation can be numerically solved by LBM which is widely used to address various fluid dynamics problems. Space is discretized by a three-dimensional cubic lattice and time is discretized as well. The generally accepted notation for lattice Boltzmann models is DdQq where D stands for space dimension and Q for the number of discrete velocities. The present model is denoted by D3Q19. Moreover, the Two Relaxation Times variant of the Multi Relaxation Times model is implemented. Bounce back boundary conditions are used on the solid-fluid interfaces. The porous medium is spatially periodic. Reconstructed media were used; they are obtained by imposing a porosity and a correlation function characterized by a correlation length. Real samples can be obtained by MicroCT. In contrast with other previous contributions, the dynamic permeability K(omega) which is a complex number, is derived by imposing an oscillating body force of pulsation omega on the unit cell and by deriving the amplitude and the phase shift of the resulting time dependent seepage velocity. The influence of two limiting parameters, namely the Knudsen number Kn and the discretization for high frequencies, on K(omega) is carefully studied for the first time. Kn is proportional to nu/(c_s H) where nu is the kinematic viscosity, c_s the speed of sound in the fluid and H a characteristic length scale of the porous medium. Several porous media such as the classical plane Poiseuille flow and the reconstructed media are used to show that it is only for small enough values of Kn that reliable results are obtained. Otherwise, the data depend on Kn and may even be totally unphysical. However, it should be noticed that the limiting value of Kn could not be derived in general since it depends very much on the structure of the medium. Problems occur
2-D Gaussian beam imaging of multicomponent seismic data in anisotropic media
NASA Astrophysics Data System (ADS)
Protasov, M. I.
2015-12-01
An approach for true-amplitude seismic beam imaging of multicomponent seismic data in 2-D anisotropic elastic media is presented and discussed. Here, the recovered true-amplitude function is a scattering potential. This approach is a migration procedure based on the weighted summation of pre-stack data. The true-amplitude weights are computed by applying Gaussian beams (GBs). We shoot a pair of properly chosen GBs with a fixed dip and opening angles from the current imaging point towards an acquisition system. This pair of beams is used to compute a true-amplitude selective image of a rapid velocity variation. The total true-amplitude image is constructed by superimposing selective images computed for a range of available dip angles. The global regularity of the GBs allows one to disregard whether a ray field is regular or irregular. P- and S-wave GBs can be used to handle raw multicomponent data without separating the waves. The use of anisotropic GBs allows one to take into account the anisotropy of the background model.
Static properties of polymer chains in porous media
NASA Astrophysics Data System (ADS)
Honeycutt, J. D.; Thirumalai, D.
1989-04-01
The static properties of a polymer molecule in a porous medium are investigated. The porous medium is simulated using a site percolation model in which the various sites are occupied (or unoccupied) randomly. A freely jointed chain is allowed to move in continuous space between the obstacles. Effects of excluded volume interactions between the links have also been studied. Using a generalized Flory theory, we have shown that, when the strength of disorder is large enough, the mean square end-to-end distance scales as N2ν, where N is the number of links in the chain, and ν takes on a value different from that for a free chain. Under these conditions, the polymer assumes a compact, globule-like conformation. For sufficiently large N, the Flory theory gives ν=1/(d+2) for freely jointed chains and ν=1/d for chains with excluded volume. Various correlation functions such as the distribution of the end-to-end distance and density profile of monomers with respect to the center of mass of the chain have been computed using Monte Carlo simulations. These results are interpreted using scaling concepts and an approximate variational theory based on replica methods. The limitations of the replica variational theory are assessed by an application to the directed polymer in a quenched random environment. We have also studied the shape fluctuations that the polymer molecule undergoes in the random environment. It is argued that these shape fluctuations are relevant to the transport mechanism of polymers in random media. The results obtained for the porous media are contrasted with those found for polymers in media where the obstacles are arranged in a regular manner.
Direct numerical simulation of inertial flows in porous media
NASA Astrophysics Data System (ADS)
Apte, S.; Finn, J.; Wood, B. D.
2010-12-01
At modest flow rates (10 ≤ Re ≤ 300) through porous media and packed beds, fluid inertia can result in complex steady and unsteady recirculation regions, dependent on the local pore geometry. Body fitted CFD is a broadly used design and analysis tool for flows in porous media and packed bed type reactors. Unfortunately, the inherent complexities of porous media make unstructured mesh generation a difficult and time consuming step in the simulation process. To accurately capture the inertial dynamics using high-fidelity direct simulations, body fitted meshes must be high quality and sufficiently refined. We present methods to parameterize and simplify mesh generation for packed beds, with an eye toward obtaining efficient mesh independence for Reynolds numbers in the inertial and unsteady regimes. The crux of mesh generation for packed beds is dealing with sphere-sphere or sphere-wall contact points, where a geometric singularity exists. To handle the sphere-sphere and sphere-wall contact points, we use a fillet bridge model, in which every pair of contacting entities are bridged by a fillet, eliminating a small fluid region near the contact point. This results in a continuous surface mesh which does not require resizing of the spheres and can accommodate prism cells for improved boundary layer resolution. A second order accurate, parallel, incompressible flow solver [Moin and Apte, AIAA J. 2006] is used to simulate flow through three different sphere packings: a periodic simple cubic packing, a wall bounded hexagonal close packing, and a randomly packed tube. Mesh independence is assessed using several measures including Ergun pressure drop coefficients, viscous and pressure components of drag force, kinetic energy, kinetic energy dissipation and interstitial velocity profiles. The results of these test cases are used to determine the feasibility of accurate and very large scale simulations of flow through a randomly packed bed of 103 pores. Preliminary results
Theoretical analysis of three phase flow experiments in porous media
Sahni, A.; Blunt, M.; Guzman, R.
1996-12-31
This paper presents a framework for analyzing three phase flow experiments in porous media. We analyze previously published data from a dynamic displacement experiment, where three phase relative permeabilities have been measured. The relative permeability of each phase is, to a good approximation, a polynomial function of the saturation of that phase. Then analytically, we calculate the saturation paths and recovery for the experiments using the method of characteristics. The predicted paths agree well with the experimental measurements, and with one-dimensional numerical solutions.
NMR imaging of immiscible displacements in porous media
Majors, P.D.; Li, P.; Peters, E.J.
1995-12-31
We introduce a rapid, quantitative nuclear magnetic resonance imaging (NMRI) technique to resolve and measure multiple fluid phases in porous media. Liquids are resolved on the basis of their NMR spin-spin (T{sub 2}) relaxation times, and their intensities are corrected via attenuation analysis. The spatially resolved and corrected NMRI intensities are normalized to yield fluid saturations. In-situ saturation measurements are presented for three immiscible (oil and water) displacements in the same Berea sandstone core. NMRI and effluent recovery methods were compared. T{sub 2} of the displacement fluids were observed to be sensitive to displacement conditions.
Rayleigh-Taylor instability of immiscible fluids in porous media
NASA Astrophysics Data System (ADS)
Kalisch, H.; Mitrovic, D.; Nordbotten, J. M.
2016-05-01
The time development of an interface separating two immiscible fluids of different densities in heterogeneous two-dimensional porous media is studied. The governing equations are simplified with the help of approximate Green's functions which allow computation of the shape of the interface directly without resolving the fluid flow in the entire domain. The new formulation is amenable to numerical approximation, and the reduction in dimension leads to a significant gain in efficiency in the numerical simulation of the interfacial dynamics. Several test cases are investigated, and the numerical solutions are compared to known exact solutions and experimental data.
Combined Evaporation and Salt Precipitation in Porous Media
NASA Astrophysics Data System (ADS)
Weisbrod, N.; Dragila, M. I.; Nachshon, U.; Or, D.; Shaharani, E.; Grader, A.
2012-12-01
The vadose zone pore water contains dissolved salts and minerals; therefore, evaporation results in high rates of salt accumulation that may change the physical and chemical properties of the porous media. Here, a series of experiments, together with a mathematical model, are presented to shed new light on these processes. Experiments included: (1) long-term column evaporation experiments to quantify changes in evaporation rates due to salt precipitation; (2) CT scans of evaporated porous media samples saturated with salt solutions, to observe salt precipitation from micro to macro scales; and (3) Infrared thermography analysis to quantify evaporation rates from porous media surfaces for homogeneous and heterogeneous conditions and constant water table, in the presence of salt precipitation. As expected, the majority of salt crystallization occurs in the upper parts of the matrix, near the evaporation front. For heterogeneous porous matrices, salt precipitation will occur mainly in the fine pore regions as preferential evaporation takes place in these locations. In addition, it was found that the precipitated NaCl salt crust diffusion coefficient for water vapor is one to two orders of magnitude lower than the vapor diffusion coefficient in free air, depending on environmental conditions and salt crystallization rates. Three new stages of evaporation were defined for saline solutions: SS1, SS2 and SS3. SS1 exhibits a low and gradual decrease in the evaporation rate due to osmotic pressure. During SS2, the evaporation rate falls progressively due to salt precipitation; SS3 is characterized by a constant low evaporation rate and determined by the diffusion rate of water vapor through the precipitated salt layer. Even though phenomenologically similar to the classical evaporation stages of pure water, these stages correspond to different mechanisms and the transition between stages can occur regardless the hydraulic conditions. As well, it was shown that matrix
A study of vapor-liquid flow in porous media
Satik, Cengiz; Yortsos, Yanis C.
1994-01-20
We study the heat transfer-driven liquid-to-vapor phase change in single-component systems in porous media by using pore network models and flow visualization experiments. Experiments using glass micromodels were conducted. The flow visualization allowed us to define the rules for the numerical pore network model. A numerical pore network model is developed for vapor-liquid displacement where fluid flow, heat transfer and capillarity are included at the pore level. We examine the growth process at two different boundary conditions.
Bacteria transport through porous media. Annual report, December 31, 1984
Yen, T.F.
1986-09-01
The following five chapters in this report have been processed separately for inclusion in the Energy Data Base: (1) theoretical model of convective diffusion of motile and non-motile bacteria toward solid surfaces; (2) interfacial electrochemistry of oxide surfaces in oil-bearing sands and sandstones; (3) effects of sodium pyrophosphate additive on the ''huff and puff''/nutrient flooding MEOR process; (4) interaction of Escherichia coli B, B/4, and bacteriophage T4D with Berea sandstone rock in relation to enhanced oil recovery; and (5) transport of bacteria in porous media and its significance in microbial enhanced oil recovery.
Flow visualization in porous media via Positron Emission Tomography
NASA Astrophysics Data System (ADS)
Khalili, A.; Basu, A. J.; Pietrzyk, U.
1998-04-01
We demonstrate here the use of a non-invasive technique based on Positron Emission Tomography (PET) in visualizing and in making quantitative measurements of scalar transport through natural opaque permeable sediments. Along with various other possibilities, this technique has the potential to help improve the understanding of processes that take place at the seabeds between the porewater and the overlying water, which result in exchange of nutrients, toxins and solute. Unlike many other methods, PET is able to produce full three-dimensional pictures of the percolation of fluid through not only "constructed" but also natural porous media.
Porous Media Approach for Modeling Closed Cell Foam
NASA Technical Reports Server (NTRS)
Ghosn, Louis J.; Sullivan, Roy M.
2006-01-01
In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is
A Porous Media Model for Blood Flow within Reticulated Foam
Ortega, J.M.
2013-01-01
A porous media model is developed for non-Newtonian blood flow through reticulated foam at Reynolds numbers ranging from 10−8 to 10. This empirical model effectively divides the pressure gradient versus flow speed curve into three regimes, in which either the non-Newtonian viscous forces, the Newtonian viscous forces, or the inertial fluid forces are most prevalent. When compared to simulation data of blood flow through two reticulated foam geometries, the model adequately captures the pressure gradient within all three regimes, especially that within the Newtonian regime where blood transitions from a power-law to a constant viscosity fluid. PMID:24031095
Focus on the physics of magnetic resonance on porous media
NASA Astrophysics Data System (ADS)
Song, Yi-Qiao
2012-05-01
Porous media are ubiquitous in our environment and their application is extremely broad. The common connection between these diverse materials is the importance of the microstructure in determining the physical, chemical and biological functions and properties. Magnetic resonance and its imaging modality have been essential for noninvasive characterization of these materials in the development of catalysts, understanding cement hydration, fluid transport in rocks and soil, geological prospecting and characterization of tissue properties for medical diagnosis. This focus issue highlights recent NMR/MRI technical development, the underlying physics associated with the new technology, as well as novel applications of the technologies.
Heat Transfer through Porous Media in the Counterflow Regime of He II
Maekawa, R.; Baudouy, B.
2004-06-23
Experimental results are presented for He II counter flow through Al2O3 porous media. Tests have been performed on three porous disks with different thicknesses, 2, 3 and 4 mm, having the same porosity of 32 %, average pore diameter of 2 {mu}m and permeability in the range of 10-14 m2. Temperature and pressure differences were measured across porous media from 1.4 K to 2.1 K in the saturated superfluid helium. The influence on the porous media thickness to the heat transfer is clearly seen on the typical linear Darcy regime and the turbulent Gorter-Mellink regime. Deviation from these classical laws is observed for large temperature difference that can be attributed to the change of He II heat transfer properties due to the complex flow paths of porous media. The effect of porous media thickness to the He II heat transfer is discussed.
NASA Astrophysics Data System (ADS)
Durant, K. A.; Duchateau, C.; Kovscek, A. R.
2012-12-01
Geological sequestration in saline aquifers is a process that reduces carbon dioxide emissions and has become increasingly significant because it may help to mitigate global climate change and enhance energy sustainability. Understanding the flow behavior and stability of partially miscible fluids - such as carbon dioxide and water - in porous media is the main goal of this research. Additionally, prediction of unstable flows in porous media presents an ongoing challenge for engineers and scientists; therefore, this research also has positive implications in areas such as enhanced oil recovery and the development of chemical reactors. In order to improve the knowledge base of this intricate problem, stability and flow characteristics must be examined experimentally in three ways: (i) immiscible flows in homogeneous porous media, (ii) effects of heterogeneity, and (iii) the impact of partial miscibility of the fluids. This study focused specifically on immiscible flows and was accomplished by saturating a 2-D homogeneous, silicon-etched micromodel with UV sensitive dyed water and subsequently injecting heptane to displace the water. Heptane is more viscous than carbon dioxide; therefore various concentrations of glycerin were added to the aqueous phase to alter the viscosity ratio so that the heptane-water displacement in the micromodel was comparable to the carbon dioxide-water displacement in an aquifer. As the heptane moved through the micromodel, a high-speed camera system imaged the gradual displacement changes, seen due to the color changes in the UV dyed water. The extent of fingering in the micromodel is the physical characteristic used to determine stability. It is found that high flow rates and glycerin concentrations cause viscous fingering while low flow rates and glycerin concentrations result in capillary fingering.
Numerical method for computing flow through partially saturated porous media
NASA Astrophysics Data System (ADS)
Eaton, R. R.
This paper discusses the development of the finite element computer code SAGUARO which calculates the two-dimensional flow of mass and energy through porous media. The media may be saturated or partially saturated. SAGUARO solves the parabolic time-dependent mass transport equation which accounts for the presence of partially saturated zones through the use of highly non-linear material characteristic curves. The energy equation accounts for the possibility of partially-saturated regions by adjusting the thermal capacitances and thermal conductivities according to the volume fraction of water present in the local pores. The code capabilities are demonstrated through the presentation of a sample problem involving the one dimensional calculation of simultaneous energy transfer and water infiltration into partially saturated hard rock.
Vorticity and upscaled dispersion in 3D heterogeneous porous media
NASA Astrophysics Data System (ADS)
Di Dato, Mariaines; Chiogna, Gabriele; de Barros, Felipe; Bellin, Alberto; Fiori, Aldo
2015-04-01
Modeling flow in porous media is relevant for many environmental, energy and industrial applications. From an environmental perspective, the relevance of porous media flow becomes evident in subsurface hydrology. In general, flow in natural porous media is creeping, yet the large variability in the hydraulic conductivity values encountered in natural aquifers leads to highly heterogeneous flow fields. This natural variability in the conductivity field will affect both dilution rates of chemical species and reactive mixing. A physical consequence of this heterogeneity is also the presence of a various localized kinematical features such as straining, shearing and vorticity in aquifers, which will influence the shape of solute clouds and its fate and transport. This work aims in fundamentally characterizing the vorticity field in spatially heterogeneous flow fields as a function of their statistical properties in order to analyze the impact on transport processes. In our study, three-dimensional porous formations are constructed with an ensemble of N independent, non-overlapping spheroidal inclusions submerged into an homogeneous matrix, of conductivity K0. The inclusions are randomly located in a domain of volume W and are fully characterized by the geometry of spheroid (oblate or prolate), their conductivity K (random and drawn from a given probability density function fκ), the centroid location ¯x, the axes ratio e, the orientation of the rotational axis (α1,α2) and the volume w. Under the assumption of diluted medium, the flow problem is solved analitically by means of only two parameters: the conductivity contrast κ = K/K0 and the volume fraction n = Nw/W . Through the variation of these parameters of the problem, it is possible to approximate the structure of natural heterogeneous porous media. Using a random distribution of the orientation of the inclusions, we create media defined by the same global anisotropy f = Iz/Ix but different micro
Towards a porous media model of the human lung
NASA Astrophysics Data System (ADS)
DeGroot, Christopher T.; Straatman, Anthony G.
2012-05-01
In this article, progress towards building a complete porous media model of the human lung is discussed. While the recent trend in computational fluid dynamics studies of airflow in the human lung has been to continually increase the size and detail of the airway tree under consideration, it is proposed in this work that simulating flow in the human lung as a coupled fluid-porous system is an effective method to simulate the flow in the whole lung. Under the proposed modeling paradigm, a truncated airway tree constitutes a fluid region which is coupled to a porous region that represents the remainder of the lung volume, containing small airways and alveoli. The first part of this work describes pore-level simulations conducted in an alveolated duct geometry, which are present in large quantities in the human lung, to determine its permeability. Next, volume-averaged simulations incorporating the results of the pore-level simulations and using a realistic lung geometry based on computed tomography images are discussed along with future directions for this work.
Microscale simulation of particle deposition in porous media.
Boccardo, Gianluca; Marchisio, Daniele L; Sethi, Rajandrea
2014-03-01
In this work several geometries, each representing a different porous medium, are considered to perform detailed computational fluid dynamics simulation for fluid flow, particle transport and deposition. Only Brownian motions and steric interception are accounted for as deposition mechanisms. Firstly pressure drop in each porous medium is analyzed in order to determine an effective grain size, by fitting the results with the Ergun law. Then grid independence is assessed. Lastly, particle transport in the system is investigated via Eulerian steady-state simulations, where particle concentration is solved for, not following explicitly particles' trajectories, but solving the corresponding advection-diffusion equation. An assumption was made in considering favorable collector-particle interactions, resulting in a "perfect sink" boundary condition for the collectors. The gathered simulation data are used to calculate the deposition efficiency due to Brownian motions and steric interception. The original Levich law for one simple circular collector is verified; subsequently porous media constituted by a packing of collectors are scrutinized. Results show that the interactions between the different collectors result in behaviors which are not in line with the theory developed by Happel and co-workers, highlighting a different dependency of the deposition efficiency on the dimensionless groups involved in the relevant correlations. PMID:24407681
NMR Studies of Biodegradation Reactions in Porous Media
NASA Astrophysics Data System (ADS)
Mitreiter, I.; Oswald, S. E.; Stallmach, F.
2007-12-01
Subsurface contamination caused by organic compounds is a widespread environmental problem. Biodegradation is the main process to reduce the mass of organic contaminants in soils and groundwater. Often the biodegradation is limited by the supply of electron acceptors at the location of the contaminants, and mixing and diffusion are therefore coupled to the degradation process. Nuclear magnetic resonance (NMR) allows non- invasive and non-destructive insight into diffusion processes and effects of microbial biomass on the pore scale, where the mixing and degradation processes actually are taking place. Studying key processes non-invasively with high temporal resolution will contribute to a better understanding of the effective biodegradation rates of organic contaminants in saturated porous systems. The application of 1H NMR relaxometry and PFG NMR diffusometry provide the possibilities to study diffusive transport at the mixing zones in porous media, where contaminants and electron acceptors show strong concentration gradients due to degradation reactions. Nuclear spin relaxation times of water reflect its interaction with the porous media as well as the biomass. To localize microbial biomass the mobility of the water molecules (relaxation times) can be classified corresponding to bound water in the biomass and bulk water. Oxygen is the most important electron acceptor that stimulate the activity and growth of aerobic microbes, and iron(III) a major one for anaerobes. In our work we applied low field (0.2 T) and high field (3 T) NMR measurements to analyze the diffusion in water and biomass and to monitor the concentration changes of electron acceptors. The experimental setup consists of small columns containing saturated porous media applying inversion recovery (IR) and pulsed field gradient (PFG) sequences. We show that both, oxygen as well as iron(III) affect the relaxation times by their paramagnetic properties, and can be determined by NMR relaxometry in
Mobility of Multi-walled Carbon Nanotubes in Porous Media
NASA Astrophysics Data System (ADS)
O'Carroll, D. M.; Liu, X.; Petersen, E.; Huang, Q.; Anderson, L.
2007-12-01
Engineered multi-walled carbon nanotubes (MWCNTs) are the subject of intense research and are expected to gain widespread usage in a broad variety of commercial products. However concerns have been raised regarding their potential environmental and health risks. The mobility of MWCNTs in porous media is examined in this study through one dimensional flow-through column experiments under conditions representative of subsurface and drinking water treatment systems. The goal of this work was to determine dominant MWCNT removal mechanisms and factors that control MWCNT transport. Results demonstrate that pore water velocity strongly influenced MWCNT transport, a result that stands in contrast to traditional colloid filtration theory, which suggests a relatively minor effect of flow velocity in comparison to Brownian diffusion. Experiments conducted at different ionic strengths indicate that both particle deposition and straining are important MWCNT removal mechanisms from the aqueous phase. Given these findings, traditional colloid filtration theory may not be appropriate for the prediction of MWCNT mobility in porous media. This may be due to the large aspect ratio of the MWCNTs and the importance of straining in MWCNT removal.
Dynamic permeability of porous media by the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Pazdniakou, A.; Adler, P. M.
2013-12-01
The lattice Boltzmann method (LBM) is applied to calculate the dynamic permeability K(ω) of porous media; an oscillating macroscopic pressure gradient is imposed in order to generate oscillating flows. The LBM simulation yields the time dependent seepage velocity of amplitude A and phase shift B which are used to calculate K(ω). The procedure is validated for plane Poiseuille flows where excellent agreement with the analytical solution is obtained. The limitations of the method are discussed. When the ratio between the kinematic viscosity and the characteristic size of the pores is high, the corresponding Knudsen number Kn is high and the numerical values of K(ω) are incorrect with a positive imaginary part; it is only when Kn is small enough that correct values are obtained. The influence of the time discretization of the oscillating body force is studied; simulation results are influenced by an insufficient discretization, i.e., it is necessary to avoid using too high frequencies. The influence of absolute errors in the seepage velocity amplitude δA and the phase shift δB on K(ω) shows that for high ω even small errors in B can cause drastic errors in ReK(ω). The dynamic permeability of reconstructed and real (sandstone) porous media is calculated for a large range of frequencies and the universal scaling behavior is verified. Very good correspondences with the theoretical predictions are observed.
Influence of biofilms on transport properties in porous media
NASA Astrophysics Data System (ADS)
Davit, Y.
2015-12-01
Microbial activity and biofilm growth in porous media can drastically modify transport properties such as permeability, longitudinal and transverse dispersion or effective reaction rates. Understanding these effects has proven to be a considerable challenge. Advances in this field have been hindered by the difficulty of modeling and visualizing these multi-phase non-linear effects across a broad range of spatial and temporal scales. To address these issues, we are developing a strategy that combines imaging techniques based on x-ray micro-tomography with homogenization of pore-scale transport equations. Here, we review recent progress in x-ray imaging of biofilms in porous media, with a particular focus on the contrast agents that are used to differentiate between the fluid and biofilm phases. We further show how the 3D distribution of the different phases can be used to extract specific information about the biofilm and how effective properties can be calculated via the resolution of closure problems. These closure problems are obtained using the method of volume averaging and must be adapted to the problem of interest. In hydrological systems, we show that a generic formulation for reactive solute transport is based on a domain decomposition approach at the micro-scale yielding macro-scale models reminiscent of multi-rate mass transfer approaches.
Evaluation of liquid aerosol transport through porous media.
Hall, R; Murdoch, L; Falta, R; Looney, B; Riha, B
2016-07-01
Application of remediation methods in contaminated vadose zones has been hindered by an inability to effectively distribute liquid- or solid-phase amendments. Injection as aerosols in a carrier gas could be a viable method for achieving useful distributions of amendments in unsaturated materials. The objectives of this work were to characterize radial transport of aerosols in unsaturated porous media, and to develop capabilities for predicting results of aerosol injection scenarios at the field-scale. Transport processes were investigated by conducting lab-scale injection experiments with radial flow geometry, and predictive capabilities were obtained by developing and validating a numerical model for simulating coupled aerosol transport, deposition, and multi-phase flow in porous media. Soybean oil was transported more than 2m through sand by injecting it as micron-scale aerosol droplets. Oil saturation in the sand increased with time to a maximum of 0.25, and decreased with radial distance in the experiments. The numerical analysis predicted the distribution of oil saturation with only minor calibration. The results indicated that evolution of oil saturation was controlled by aerosol deposition and subsequent flow of the liquid oil, and simulation requires including these two coupled processes. The calibrated model was used to evaluate field applications. The results suggest that amendments can be delivered to the vadose zone as aerosols, and that gas injection rate and aerosol particle size will be important controls on the process. PMID:27149690
Multicomponent, multiphase flow in porous media with temperature variation
Wingard, J.S.; Orr, F.M. Jr.
1990-10-01
Recovery of hydrocarbons from porous media is an ongoing concern. Advanced techniques augment conventional recovery methods by injecting fluids that favorably interact with the oil. These fluids interact with the oil by energy transfer, in the case of steam injection, or by mass transfer, as in a miscible gas flood. Often both thermal and compositional considerations are important. An understanding of these injection methods requires knowledge of how temperature variations, phase equilibrium and multiphase flow in porous media interact. The material balance for each component and energy balance are cast as a system of non-strictly hyperbolic partial differential equations. This system of equations is solved using the method of characteristics. The model takes into account the phase behavior by using the Peng-Robinson equation of state to partition the individual components into different phases. Temperature effects are accounted for by the energy balance. Flow effects are modelled by using fractional flow curves and a Stone's three phase relative permeability model. Three problems are discussed. The first problem eliminates the phase behavior aspect of the problem by studying the flow of a single component as it undergoes an isothermal phase change. The second couples the effects of temperature and flow behavior by including a second component that is immiscible with the original component. Phase behavior is added by using a set of three partially miscible components that partition into two or three separate phases. 66 refs., 54 figs., 14 tabs.
Electrokinetic induced solute dispersion in porous media; pore network modeling
NASA Astrophysics Data System (ADS)
Li, Shuai; Schotting, Ruud; Raoof, Amir
2013-04-01
Electrokinetic flow plays an important role in remediation process, separation technique, and chromatography. The solute dispersion is a key parameter to determine transport efficiency. In this study, we present the electrokinetic effects on solute dispersion in porous media at the pore scale, using a pore network model. The analytical solution of the electrokinetic coupling coefficient was obtained to quantity the fluid flow velocity in a cylinder capillary. The effect of electrical double layer on the electrokinetic coupling coefficient was investigated by applying different ionic concentration. By averaging the velocity over cross section within a single pore, the average flux was obtained. Applying such single pore relationships, in the thin electrical double layer limit, to each and every pore within the pore network, potential distribution and the induced fluid flow was calculated for the whole domain. The resulting pore velocities were used to simulate solute transport within the pore network. By averaging the results, we obtained the breakthrough curve (BTC) of the average concentration at the outlet of the pore network. Optimizing the solution of continuum scale advection-dispersion equation to such a BTC, solute dispersion coefficient was estimated. We have compared the dispersion caused by electrokinetic flow and pure pressure driven flow under different Peclet number values. In addition, the effect of microstructure and topological properties of porous media on fluid flow and solute dispersion is presented, mainly based on different pore coordination numbers.
Pore-network modeling of biofilm evolution in porous media.
Ezeuko, C C; Sen, A; Grigoryan, A; Gates, I D
2011-10-01
The influence of bacterial biomass on hydraulic properties of porous media (bioclogging) has been explored as a viable means for optimizing subsurface bioremediation and microbial enhanced oil recovery. In this study, we present a pore network simulator for modeling biofilm evolution in porous media including hydrodynamics and nutrient transport based on coupling of advection transport with Fickian diffusion and a reaction term to account for nutrient consumption. Biofilm has non-zero permeability permitting liquid flow and transport through the biofilm itself. To handle simultaneous mass transfer in both liquid and biofilm in a pore element, a dual-diffusion mass transfer model is introduced. The influence of nutrient limitation on predicted results is explored. Nutrient concentration in the network is affected by diffusion coefficient for nutrient transfer across biofilm (compared to water/water diffusion coefficient) under advection dominated transport, represented by mass transport Péclet number >1. The model correctly predicts a dependence of rate of biomass accumulation on inlet concentration. Poor network connectivity shows a significantly large reduction of permeability, for a small biomass pore volume. PMID:21520022