Disorder-induced stiffness degradation of highly disordered porous materials

NASA Astrophysics Data System (ADS)

Laubie, Hadrien; Monfared, Siavash; Radjaï, Farhang; Pellenq, Roland; Ulm, Franz-Josef

2017-09-01

The effective mechanical behavior of multiphase solid materials is generally modeled by means of homogenization techniques that account for phase volume fractions and elastic moduli without considering the spatial distribution of the different phases. By means of extensive numerical simulations of randomly generated porous materials using the lattice element method, the role of local textural properties on the effective elastic properties of disordered porous materials is investigated and compared with different continuum micromechanics-based models. It is found that the pronounced disorder-induced stiffness degradation originates from stress concentrations around pore clusters in highly disordered porous materials. We identify a single disorder parameter, φsa, which combines a measure of the spatial disorder of pores (the clustering index, sa) with the pore volume fraction (the porosity, φ) to scale the disorder-induced stiffness degradation. Thus, we conclude that the classical continuum micromechanics models with one spherical pore phase, due to their underlying homogeneity assumption fall short of addressing the clustering effect, unless additional texture information is introduced, e.g. in form of the shift of the percolation threshold with disorder, or other functional relations between volume fractions and spatial disorder; as illustrated herein for a differential scheme model representative of a two-phase (solid-pore) composite model material.

Modeling of batch sorber system: kinetic, mechanistic, and thermodynamic modeling

NASA Astrophysics Data System (ADS)

Mishra, Vishal

2017-10-01

The present investigation has dealt with the biosorption of copper and zinc ions on the surface of egg-shell particles in the liquid phase. Various rate models were evaluated to elucidate the kinetics of copper and zinc biosorptions, and the results indicated that the pseudo-second-order model was more appropriate than the pseudo-first-order model. The curve of the initial sorption rate versus the initial concentration of copper and zinc ions also complemented the results of the pseudo-second-order model. Models used for the mechanistic modeling were the intra-particle model of pore diffusion and Bangham's model of film diffusion. The results of the mechanistic modeling together with the values of pore and film diffusivities indicated that the preferential mode of the biosorption of copper and zinc ions on the surface of egg-shell particles in the liquid phase was film diffusion. The results of the intra-particle model showed that the biosorption of the copper and zinc ions was not dominated by the pore diffusion, which was due to macro-pores with open-void spaces present on the surface of egg-shell particles. The thermodynamic modeling reproduced the fact that the sorption of copper and zinc was spontaneous, exothermic with the increased order of the randomness at the solid-liquid interface.

Fabrication of scalable tissue engineering scaffolds with dual-pore microarchitecture by combining 3D printing and particle leaching.

PubMed

Mohanty, Soumyaranjan; Sanger, Kuldeep; Heiskanen, Arto; Trifol, Jon; Szabo, Peter; Dufva, Marin; Emnéus, Jenny; Wolff, Anders

2016-04-01

Limitations in controlling scaffold architecture using traditional fabrication techniques are a problem when constructing engineered tissues/organs. Recently, integration of two pore architectures to generate dual-pore scaffolds with tailored physical properties has attracted wide attention in tissue engineering community. Such scaffolds features primary structured pores which can efficiently enhance nutrient/oxygen supply to the surrounding, in combination with secondary random pores, which give high surface area for cell adhesion and proliferation. Here, we present a new technique to fabricate dual-pore scaffolds for various tissue engineering applications where 3D printing of poly(vinyl alcohol) (PVA) mould is combined with salt leaching process. In this technique the sacrificial PVA mould, determining the structured pore architecture, was filled with salt crystals to define the random pore regions of the scaffold. After crosslinking the casted polymer the combined PVA-salt mould was dissolved in water. The technique has advantages over previously reported ones, such as automated assembly of the sacrificial mould, and precise control over pore architecture/dimensions by 3D printing parameters. In this study, polydimethylsiloxane and biodegradable poly(ϵ-caprolactone) were used for fabrication. However, we show that this technique is also suitable for other biocompatible/biodegradable polymers. Various physical and mechanical properties of the dual-pore scaffolds were compared with control scaffolds with either only structured or only random pores, fabricated using previously reported methods. The fabricated dual-pore scaffolds supported high cell density, due to the random pores, in combination with uniform cell distribution throughout the scaffold, and higher cell proliferation and viability due to efficient nutrient/oxygen transport through the structured pores. In conclusion, the described fabrication technique is rapid, inexpensive, scalable, and compatible with different polymers, making it suitable for engineering various large scale organs/tissues. Copyright © 2015. Published by Elsevier B.V.

Solvable continuous-time random walk model of the motion of tracer particles through porous media.

PubMed

Fouxon, Itzhak; Holzner, Markus

2016-08-01

We consider the continuous-time random walk (CTRW) model of tracer motion in porous medium flows based on the experimentally determined distributions of pore velocity and pore size reported by Holzner et al. [M. Holzner et al., Phys. Rev. E 92, 013015 (2015)PLEEE81539-375510.1103/PhysRevE.92.013015]. The particle's passing through one channel is modeled as one step of the walk. The step (channel) length is random and the walker's velocity at consecutive steps of the walk is conserved with finite probability, mimicking that at the turning point there could be no abrupt change of velocity. We provide the Laplace transform of the characteristic function of the walker's position and reductions for different cases of independence of the CTRW's step duration τ, length l, and velocity v. We solve our model with independent l and v. The model incorporates different forms of the tail of the probability density of small velocities that vary with the model parameter α. Depending on that parameter, all types of anomalous diffusion can hold, from super- to subdiffusion. In a finite interval of α, ballistic behavior with logarithmic corrections holds, which was observed in a previously introduced CTRW model with independent l and τ. Universality of tracer diffusion in the porous medium is considered.

Porous media flux sensitivity to pore-scale geostatistics: A bottom-up approach

NASA Astrophysics Data System (ADS)

Di Palma, P. R.; Guyennon, N.; Heße, F.; Romano, E.

2017-04-01

Macroscopic properties of flow through porous media can be directly computed by solving the Navier-Stokes equations at the scales related to the actual flow processes, while considering the porous structures in an explicit way. The aim of this paper is to investigate the effects of the pore-scale spatial distribution on seepage velocity through numerical simulations of 3D fluid flow performed by the lattice Boltzmann method. To this end, we generate multiple random Gaussian fields whose spatial correlation follows an assigned semi-variogram function. The Exponential and Gaussian semi-variograms are chosen as extreme-cases of correlation for short distances and statistical properties of the resulting porous media (indicator field) are described using the Matèrn covariance model, with characteristic lengths of spatial autocorrelation (pore size) varying from 2% to 13% of the linear domain. To consider the sensitivity of the modeling results to the geostatistical representativeness of the domain as well as to the adopted resolution, porous media have been generated repetitively with re-initialized random seeds and three different resolutions have been tested for each resulting realization. The main difference among results is observed between the two adopted semi-variograms, indicating that the roughness (short distances autocorrelation) is the property mainly affecting the flux. However, computed seepage velocities show additionally a wide variability (about three orders of magnitude) for each semi-variogram model in relation to the assigned correlation length, corresponding to pore sizes. The spatial resolution affects more the results for short correlation lengths (i.e., small pore sizes), resulting in an increasing underestimation of the seepage velocity with the decreasing correlation length. On the other hand, results show an increasing uncertainty as the correlation length approaches the domain size.

A kinetic Monte Carlo approach to study fluid transport in pore networks

NASA Astrophysics Data System (ADS)

Apostolopoulou, M.; Day, R.; Hull, R.; Stamatakis, M.; Striolo, A.

2017-10-01

The mechanism of fluid migration in porous networks continues to attract great interest. Darcy's law (phenomenological continuum theory), which is often used to describe macroscopically fluid flow through a porous material, is thought to fail in nano-channels. Transport through heterogeneous and anisotropic systems, characterized by a broad distribution of pores, occurs via a contribution of different transport mechanisms, all of which need to be accounted for. The situation is likely more complicated when immiscible fluid mixtures are present. To generalize the study of fluid transport through a porous network, we developed a stochastic kinetic Monte Carlo (KMC) model. In our lattice model, the pore network is represented as a set of connected finite volumes (voxels), and transport is simulated as a random walk of molecules, which "hop" from voxel to voxel. We simulated fluid transport along an effectively 1D pore and we compared the results to those expected by solving analytically the diffusion equation. The KMC model was then implemented to quantify the transport of methane through hydrated micropores, in which case atomistic molecular dynamic simulation results were reproduced. The model was then used to study flow through pore networks, where it was able to quantify the effect of the pore length and the effect of the network's connectivity. The results are consistent with experiments but also provide additional physical insights. Extension of the model will be useful to better understand fluid transport in shale rocks.

Porous Electrodes I: Numerical Simulation Using Random Network and Single Pore Models.

DTIC Science & Technology

1984-01-31

characteristic of Zn and ZnO ) and scaling them down to the magnitude of a unit pore size- approximately 10i in diameter. We define a characteristic...supported by the Office of Naval Research under contract N00014-81-K-0339. 11 - 15 - REFERENCES 1. R. de Levis in Advances in Electrochemistry and...P. 3. Hendra Dr. C. E. Mueller Department of Chemistry The Electrochemistry Branch University of Southampton Naval Surface Weapons Center

Characterization of Pore Defects and Fatigue Cracks in Die Cast AM60 Using 3D X-ray Computed Tomography

NASA Astrophysics Data System (ADS)

Yang, Zhuofei; Kang, Jidong; Wilkinson, David S.

2015-08-01

AM60 high pressure die castings have been used in automobile applications to reduce the weight of vehicles. However, the pore defects that are inherent in die casting may negatively affect mechanical properties, especially the fatigue properties. Here we have studied damage ( e.g., pore defects, fatigue cracks) during strained-controlled fatigue using 3-dimensional X-ray computed tomography (XCT). The fatigue test was interrupted every 2000 cycles and the specimen was removed to be scanned using a desktop micro-CT system. XCT reveals pore defects, cracks, and fracture surfaces. The results show that pores can be accurately measured and modeled in 3D. Defect bands are found to be made of pores under 50 µm (based on volume-equivalent sphere diameter). Larger pores are randomly distributed in the region between the defect bands. Observation of fatigue cracks by XCT is performed in three ways such that the 3D model gives the best illustration of crack-porosity interaction while the other two methods, with the cracks being viewed on transverse or longitudinal cross sections, have better detectability on crack initiation and crack tip observation. XCT is also of value in failure analysis on fracture surfaces. By assessing XCT data during fatigue testing and observing fracture surfaces on a 3D model, a better understanding on the crack initiation, crack-porosity interaction, and the morphology of fracture surface is achieved.

Pore-scale simulation of liquid CO2 displacement of water using a two-phase lattice Boltzmann model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Haihu; Valocchi, Albert J.; Werth, Charles J.

A lattice Boltzmann color-fluid model, which was recently proposed by Liu et al. [H. Liu, A.J. Valocchi, and Q. Kang. Three-dimensional lattice Boltzmann model for immiscible two-phase flow simulations. Phys. Rev. E, 85:046309, 2012.] based on a concept of continuum surface force, is improved to simulate immiscible two-phase flows in porous media. The new improvements allow the model to account for different kinematic viscosities of both fluids and to model fluid-solid interactions. The capability and accuracy of this model is first validated by two benchmark tests: a layered two-phase flow with a viscosity ratio, and a dynamic capillary intrusion. Thismore » model is then used to simulate liquid CO2 (LCO2) displacing water in a dual-permeability pore network. The extent and behavior of LCO2 preferential flow (i.e., fingering) is found to depend on the capillary number (Ca), and three different displacement patterns observed in previous micromodel experiments are reproduced. The predicted variation of LCO2 saturation with Ca, as well as variation of specific interfacial length with LCO2 saturation, are both in good agreement with the experimental observations. To understand the effect of heterogeneity on pore-scale displacement, we also simulate LCO2 displacing water in a randomly heterogeneous pore network, which has the same size and porosity as the dual-permeability pore network. In comparison to the dual-permeability case, the transition from capillary fingering to viscous fingering occurs at a higher Ca, and LCO2 saturation is higher at low Ca but lower at high Ca. In either pore network, the LCO2-water specific interfacial length is found to obey a power-law dependence on LCO2 saturation.« less

Prediction of the filtrate particle size distribution from the pore size distribution in membrane filtration: Numerical correlations from computer simulations

NASA Astrophysics Data System (ADS)

Marrufo-Hernández, Norma Alejandra; Hernández-Guerrero, Maribel; Nápoles-Duarte, José Manuel; Palomares-Báez, Juan Pedro; Chávez-Rojo, Marco Antonio

2018-03-01

We present a computational model that describes the diffusion of a hard spheres colloidal fluid through a membrane. The membrane matrix is modeled as a series of flat parallel planes with circular pores of different sizes and random spatial distribution. This model was employed to determine how the size distribution of the colloidal filtrate depends on the size distributions of both, the particles in the feed and the pores of the membrane, as well as to describe the filtration kinetics. A Brownian dynamics simulation study considering normal distributions was developed in order to determine empirical correlations between the parameters that characterize these distributions. The model can also be extended to other distributions such as log-normal. This study could, therefore, facilitate the selection of membranes for industrial or scientific filtration processes once the size distribution of the feed is known and the expected characteristics in the filtrate have been defined.

Correlation effects during liquid infiltration into hydrophobic nanoporous media

DOE Office of Scientific and Technical Information (OSTI.GOV)

Borman, V. D., E-mail: vdborman@mephi.ru; Belogorlov, A. A.; Byrkin, V. A.

To explain the thermal effects observed during the infiltration of a nonwetting liquid into a disordered nanoporous medium, we have constructed a model that includes correlation effects in a disordered medium. It is based on analytical methods of the percolation theory. The infiltration of a porous medium is considered as the infiltration of pores in an infinite cluster of interconnected pores. Using the model of randomly situated spheres (RSS), we have been able to take into account the correlation effect of the spatial arrangement and connectivity of pores in the medium. The other correlation effect of the mutual arrangement ofmore » filled and empty pores on the shell of an infinite percolation cluster of filled pores determines the infiltration fluctuation probability. This probability has been calculated analytically. Allowance for these correlation effects during infiltration and defiltration makes it possible to suggest a physical mechanism of the contact angle hysteresis and to calculate the dependences of the contact angles on the degree of infiltration, porosity of the medium, and temperature. Based on the suggested model, we have managed to describe the temperature dependences of the infiltration and defiltration pressures and the thermal effects that accompany the absorption of energy by disordered porous medium-nonwetting liquid systems with various porosities in a unified way.« less

Simulating Pre-Asymptotic, Non-Fickian Transport Although Doing Simple Random Walks - Supported By Empirical Pore-Scale Velocity Distributions and Memory Effects

NASA Astrophysics Data System (ADS)

Most, S.; Jia, N.; Bijeljic, B.; Nowak, W.

2016-12-01

Pre-asymptotic characteristics are almost ubiquitous when analyzing solute transport processes in porous media. These pre-asymptotic aspects are caused by spatial coherence in the velocity field and by its heterogeneity. For the Lagrangian perspective of particle displacements, the causes of pre-asymptotic, non-Fickian transport are skewed velocity distribution, statistical dependencies between subsequent increments of particle positions (memory) and dependence between the x, y and z-components of particle increments. Valid simulation frameworks should account for these factors. We propose a particle tracking random walk (PTRW) simulation technique that can use empirical pore-space velocity distributions as input, enforces memory between subsequent random walk steps, and considers cross dependence. Thus, it is able to simulate pre-asymptotic non-Fickian transport phenomena. Our PTRW framework contains an advection/dispersion term plus a diffusion term. The advection/dispersion term produces time-series of particle increments from the velocity CDFs. These time series are equipped with memory by enforcing that the CDF values of subsequent velocities change only slightly. The latter is achieved through a random walk on the axis of CDF values between 0 and 1. The virtual diffusion coefficient for that random walk is our only fitting parameter. Cross-dependence can be enforced by constraining the random walk to certain combinations of CDF values between the three velocity components in x, y and z. We will show that this modelling framework is capable of simulating non-Fickian transport by comparison with a pore-scale transport simulation and we analyze the approach to asymptotic behavior.

Modeling and simulation of Cu diffusion and drift in porous CMOS backend dielectrics

NASA Astrophysics Data System (ADS)

Ali, R.; Fan, Y.; King, S.; Orlowski, M.

2018-06-01

With the advent of porous dielectrics, Cu drift-diffusion reliability issues in CMOS backend have only been exacerbated. In this regard, a modeling and simulation study of Cu atom/ion drift-diffusion in porous dielectrics is presented to assess the backend reliability and to explore conditions for a reliable Resistive Random Access Memory (RRAM) operation. The numerical computation, using elementary jump frequencies for a random walk in 2D and 3D, is based on an extended adjacency tensor concept. It is shown that Cu diffusion and drift transport are affected as much by the level of porosity as by the pore morphology. Allowance is made for different rates of Cu dissolution into the dielectric and for Cu absorption and transport at and on the inner walls of the pores. Most of the complex phenomena of the drift-diffusion transport in porous media can be understood in terms of local lateral and vertical gradients and the degree of their perturbation caused by the presence of pores in the transport domain. The impact of pore morphology, related to the concept of tortuosity, is discussed in terms of "channeling" and "trapping" effects. The simulations are calibrated to experimental results of porous SiCOH layers of 25 nm thickness, sandwiched between Cu and Pt(W) electrodes with experimental porosity levels of 0%, 8%, 12%, and 25%. We find that porous SICOH is more immune to Cu+ drift at 300 K than non-porous SICOH.

Application of real rock pore-threat statistics to a regular pore network model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rakibul, M.; Sarker, H.; McIntyre, D.

2011-01-01

This work reports the application of real rock statistical data to a previously developed regular pore network model in an attempt to produce an accurate simulation tool with low computational overhead. A core plug from the St. Peter Sandstone formation in Indiana was scanned with a high resolution micro CT scanner. The pore-throat statistics of the three-dimensional reconstructed rock were extracted and the distribution of the pore-throat sizes was applied to the regular pore network model. In order to keep the equivalent model regular, only the throat area or the throat radius was varied. Ten realizations of randomly distributed throatmore » sizes were generated to simulate the drainage process and relative permeability was calculated and compared with the experimentally determined values of the original rock sample. The numerical and experimental procedures are explained in detail and the performance of the model in relation to the experimental data is discussed and analyzed. Petrophysical properties such as relative permeability are important in many applied fields such as production of petroleum fluids, enhanced oil recovery, carbon dioxide sequestration, ground water flow, etc. Relative permeability data are used for a wide range of conventional reservoir engineering calculations and in numerical reservoir simulation. Two-phase oil water relative permeability data are generated on the same core plug from both pore network model and experimental procedure. The shape and size of the relative permeability curves were compared and analyzed and good match has been observed for wetting phase relative permeability but for non-wetting phase, simulation results were found to be deviated from the experimental ones. Efforts to determine petrophysical properties of rocks using numerical techniques are to eliminate the necessity of regular core analysis, which can be time consuming and expensive. So a numerical technique is expected to be fast and to produce reliable results. In applied engineering, sometimes quick result with reasonable accuracy is acceptable than the more time consuming results. Present work is an effort to check the accuracy and validity of a previously developed pore network model for obtaining important petrophysical properties of rocks based on cutting-sized sample data.« less

Application of real rock pore-throat statistics to a regular pore network model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sarker, M.R.; McIntyre, D.; Ferer, M.

2011-01-01

This work reports the application of real rock statistical data to a previously developed regular pore network model in an attempt to produce an accurate simulation tool with low computational overhead. A core plug from the St. Peter Sandstone formation in Indiana was scanned with a high resolution micro CT scanner. The pore-throat statistics of the three-dimensional reconstructed rock were extracted and the distribution of the pore-throat sizes was applied to the regular pore network model. In order to keep the equivalent model regular, only the throat area or the throat radius was varied. Ten realizations of randomly distributed throatmore » sizes were generated to simulate the drainage process and relative permeability was calculated and compared with the experimentally determined values of the original rock sample. The numerical and experimental procedures are explained in detail and the performance of the model in relation to the experimental data is discussed and analyzed. Petrophysical properties such as relative permeability are important in many applied fields such as production of petroleum fluids, enhanced oil recovery, carbon dioxide sequestration, ground water flow, etc. Relative permeability data are used for a wide range of conventional reservoir engineering calculations and in numerical reservoir simulation. Two-phase oil water relative permeability data are generated on the same core plug from both pore network model and experimental procedure. The shape and size of the relative permeability curves were compared and analyzed and good match has been observed for wetting phase relative permeability but for non-wetting phase, simulation results were found to be deviated from the experimental ones. Efforts to determine petrophysical properties of rocks using numerical techniques are to eliminate the necessity of regular core analysis, which can be time consuming and expensive. So a numerical technique is expected to be fast and to produce reliable results. In applied engineering, sometimes quick result with reasonable accuracy is acceptable than the more time consuming results. Present work is an effort to check the accuracy and validity of a previously developed pore network model for obtaining important petrophysical properties of rocks based on cutting-sized sample data. Introduction« less

Confinement-Driven Phase Separation of Quantum Liquid Mixtures

NASA Astrophysics Data System (ADS)

Prisk, T. R.; Pantalei, C.; Kaiser, H.; Sokol, P. E.

2012-08-01

We report small-angle neutron scattering studies of liquid helium mixtures confined in Mobil Crystalline Material-41 (MCM-41), a porous silica glass with narrow cylindrical nanopores (d=3.4nm). MCM-41 is an ideal model adsorbent for fundamental studies of gas sorption in porous media because its monodisperse pores are arranged in a 2D triangular lattice. The small-angle scattering consists of a series of diffraction peaks whose intensities are determined by how the imbibed liquid fills the pores. Pure He4 adsorbed in the pores show classic, layer-by-layer film growth as a function of pore filling, leaving the long range symmetry of the system intact. In contrast, the adsorption of He3-He4 mixtures produces a structure incommensurate with the pore lattice. Neither capillary condensation nor preferential adsorption of one helium isotope to the pore walls can provide the symmetry-breaking mechanism. The scattering is consistent with the formation of randomly distributed liquid-liquid microdomains ˜2.3nm in size, providing evidence that confinement in a nanometer scale capillary can drive local phase separation in quantum liquid mixtures.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sadovskii, V. M., E-mail: sadov@icm.krasn.ru; Sadovskaya, O. V., E-mail: o-sadov@icm.krasn.ru

Based on the generalized rheological method, the mathematical model describing small deformations of a single-phase porous medium without regard to the effects of a fluid or gas in pores is constructed. The change in resistance of a material to the external mechanical impacts at the moment of pore collapse is taken into account by means of the von Mises–Schleicher strength condition. In order to consider irreversible deformations, alongside with the classical yield conditions by von Mises and Tresca– Saint-Venant, the special condition modeling the plastic loss of stability of a porous skeleton is used. The random nature of the poremore » size distribution is taken into account. It is shown that the proposed mathematical model satisfies the principles of thermodynamics of irreversible processes. Phenomenological parameters of the model are determined on the basis of the approximate calculation of the problem on quasi-static loading of a cubic periodicity cell with spherical voids. In the framework of the obtained model, the process of propagation of plane longitudinal waves of the compression in a homogenous porous medium, accompanied by the plastic deformation of a skeleton and the collapse of pores, is analyzed.« less

One Hundred Ways to be Non-Fickian - A Rigorous Multi-Variate Statistical Analysis of Pore-Scale Transport

NASA Astrophysics Data System (ADS)

Most, Sebastian; Nowak, Wolfgang; Bijeljic, Branko

2015-04-01

Fickian transport in groundwater flow is the exception rather than the rule. Transport in porous media is frequently simulated via particle methods (i.e. particle tracking random walk (PTRW) or continuous time random walk (CTRW)). These methods formulate transport as a stochastic process of particle position increments. At the pore scale, geometry and micro-heterogeneities prohibit the commonly made assumption of independent and normally distributed increments to represent dispersion. Many recent particle methods seek to loosen this assumption. Hence, it is important to get a better understanding of the processes at pore scale. For our analysis we track the positions of 10.000 particles migrating through the pore space over time. The data we use come from micro CT scans of a homogeneous sandstone and encompass about 10 grain sizes. Based on those images we discretize the pore structure and simulate flow at the pore scale based on the Navier-Stokes equation. This flow field realistically describes flow inside the pore space and we do not need to add artificial dispersion during the transport simulation. Next, we use particle tracking random walk and simulate pore-scale transport. Finally, we use the obtained particle trajectories to do a multivariate statistical analysis of the particle motion at the pore scale. Our analysis is based on copulas. Every multivariate joint distribution is a combination of its univariate marginal distributions. The copula represents the dependence structure of those univariate marginals and is therefore useful to observe correlation and non-Gaussian interactions (i.e. non-Fickian transport). The first goal of this analysis is to better understand the validity regions of commonly made assumptions. We are investigating three different transport distances: 1) The distance where the statistical dependence between particle increments can be modelled as an order-one Markov process. This would be the Markovian distance for the process, where the validity of yet-unexplored non-Gaussian-but-Markovian random walks start. 2) The distance where bivariate statistical dependence simplifies to a multi-Gaussian dependence based on simple linear correlation (validity of correlated PTRW/CTRW). 3) The distance of complete statistical independence (validity of classical PTRW/CTRW). The second objective is to reveal characteristic dependencies influencing transport the most. Those dependencies can be very complex. Copulas are highly capable of representing linear dependence as well as non-linear dependence. With that tool we are able to detect persistent characteristics dominating transport even across different scales. The results derived from our experimental data set suggest that there are many more non-Fickian aspects of pore-scale transport than the univariate statistics of longitudinal displacements. Non-Fickianity can also be found in transverse displacements, and in the relations between increments at different time steps. Also, the found dependence is non-linear (i.e. beyond simple correlation) and persists over long distances. Thus, our results strongly support the further refinement of techniques like correlated PTRW or correlated CTRW towards non-linear statistical relations.

An effective medium approach to modelling the pressure-dependent electrical properties of porous rocks

NASA Astrophysics Data System (ADS)

Han, Tongcheng

2018-07-01

Understanding the electrical properties of rocks under varying pressure is important for a variety of geophysical applications. This study proposes an approach to modelling the pressure-dependent electrical properties of porous rocks based on an effective medium model. The so-named Textural model uses the aspect ratios and pressure-dependent volume fractions of the pores and the aspect ratio and electrical conductivity of the matrix grains. The pores were represented by randomly oriented stiff and compliant spheroidal shapes with constant aspect ratios, and their pressure-dependent volume fractions were inverted from the measured variation of total porosity with differential pressure using a dual porosity model. The unknown constant stiff and compliant pore aspect ratios and the aspect ratio and electrical conductivity of the matrix grains were inverted by best fitting the modelled electrical formation factor to the measured data. Application of the approach to three sandstone samples covering a broad porosity range showed that the pressure-dependent electrical properties can be satisfactorily modelled by the proposed approach. The results demonstrate that the dual porosity concept is sufficient to explain the electrical properties of porous rocks under pressure through the effective medium model scheme.

Topological characteristics underpin intermittency and anomalous transport behavior in soil-like porous media

NASA Astrophysics Data System (ADS)

Holzner, M.; Morales, V.; Willmann, M.; Jerjen, I.; Kaufmann, R.; Dentz, M.

2016-12-01

Continuum models of porous media are based on the validity of the Darcy equation for fluid and Fick's law for scalar fluxes on a representative elementary volume. Fluctuations of pore-scale flow and scalar transport are averaged out and represented in terms of effective parameters such as hydrodynamic dispersion. However, the intermittent behavior of pore-scale flow impacts on the nature of particle and scalar transport, and it determines the way dissolved substances mix and react. The understanding of the origin of these processes is of both fundamental and practical importance in applications ranging from reactive transport in groundwater flow to diffusion in fuel cells or biological systems. A central issue in porous medium flow is therefore to relate intermittent behavior of Lagrangian velocity at pore scale imposed by the complex pore network geometry to transport properties at larger scales. Lagrangian measurements in porous systems are nonetheless scarce and most experimental techniques do not provide access to all three velocity components. In this contribution we report 3D measurements of Lagrangian velocity in soil-like porous media. We complement these measurements with detailed X-ray scans of the pore network. We find sharp velocity transitions close to pore throats, and low flow variability in the pore bodies, which gives rise to stretched exponential Lagrangian velocity and acceleration distributions characterized by a sharp peak at low velocity and a superlinear evolution of particle dispersion. We demonstrate that porosity and pore size distribution alone cannot explain the observed features of the flow. Rather, anomalous transport is better interpreted in terms of how pores of various geometries are interconnected. We reproduce the main observations using a continuous-time random walk (CTRW) model revealing the main features that control the system and showing the potential of this simple model to capture transport in complex geometries.

Transport upscaling from pore- to Darcy-scale: Incorporating pore-scale Berea sandstone Lagrangian velocity statistics into a Darcy-scale transport CTRW model

NASA Astrophysics Data System (ADS)

Puyguiraud, Alexandre; Dentz, Marco; Gouze, Philippe

2017-04-01

For the past several years a lot of attention has been given to pore-scale flow in order to understand and model transport, mixing and reaction in porous media. Nevertheless we believe that an accurate study of spatial and temporal evolution of velocities could bring important additional information for the upscaling from pore to higher scales. To gather these pieces of information, we perform Stokes flow simulations on pore-scale digitized images of a Berea sandstone core. First, micro-tomography (XRMT) imaging and segmentation processes allow us to obtain 3D black and white images of the sample [1]. Then we used an OpenFoam solver to perform the Stokes flow simulations mentioned above, which gives us the velocities at the interfaces of a cubic mesh. Subsequently, we use a particle streamline reconstruction technique which uses the Eulerian velocity field previously obtained. This technique, based on a modified Pollock algorithm [2], enables us to make particle tracking simulations on the digitized sample. In order to build a stochastic pore-scale transport model, we analyze the Lagrangian velocity series in two different ways. First we investigate the velocity evolution by sampling isochronically (t-Lagrangian), and by studying its statistical properties in terms of one- and two-points statistics. Intermittent patterns can be observed. These are due to the persistance of low velocities over a characteristic space length. Other results are investigated, such as correlation functions and velocity PDFs, which permit us to study more deeply this persistence in the velocities and to compute the correlation times. However, with the second approach, doing these same analysis in space by computing the velocities equidistantly, enables us to remove the intermittency shown in the temporal evolution and to model these velocity series as a Markov process. This renders the stochastic particle dynamics into a CTRW [3]. [1] Gjetvaj, F., A. Russian, P. Gouze, and M. Dentz (2015), Dual control of flow field heterogeneity and immobile porosity on non-Fickian transport in Berea sandstone, Water Resour. Res., 51, 8273-8293, doi:10.1002/2015WR017645. [2] Mostaghimi, P., Bijeljic, B., Blunt, M. (2012). Simulation of Flow and Dispersion on Pore-Space Images. Society of Petroleum Engineers. doi:10.2118/135261-PA. [3] Dentz, M., P. K. Kang, A. Comolli, T. Le Borgne, and D. R. Lester, Continuous time random walks for the evolution of Lagrangian velocities, Phys. Rev. Fluids, 2016. Keywords: Porescale, particle tracking, transport, Lagrangian velocity, ergodicity, Markovianity, continuous time random walks, upscaling.

Pore Pressure and Stress Distributions Around a Hydraulic Fracture in Heterogeneous Rock

NASA Astrophysics Data System (ADS)

Gao, Qian; Ghassemi, Ahmad

2017-12-01

One of the most significant characteristics of unconventional petroleum bearing formations is their heterogeneity, which affects the stress distribution, hydraulic fracture propagation and also fluid flow. This study focuses on the stress and pore pressure redistributions during hydraulic stimulation in a heterogeneous poroelastic rock. Lognormal random distributions of Young's modulus and permeability are generated to simulate the heterogeneous distributions of material properties. A 3D fully coupled poroelastic model based on the finite element method is presented utilizing a displacement-pressure formulation. In order to verify the model, numerical results are compared with analytical solutions showing excellent agreements. The effects of heterogeneities on stress and pore pressure distributions around a penny-shaped fracture in poroelastic rock are then analyzed. Results indicate that the stress and pore pressure distributions are more complex in a heterogeneous reservoir than in a homogeneous one. The spatial extent of stress reorientation during hydraulic stimulations is a function of time and is continuously changing due to the diffusion of pore pressure in the heterogeneous system. In contrast to the stress distributions in homogeneous media, irregular distributions of stresses and pore pressure are observed. Due to the change of material properties, shear stresses and nonuniform deformations are generated. The induced shear stresses in heterogeneous rock cause the initial horizontal principal stresses to rotate out of horizontal planes.

A phase-field approach to model multi-axial and microstructure dependent fracture in nuclear grade graphite

DOE PAGES

Chakraborty, Pritam; Sabharwall, Piyush; Carroll, Mark C.

2016-04-07

The fracture behavior of nuclear grade graphites is strongly influenced by underlying microstructural features such as the character of filler particles, and the distribution of pores and voids. These microstructural features influence the crack nucleation and propagation behavior, resulting in quasi-brittle fracture with a tortuous crack path and significant scatter in measured bulk strength. This paper uses a phase-field method to model the microstructural and multi-axial fracture in H-451, a historic variant of nuclear graphite that provides the basis for an idealized study on a legacy grade. The representative volume elements are constructed from randomly located pores with random sizemore » obtained from experimentally determined log-normal distribution. The representative volume elements are then subjected to simulated multi-axial loading, and a reasonable agreement of the resulting fracture stress with experiments is obtained. Finally, quasi-brittle stress-strain evolution with a tortuous crack path is also observed from the simulations and is consistent with experimental results.« less

Interaction and its induced inhibiting or synergistic effects during co-gasification of coal char and biomass char.

PubMed

Ding, Liang; Zhang, Yongqi; Wang, Zhiqing; Huang, Jiejie; Fang, Yitian

2014-12-01

Co-gasification of coal char and biomass char was conducted to investigate the interactions between them. And random pore model (RPM) and modified random pore model (MRPM) were applied to describe the gasification behaviors of the samples. The results show that inhibiting effect was observed during co-gasification of corn stalk char with Hulunbeier lignite coal char, while synergistic effects were observed during co-gasification of corn stalk char with Shenmu bituminous coal char and Jincheng anthracite coal char. The inhibiting effect was attributed to the intimate contact and comparable gasification rate between biomass char and coal char, and the loss of the active form of potassium caused by the formation of KAlSiO4, which was proved to be inactive during gasification. While the synergistic effect was caused by the high potassium content of biomass char and the significant difference of reaction rate between coal char and biomass char during gasification. Copyright © 2014 Elsevier Ltd. All rights reserved.

Characterization of porosity via secondary reactions. Final technical report, 1 September 1991--30 November 1995

DOE Office of Scientific and Technical Information (OSTI.GOV)

Calo, J.M.; Zhang, L.; Hall, P.J.

1997-09-01

A new approach to the study of porosity and porosity development in coal chars during gasification was investigated. This approach involves the establishment of the relationships between the amount and type of surface complexes evolved during post-activation temperature programmed desorption (TPD), and the porosity, as measured by gas adsorption and small angle neutron scattering (SANS) techniques. With this new method, the total surface area and micropore volume can be determined by the interpretation of post-activation TPD spectra. The primary conclusion of this work is that it is possible to predict total surface area and micropore volume from TPD spectra. Frommore » the extended random pore model, additional information about the micropore surface area, the nonmicroporous surface area, and the mean micropore size development as a function of reaction time (or burn-off) can also be predicted. Therefore, combining the TPD technique and the extended random pore model provides a new method for the characterization of char porosity.« less

Asymptotic limits of some models for sound propagation in porous media and the assignment of the pore characteristic lengths.

PubMed

Horoshenkov, Kirill V; Groby, Jean-Philippe; Dazel, Olivier

2016-05-01

Modeling of sound propagation in porous media requires the knowledge of several intrinsic material parameters, some of which are difficult or impossible to measure directly, particularly in the case of a porous medium which is composed of pores with a wide range of scales and random interconnections. Four particular parameters which are rarely measured non-acoustically, but used extensively in a number of acoustical models, are the viscous and thermal characteristic lengths, thermal permeability, and Pride parameter. The main purpose of this work is to show how these parameters relate to the pore size distribution which is a routine characteristic measured non-acoustically. This is achieved through the analysis of the asymptotic behavior of four analytical models which have been developed previously to predict the dynamic density and/or compressibility of the equivalent fluid in a porous medium. In this work the models proposed by Johnson, Koplik, and Dashn [J. Fluid Mech. 176, 379-402 (1987)], Champoux and Allard [J. Appl. Phys. 70(4), 1975-1979 (1991)], Pride, Morgan, and Gangi [Phys. Rev. B 47, 4964-4978 (1993)], and Horoshenkov, Attenborough, and Chandler-Wilde [J. Acoust. Soc. Am. 104, 1198-1209 (1998)] are compared. The findings are then used to compare the behavior of the complex dynamic density and compressibility of the fluid in a material pore with uniform and variable cross-sections.

Intermittent Lagrangian velocities and accelerations in three-dimensional porous medium flow.

PubMed

Holzner, M; Morales, V L; Willmann, M; Dentz, M

2015-07-01

Intermittency of Lagrangian velocity and acceleration is a key to understanding transport in complex systems ranging from fluid turbulence to flow in porous media. High-resolution optical particle tracking in a three-dimensional (3D) porous medium provides detailed 3D information on Lagrangian velocities and accelerations. We find sharp transitions close to pore throats, and low flow variability in the pore bodies, which gives rise to stretched exponential Lagrangian velocity and acceleration distributions characterized by a sharp peak at low velocity, superlinear evolution of particle dispersion, and double-peak behavior in the propagators. The velocity distribution is quantified in terms of pore geometry and flow connectivity, which forms the basis for a continuous-time random-walk model that sheds light on the observed Lagrangian flow and transport behaviors.

Porous medium acoustics of wave-induced vorticity diffusion

NASA Astrophysics Data System (ADS)

Müller, T. M.; Sahay, P. N.

2011-02-01

A theory for attenuation and dispersion of elastic waves due to wave-induced generation of vorticity at pore-scale heterogeneities in a macroscopically homogeneous porous medium is developed. The diffusive part of the vorticity field associated with a viscous wave in the pore space—the so-called slow shear wave—is linked to the porous medium acoustics through incorporation of the fluid strain rate tensor of a Newtonian fluid in the poroelastic constitutive relations. The method of statistical smoothing is then used to derive dynamic-equivalent elastic wave velocities accounting for the conversion scattering process into the diffusive slow shear wave in the presence of randomly distributed pore-scale heterogeneities. The result is a simple model for wave attenuation and dispersion associated with the transition from viscosity- to inertia-dominated flow regime.

Forecasting the brittle failure of heterogeneous, porous geomaterials

NASA Astrophysics Data System (ADS)

Vasseur, Jérémie; Wadsworth, Fabian; Heap, Michael; Main, Ian; Lavallée, Yan; Dingwell, Donald

2017-04-01

Heterogeneity develops in magmas during ascent and is dominated by the development of crystal and importantly, bubble populations or pore-network clusters which grow, interact, localize, coalesce, outgas and resorb. Pore-scale heterogeneity is also ubiquitous in sedimentary basin fill during diagenesis. As a first step, we construct numerical simulations in 3D in which randomly generated heterogeneous and polydisperse spheres are placed in volumes and which are permitted to overlap with one another, designed to represent the random growth and interaction of bubbles in a liquid volume. We use these simulated geometries to show that statistical predictions of the inter-bubble lengthscales and evolving bubble surface area or cluster densities can be made based on fundamental percolation theory. As a second step, we take a range of well constrained random heterogeneous rock samples including sandstones, andesites, synthetic partially sintered glass bead samples, and intact glass samples and subject them to a variety of stress loading conditions at a range of temperatures until failure. We record in real time the evolution of the number of acoustic events that precede failure and show that in all scenarios, the acoustic event rate accelerates toward failure, consistent with previous findings. Applying tools designed to forecast the failure time based on these precursory signals, we constrain the absolute error on the forecast time. We find that for all sample types, the error associated with an accurate forecast of failure scales non-linearly with the lengthscale between the pore clusters in the material. Moreover, using a simple micromechanical model for the deformation of porous elastic bodies, we show that the ratio between the equilibrium sub-critical crack length emanating from the pore clusters relative to the inter-pore lengthscale, provides a scaling for the error on forecast accuracy. Thus for the first time we provide a potential quantitative correction for forecasting the failure of porous brittle solids that build the Earth's crust.

A Binomial Modeling Approach for Upscaling Colloid Transport Under Unfavorable Attachment Conditions: Emergent Prediction of Nonmonotonic Retention Profiles

NASA Astrophysics Data System (ADS)

Hilpert, Markus; Johnson, William P.

2018-01-01

We used a recently developed simple mathematical network model to upscale pore-scale colloid transport information determined under unfavorable attachment conditions. Classical log-linear and nonmonotonic retention profiles, both well-reported under favorable and unfavorable attachment conditions, respectively, emerged from our upscaling. The primary attribute of the network is colloid transfer between bulk pore fluid, the near-surface fluid domain (NSFD), and attachment (treated as irreversible). The network model accounts for colloid transfer to the NSFD of downgradient grains and for reentrainment to bulk pore fluid via diffusion or via expulsion at rear flow stagnation zones (RFSZs). The model describes colloid transport by a sequence of random trials in a one-dimensional (1-D) network of Happel cells, which contain a grain and a pore. Using combinatorial analysis that capitalizes on the binomial coefficient, we derived from the pore-scale information the theoretical residence time distribution of colloids in the network. The transition from log-linear to nonmonotonic retention profiles occurs when the conditions underlying classical filtration theory are not fulfilled, i.e., when an NSFD colloid population is maintained. Then, nonmonotonic retention profiles result potentially both for attached and NSFD colloids. The concentration maxima shift downgradient depending on specific parameter choice. The concentration maxima were also shown to shift downgradient temporally (with continued elution) under conditions where attachment is negligible, explaining experimentally observed downgradient transport of retained concentration maxima of adhesion-deficient bacteria. For the case of zero reentrainment, we develop closed-form, analytical expressions for the shape, and the maximum of the colloid retention profile.

X-ray microanalysis of porous materials using Monte Carlo simulations.

PubMed

Poirier, Dominique; Gauvin, Raynald

2011-01-01

Quantitative X-ray microanalysis models, such as ZAF or φ(ρz) methods, are normally based on solid, flat-polished specimens. This limits their use in various domains where porous materials are studied, such as powder metallurgy, catalysts, foams, etc. Previous experimental studies have shown that an increase in porosity leads to a deficit in X-ray emission for various materials, such as graphite, Cr(2) O(3) , CuO, ZnS (Ichinokawa et al., '69), Al(2) O(3) , and Ag (Lakis et al., '92). However, the mechanisms responsible for this decrease are unclear. The porosity by itself does not explain the loss in intensity, other mechanisms have therefore been proposed, such as extra energy loss by the diffusion of electrons by surface plasmons generated at the pores-solid interfaces, surface roughness, extra charging at the pores-solid interface, or carbon diffusion in the pores. However, the exact mechanism is still unclear. In order to better understand the effects of porosity on quantitative microanalysis, a new approach using Monte Carlo simulations was developed by Gauvin (2005) using a constant pore size. In this new study, the X-ray emissions model was modified to include a random log normal distribution of pores size in the simulated materials. This article presents, after a literature review of the previous works performed about X-ray microanalysis of porous materials, some of the results obtained with Gauvin's modified model. They are then compared with experimental results. Copyright © 2011 Wiley Periodicals, Inc.

Microbial dispersal in unsaturated porous media: Characteristics of motile bacterial cell motions in unsaturated angular pore networks

NASA Astrophysics Data System (ADS)

Ebrahimi, Ali N.; Or, Dani

2014-09-01

The dispersal rates of self-propelled microorganisms affect their spatial interactions and the ecological functioning of microbial communities. Microbial dispersal rates affect risk of contamination of water resources by soil-borne pathogens, the inoculation of plant roots, or the rates of spoilage of food products. In contrast with the wealth of information on microbial dispersal in water replete systems, very little is known about their dispersal rates in unsaturated porous media. The fragmented aqueous phase occupying complex soil pore spaces suppress motility and limits dispersal ranges in unsaturated soil. The primary objective of this study was to systematically evaluate key factors that shape microbial dispersal in model unsaturated porous media to quantify effects of saturation, pore space geometry, and chemotaxis on characteristics of principles that govern motile microbial dispersion in unsaturated soil. We constructed a novel 3-D angular pore network model (PNM) to mimic aqueous pathways in soil for different hydration conditions; within the PNM, we employed an individual-based model that considers physiological and biophysical properties of motile and chemotactic bacteria. The effects of hydration conditions on first passage times in different pore networks were studied showing that fragmentation of aquatic habitats under dry conditions sharply suppresses nutrient transport and microbial dispersal rates in good agreement with limited experimental data. Chemotactically biased mean travel speed of microbial cells across 9 mm saturated PNM was ˜3 mm/h decreasing exponentially to 0.45 mm/h for the PNM at matric potential of -15 kPa (for -35 kPa, dispersal practically ceases and the mean travel time to traverse the 9 mm PNM exceeds 1 year). Results indicate that chemotaxis enhances dispersal rates by orders of magnitude relative to random (diffusive) motions. Model predictions considering microbial cell sizes relative to available liquid pathways sizes were in good agreement with experimental results for unsaturated soils. The new modeling platform enables quantitative consideration of key biophysical factors (e.g., pore space heterogeneities and hydration conditions) governing microbial interactions in 3-D soil pore spaces.

Gravitational Effects on Closed-Cellular-Foam Microstructure

NASA Technical Reports Server (NTRS)

Noever, David A.; Cronise, Raymond J.; Wessling, Francis C.; McMannus, Samuel P.; Mathews, John; Patel, Darayas

1996-01-01

Polyurethane foam has been produced in low gravity for the first time. The cause and distribution of different void or pore sizes are elucidated from direct comparison of unit-gravity and low-gravity samples. Low gravity is found to increase the pore roundness by 17% and reduce the void size by 50%. The standard deviation for pores becomes narrower (a more homogeneous foam is produced) in low gravity. Both a Gaussian and a Weibull model fail to describe the statistical distribution of void areas, and hence the governing dynamics do not combine small voids in either a uniform or a dependent fashion to make larger voids. Instead, the void areas follow an exponential law, which effectively randomizes the production of void sizes in a nondependent fashion consistent more with single nucleation than with multiple or combining events.

Percolation Pore Network Study on the Residue Gas Saturation of Dry Reservoir Rocks

NASA Astrophysics Data System (ADS)

Cheng, T.; Tang, Y. B.; Zou, G. Y.; Jiang, K.; Li, M.

2014-12-01

We tried to model the effect of pore size heterogeneity and pore connectivity on the residue gas saturation for dry gas reservoir rocks. If we consider that snap-off does not exist and only piston displacement takes place in all pores with the same size during imbibition process, in the extreme case, the residue gas saturation will be equal to zero. Thus we can suppose that the residue gas saturation of dry rocks is mainly controlled by the pore size distribution. To verify the assumption, percolation pore networks (i.e., three-dimensional simple cubic (SC) and body-center cubic (BCC)) were used in the study. The connectivity and the pore size distribution in percolation pore network could be changed randomly. The concept of water phase connectivity zw(i.e., water coordination number) and gas phase connectivity zg (i.e., gas coordination number) was introduced here. zw and zg will change during simulation and can be estimated numerically from the results of simulations through gradually saturated networks by water. The Simulation results show that when zg less than or equal to 1.5 during water quasi - static imbibition, the gas will be trapped in rock pores. Network simulation results also shows that the residue gas saturation Srg follows a power law relationship (i.e.,Srg∝σrα, where σr is normalized standard deviation of the pore radius distribution, and exponent α is a function of coordination number). This indicates that the residue gas saturation has no explicit relationship with porosity and permeability as it should have in light of previous study, pore radius distribution is the principal factor in determining the residue gas saturation of dry reservoir rocks.

Characterization and Upscaling of Pore Scale Hydrodynamic Mass Transfer

NASA Astrophysics Data System (ADS)

Gouze, P.; Roubinet, D.; Dentz, M.; Planes, V.; Russian, A.

2017-12-01

Imaging reservoir rocks in 3D using X-ray microtomography with spatial resolution ranging from about 1 to 10 mm provides us a unique opportunity not only to characterize pore space geometry but also for simulating hydrodynamical processes. Yet, pores and throats displaying sizes smaller than the resolution cannot be distinguished on the images and must be assigned to a so called microporous phase during the process of image segmentation. Accordingly one simulated mass transfers caused by advection and diffusion in the connected pores (mobile domain) and diffusion in the microporous clusters (immobile domain) using Time Domain Random Walk (TDRW) and developed a set of metrics that can be used to monitor the different mechanisms of transport in the sample, the final objective being of proposing a simple but accurate upscaled 1D model in which the particle travel times in the mobile and immobile domain and the number of mobile-immobile transfer events (called trapping events) are independently distributed random variables characterized by PDFs. For TDRW the solute concentration is represented by the density distribution of non-interacting point-like solute particles which move due to advection and dispersion. The set of metrics derives from different spatial and temporal statistical analyses of the particle motion, and is used for characterizing the particles transport (i) in the mobile domain in relation with the velocity field properties, (ii) in the immobile domain in relation with the structure and the properties of microporous phase and at the mobile-immobile interface. We specifically focused on how to model the trapping frequency and rate into the immobile domain in relation with the structure and the spatial distribution of the mobile-immobile domain interface. This thorough analysis of the particle motion for both simple artificial structures and real rock images allowed us to derive the parametrization of the upscaled 1D model.

Dispersion upscaling from a pore scale characterization of Lagrangian velocities

NASA Astrophysics Data System (ADS)

Turuban, Régis; de Anna, Pietro; Jiménez-Martínez, Joaquín; Tabuteau, Hervé; Méheust, Yves; Le Borgne, Tanguy

2013-04-01

Mixing and reactive transport are primarily controlled by the interplay between diffusion, advection and reaction at pore scale. Yet, how the distribution and spatial correlation of the velocity field at pore scale impact these processes is still an open question. Here we present an experimental investigation of the distribution and correlation of pore scale velocities and its relation with upscaled dispersion. We use a quasi two-dimensional (2D) horizontal set up, consisting of two glass plates filled with cylinders representing the grains of the porous medium : the cell is built by soft lithography technique, wich allows for full control of the system geometry. The local velocity field is quantified from particle tracking velocimetry using microspheres that are advected with the pore scale flow. Their displacement is purely advective, as the particle size is chosen large enough to avoid diffusion. We thus obtain particle trajectories as well as lagrangian velocities in the entire system. The measured velocity field shows the existence of a network of preferential flow paths in channels with high velocities, as well as very low velocity in stagnation zones, with a non Gaussian distribution. Lagrangian velocities are long range correlated in time, which implies a non-fickian scaling of the longitudinal variance of particle positions. To upscale this process we develop an effective transport model, based on correlated continous time random walk, which is entirely parametrized by the pore scale velocity distribution and correlation. The model predictions are compared with conservative tracer test data for different Peclet numbers. Furthermore, we investigate the impact of different pore geometries on the distribution and correlation of Lagrangian velocities and we discuss the link between these properties and the effective dispersion behavior.

Characteristics of nuclepore filters with large pore size—I. Physical properties

NASA Astrophysics Data System (ADS)

John, W.; Hering, S.; Reischl, G.; Sasaki, G.; Goren, S.

Measurements of pore diameter, pore density and filter thickness have been made on Nuclepore filters of 5, 8 and 12 μm pore size. The areal distribution of the pores is random, as verified by total hole counts and by counts of overlapping holes. Filter thicknesses decrease with increasing pore diameter. The Hagen-Poiseuille formula accounts for less than half of the measured pressure drop across 12 μm pore size filters. A new calculation, including a term for the pressure drop external to the filter, accounts quantitatively for the observations. There are sufficient variations among filter batches to require knowledge of the filter parameters for each batch to ensure accurate measurements using these filters.

Modeling of time dependent localized flow shear stress and its impact on cellular growth within additive manufactured titanium implants

PubMed Central

Zhang, Ziyu; Yuan, Lang; Lee, Peter D; Jones, Eric; Jones, Julian R

2014-01-01

Bone augmentation implants are porous to allow cellular growth, bone formation and fixation. However, the design of the pores is currently based on simple empirical rules, such as minimum pore and interconnects sizes. We present a three-dimensional (3D) transient model of cellular growth based on the Navier–Stokes equations that simulates the body fluid flow and stimulation of bone precursor cellular growth, attachment, and proliferation as a function of local flow shear stress. The model's effectiveness is demonstrated for two additive manufactured (AM) titanium scaffold architectures. The results demonstrate that there is a complex interaction of flow rate and strut architecture, resulting in partially randomized structures having a preferential impact on stimulating cell migration in 3D porous structures for higher flow rates. This novel result demonstrates the potential new insights that can be gained via the modeling tool developed, and how the model can be used to perform what-if simulations to design AM structures to specific functional requirements. PMID:24664988

Hysteresis of liquid adsorption in porous media by coarse-grained Monte Carlo with direct experimental validation

NASA Astrophysics Data System (ADS)

Zeidman, Benjamin D.; Lu, Ning; Wu, David T.

2016-05-01

The effects of path-dependent wetting and drying manifest themselves in many types of physical systems, including nanomaterials, biological systems, and porous media such as soil. It is desirable to better understand how these hysteretic macroscopic properties result from a complex interplay between gasses, liquids, and solids at the pore scale. Coarse-Grained Monte Carlo (CGMC) is an appealing approach to model these phenomena in complex pore spaces, including ones determined experimentally. We present two-dimensional CGMC simulations of wetting and drying in two systems with pore spaces determined by sections from micro X-ray computed tomography: a system of randomly distributed spheres and a system of Ottawa sand. Results for the phase distribution, water uptake, and matric suction when corrected for extending to three dimensions show excellent agreement with experimental measurements on the same systems. This supports the hypothesis that CGMC can generate metastable configurations representative of experimental hysteresis and can also be used to predict hysteretic constitutive properties of particular experimental systems, given pore space images.

Hysteresis of liquid adsorption in porous media by coarse-grained Monte Carlo with direct experimental validation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zeidman, Benjamin D.; Lu, Ning; Wu, David T., E-mail: dwu@mines.edu

2016-05-07

The effects of path-dependent wetting and drying manifest themselves in many types of physical systems, including nanomaterials, biological systems, and porous media such as soil. It is desirable to better understand how these hysteretic macroscopic properties result from a complex interplay between gasses, liquids, and solids at the pore scale. Coarse-Grained Monte Carlo (CGMC) is an appealing approach to model these phenomena in complex pore spaces, including ones determined experimentally. We present two-dimensional CGMC simulations of wetting and drying in two systems with pore spaces determined by sections from micro X-ray computed tomography: a system of randomly distributed spheres andmore » a system of Ottawa sand. Results for the phase distribution, water uptake, and matric suction when corrected for extending to three dimensions show excellent agreement with experimental measurements on the same systems. This supports the hypothesis that CGMC can generate metastable configurations representative of experimental hysteresis and can also be used to predict hysteretic constitutive properties of particular experimental systems, given pore space images.« less

An Amorphous Network Model for Capillary Flow and Dispersion in a Partially Saturated Porous Medium

NASA Astrophysics Data System (ADS)

Simmons, C. S.; Rockhold, M. L.

2013-12-01

Network models of capillary flow are commonly used to represent conduction of fluids at pore scales. Typically, a flow system is described by a regular geometric lattice of interconnected tubes. Tubes constitute the pore throats, while connection junctions (nodes) are pore bodies. Such conceptualization of the geometry, however, is questionable for the pore scale, where irregularity clearly prevails, although prior published models using a regular lattice have demonstrated successful descriptions of the flow in the bulk medium. Here a network is allowed to be amorphous, and is not subject to any particular lattice structure. Few network flow models have treated partially saturated or even multiphase conditions. The research trend is toward using capillary tubes with triangular or square cross sections that have corners and always retain some fluid by capillarity when drained. In contrast, this model uses only circular capillaries, whose filled state is controlled by a capillary pressure rule for the junctions. The rule determines which capillary participate in the flow under an imposed matric potential gradient during steady flow conditions. Poiseuille's Law and Laplace equation are used to describe flow and water retention in the capillary units of the model. A modified conjugate gradient solution for steady flow that tracks which capillary in an amorphous network contribute to fluid conduction was devised for partially saturated conditions. The model thus retains the features of classical capillary models for determining hydraulic flow properties under unsaturated conditions based on distribution of non-interacting tubes, but now accounts for flow exchange at junctions. Continuity of the flow balance at every junction is solved simultaneously. The effective water retention relationship and unsaturated permeability are evaluated for an extensive enough network to represent a small bulk sample of porous medium. The model is applied for both a hypothetically randomly generate network and for a directly measured porous medium structure, by means of xray-CT scan. A randomly generated network has the benefit of providing ensemble averages for sample replicates of a medium's properties, whereas network structure measurements are expected to be more predictive. Dispersion of solute in a network flow is calculate by using particle tracking to determine the travel time breakthrough between inflow and outflow boundaries. The travel time distribution can exhibit substantial skewness that reflects both network velocity variability and mixing dilution at junctions. When local diffusion is not included, and transport is strictly advective, then the skew breakthrough is not due to mobile-immobile flow region behavior. The approach of dispersivity to its asymptotic value with sample size is examined, and may be only an indicator of particular stochastic flow variation. It is not proven that a simplified network flow model can accurately predict the hydraulic properties of a sufficiently large-size medium sample, but such a model can at least demonstrate macroscopic flow resulting from the interaction of physical processes at pore scales.

Upscaling of spectral induced polarization response using random tube networks

NASA Astrophysics Data System (ADS)

Maineult, Alexis; Revil, André; Camerlynck, Christian; Florsch, Nicolas; Titov, Konstantin

2017-05-01

In order to upscale the induced polarization (IP) response of porous media, from the pore scale to the sample scale, we implement a procedure to compute the macroscopic complex resistivity response of random tube networks. A network is made of a 2-D square-meshed grid of connected tubes, which obey to a given tube radius distribution. In a simplified approach, the electrical impedance of each tube follows a local Pelton resistivity model, with identical resistivity, chargeability and Cole-Cole exponent values for all the tubes-only the time constant varies, as it depends on the radius of each tube and on a diffusion coefficient also identical for all the tubes. By solving the conservation law for the electrical charge, the macroscopic IP response of the network is obtained. We fit successfully the macroscopic complex resistivity also by a Pelton resistivity model. Simulations on uncorrelated and correlated networks, for which the tube radius distribution is so that the decimal logarithm of the radius is normally distributed, evidence that the local and macroscopic model parameters are the same, except the Cole-Cole exponent: its macroscopic value diminishes with increasing heterogeneity (i.e. with increasing standard deviation of the radius distribution), compared to its local value. The methodology is also applied to six siliciclastic rock samples, for which the pore radius distributions from mercury porosimetry are available. These samples exhibit the same behaviour as synthetic media, that is, the macroscopic Cole-Cole exponent is always lower than the local one. As a conclusion, the pore network method seems to be a promising tool for studying the upscaling of the IP response of porous media.

Coupled flow and deformations in granular systems beyond the pendular regime

NASA Astrophysics Data System (ADS)

Yuan, Chao; Chareyre, Bruno; Darve, Felix

2017-06-01

A pore-scale numerical model is proposed for simulating the quasi-static primary drainage and the hydro-mechanical couplings in multiphase granular systems. The solid skeleton is idealized to a dense random packing of polydisperse spheres by DEM. The fluids (nonwetting and wetting phases) space is decomposed to a network of tetrahedral pores based on the Regular Triangulation method. The local drainage rules and invasion logic are defined. The fluid forces acting on solid grains are formulated. The model can simulate the hydraulic evolution from a fully saturated state to a low level of saturation but beyond the pendular regime. The features of wetting phase entrapments and capillary fingering can also be reproduced. Finally, a primary drainage test is performed on a 40,000 spheres of sample. The water retention curve is obtained. The solid skeleton first shrinks then swells.

Unifying Pore Network Modeling, Continuous Time Random Walk Theory and Experiment - Accomplishments and Future Directions

NASA Astrophysics Data System (ADS)

Bijeljic, B.

2008-05-01

This talk will describe and highlight the advantages offered by a methodology that unifies pore network modeling, CTRW theory and experiment in description of solute dispersion in porous media. Solute transport in a porous medium is characterized by the interplay of advection and diffusion (described by Peclet number, Pe) that cause spreading of solute particles. This spreading is traditionally described by dispersion coefficients, D, defined by σ 2 = 2Dt, where σ 2 is the variance of the solute position and t is the time. Using a pore-scale network model based on particle tracking, the rich Peclet- number dependence of dispersion coefficient is predicted from first principles and is shown to compare well with experimental data for restricted diffusion, transition, power-law and mechanical dispersion regimes in the asymptotic limit. In the asymptotic limit D is constant and can be used in an averaged advection-dispersion equation. However, it is highly important to recognize that, until the velocity field is fully sampled, the particle transport is non-Gaussian and D possesses temporal or spatial variation. Furthermore, temporal probability density functions (PDF) of tracer particles are studied in pore networks and an excellent agreement for the spectrum of transition times for particles from pore to pore is obtained between network model results and CTRW theory. Based on the truncated power-law interpretation of PDF-s, the physical origin of the power-law scaling of dispersion coefficient vs. Peclet number has been explained for unconsolidated porous media, sands and a number of sandstones, arriving at the same conclusion from numerical network modelling, analytic CTRW theory and experiment. Future directions for further applications of the methodology presented are discussed in relation to the scale- dependent solute dispersion and reactive transport. Significance of pre-asymptotic dispersion in porous media is addressed from pore-scale upwards and the impact of heterogeneity is discussed. The length traveled by solute plumes before Gaussian behaviour is reached increases with an increase in heterogeneity and/or Pe. This opens up the question on the nature of dispersion in natural systems where the heterogeneities at the larger scales will profoundly increase the range of velocities in the aquifer, thus considerably delaying the asymptotic approach to Gaussian behaviour. As a consequence, the asymptotic behaviour might not be reached at the field scale.

The effect of pore-scale geometry and wettability on two-phase relative permeabilities within elementary cells

NASA Astrophysics Data System (ADS)

Bianchi Janetti, Emanuela; Riva, Monica; Guadagnini, Alberto

2017-04-01

We study the relative role of the complex pore space geometry and wettability of the solid matrix on the quantification of relative permeabilities characterizing steady state immiscible two-phase flow in porous media. We do so by considering elementary cells, which are typically employed in upscaling frameworks based on, e.g., homogenization or volume averaging. In this context one typically relies on the solution of pore-scale physics at a scale which is much smaller than that of an investigated porous system. Pressure-driven two-phase flow following simultaneous co-current injection of water and oil is numerically solved for a suite of regular and stochastically generated two-dimensional explicit elementary cells with fixed porosity and sharing main topological/morphological features. We show that relative permeabilities of the randomly generated elementary cells are significantly influenced by the formation of preferential percolation paths (principal pathways), giving rise to a strongly nonuniform distribution of fluid fluxes. These pathways are a result of the spatially variable resistance that the random pore structures exert on the fluid. The overall effect on relative permeabilities of the diverse organization of principal pathways, as driven by a given random realization at the scale of the unit cell, is significantly larger than that of the wettability of the host rock. In contrast to what can be observed for the random cells analyzed, relative permeabilities of regular cells display a clear trend with contact angle at the investigated scale. Our findings suggest the need to perform systematic upscaling studies in a stochastic context, to propagate the effects of uncertain pore space geometries to a probabilistic description of relative permeability curves at the continuum scale.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tammas-Williams, S., E-mail: Samuel.tammas-wiliams@manchester.ac.uk; Department of Materials Science and Engineering, University of Sheffield, Sheffield S1 3JD; Zhao, H.

Selective Electron Beam Melting (SEBM) is a promising powder bed Additive Manufacturing technique for near-net-shape manufacture of high-value titanium components. However without post-manufacture HIPing the fatigue life of SEBM parts is currently dominated by the presence of porosity. In this study, the size, volume fraction, and spatial distribution of the pores in model samples have been characterised in 3D, using X-ray Computed Tomography, and correlated to the process variables. The average volume fraction of the pores (< 0.2%) was measured to be lower than that usually observed in competing processes, such as selective laser melting, but a strong relationship wasmore » found with the different beam strategies used to contour, and infill by hatching, a part section. The majority of pores were found to be small spherical gas pores, concentrated in the infill hatched region; this was attributed to the lower energy density and less focused beam used in the infill strategy allowing less opportunity for gas bubbles to escape the melt pool. Overall, increasing the energy density or focus of the beam was found to correlate strongly to a reduction in the level of gas porosity. Rarer irregular shaped pores were mostly located in the contour region and have been attributed to a lack of fusion between powder particles. - Graphical abstract: Display Omitted - Highlights: • Vast majority of defects detected were small spherical gas pores. • Gas bubbles trapped in the powder granules expand and coalesce in the melt pool. • Pores have been shown not to be randomly distributed. • Larger and deeper melt pools give more opportunity for gas to escape. • Minor changes to melt strategy result in significant reductions in pore population.« less

Observations of the Dynamic Connectivity of the Non-Wetting Phase During Steady State Flow at the Pore Scale Using 3D X-ray Microtomography

NASA Astrophysics Data System (ADS)

Reynolds, C. A.; Menke, H. P.; Blunt, M. J.; Krevor, S. C.

2015-12-01

We observe a new type of non-wetting phase flow using time-resolved pore scale imaging. The traditional conceptual model of drainage involves a non-wetting phase invading a porous medium saturated with a wetting phase as either a fixed, connected flow path through the centres of pores or as discrete ganglia which move individually through the pore space, depending on the capillary number. We observe a new type of flow behaviour at low capillary number in which the flow of the non-wetting phase occurs through networks of persistent ganglia that occupy the large pores but continuously rearrange their connectivity (Figure 1). Disconnections and reconnections occur randomly to provide short-lived pseudo-steady state flow paths between pores. This process is distinctly different to the notion of flowing ganglia which coalesce and break-up. The size distribution of ganglia is dependent on capillary number. Experiments were performed by co-injecting N2and 25 wt% KI brine into a Bentheimer sandstone core (4mm diameter, 35mm length) at 50°C and 10 MPa. Drainage was performed at three flow rates (0.04, 0.3 and 1 ml/min) at a constant fractional flow of 0.5 and the variation in ganglia populations and connectivity observed. We obtained images of the pore space during steady state flow with a time resolution of 43 s over 1-2 hours. Experiments were performed at the Diamond Light Source synchrotron. Figure 1. The position of N2 in the pore space during steady state flow is summed over 40 time steps. White indicates that N2 occupies the space over >38 time steps and red <5 time steps.

Cycle flux algebra for ion and water flux through the KcsA channel single-file pore links microscopic trajectories and macroscopic observables.

PubMed

Oiki, Shigetoshi; Iwamoto, Masayuki; Sumikama, Takashi

2011-01-31

In narrow pore ion channels, ions and water molecules diffuse in a single-file manner and cannot pass each other. Under such constraints, ion and water fluxes are coupled, leading to experimentally observable phenomena such as the streaming potential. Analysis of this coupled flux would provide unprecedented insights into the mechanism of permeation. In this study, ion and water permeation through the KcsA potassium channel was the focus, for which an eight-state discrete-state Markov model has been proposed based on the crystal structure, exhibiting four ion-binding sites. Random transitions on the model lead to the generation of the net flux. Here we introduced the concept of cycle flux to derive exact solutions of experimental observables from the permeation model. There are multiple cyclic paths on the model, and random transitions complete the cycles. The rate of cycle completion is called the cycle flux. The net flux is generated by a combination of cyclic paths with their own cycle flux. T.L. Hill developed a graphical method of exact solutions for the cycle flux. This method was extended to calculate one-way cycle fluxes of the KcsA channel. By assigning the stoichiometric numbers for ion and water transfer to each cycle, we established a method to calculate the water-ion coupling ratio (CR(w-i)) through cycle flux algebra. These calculations predicted that CR(w-i) would increase at low potassium concentrations. One envisions an intuitive picture of permeation as random transitions among cyclic paths, and the relative contributions of the cycle fluxes afford experimental observables.

Cycle Flux Algebra for Ion and Water Flux through the KcsA Channel Single-File Pore Links Microscopic Trajectories and Macroscopic Observables

PubMed Central

Oiki, Shigetoshi; Iwamoto, Masayuki; Sumikama, Takashi

2011-01-01

In narrow pore ion channels, ions and water molecules diffuse in a single-file manner and cannot pass each other. Under such constraints, ion and water fluxes are coupled, leading to experimentally observable phenomena such as the streaming potential. Analysis of this coupled flux would provide unprecedented insights into the mechanism of permeation. In this study, ion and water permeation through the KcsA potassium channel was the focus, for which an eight-state discrete-state Markov model has been proposed based on the crystal structure, exhibiting four ion-binding sites. Random transitions on the model lead to the generation of the net flux. Here we introduced the concept of cycle flux to derive exact solutions of experimental observables from the permeation model. There are multiple cyclic paths on the model, and random transitions complete the cycles. The rate of cycle completion is called the cycle flux. The net flux is generated by a combination of cyclic paths with their own cycle flux. T.L. Hill developed a graphical method of exact solutions for the cycle flux. This method was extended to calculate one-way cycle fluxes of the KcsA channel. By assigning the stoichiometric numbers for ion and water transfer to each cycle, we established a method to calculate the water-ion coupling ratio (CR w-i) through cycle flux algebra. These calculations predicted that CR w-i would increase at low potassium concentrations. One envisions an intuitive picture of permeation as random transitions among cyclic paths, and the relative contributions of the cycle fluxes afford experimental observables. PMID:21304994

How realistic is the pore size distribution calculated from adsorption isotherms if activated carbon is composed of fullerene-like fragments?

PubMed

Terzyk, Artur P; Furmaniak, Sylwester; Harris, Peter J F; Gauden, Piotr A; Włoch, Jerzy; Kowalczyk, Piotr; Rychlicki, Gerhard

2007-11-28

A plausible model for the structure of non-graphitizing carbon is one which consists of curved, fullerene-like fragments grouped together in a random arrangement. Although this model was proposed several years ago, there have been no attempts to calculate the properties of such a structure. Here, we determine the density, pore size distribution and adsorption properties of a model porous carbon constructed from fullerene-like elements. Using the method proposed recently by Bhattacharya and Gubbins (BG), which was tested in this study for ideal and defective carbon slits, the pore size distributions (PSDs) of the initial model and two related carbon models are calculated. The obtained PSD curves show that two structures are micro-mesoporous (with different ratio of micro/mesopores) and the third is strictly microporous. Using the grand canonical Monte Carlo (GCMC) method, adsorption isotherms of Ar (87 K) are simulated for all the structures. Finally PSD curves are calculated using the Horvath-Kawazoe, non-local density functional theory (NLDFT), Nguyen and Do, and Barrett-Joyner-Halenda (BJH) approaches, and compared with those predicted by the BG method. This is the first study in which different methods of calculation of PSDs for carbons from adsorption data can be really verified, since absolute (i.e. true) PSDs are obtained using the BG method. This is also the first study reporting the results of computer simulations of adsorption on fullerene-like carbon models.

Synthesis, characterization, and evaluation of a superficially porous particle with unique, elongated pore channels normal to the surface.

PubMed

Wei, Ta-Chen; Mack, Anne; Chen, Wu; Liu, Jia; Dittmann, Monika; Wang, Xiaoli; Barber, William E

2016-04-01

In recent years, superficially porous particles (SPPs) have drawn great interest because of their special particle characteristics and improvement in separation efficiency. Superficially porous particles are currently manufactured by adding silica nanoparticles onto solid cores using either a multistep multilayer process or one-step coacervation process. The pore size is mainly controlled by the size of the silica nanoparticles and the tortuous pore channel geometry is determined by how those nanoparticles randomly aggregate. Such tortuous pore structure is also similar to that of all totally porous particles used in HPLC today. In this article, we report on the development of a next generation superficially porous particle with a unique pore structure that includes a thinner shell thickness and ordered pore channels oriented normal to the particle surface. The method of making the new superficially porous particles is a process called pseudomorphic transformation (PMT), which is a form of micelle templating. Porosity is no longer controlled by randomly aggregated nanoparticles but rather by micelles that have an ordered liquid crystal structure. The new particle possesses many advantages such as a narrower particle size distribution, thinner porous layer with high surface area and, most importantly, highly ordered, non-tortuous pore channels oriented normal to the particle surface. This PMT process has been applied to make 1.8-5.1μm SPPs with pore size controlled around 75Å and surface area around 100m(2)/g. All particles with different sizes show the same unique pore structure with tunable pore size and shell thickness. The impact of the novel pore structure on the performance of these particles is characterized by measuring van Deemter curves and constructing kinetic plots. Reduced plate heights as low as 1.0 have been achieved on conventional LC instruments. This indicates higher efficiency of such particles compared to conventional totally porous and superficially porous particles. Copyright © 2016 Elsevier B.V. All rights reserved.

Modeling of time dependent localized flow shear stress and its impact on cellular growth within additive manufactured titanium implants.

PubMed

Zhang, Ziyu; Yuan, Lang; Lee, Peter D; Jones, Eric; Jones, Julian R

2014-11-01

Bone augmentation implants are porous to allow cellular growth, bone formation and fixation. However, the design of the pores is currently based on simple empirical rules, such as minimum pore and interconnects sizes. We present a three-dimensional (3D) transient model of cellular growth based on the Navier-Stokes equations that simulates the body fluid flow and stimulation of bone precursor cellular growth, attachment, and proliferation as a function of local flow shear stress. The model's effectiveness is demonstrated for two additive manufactured (AM) titanium scaffold architectures. The results demonstrate that there is a complex interaction of flow rate and strut architecture, resulting in partially randomized structures having a preferential impact on stimulating cell migration in 3D porous structures for higher flow rates. This novel result demonstrates the potential new insights that can be gained via the modeling tool developed, and how the model can be used to perform what-if simulations to design AM structures to specific functional requirements. © 2014 Wiley Periodicals, Inc.

Upscaling of dilution and mixing using a trajectory based Spatial Markov random walk model in a periodic flow domain

NASA Astrophysics Data System (ADS)

Sund, Nicole L.; Porta, Giovanni M.; Bolster, Diogo

2017-05-01

The Spatial Markov Model (SMM) is an upscaled model that has been used successfully to predict effective mean transport across a broad range of hydrologic settings. Here we propose a novel variant of the SMM, applicable to spatially periodic systems. This SMM is built using particle trajectories, rather than travel times. By applying the proposed SMM to a simple benchmark problem we demonstrate that it can predict mean effective transport, when compared to data from fully resolved direct numerical simulations. Next we propose a methodology for using this SMM framework to predict measures of mixing and dilution, that do not just depend on mean concentrations, but are strongly impacted by pore-scale concentration fluctuations. We use information from trajectories of particles to downscale and reconstruct pore-scale approximate concentration fields from which mixing and dilution measures are then calculated. The comparison between measurements from fully resolved simulations and predictions with the SMM agree very favorably.

A multi-scale homogenization model for fine-grained porous viscoplastic polycrystals: I - Finite-strain theory

NASA Astrophysics Data System (ADS)

Song, Dawei; Ponte Castañeda, P.

2018-06-01

We make use of the recently developed iterated second-order homogenization method to obtain finite-strain constitutive models for the macroscopic response of porous polycrystals consisting of large pores randomly distributed in a fine-grained polycrystalline matrix. The porous polycrystal is modeled as a three-scale composite, where the grains are described by single-crystal viscoplasticity and the pores are assumed to be large compared to the grain size. The method makes use of a linear comparison composite (LCC) with the same substructure as the actual nonlinear composite, but whose local properties are chosen optimally via a suitably designed variational statement. In turn, the effective properties of the resulting three-scale LCC are determined by means of a sequential homogenization procedure, utilizing the self-consistent estimates for the effective behavior of the polycrystalline matrix, and the Willis estimates for the effective behavior of the porous composite. The iterated homogenization procedure allows for a more accurate characterization of the properties of the matrix by means of a finer "discretization" of the properties of the LCC to obtain improved estimates, especially at low porosities, high nonlinearties and high triaxialities. In addition, consistent homogenization estimates for the average strain rate and spin fields in the pores and grains are used to develop evolution laws for the substructural variables, including the porosity, pore shape and orientation, as well as the "crystallographic" and "morphological" textures of the underlying matrix. In Part II of this work has appeared in Song and Ponte Castañeda (2018b), the model will be used to generate estimates for both the instantaneous effective response and the evolution of the microstructure for porous FCC and HCP polycrystals under various loading conditions.

Modeling postshock evolution of large electropores

NASA Astrophysics Data System (ADS)

Neu, John C.; Krassowska, Wanda

2003-02-01

The Smoluchowski equation (SE), which describes the evolution of pores created by electric shocks, cannot be applied to modeling large and long-lived pores for two reasons: (1) it does not predict pores of radius above 20 nm without also predicting membrane rupture; (2) it does not predict postshock growth of pores. This study proposes a model in which pores are coupled by membrane tension, resulting in a nonlinear generalization of SE. The predictions of the model are explored using examples of homogeneous (all pore radii r are equal) and heterogeneous (0⩽r⩽rmax) distributions of pores. Pores in a homogeneous population either shrink to zero or assume a stable radius corresponding to the minimum of the bilayer energy. For a heterogeneous population, such a stable radius does not exist. All pores, except rmax, shrink to zero and rmax grows to infinity. However, the unbounded growth of rmax is not physical because the number of pores per cell decreases in time and the continuum model loses validity. When the continuum formulation is replaced by the discrete one, the model predicts the coarsening process: all pores, except rmax, shrink to zero and rmax assumes a stable radius. Thus, the model with tension-coupled pores does not predict membrane rupture and the predicted postshock growth of pores is consistent with experimental evidence.

Improving the prospects of cleavage-based nanopore sequencing engines

NASA Astrophysics Data System (ADS)

Brady, Kyle T.; Reiner, Joseph E.

2015-08-01

Recently proposed methods for DNA sequencing involve the use of cleavage-based enzymes attached to the opening of a nanopore. The idea is that DNA interacting with either an exonuclease or polymerase protein will lead to a small molecule being cleaved near the mouth of the nanopore, and subsequent entry into the pore will yield information about the DNA sequence. The prospects for this approach seem promising, but it has been shown that diffusion related effects impose a limit on the capture probability of molecules by the pore, which limits the efficacy of the technique. Here, we revisit the problem with the goal of optimizing the capture probability via a step decrease in the nucleotide diffusion coefficient between the pore and bulk solutions. It is shown through random walk simulations and a simplified analytical model that decreasing the molecule's diffusion coefficient in the bulk relative to its value in the pore increases the nucleotide capture probability. Specifically, we show that at sufficiently high applied transmembrane potentials (≥100 mV), increasing the potential by a factor f is equivalent to decreasing the diffusion coefficient ratio Dbulk/Dpore by the same factor f. This suggests a promising route toward implementation of cleavage-based sequencing protocols. We also discuss the feasibility of forming a step function in the diffusion coefficient across the pore-bulk interface.

Modeling non-Fickian dispersion by use of the velocity PDF on the pore scale

NASA Astrophysics Data System (ADS)

Kooshapur, Sheema; Manhart, Michael

2015-04-01

For obtaining a description of reactive flows in porous media, apart from the geometrical complications of resolving the velocities and scalar values, one has to deal with the additional reactive term in the transport equation. An accurate description of the interface of the reacting fluids - which is strongly influenced by dispersion- is essential for resolving this term. In REV-based simulations the reactive term needs to be modeled taking sub-REV fluctuations and possibly non-Fickian dispersion into account. Non-Fickian dispersion has been observed in strongly heterogeneous domains and in early phases of transport. A fully resolved solution of the Navier-Stokes and transport equations which yields a detailed description of the flow properties, dispersion, interfaces of fluids, etc. however, is not practical for domains containing more than a few thousand grains, due to the huge computational effort required. Through Probability Density Function (PDF) based methods, the velocity distribution in the pore space can facilitate the understanding and modelling of non-Fickian dispersion [1,2]. Our aim is to model the transition between non-Fickian and Fickian dispersion in a random sphere pack within the framework of a PDF based transport model proposed by Meyer and Tchelepi [1,3]. They proposed a stochastic transport model where velocity components of tracer particles are represented by a continuous Markovian stochastic process. In addition to [3], we consider the effects of pore scale diffusion and formulate a different stochastic equation for the increments in velocity space from first principles. To assess the terms in this equation, we performed Direct Numerical Simulations (DNS) for solving the Navier-Stokes equation on a random sphere pack. We extracted the PDFs and statistical moments (up to the 4th moment) of the stream-wise velocity, u, and first and second order velocity derivatives both independent and conditioned on velocity. By using this data and combining the Taylor expansion of velocity increments, du, and the Langevin equation for point particles we obtained the components of velocity fluxes which point to a drift and diffusion behavior in the velocity space. Thus a partial differential equation for the velocity PDF has been formulated that constitutes an advection-diffusion equation in velocity space (a Fokker-Planck equation) in which the drift and diffusion coefficients are obtained using the velocity conditioned statistics of the derivatives of the pore scale velocity field. This has been solved by both a Random Walk (RW) model and a Finite Volume method. We conclude that both, these methods are able to simulate the velocity PDF obtained by DNS. References [1] D. W. Meyer, P. Jenny, H.A.Tschelepi, A joint velocity-concentration PDF method for traqcer flow in heterogeneous porous media, Water Resour.Res., 46, W12522, (2010). [2] Nowak, W., R. L. Schwede, O. A. Cirpka, and I. Neuweiler, Probability density functions of hydraulic head and velocity in three-dimensional heterogeneous porous media, Water Resour.Res., 44, W08452, (2008) [3] D. W. Meyer, H. A. Tchelepi, Particle-based transport model with Markovian velocity processes for tracer dispersion in highly heterogeneous porous media, Water Resour. Res., 46, W11552, (2010)

Controlled nucleation in freeze-drying: effects on pore size in the dried product layer, mass transfer resistance, and primary drying rate.

PubMed

Konstantinidis, Alex K; Kuu, Wei; Otten, Lori; Nail, Steven L; Sever, Robert R

2011-08-01

A novel and scalable method has been developed to enable control of the ice nucleation step for the freezing process during lyophilization. This method manipulates the chamber pressure of the freeze dryer to simultaneously induce nucleation in all product vials at a desired temperature. The effects of controlled nucleation on the drying rate of various formulations including 5% (w/w) mannitol, 5% (w/w) sucrose, and a mixture of 3% (w/w) mannitol and 2% (w/w) sucrose were studied. For a 5% (w/w) mannitol, uncontrolled ice nucleation occurred randomly at product temperatures between -8.0°C and -15.9°C as the vials were cooled to -40°C. Controlled ice nucleation was achieved at product temperatures between -2.3°C and -3.7°C. The effect of nucleation control on the effective pore radius (r(e) ) of the cake was determined from the product temperature profiles using a pore diffusion model in combination with a nonlinear parameter estimation approach reported earlier. Results show that the value of r(e) for 5% (w/w) mannitol was enlarged from 13 to 27 μm by uniformly inducing nucleation at higher temperatures. Applying the resistance parameters obtained from the pore diffusion model for 5% (w/w) mannitol, optimized cycles were theoretically generated and experimentally tested, resulting in a 41% reduction in primary drying time. Copyright © 2011 Wiley-Liss, Inc.

Towards a rigorous mesoscale modeling of reactive flow and transport in an evolving porous medium and its applications to soil science

NASA Astrophysics Data System (ADS)

Ray, Nadja; Rupp, Andreas; Knabner, Peter

2016-04-01

Soil is arguably the most prominent example of a natural porous medium that is composed of a porous matrix and a pore space. Within this framework and in terms of soil's heterogeneity, we first consider transport and fluid flow at the pore scale. From there, we develop a mechanistic model and upscale it mathematically to transfer our model from the small scale to that of the mesoscale (laboratory scale). The mathematical framework of (periodic) homogenization (in principal) rigorously facilitates such processes by exactly computing the effective coefficients/parameters by means of the pore geometry and processes. In our model, various small-scale soil processes may be taken into account: molecular diffusion, convection, drift emerging from electric forces, and homogeneous reactions of chemical species in a solvent. Additionally, our model may consider heterogeneous reactions at the porous matrix, thus altering both the porosity and the matrix. Moreover, our model may additionally address biophysical processes, such as the growth of biofilms and how this affects the shape of the pore space. Both of the latter processes result in an intrinsically variable soil structure in space and time. Upscaling such models under the assumption of a locally periodic setting must be performed meticulously to preserve information regarding the complex coupling of processes in the evolving heterogeneous medium. Generally, a micro-macro model emerges that is then comprised of several levels of couplings: Macroscopic equations that describe the transport and fluid flow at the scale of the porous medium (mesoscale) include averaged time- and space-dependent coefficient functions. These functions may be explicitly computed by means of auxiliary cell problems (microscale). Finally, the pore space in which the cell problems are defined is time- and space dependent and its geometry inherits information from the transport equation's solutions. Numerical computations using mixed finite elements and potentially random initial data, e.g. that of porosity, complement our theoretical results. Our investigations contribute to the theoretical understanding of the link between soil formation and soil functions. This general framework may be applied to various problems in soil science for a range of scales, such as the formation and turnover of microaggregates or soil remediation.

Impact of geometrical properties on permeability and fluid phase distribution in porous media

NASA Astrophysics Data System (ADS)

Lehmann, P.; Berchtold, M.; Ahrenholz, B.; Tölke, J.; Kaestner, A.; Krafczyk, M.; Flühler, H.; Künsch, H. R.

2008-09-01

To predict fluid phase distribution in porous media, the effect of geometric properties on flow processes must be understood. In this study, we analyze the effect of volume, surface, curvature and connectivity (the four Minkowski functionals) on the hydraulic conductivity and the water retention curve. For that purpose, we generated 12 artificial structures with 800 3 voxels (the units of a 3D image) and compared them with a scanned sand sample of the same size. The structures were generated with a Boolean model based on a random distribution of overlapping ellipsoids whose size and shape were chosen to fulfill the criteria of the measured functionals. The pore structure of sand material was mapped with X-rays from synchrotrons. To analyze the effect of geometry on water flow and fluid distribution we carried out three types of analysis: Firstly, we computed geometrical properties like chord length, distance from the solids, pore size distribution and the Minkowski functionals as a function of pore size. Secondly, the fluid phase distribution as a function of the applied pressure was calculated with a morphological pore network model. Thirdly, the permeability was determined using a state-of-the-art lattice-Boltzmann method. For the simulated structure with the true Minkowski functionals the pores were larger and the computed air-entry value of the artificial medium was reduced to 85% of the value obtained from the scanned sample. The computed permeability for the geometry with the four fitted Minkowski functionals was equal to the permeability of the scanned image. The permeability was much more sensitive to the volume and surface than to curvature and connectivity of the medium. We conclude that the Minkowski functionals are not sufficient to characterize the geometrical properties of a porous structure that are relevant for the distribution of two fluid phases. Depending on the procedure to generate artificial structures with predefined Minkowski functionals, structures differing in pore size distribution can be obtained.

A simple shape-free model for pore-size estimation with positron annihilation lifetime spectroscopy

NASA Astrophysics Data System (ADS)

Wada, Ken; Hyodo, Toshio

2013-06-01

Positron annihilation lifetime spectroscopy is one of the methods for estimating pore size in insulating materials. We present a shape-free model to be used conveniently for such analysis. A basic model in classical picture is modified by introducing a parameter corresponding to an effective size of the positronium (Ps). This parameter is adjusted so that its Ps-lifetime to pore-size relation merges smoothly with that of the well-established Tao-Eldrup model (with modification involving the intrinsic Ps annihilation rate) applicable to very small pores. The combined model, i.e., modified Tao-Eldrup model for smaller pores and the modified classical model for larger pores, agrees surprisingly well with the quantum-mechanics based extended Tao-Eldrup model, which deals with Ps trapped in and thermally equilibrium with a rectangular pore.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ahmed, K.; Tonks, M.; Zhang, Y.

A detailed phase field model for the effect of pore drag on grain growth kinetics was implemented in MARMOT. The model takes into consideration both the curvature-driven grain boundary motion and pore migration by surface diffusion. As such, the model accounts for the interaction between pore and grain boundary kinetics, which tends to retard the grain growth process. Our 2D and 3D simulations demonstrate that the model capture all possible pore-grain boundary interactions proposed in theoretical models. For high enough surface mobility, the pores move along with the migrating boundary as a quasi-rigid-body, albeit hindering its migration rate compared tomore » the pore-free case. For less mobile pores, the migrating boundary can separate from the pores. For the pore-controlled grain growth kinetics, the model predicts a strong dependence of the growth rate on the number of pores, pore size, and surface diffusivity in agreement with theroretical models. An evolution equation for the grain size that includes these parameters was derived and showed to agree well with numerical solution. It shows a smooth transition from boundary-controlled kinetics to pore-controlled kinetics as the surface diffusivity decreases or the number of pores or their size increases. This equation can be utilized in BISON to give accurate estimate for the grain size evolution. This will be accomplished in the near future. The effect of solute drag and anisotropy of grain boundary on grain growth will be investigated in future studies.« less

Modeling for stress-strain curve of a porous NiTi under compressive loading

NASA Astrophysics Data System (ADS)

Zhao, Ying; Taya, Minoru

2005-05-01

Two models for predicting the stress-strain curve of porous NiTi under compressive loading are presented in this paper. Porous NiTi shape memory alloy is investigated as a composite composed of solid NiTi as matrix and pores as inclusions. Eshelby"s equivalent inclusion method and Mori-Tanaka"s mean-field theory are employed in both models. In the first model, the geometry of the pores is assumed as sphere. The composite is with close-cells. While in the second model, two geometries of the pores, sphere and ellipsoid, are investigated. The pores are interconnected to each other forming an open-cell microstructure. The two adjacent pores connected along equator ring are investigated as a unit. Two pores interact with each other as they are connected. The average eigenstrain of each unit is obtained by taking the average of each pore"s eigenstrain. The stress-strain curves of porous shape memory alloy with spherical pores and ellipsoidal pores are compared, it is found that the shape of the pores has a nonignorable influence on the mechanical property of the porous NiTi. Comparison of the stress-strain curves of the two models shows that introducing of the average eigenstrains in the second model makes the predictions more agreeable to the experimental results.

Data Assimilation by Ensemble Kalman Filter during One-Dimensional Nonlinear Consolidation in Randomly Heterogeneous Highly Compressible Aquitards

NASA Astrophysics Data System (ADS)

Zapata Norberto, B.; Morales-Casique, E.; Herrera, G. S.

2017-12-01

Severe land subsidence due to groundwater extraction may occur in multiaquifer systems where highly compressible aquitards are present. The highly compressible nature of the aquitards leads to nonlinear consolidation where the groundwater flow parameters are stress-dependent. The case is further complicated by the heterogeneity of the hydrogeologic and geotechnical properties of the aquitards. We explore the effect of realistic vertical heterogeneity of hydrogeologic and geotechnical parameters on the consolidation of highly compressible aquitards by means of 1-D Monte Carlo numerical simulations. 2000 realizations are generated for each of the following parameters: hydraulic conductivity (K), compression index (Cc) and void ratio (e). The correlation structure, the mean and the variance for each parameter were obtained from a literature review about field studies in the lacustrine sediments of Mexico City. The results indicate that among the parameters considered, random K has the largest effect on the ensemble average behavior of the system. Random K leads to the largest variance (and therefore largest uncertainty) of total settlement, groundwater flux and time to reach steady state conditions. We further propose a data assimilation scheme by means of ensemble Kalman filter to estimate the ensemble mean distribution of K, pore-pressure and total settlement. We consider the case where pore-pressure measurements are available at given time intervals. We test our approach by generating a 1-D realization of K with exponential spatial correlation, and solving the nonlinear flow and consolidation problem. These results are taken as our "true" solution. We take pore-pressure "measurements" at different times from this "true" solution. The ensemble Kalman filter method is then employed to estimate ensemble mean distribution of K, pore-pressure and total settlement based on the sequential assimilation of these pore-pressure measurements. The ensemble-mean estimates from this procedure closely approximate those from the "true" solution. This procedure can be easily extended to other random variables such as compression index and void ratio.

Pore-scale modeling of capillary trapping in water-wet porous media: A new cooperative pore-body filling model

NASA Astrophysics Data System (ADS)

Ruspini, L. C.; Farokhpoor, R.; Øren, P. E.

2017-10-01

We present a pore-network model study of capillary trapping in water-wet porous media. The amount and distribution of trapped non-wetting phase is determined by the competition between two trapping mechanisms - snap-off and cooperative pore-body filling. We develop a new model to describe the pore-body filling mechanism in geologically realistic pore-networks. The model accounts for the geometrical characteristics of the pore, the spatial location of the connecting throats and the local fluid topology at the time of the displacement. We validate the model by comparing computed capillary trapping curves with published data for four different water-wet rocks. Computations are performed on pore-networks extracted from micro-CT images and process-based reconstructions of the actual rocks used in the experiments. Compared with commonly used stochastic models, the new model describes more accurately the experimental measurements, especially for well connected porous systems where trapping is controlled by subtleties of the pore structure. The new model successfully predicts relative permeabilities and residual saturation for Bentheimer sandstone using in-situ measured contact angles as input to the simulations. The simulated trapped cluster size distributions are compared with predictions from percolation theory.

Coverage dependent molecular assembly of anthraquinone on Au(111)

NASA Astrophysics Data System (ADS)

DeLoach, Andrew S.; Conrad, Brad R.; Einstein, T. L.; Dougherty, Daniel B.

2017-11-01

A scanning tunneling microscopy study of anthraquinone (AQ) on the Au(111) surface shows that the molecules self-assemble into several structures depending on the local surface coverage. At high coverages, a close-packed saturated monolayer is observed, while at low coverages, mobile surface molecules coexist with stable chiral hexamer clusters. At intermediate coverages, a disordered 2D porous network interlinking close-packed islands is observed in contrast to the giant honeycomb networks observed for the same molecule on Cu(111). This difference verifies the predicted extreme sensitivity [J. Wyrick et al., Nano Lett. 11, 2944 (2011)] of the pore network to small changes in the surface electronic structure. Quantitative analysis of the 2D pore network reveals that the areas of the vacancy islands are distributed log-normally. Log-normal distributions are typically associated with the product of random variables (multiplicative noise), and we propose that the distribution of pore sizes for AQ on Au(111) originates from random linear rate constants for molecules to either desorb from the surface or detach from the region of a nucleated pore.

Coverage dependent molecular assembly of anthraquinone on Au(111).

PubMed

DeLoach, Andrew S; Conrad, Brad R; Einstein, T L; Dougherty, Daniel B

2017-11-14

A scanning tunneling microscopy study of anthraquinone (AQ) on the Au(111) surface shows that the molecules self-assemble into several structures depending on the local surface coverage. At high coverages, a close-packed saturated monolayer is observed, while at low coverages, mobile surface molecules coexist with stable chiral hexamer clusters. At intermediate coverages, a disordered 2D porous network interlinking close-packed islands is observed in contrast to the giant honeycomb networks observed for the same molecule on Cu(111). This difference verifies the predicted extreme sensitivity [J. Wyrick et al., Nano Lett. 11, 2944 (2011)] of the pore network to small changes in the surface electronic structure. Quantitative analysis of the 2D pore network reveals that the areas of the vacancy islands are distributed log-normally. Log-normal distributions are typically associated with the product of random variables (multiplicative noise), and we propose that the distribution of pore sizes for AQ on Au(111) originates from random linear rate constants for molecules to either desorb from the surface or detach from the region of a nucleated pore.

Flow and Transport in Complex Microporous Carbonates as a Consequence of Separation of Scales

NASA Astrophysics Data System (ADS)

Bijeljic, B.; Raeini, A. Q.; Lin, Q.; Blunt, M. J.

2017-12-01

Some of the most important examples of flow and transport in complex pore structures are found in subsurface applications such as contaminant hydrology, carbon storage and enhanced oil recovery. Carbonate rock structures contain most of the world's oil reserves, considerable amount of water reserves, and potentially hold a storage capacity for carbon dioxide. However, this type of pore space is difficult to represent due to complexities associated with a wide range of pore sizes and variation in connectivity which poses a considerable challenge for quantitative predictions of transport across multiple scales.A new concept unifying X-ray tomography experiment and direct numerical simulation has been developed that relies on full description flow and solute transport at the pore scale. Differential imaging method (Lin et al. 2016) provides rich information in microporous space, while advective and diffusive mass transport are simulated on micro-CT images of pore-space: Navier-Stokes equations are solved for flow in the image voxels comprising the pore space, streamline-based simulation is used to account for advection, and diffusion is superimposed by random walk.Quantitative validation has been done on analytical solutions for diffusion and by comparing the model predictions versus the experimental NMR measurements in the dual porosity beadpack. Furthermore, we discriminate signatures of multi-scale transport behaviour for a range of carbonate rock (Figure 1), dependent on the heterogeneity of the inter- and intra-grain pore space, heterogeneity in the flow field, and the mass transfer characteristics of the porous media. Finally, we demonstrate the predictive capabilities of the model through an analysis that includes a number of probability density functions flow and transport (PDFs) measures of non-Fickian transport on the micro-CT i935mages. In complex porous media separation of scales exists, leading to flow and transport signatures that need to be described by multiple functions with distinct flow field and transport characteristics. Reference: Lin, Q., Al-Khulaifi Y., Blunt, M.J. and Bijeljic B. (2016). Advances in Water Resources, 96, 306-322, doi:10.1016/j.advwatres.2016.08.002.

Pore scale modelling of electrical and hydraulic properties of a semi-consolidated sandstone under unsaturated conditions

NASA Astrophysics Data System (ADS)

Cassiani, G.; dalla, E.; Brovelli, A.; Pitea, D.; Binley, A. M.

2003-04-01

The development of reliable constitutive laws to translate geophysical properties into hydrological ones is the fundamental step for successful applications of hydrogeophysical techniques. Many such laws have been proposed and applied, particularly with regard to two types of relationships: (a) between moisture content and dielectric properties, and (b) between electrical resistivity, rock structure and water saturation. The classical Archie's law belongs to this latter category. Archie's relationship has been widely used, starting from borehole logs applications, to translate geoelectrical measurements into estimates of saturation. However, in spite of its popularity, it remains an empirical relationship, the parameters of which must be calibrated case by case, e.g. on laboratory data. Pore-scale models have been recently recognized and used as powerful tools to investigate the constitutive relations of multiphase soils from a pore-scale point of view, because they bridge the microscopic and macroscopic scales. In this project, we develop and validate a three-dimensional pore-scale method to compute electrical properties of unsaturated and saturated porous media. First we simulate a random packing of spheres [1] that obeys the grain-size distribution and porosity of an experimental porous medium system; then we simulate primary drainage with a morphological approach [2]; finally, for each state of saturation during the drainage process, we solve the electrical conduction equation within the grain structure with a new numerical model and compute the apparent electrical resistivity of the porous medium. We apply the new method to a semi-consolidated Permo-Triassic Sandstone from the UK (Sherwood Sandstone) for which both pressure-saturation (Van Genuchten) and Archie's law parameters have been measured on laboratory samples. A comparison between simulated and measured relationships has been performed.

Pore-Scale Hydrodynamics in a Progressively Bioclogged Three-Dimensional Porous Medium: 3-D Particle Tracking Experiments and Stochastic Transport Modeling

NASA Astrophysics Data System (ADS)

Carrel, M.; Morales, V. L.; Dentz, M.; Derlon, N.; Morgenroth, E.; Holzner, M.

2018-03-01

Biofilms are ubiquitous bacterial communities that grow in various porous media including soils, trickling, and sand filters. In these environments, they play a central role in services ranging from degradation of pollutants to water purification. Biofilms dynamically change the pore structure of the medium through selective clogging of pores, a process known as bioclogging. This affects how solutes are transported and spread through the porous matrix, but the temporal changes to transport behavior during bioclogging are not well understood. To address this uncertainty, we experimentally study the hydrodynamic changes of a transparent 3-D porous medium as it experiences progressive bioclogging. Statistical analyses of the system's hydrodynamics at four time points of bioclogging (0, 24, 36, and 48 h in the exponential growth phase) reveal exponential increases in both average and variance of the flow velocity, as well as its correlation length. Measurements for spreading, as mean-squared displacements, are found to be non-Fickian and more intensely superdiffusive with progressive bioclogging, indicating the formation of preferential flow pathways and stagnation zones. A gamma distribution describes well the Lagrangian velocity distributions and provides parameters that quantify changes to the flow, which evolves from a parallel pore arrangement under unclogged conditions, toward a more serial arrangement with increasing clogging. Exponentially evolving hydrodynamic metrics agree with an exponential bacterial growth phase and are used to parameterize a correlated continuous time random walk model with a stochastic velocity relaxation. The model accurately reproduces transport observations and can be used to resolve transport behavior at intermediate time points within the exponential growth phase considered.

Unifying Pore Network Modeling, Continuous Time Random Walk (CTRW) Theory and Experiment to Describe Impact of Spatial Heterogeneities on Solute Dispersion at Multiple Length-scales

NASA Astrophysics Data System (ADS)

Bijeljic, B.; Blunt, M. J.; Rhodes, M. E.

2009-04-01

This talk will describe and highlight the advantages offered by a novel methodology that unifies pore network modeling, CTRW theory and experiment in description of solute dispersion in porous media. Solute transport in a porous medium is characterized by the interplay of advection and diffusion (described by Peclet number, Pe) that cause dispersion of solute particles. Dispersion is traditionally described by dispersion coefficients, D, that are commonly calculated from the spatial moments of the plume. Using a pore-scale network model based on particle tracking, the rich Peclet-number dependence of dispersion coefficient is predicted from first principles and is shown to compare well with experimental data for restricted diffusion, transition, power-law and mechanical dispersion regimes in the asymptotic limit. In the asymptotic limit D is constant and can be used in an averaged advection-dispersion equation. However, it is highly important to recognize that, until the velocity field is fully sampled, the particle transport is non-Gaussian and D possesses temporal or spatial variation. Furthermore, temporal probability density functions (PDF) of tracer particles are studied in pore networks and an excellent agreement for the spectrum of transition times for particles from pore to pore is obtained between network model results and CTRW theory. Based on the truncated power-law interpretation of PDF-s, the physical origin of the power-law scaling of dispersion coefficient vs. Peclet number has been explained for unconsolidated porous media, sands and a number of sandstones, arriving at the same conclusion from numerical network modelling, analytic CTRW theory and experiment. The length traveled by solute plumes before Gaussian behaviour is reached increases with an increase in heterogeneity and/or Pe. This opens up the question on the nature of dispersion in natural systems where the heterogeneities at the larger scales will significantly increase the range of velocities in the reservoir, thus significantly delaying the asymptotic approach to Gaussian behaviour. As a consequence, the asymptotic behaviour might not be reached at the field scale. This is illustrated by the multi-scale approach in which transport at core, gridblock and field scale is viewed as a series of particle transitions between discrete nodes governed by probability distributions. At each scale of interest a distribution that represents transport physics (and the heterogeneity) is used as an input to model a subsequent reservoir scale. The extensions to reactive transport are discussed.

Water adsorption isotherms on porous onionlike carbonaceous particles. Simulations with the grand canonical Monte Carlo method

NASA Astrophysics Data System (ADS)

Hantal, György; Picaud, Sylvain; Hoang, Paul N. M.; Voloshin, Vladimir P.; Medvedev, Nikolai N.; Jedlovszky, Pál

2010-10-01

The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water molecules on different types of model soot particles. These soot models are constructed by first removing atoms from onion-fullerene structures in order to create randomly distributed pores inside the soot, and then performing molecular dynamics simulations, based on the reactive adaptive intermolecular reactive empirical bond order (AIREBO) description of the interaction between carbon atoms, to optimize the resulting structures. The obtained results clearly show that the main driving force of water adsorption on soot is the possibility of the formation of new water-water hydrogen bonds with the already adsorbed water molecules. The shape of the calculated water adsorption isotherms at 298 K strongly depends on the possible confinement of the water molecules in pores of the carbonaceous structure. We found that there are two important factors influencing the adsorption ability of soot. The first of these factors, dominating at low pressures, is the ability of the soot of accommodating the first adsorbed water molecules at strongly hydrophilic sites. The second factor concerns the size and shape of the pores, which should be such that the hydrogen bonding network of the water molecules filling them should be optimal. This second factor determines the adsorption properties at higher pressures.

Simulation of solute transport through heterogeneous networks: analysis using the method of moments and the statistics of local transport characteristics.

PubMed

Li, Min; Qi, Tao; Bernabé, Yves; Zhao, Jinzhou; Wang, Ying; Wang, Dong; Wang, Zheming

2018-02-28

We used a time domain random walk approach to simulate passive solute transport in networks. In individual pores, solute transport was modeled as a combination of Poiseuille flow and Taylor dispersion. The solute plume data were interpreted via the method of moments. Analysis of the first and second moments showed that the longitudinal dispersivity increased with increasing coefficient of variation of the pore radii CV and decreasing pore coordination number Z. The third moment was negative and its magnitude grew linearly with time, meaning that the simulated dispersion was intrinsically non-Fickian. The statistics of the Eulerian mean fluid velocities [Formula: see text], the Taylor dispersion coefficients [Formula: see text] and the transit times [Formula: see text] were very complex and strongly affected by CV and Z. In particular, the probability of occurrence of negative velocities grew with increasing CV and decreasing Z. Hence, backward and forward transit times had to be distinguished. The high-τ branch of the transit-times probability curves had a power law form associated to non-Fickian behavior. However, the exponent was insensitive to pore connectivity, although variations of Z affected the third moment growth. Thus, we conclude that both the low- and high-τ branches played a role in generating the observed non-Fickian behavior.

Uncertainty Quantification of Nonlinear Electrokinetic Response in a Microchannel-Membrane Junction

NASA Astrophysics Data System (ADS)

Alizadeh, Shima; Iaccarino, Gianluca; Mani, Ali

2015-11-01

We have conducted uncertainty quantification (UQ) for electrokinetic transport of ionic species through a hybrid microfluidic system using different probabilistic techniques. The system of interest is an H-configuration consisting of two parallel microchannels that are connected via a nafion junction. This system is commonly used for ion preconcentration and stacking by utilizing a nonlinear response at the channel-nafion junction that leads to deionization shocks. In this work, the nafion medium is modeled as many parallel nano-pores where, the nano-pore diameter, nafion porosity, and surface charge density are independent random variables. We evaluated the resulting uncertainty on the ion concentration fields as well as the deionization shock location. The UQ methods predicted consistent statistics for the outputs and the results revealed that the shock location is weakly sensitive to the nano-pore surface charge and primarily driven by nano-pore diameters. The present study can inform the design of electrokinetic networks with increased robustness to natural manufacturing variability. Applications include water desalination and lab-on-a-chip systems. Shima is a graduate student in the department of Mechanical Engineering at Stanford University. She received her Master's degree from Stanford in 2011. Her research interests include Electrokinetics in porous structures and high performance computing.

Anomalously slow relaxation of the system of liquid clusters in a disordered nanoporous medium according to the self-organized criticality scenario

NASA Astrophysics Data System (ADS)

Borman, V. D.; Tronin, V. N.; Byrkin, V. A.

2016-04-01

We propose a physical model of a relaxation of states of clusters of nonwetting liquid confined in a random nanoporous medium. The relaxation is occurred by the self-organized criticality (SOC) scenario. Process is characterized by waiting for fluctuation necessary for overcoming of a local energy barrier with the subsequent avalanche hydrodynamic extrusion of the liquid by surface forces of the nonwetting frame. The dependence of the interaction between local configurations on the number of filled pores belonging to the infinite percolation cluster of filled pores serves as an internal feedback initiating the SOC process. The calculations give a power-law time dependence of the relative volume θ of the confined liquid θ ∼t-ν (ν ∼ 0.2) as in the picture of relaxation in the mean field approximation. The model of the relaxation of the porous medium with the nonwetting liquid demonstrates possible mechanisms and scenarios of SOC for relaxation of other disordered systems.

Lattice Boltzmann simulation of the gas-solid adsorption process in reconstructed random porous media.

PubMed

Zhou, L; Qu, Z G; Ding, T; Miao, J Y

2016-04-01

The gas-solid adsorption process in reconstructed random porous media is numerically studied with the lattice Boltzmann (LB) method at the pore scale with consideration of interparticle, interfacial, and intraparticle mass transfer performances. Adsorbent structures are reconstructed in two dimensions by employing the quartet structure generation set approach. To implement boundary conditions accurately, all the porous interfacial nodes are recognized and classified into 14 types using a proposed universal program called the boundary recognition and classification program. The multiple-relaxation-time LB model and single-relaxation-time LB model are adopted to simulate flow and mass transport, respectively. The interparticle, interfacial, and intraparticle mass transfer capacities are evaluated with the permeability factor and interparticle transfer coefficient, Langmuir adsorption kinetics, and the solid diffusion model, respectively. Adsorption processes are performed in two groups of adsorbent media with different porosities and particle sizes. External and internal mass transfer resistances govern the adsorption system. A large porosity leads to an early time for adsorption equilibrium because of the controlling factor of external resistance. External and internal resistances are dominant at small and large particle sizes, respectively. Particle size, under which the total resistance is minimum, ranges from 3 to 7 μm with the preset parameters. Pore-scale simulation clearly explains the effect of both external and internal mass transfer resistances. The present paper provides both theoretical and practical guidance for the design and optimization of adsorption systems.

Lattice Boltzmann simulation of the gas-solid adsorption process in reconstructed random porous media

NASA Astrophysics Data System (ADS)

Zhou, L.; Qu, Z. G.; Ding, T.; Miao, J. Y.

2016-04-01

The gas-solid adsorption process in reconstructed random porous media is numerically studied with the lattice Boltzmann (LB) method at the pore scale with consideration of interparticle, interfacial, and intraparticle mass transfer performances. Adsorbent structures are reconstructed in two dimensions by employing the quartet structure generation set approach. To implement boundary conditions accurately, all the porous interfacial nodes are recognized and classified into 14 types using a proposed universal program called the boundary recognition and classification program. The multiple-relaxation-time LB model and single-relaxation-time LB model are adopted to simulate flow and mass transport, respectively. The interparticle, interfacial, and intraparticle mass transfer capacities are evaluated with the permeability factor and interparticle transfer coefficient, Langmuir adsorption kinetics, and the solid diffusion model, respectively. Adsorption processes are performed in two groups of adsorbent media with different porosities and particle sizes. External and internal mass transfer resistances govern the adsorption system. A large porosity leads to an early time for adsorption equilibrium because of the controlling factor of external resistance. External and internal resistances are dominant at small and large particle sizes, respectively. Particle size, under which the total resistance is minimum, ranges from 3 to 7 μm with the preset parameters. Pore-scale simulation clearly explains the effect of both external and internal mass transfer resistances. The present paper provides both theoretical and practical guidance for the design and optimization of adsorption systems.

Percolation Laws of a Fractal Fracture-Pore Double Medium

NASA Astrophysics Data System (ADS)

Zhao, Yangsheng; Feng, Zengchao; Lv, Zhaoxing; Zhao, Dong; Liang, Weiguo

2016-12-01

The fracture-pore double porosity medium is one of the most common media in nature, for example, rock mass in strata. Fracture has a more significant effect on fluid flow than a pore in a fracture-pore double porosity medium. Hence, the fracture effect on percolation should be considered when studying the percolation phenomenon in porous media. In this paper, based on the fractal distribution law, three-dimensional (3D) fracture surfaces, and two-dimensional (2D) fracture traces in rock mass, the locations of fracture surfaces or traces are determined using a random function of uniform distribution. Pores are superimposed to build a fractal fracture-pore double medium. Numerical experiments were performed to show percolation phenomena in the fracture-pore double medium. The percolation threshold can be determined from three independent variables (porosity n, fracture fractal dimension D, and initial value of fracture number N0). Once any two are determined, the percolation probability exists at a critical point with the remaining parameter changing. When the initial value of the fracture number is greater than zero, the percolation threshold in the fracture-pore medium is much smaller than that in a pore medium. When the fracture number equals zero, the fracture-pore medium degenerates to a pore medium, and both percolation thresholds are the same.

Three-dimensional Mesoscale Simulations of Detonation Initiation in Energetic Materials with Density-based Kinetics

NASA Astrophysics Data System (ADS)

Jackson, Thomas; Jost, A. M.; Zhang, Ju; Sridharan, P.; Amadio, G.

2017-06-01

In this work we present three-dimensional mesoscale simulations of detonation initiation in energetic materials. We solve the reactive Euler equations, with the energy equation augmented by a power deposition term. The reaction rate at the mesoscale is modelled using a density-based kinetics scheme, adapted from standard Ignition and Growth models. The deposition term is based on previous results of simulations of pore collapse at the microscale, modelled at the mesoscale as hot-spots. We carry out three-dimensional mesoscale simulations of random packs of HMX crystals in a binder, and show that the transition between no-detonation and detonation depends on the number density of the hot-spots, the initial radius of the hot-spot, the post-shock pressure of an imposed shock, and the amplitude of the power deposition term. The trends of transition at lower pressure of the imposed shock for larger number density of pore observed in experiments is reproduced. Initial attempts to improve the agreement between the simulation and experiments through calibration of various parameters will also be made.

Upscaling transport of a reacting solute through a peridocially converging-diverging channel at pre-asymptotic times

NASA Astrophysics Data System (ADS)

Sund, Nicole L.; Bolster, Diogo; Dawson, Clint

2015-11-01

In this study we extend the Spatial Markov model, which has been successfully used to upscale conservative transport across a diverse range of porous media flows, to test if it can accurately upscale reactive transport, defined by a spatially heterogeneous first order degradation rate. We test the model in a well known highly simplified geometry, commonly considered as an idealized pore or fracture structure, a periodic channel with wavy boundaries. The edges of the flow domain have a layer through which there is no flow, but in which diffusion of a solute still occurs. Reactions are confined to this region. We demonstrate that the Spatial Markov model, an upscaled random walk model that enforces correlation between successive jumps, can reproduce breakthrough curves measured from microscale simulations that explicitly resolve all pertinent processes. We also demonstrate that a similar random walk model that does not enforce successive correlations is unable to reproduce all features of the measured breakthrough curves.

Validation of model predictions of pore-scale fluid distributions during two-phase flow

NASA Astrophysics Data System (ADS)

Bultreys, Tom; Lin, Qingyang; Gao, Ying; Raeini, Ali Q.; AlRatrout, Ahmed; Bijeljic, Branko; Blunt, Martin J.

2018-05-01

Pore-scale two-phase flow modeling is an important technology to study a rock's relative permeability behavior. To investigate if these models are predictive, the calculated pore-scale fluid distributions which determine the relative permeability need to be validated. In this work, we introduce a methodology to quantitatively compare models to experimental fluid distributions in flow experiments visualized with microcomputed tomography. First, we analyzed five repeated drainage-imbibition experiments on a single sample. In these experiments, the exact fluid distributions were not fully repeatable on a pore-by-pore basis, while the global properties of the fluid distribution were. Then two fractional flow experiments were used to validate a quasistatic pore network model. The model correctly predicted the fluid present in more than 75% of pores and throats in drainage and imbibition. To quantify what this means for the relevant global properties of the fluid distribution, we compare the main flow paths and the connectivity across the different pore sizes in the modeled and experimental fluid distributions. These essential topology characteristics matched well for drainage simulations, but not for imbibition. This suggests that the pore-filling rules in the network model we used need to be improved to make reliable predictions of imbibition. The presented analysis illustrates the potential of our methodology to systematically and robustly test two-phase flow models to aid in model development and calibration.

Pore-Scale Modeling of Pore Structure Effects on P-Wave Scattering Attenuation in Dry Rocks

PubMed Central

Li, Tianyang; Qiu, Hao; Wang, Feifei

2015-01-01

Underground rocks usually have complex pore system with a variety of pore types and a wide range of pore size. The effects of pore structure on elastic wave attenuation cannot be neglected. We investigated the pore structure effects on P-wave scattering attenuation in dry rocks by pore-scale modeling based on the wave theory and the similarity principle. Our modeling results indicate that pore size, pore shape (such as aspect ratio), and pore density are important factors influencing P-wave scattering attenuation in porous rocks, and can explain the variation of scattering attenuation at the same porosity. From the perspective of scattering attenuation, porous rocks can safely suit to the long wavelength assumption when the ratio of wavelength to pore size is larger than 15. Under the long wavelength condition, the scattering attenuation coefficient increases as a power function as the pore density increases, and it increases exponentially with the increase in aspect ratio. For a certain porosity, rocks with smaller aspect ratio and/or larger pore size have stronger scattering attenuation. When the pore aspect ratio is larger than 0.5, the variation of scattering attenuation at the same porosity is dominantly caused by pore size and almost independent of the pore aspect ratio. These results lay a foundation for pore structure inversion from elastic wave responses in porous rocks. PMID:25961729

Numerical simulation on hydromechanical coupling in porous media adopting three-dimensional pore-scale model.

PubMed

Liu, Jianjun; Song, Rui; Cui, Mengmeng

2014-01-01

A novel approach of simulating hydromechanical coupling in pore-scale models of porous media is presented in this paper. Parameters of the sandstone samples, such as the stress-strain curve, Poisson's ratio, and permeability under different pore pressure and confining pressure, are tested in laboratory scale. The micro-CT scanner is employed to scan the samples for three-dimensional images, as input to construct the model. Accordingly, four physical models possessing the same pore and rock matrix characteristics as the natural sandstones are developed. Based on the micro-CT images, the three-dimensional finite element models of both rock matrix and pore space are established by MIMICS and ICEM software platform. Navier-Stokes equation and elastic constitutive equation are used as the mathematical model for simulation. A hydromechanical coupling analysis in pore-scale finite element model of porous media is simulated by ANSYS and CFX software. Hereby, permeability of sandstone samples under different pore pressure and confining pressure has been predicted. The simulation results agree well with the benchmark data. Through reproducing its stress state underground, the prediction accuracy of the porous rock permeability in pore-scale simulation is promoted. Consequently, the effects of pore pressure and confining pressure on permeability are revealed from the microscopic view.

Numerical Simulation on Hydromechanical Coupling in Porous Media Adopting Three-Dimensional Pore-Scale Model

PubMed Central

Liu, Jianjun; Song, Rui; Cui, Mengmeng

2014-01-01

A novel approach of simulating hydromechanical coupling in pore-scale models of porous media is presented in this paper. Parameters of the sandstone samples, such as the stress-strain curve, Poisson's ratio, and permeability under different pore pressure and confining pressure, are tested in laboratory scale. The micro-CT scanner is employed to scan the samples for three-dimensional images, as input to construct the model. Accordingly, four physical models possessing the same pore and rock matrix characteristics as the natural sandstones are developed. Based on the micro-CT images, the three-dimensional finite element models of both rock matrix and pore space are established by MIMICS and ICEM software platform. Navier-Stokes equation and elastic constitutive equation are used as the mathematical model for simulation. A hydromechanical coupling analysis in pore-scale finite element model of porous media is simulated by ANSYS and CFX software. Hereby, permeability of sandstone samples under different pore pressure and confining pressure has been predicted. The simulation results agree well with the benchmark data. Through reproducing its stress state underground, the prediction accuracy of the porous rock permeability in pore-scale simulation is promoted. Consequently, the effects of pore pressure and confining pressure on permeability are revealed from the microscopic view. PMID:24955384

Upscaling pore pressure-dependent gas permeability in shales

NASA Astrophysics Data System (ADS)

Ghanbarian, Behzad; Javadpour, Farzam

2017-04-01

Upscaling pore pressure dependence of shale gas permeability is of great importance and interest in the investigation of gas production in unconventional reservoirs. In this study, we apply the Effective Medium Approximation, an upscaling technique from statistical physics, and modify the Doyen model for unconventional rocks. We develop an upscaling model to estimate the pore pressure-dependent gas permeability from pore throat size distribution, pore connectivity, tortuosity, porosity, and gas characteristics. We compare our adapted model with six data sets: three experiments, one pore-network model, and two lattice-Boltzmann simulations. Results showed that the proposed model estimated the gas permeability within a factor of 3 of the measurements/simulations in all data sets except the Eagle Ford experiment for which we discuss plausible sources of discrepancies.

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Xiaofan; Mehmani, Yashar; Perkins, William A.

2016-09-01

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include methods that 1) explicitly model the three-dimensional geometry of pore spaces and 2) those that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of class 1, based on direct numerical simulation using computational fluid dynamics (CFD) codes, against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of class 1 based on the immersed-boundary method (IMB),more » lattice Boltzmann method (LBM), smoothed particle hydrodynamics (SPH), as well as a model of class 2 (a pore-network model or PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results with previously reported experimental observations. Experimental observations are limited to measured pore-scale velocities, so solute transport comparisons are made only among the various models. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations).« less

Can the three pore model correctly describe peritoneal transport of protein?

PubMed

Waniewski, Jacek; Poleszczuk, Jan; Antosiewicz, Stefan; Baczynński, Daniel; Gałach, Magda; Pietribiasi, Mauro; Wanńkowicz, Zofia

2014-01-01

The three pore model (3PM) includes large pores for the description of protein leak to the peritoneal cavity during peritoneal dialysis. However, the reliability of this description has been not fully tested against clinical data yet. Peritoneal transport parameters were estimated using 3PM, extended 3p model (with estimation of fraction of large pores, ext3PM), ext3PM with modified size of pores and proteins (mext3PM), and simplified two pore (2PM, small and ultrasmall pores) models for 32 patients on peritoneal dialysis investigated using the sequential peritoneal equilibration test (consecutive peritoneal equilibration test [PET]: glucose 2.27%, 4 h, and miniPET: glucose 3.86%, 1 h). Urea, creatinine, glucose, sodium, phosphate, albumin, and IgM concentrations were measured in dialysis fluid and plasma. Ext3PM and mext3PM, with large pore fraction of about 0.14, provided a good description of fluid and small solute kinetics, but their predictions for albumin transport were less accurate. Two pore model precisely described the data on fluid and small solute transport. The 3p models could not describe the diffusive-convective transport of albumin as precisely as the transport of fluid, small solutes, and IgM. The 2p model (not applicable for proteins) was an efficient tool for modeling fluid and small solute transport.

Parallel capillary-tube-based extension of thermoacoustic theory for random porous media.

PubMed

Roh, Heui-Seol; Raspet, Richard; Bass, Henry E

2007-03-01

Thermoacoustic theory is extended to stacks made of random bulk media. Characteristics of the porous stack such as the tortuosity and dynamic shape factors are introduced into the thermoacoustic wave equation in the low reduced frequency approximation. Basic thermoacoustic equations for a bulk porous medium are formulated analogously to the equations for a single pore. Use of different dynamic shape factors for the viscous and thermal effects is adopted and scaling using the dynamic shape factors and tortuosity is demonstrated. Comparisons of the calculated and experimentally derived thermoacoustic properties of reticulated vitreous carbon and aluminum foam show good agreement. A consistent mathematical model of sound propagation in a random porous medium with an imposed temperature is developed. This treatment leads to an expression for the coefficient of the temperature gradient in terms of scaled cylindrical thermoviscous functions.

Review of pore network modelling of porous media: Experimental characterisations, network constructions and applications to reactive transport.

PubMed

Xiong, Qingrong; Baychev, Todor G; Jivkov, Andrey P

2016-09-01

Pore network models have been applied widely for simulating a variety of different physical and chemical processes, including phase exchange, non-Newtonian displacement, non-Darcy flow, reactive transport and thermodynamically consistent oil layers. The realism of such modelling, i.e. the credibility of their predictions, depends to a large extent on the quality of the correspondence between the pore space of a given medium and the pore network constructed as its representation. The main experimental techniques for pore space characterisation, including direct imaging, mercury intrusion porosimetry and gas adsorption, are firstly summarised. A review of the main pore network construction techniques is then presented. Particular focus is given on how such constructions are adapted to the data from experimentally characterised pore systems. Current applications of pore network models are considered, with special emphasis on the effects of adsorption, dissolution and precipitation, as well as biomass growth, on transport coefficients. Pore network models are found to be a valuable tool for understanding and predicting meso-scale phenomena, linking single pore processes, where other techniques are more accurate, and the homogenised continuum porous media, used by engineering community. Copyright © 2016 The Authors. Published by Elsevier B.V. All rights reserved.

Comparison of Effective Medium Schemes For Seismic Velocities in Cracked Anisotropic Rock

NASA Astrophysics Data System (ADS)

Morshed, S.; Chesnokov, E.

2017-12-01

Understanding of elastic properties of reservoir rock is necessary for meaningful interpretation and analysis of seismic measurements. The elastic properties of a rock are controlled by the microstructural properties such as mineralogical composition, pore and crack distribution, texture and pore connectivity. However, seismic scale is much larger than microstructure scale. Understanding of macroscopic properties at relevant seismic scale (e.g. borehole sonic data) comes from effective medium theory (EMT). However, most of the effective medium theories fail at high crack density as the interactions of strain fields of the cracks can't be ignored. We compare major EMT schemes from low to high crack density. While at low crack density all method gives similar results, at high crack density they differ significantly. Then, we focus on generalized singular approximation (GSA) and effective field (EF) method as they allow cracks beyond the limit of dilute concentrations. Additionally, we use grain contact (GC) method to examine the stiffness constants of the rock matrix. We prepare simple models of a multiphase media containing low to high concentrations of isolated pores. Randomly oriented spherical pores and horizontally oriented ellipsoidal (aspect ratio =0.1) pores have been considered. For isolated spherical pores, all the three methods show exactly same or similar results. However, inclusion interactions are different in different directions in case of horizontal ellipsoidal pores and individual stiffness constants differ greatly from one method to another at different crack density. Stiffness constants remain consistent in GSA method whereas some components become unusual in EF method at a higher crack density (>0.15). Finally, we applied GSA method to interpret ultrasonic velocities of core samples. Mineralogical composition from X-ray diffraction (XRD) data and lab measured porosity data have been utilized. Both compressional and shear wave velocities from GSA method show good fit with the lab measured velocities.

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE PAGES

Mehmani, Yashar; Schoenherr, Martin; Pasquali, Andrea; ...

2015-09-28

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This paper provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.« less

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Xiaofan; Mehmani, Yashar; Perkins, William A.

2016-09-01

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.« less

Permeability-Selectivity Analysis of Microfiltration and Ultrafiltration Membranes: Effect of Pore Size and Shape Distribution and Membrane Stretching.

PubMed

Siddiqui, Muhammad Usama; Arif, Abul Fazal Muhammad; Bashmal, Salem

2016-08-06

We present a modeling approach to determine the permeability-selectivity tradeoff for microfiltration and ultrafiltration membranes with a distribution of pore sizes and pore shapes. Using the formulated permeability-selectivity model, the effect of pore aspect ratio and pore size distribution on the permeability-selectivity tradeoff of the membrane is analyzed. A finite element model is developed to study the effect of membrane stretching on the distribution of pore sizes and shapes in the stretched membrane. The effect of membrane stretching on the permeability-selectivity tradeoff of membranes is also analyzed. The results show that increasing pore aspect ratio improves membrane performance while increasing the width of pore size distribution deteriorates the performance. It was also found that the effect of membrane stretching on the permeability-selectivity tradeoff is greatly affected by the uniformity of pore distribution in the membrane. Stretching showed a positive shift in the permeability-selectivity tradeoff curve of membranes with well-dispersed pores while in the case of pore clustering, a negative shift in the permeability-selectivity tradeoff curve was observed.

The Application of Fractal and Multifractal Theory in Hydraulic-Flow-Unit Characterization and Permeability Estimation

NASA Astrophysics Data System (ADS)

Chen, X.; Yao, G.; Cai, J.

2017-12-01

Pore structure characteristics are important factors in influencing the fluid transport behavior of porous media, such as pore-throat ratio, pore connectivity and size distribution, moreover, wettability. To accurately characterize the diversity of pore structure among HFUs, five samples selected from different HFUs (porosities are approximately equal, however permeability varies widely) were chosen to conduct micro-computerized tomography test to acquire direct 3D images of pore geometries and to perform mercury injection experiments to obtain the pore volume-radii distribution. To characterize complex and high nonlinear pore structure of all samples, three classic fractal geometry models were applied. Results showed that each HFU has similar box-counting fractal dimension and generalized fractal dimension in the number-area model, but there are significant differences in multifractal spectrums. In the radius-volume model, there are three obvious linear segments, corresponding to three fractal dimension values, and the middle one is proved as the actual fractal dimension according to the maximum radius. In the number-radius model, the spherical-pore size distribution extracted by maximum ball algorithm exist a decrease in the number of small pores compared with the fractal power rate rather than the traditional linear law. Among the three models, only multifractal analysis can classify the HFUs accurately. Additionally, due to the tightness and low-permeability in reservoir rocks, connate water film existing in the inner surface of pore channels commonly forms bound water. The conventional model which is known as Yu-Cheng's model has been proved to be typically not applicable. Considering the effect of irreducible water saturation, an improved fractal permeability model was also deduced theoretically. The comparison results showed that the improved model can be applied to calculate permeability directly and accurately in such unconventional rocks.

Long-Term Transport of Cryptosporidium Parvum

NASA Astrophysics Data System (ADS)

Andrea, C.; Harter, T.; Hou, L.; Atwill, E. R.; Packman, A.; Woodrow-Mumford, K.; Maldonado, S.

2005-12-01

The protozoan pathogen Cryptosporidium parvum is a leading cause of waterborne disease. Subsurface transport and filtration in natural and artificial porous media are important components of the environmental pathway of this pathogen. It has been shown that the oocysts of C. parvum show distinct colloidal properties. We conducted a series of laboratory studies on sand columns (column length: 10 cm - 60 cm, flow rates: 0.7 m/d - 30 m/d, ionic strength: 0.01 - 100 mM, filter grain size: 0.2 - 2 mm, various solution chemistry). Breakthrough curves were measured over relatively long time-periods (hundreds to thousands of pore volumes). We show that classic colloid filtration theory is a reasonable tool for predicting the initial breakthrough, but it is inadequate to explain the significant tailing observed in the breakthrough of C. parvum oocyst through sand columns. We discuss the application of the Continuous Time Random Walk approach to account for the strong tailing that was observed in our experiments. The CTRW is generalized transport modeling framework, which includes the classic advection-dispersion equation (ADE), the fractional ADE, and the multi-rate mass transfer model as special cases. Within this conceptual framework, it is possible to distinguish between the contributions of pore-scale geometrical (physical) disorder and of pore-scale physico-chemical heterogeneities (e.g., of the filtration, sorption, desorption processes) to the transport of C. parvum oocysts.

Modeling of Processing-Induced Pore Morphology in an Additively-Manufactured Ti-6Al-4V Alloy

PubMed Central

Kabir, Mohammad Rizviul; Richter, Henning

2017-01-01

A selective laser melting (SLM)-based, additively-manufactured Ti-6Al-4V alloy is prone to the accumulation of undesirable defects during layer-by-layer material build-up. Defects in the form of complex-shaped pores are one of the critical issues that need to be considered during the processing of this alloy. Depending on the process parameters, pores with concave or convex boundaries may occur. To exploit the full potential of additively-manufactured Ti-6Al-4V, the interdependency between the process parameters, pore morphology, and resultant mechanical properties, needs to be understood. By incorporating morphological details into numerical models for micromechanical analyses, an in-depth understanding of how these pores interact with the Ti-6Al-4V microstructure can be gained. However, available models for pore analysis lack a realistic description of both the Ti-6Al-4V grain microstructure, and the pore geometry. To overcome this, we propose a comprehensive approach for modeling and discretizing pores with complex geometry, situated in a polycrystalline microstructure. In this approach, the polycrystalline microstructure is modeled by means of Voronoi tessellations, and the complex pore geometry is approximated by strategically combining overlapping spheres of varied sizes. The proposed approach provides an elegant way to model the microstructure of SLM-processed Ti-6Al-4V containing pores or crack-like voids, and makes it possible to investigate the relationship between process parameters, pore morphology, and resultant mechanical properties in a finite-element-based simulation framework. PMID:28772504

Modeling of Processing-Induced Pore Morphology in an Additively-Manufactured Ti-6Al-4V Alloy.

PubMed

Kabir, Mohammad Rizviul; Richter, Henning

2017-02-08

A selective laser melting (SLM)-based, additively-manufactured Ti-6Al-4V alloy is prone to the accumulation of undesirable defects during layer-by-layer material build-up. Defects in the form of complex-shaped pores are one of the critical issues that need to be considered during the processing of this alloy. Depending on the process parameters, pores with concave or convex boundaries may occur. To exploit the full potential of additively-manufactured Ti-6Al-4V, the interdependency between the process parameters, pore morphology, and resultant mechanical properties, needs to be understood. By incorporating morphological details into numerical models for micromechanical analyses, an in-depth understanding of how these pores interact with the Ti-6Al-4V microstructure can be gained. However, available models for pore analysis lack a realistic description of both the Ti-6Al-4V grain microstructure, and the pore geometry. To overcome this, we propose a comprehensive approach for modeling and discretizing pores with complex geometry, situated in a polycrystalline microstructure. In this approach, the polycrystalline microstructure is modeled by means of Voronoi tessellations, and the complex pore geometry is approximated by strategically combining overlapping spheres of varied sizes. The proposed approach provides an elegant way to model the microstructure of SLM-processed Ti-6Al-4V containing pores or crack-like voids, and makes it possible to investigate the relationship between process parameters, pore morphology, and resultant mechanical properties in a finite-element-based simulation framework.

New general pore size distribution model by classical thermodynamics application: Activated carbon

USGS Publications Warehouse

Lordgooei, M.; Rood, M.J.; Rostam-Abadi, M.

2001-01-01

A model is developed using classical thermodynamics to characterize pore size distributions (PSDs) of materials containing micropores and mesopores. The thermal equation of equilibrium adsorption (TEEA) is used to provide thermodynamic properties and relate the relative pore filling pressure of vapors to the characteristic pore energies of the adsorbent/adsorbate system for micropore sizes. Pore characteristic energies are calculated by averaging of interaction energies between adsorbate molecules and adsorbent pore walls as well as considering adsorbate-adsorbate interactions. A modified Kelvin equation is used to characterize mesopore sizes by considering variation of the adsorbate surface tension and by excluding the adsorbed film layer for the pore size. The modified-Kelvin equation provides similar pore filling pressures as predicted by density functional theory. Combination of these models provides a complete PSD of the adsorbent for the micropores and mesopores. The resulting PSD is compared with the PSDs from Jaroniec and Choma and Horvath and Kawazoe models as well as a first-order approximation model using Polanyi theory. The major importance of this model is its basis on classical thermodynamic properties, less simplifying assumptions in its derivation compared to other methods, and ease of use.

Comparison of NMR simulations of porous media derived from analytical and voxelized representations.

PubMed

Jin, Guodong; Torres-Verdín, Carlos; Toumelin, Emmanuel

2009-10-01

We develop and compare two formulations of the random-walk method, grain-based and voxel-based, to simulate the nuclear-magnetic-resonance (NMR) response of fluids contained in various models of porous media. The grain-based approach uses a spherical grain pack as input, where the solid surface is analytically defined without an approximation. In the voxel-based approach, the input is a computer-tomography or computer-generated image of reconstructed porous media. Implementation of the two approaches is largely the same, except for the representation of porous media. For comparison, both approaches are applied to various analytical and digitized models of porous media: isolated spherical pore, simple cubic packing of spheres, and random packings of monodisperse and polydisperse spheres. We find that spin magnetization decays much faster in the digitized models than in their analytical counterparts. The difference in decay rate relates to the overestimation of surface area due to the discretization of the sample; it cannot be eliminated even if the voxel size decreases. However, once considering the effect of surface-area increase in the simulation of surface relaxation, good quantitative agreement is found between the two approaches. Different grain or pore shapes entail different rates of increase of surface area, whereupon we emphasize that the value of the "surface-area-corrected" coefficient may not be universal. Using an example of X-ray-CT image of Fontainebleau rock sample, we show that voxel size has a significant effect on the calculated surface area and, therefore, on the numerically simulated magnetization response.

3D Nanoporous Anodic Alumina Structures for Sustained Drug Release

PubMed Central

Xifré-Pérez, Elisabet; Eckstein, Chris; Ferré-Borrull, Josep

2017-01-01

The use of nanoporous anodic alumina (NAA) for the development of drug delivery systems has gained much attention in recent years. The release of drugs loaded inside NAA pores is complex and depends on the morphology of the pores. In this study, NAA, with different three-dimensional (3D) pore structures (cylindrical pores with several pore diameters, multilayered nanofunnels, and multilayered inverted funnels) were fabricated, and their respective drug delivery rates were studied and modeled using doxorubicin as a model drug. The obtained results reveal optimal modeling of all 3D pore structures, differentiating two drug release stages. Thus, an initial short-term and a sustained long-term release were successfully modeled by the Higuchi and the Korsmeyer–Peppas equations, respectively. This study demonstrates the influence of pore geometries on drug release rates, and further presents a sustained long-term drug release that exceeds 60 days without an undesired initial burst. PMID:28825654

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Oostrom, Martinus; Mehmani, Yashar; Romero Gomez, Pedro DJ

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing.more » Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.« less

Preventing Mesh Pore Collapse by Designing Mesh Pores With Auxetic Geometries: A Comprehensive Evaluation Via Computational Modeling.

PubMed

Knight, Katrina M; Moalli, Pamela A; Abramowitch, Steven D

2018-05-01

Pelvic organ prolapse (POP) meshes are exposed to predominately tensile loading conditions in vivo that can lead to pore collapse by 70-90%, decreasing overall porosity and providing a plausible mechanism for the contraction/shrinkage of mesh observed following implantation. To prevent pore collapse, we proposed to design synthetic meshes with a macrostructure that results in auxetic behavior, the pores expand laterally, instead of contracting when loaded. Such behavior can be achieved with a range of auxetic structures/geometries. This study utilized finite element analysis (FEA) to assess the behavior of mesh models with eight auxetic pore geometries subjected to uniaxial loading to evaluate their potential to allow for pore expansion while simultaneously providing resistance to tensile loading. Overall, substituting auxetic geometries for standard pore geometries yielded more pore expansion, but often at the expense of increased model elongation, with two of the eight auxetics not able to maintain pore expansion at higher levels of tension. Meshes with stable pore geometries that remain open with loading will afford the ingrowth of host tissue into the pores and improved integration of the mesh. Given the demonstrated ability of auxetic geometries to allow for pore size maintenance (and pore expansion), auxetically designed meshes have the potential to significantly impact surgical outcomes and decrease the likelihood of major mesh-related complications.

Computational approach to integrate 3D X-ray microtomography and NMR data

NASA Astrophysics Data System (ADS)

Lucas-Oliveira, Everton; Araujo-Ferreira, Arthur G.; Trevizan, Willian A.; Fortulan, Carlos A.; Bonagamba, Tito J.

2018-07-01

Nowadays, most of the efforts in NMR applied to porous media are dedicated to studying the molecular fluid dynamics within and among the pores. These analyses have a higher complexity due to morphology and chemical composition of rocks, besides dynamic effects as restricted diffusion, diffusional coupling, and exchange processes. Since the translational nuclear spin diffusion in a confined geometry (e.g. pores and fractures) requires specific boundary conditions, the theoretical solutions are restricted to some special problems and, in many cases, computational methods are required. The Random Walk Method is a classic way to simulate self-diffusion along a Digital Porous Medium. Bergman model considers the magnetic relaxation process of the fluid molecules by including a probability rate of magnetization survival under surface interactions. Here we propose a statistical approach to correlate surface magnetic relaxivity with the computational method applied to the NMR relaxation in order to elucidate the relationship between simulated relaxation time and pore size of the Digital Porous Medium. The proposed computational method simulates one- and two-dimensional NMR techniques reproducing, for example, longitudinal and transverse relaxation times (T1 and T2, respectively), diffusion coefficients (D), as well as their correlations. For a good approximation between the numerical and experimental results, it is necessary to preserve the complexity of translational diffusion through the microstructures in the digital rocks. Therefore, we use Digital Porous Media obtained by 3D X-ray microtomography. To validate the method, relaxation times of ideal spherical pores were obtained and compared with the previous determinations by the Brownstein-Tarr model, as well as the computational approach proposed by Bergman. Furthermore, simulated and experimental results of synthetic porous media are compared. These results make evident the potential of computational physics in the analysis of the NMR data for complex porous materials.

A streamline splitting pore-network approach for computationally inexpensive and accurate simulation of transport in porous media

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mehmani, Yashar; Oostrom, Martinus; Balhoff, Matthew

2014-03-20

Several approaches have been developed in the literature for solving flow and transport at the pore-scale. Some authors use a direct modeling approach where the fundamental flow and transport equations are solved on the actual pore-space geometry. Such direct modeling, while very accurate, comes at a great computational cost. Network models are computationally more efficient because the pore-space morphology is approximated. Typically, a mixed cell method (MCM) is employed for solving the flow and transport system which assumes pore-level perfect mixing. This assumption is invalid at moderate to high Peclet regimes. In this work, a novel Eulerian perspective on modelingmore » flow and transport at the pore-scale is developed. The new streamline splitting method (SSM) allows for circumventing the pore-level perfect mixing assumption, while maintaining the computational efficiency of pore-network models. SSM was verified with direct simulations and excellent matches were obtained against micromodel experiments across a wide range of pore-structure and fluid-flow parameters. The increase in the computational cost from MCM to SSM is shown to be minimal, while the accuracy of SSM is much higher than that of MCM and comparable to direct modeling approaches. Therefore, SSM can be regarded as an appropriate balance between incorporating detailed physics and controlling computational cost. The truly predictive capability of the model allows for the study of pore-level interactions of fluid flow and transport in different porous materials. In this paper, we apply SSM and MCM to study the effects of pore-level mixing on transverse dispersion in 3D disordered granular media.« less

Reduction behavior and kinetics of vanadium-titanium sinters under high potential oxygen enriched pulverized coal injection

NASA Astrophysics Data System (ADS)

Ma, Jin-fang; Wang, Guang-wei; Zhang, Jian-liang; Li, Xin-yu; Liu, Zheng-jian; Jiao, Ke-xin; Guo, Jian

2017-05-01

In this work, the reduction behavior of vanadium-titanium sinters was studied under five different sets of conditions of pulverized coal injection with oxygen enrichment. The modified random pore model was established to analyze the reduction kinetics. The results show that the reduction rate of sinters was accelerated by an increase of CO and H2 contents. Meanwhile, with the increase in CO and H2 contents, the increasing range of the medium reduction index (MRE) of sinters decreased. The increasing oxygen enrichment ratio played a diminishing role in improving the reduction behavior of the sinters. The reducing process kinetic parameters were solved using the modified random role model. The results indicated that, with increasing oxygen enrichment, the contents of CO and H2 in the reducing gas increased. The reduction activation energy of the sinters decreased to between 20.4 and 23.2 kJ/mol.

Permeability-porosity relationships in sedimentary rocks

USGS Publications Warehouse

Nelson, Philip H.

1994-01-01

In many consolidated sandstone and carbonate formations, plots of core data show that the logarithm of permeability (k) is often linearly proportional to porosity (??). The slope, intercept, and degree of scatter of these log(k)-?? trends vary from formation to formation, and these variations are attributed to differences in initial grain size and sorting, diagenetic history, and compaction history. In unconsolidated sands, better sorting systematically increases both permeability and porosity. In sands and sandstones, an increase in gravel and coarse grain size content causes k to increase even while decreasing ??. Diagenetic minerals in the pore space of sandstones, such as cement and some clay types, tend to decrease log(k) proportionately as ?? decreases. Models to predict permeability from porosity and other measurable rock parameters fall into three classes based on either grain, surface area, or pore dimension considerations. (Models that directly incorporate well log measurements but have no particular theoretical underpinnings from a fourth class.) Grain-based models show permeability proportional to the square of grain size times porosity raised to (roughly) the fifth power, with grain sorting as an additional parameter. Surface-area models show permeability proportional to the inverse square of pore surface area times porosity raised to (roughly) the fourth power; measures of surface area include irreducible water saturation and nuclear magnetic resonance. Pore-dimension models show permeability proportional to the square of a pore dimension times porosity raised to a power of (roughly) two and produce curves of constant pore size that transgress the linear data trends on a log(k)-?? plot. The pore dimension is obtained from mercury injection measurements and is interpreted as the pore opening size of some interconnected fraction of the pore system. The linear log(k)-?? data trends cut the curves of constant pore size from the pore-dimension models, which shows that porosity reduction is always accompanied by a reduction in characteristic pore size. The high powers of porosity of the grain-based and surface-area models are required to compensate for the inclusion of the small end of the pore size spectrum.

Stochastic dynamics of intermittent pore-scale particle motion in three-dimensional porous media

NASA Astrophysics Data System (ADS)

Morales, V. L.; Dentz, M.; Willmann, M.; Holzner, M.

2017-12-01

A proper understanding of velocity dynamics is key for making transport predictions through porous media at any scale. We study the velocity evolution process from particle dynamics at the pore-scale with particular interest in preasymptotic (non-Fickian) behavior. Experimental measurements from 3-dimensional particle tracking velocimetry are used to obtain Lagrangian velocity statistics for three different types of media heterogeneity. Particle velocities are found to be intermittent in nature, log-normally distributed and non-stationary. We show that these velocity characteristics can be captured with a correlated Ornstein-Uhlenbeck process for a random walk in space that is parameterized from velocity distributions. Our simple model is rigorously tested for accurate reproduction of velocity variability in magnitude and frequency. We further show that it captures exceptionally well the preasymptotic mean and mean squared displacement in the ballistic and superdiffusive regimes, and can be extended to determine if and when Fickian behavior will be reached. Our approach reproduces both preasymptotic and asymptotic transport behavior with a single transport model, demonstrating correct description of the fundamental controls of anomalous transport.

Numerical Simulation of Multiphase Flow in Nanoporous Organic Matter With Application to Coal and Gas Shale Systems

NASA Astrophysics Data System (ADS)

Song, Wenhui; Yao, Jun; Ma, Jingsheng; Sun, Hai; Li, Yang; Yang, Yongfei; Zhang, Lei

2018-02-01

Fluid flow in nanoscale organic pores is known to be affected by fluid transport mechanisms and properties within confined pore space. The flow of gas and water shows notably different characteristics compared with conventional continuum modeling approach. A pore network flow model is developed and implemented in this work. A 3-D organic pore network model is constructed from 3-D image that is reconstructed from 2-D shale SEM image of organic-rich sample. The 3-D pore network model is assumed to be gas-wet and to contain initially gas-filled pores only, and the flow model is concerned with drainage process. Gas flow considers a full range of gas transport mechanisms, including viscous flow, Knudsen diffusion, surface diffusion, ad/desorption, and gas PVT and viscosity using a modified van der Waals' EoS and a correlation for natural gas, respectively. The influences of slip length, contact angle, and gas adsorption layer on water flow are considered. Surface tension considers the pore size and temperature effects. Invasion percolation is applied to calculate gas-water relative permeability. The results indicate that the influences of pore pressure and temperature on water phase relative permeabilities are negligible while gas phase relative permeabilities are relatively larger in higher temperatures and lower pore pressures. Gas phase relative permeability increases while water phase relative permeability decreases with the shrinkage of pore size. This can be attributed to the fact that gas adsorption layer decreases the effective flow area of the water phase and surface diffusion capacity for adsorbed gas is enhanced in small pore size.

A Comparison of Coarse-Grained and Continuum Models for Membrane Bending in Lipid Bilayer Fusion Pores

PubMed Central

Yoo, Jejoong; Jackson, Meyer B.; Cui, Qiang

2013-01-01

To establish the validity of continuum mechanics models quantitatively for the analysis of membrane remodeling processes, we compare the shape and energies of the membrane fusion pore predicted by coarse-grained (MARTINI) and continuum mechanics models. The results at these distinct levels of resolution give surprisingly consistent descriptions for the shape of the fusion pore, and the deviation between the continuum and coarse-grained models becomes notable only when the radius of curvature approaches the thickness of a monolayer. Although slow relaxation beyond microseconds is observed in different perturbative simulations, the key structural features (e.g., dimension and shape of the fusion pore near the pore center) are consistent among independent simulations. These observations provide solid support for the use of coarse-grained and continuum models in the analysis of membrane remodeling. The combined coarse-grained and continuum analysis confirms the recent prediction of continuum models that the fusion pore is a metastable structure and that its optimal shape is neither toroidal nor catenoidal. Moreover, our results help reveal a new, to our knowledge, bowing feature in which the bilayers close to the pore axis separate more from one another than those at greater distances from the pore axis; bowing helps reduce the curvature and therefore stabilizes the fusion pore structure. The spread of the bilayer deformations over distances of hundreds of nanometers and the substantial reduction in energy of fusion pore formation provided by this spread indicate that membrane fusion can be enhanced by allowing a larger area of membrane to participate and be deformed. PMID:23442963

Developing an Effective Model for Shale Gas Flow in Nano-scale Pore Clusters based on FIB-SEM Images

NASA Astrophysics Data System (ADS)

Jiang, W. B.; Lin, M.; Yi, Z. X.; Li, H. S.

2016-12-01

Nano-scale pores existed in the form of clusters are the controlling void space in shale gas reservoir. Gas transport in nanopores which has a significant influence on shale gas' recoverability displays multiple transport regimes, including viscous, slippage flow and Knudsen diffusion. In addition, it is also influenced by pore space characteristics. For convenience and efficiency consideration, it is necessary to develop an upscaling model from nano pore to pore cluster scale. Existing models are more like framework functions that provide a format, because the parameters that represent pore space characteristics are underdetermined and may have multiple possibilities. Therefore, it is urgent to make them clear and obtained a model that is closer to reality. FIB-SEM imaging technology is able to acquire three dimensional images with nanometer resolution that nano pores can be visible. Based on the images of two shale samples, we used a high-precision pore network extraction algorithm to generate equivalent pore networks and simulate multiple regime (non-Darcy) flow in it. Several structural parameters can be obtained through pore network modelling. It is found that although the throat-radius distributions are very close, throat flux-radius distributions of different samples can be divided into two categories. The variation of tortuosity with pressure and the overall trend of throat-flux distribution changes with pressure are disclosed. A deeper understanding of shale gas flow in nano-scale pore clusters is obtained. After all, an upscaling model that connects absolute permeability, apparent permeability and other characteristic parameters is proposed, and the best parameter scheme considering throat number-radius distribution and flowing porosity for this model is selected out of three schemes based on pore scale results, and it can avoid multiple-solution problem and is useful in reservoir modelling and experiment result analysis, etc. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the National Natural Science Foundation of China (41574129), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the 973 Program (2014CB239004).

Mechanical instability and percolation of deformable particles through porous networks

NASA Astrophysics Data System (ADS)

Benet, Eduard; Lostec, Guillaume; Pellegrino, John; Vernerey, Franck

2018-04-01

The transport of micron-sized particles such as bacteria, cells, or synthetic lipid vesicles through porous spaces is a process relevant to drug delivery, separation systems, or sensors, to cite a few examples. Often, the motion of these particles depends on their ability to squeeze through small constrictions, making their capacity to deform an important factor for their permeation. However, it is still unclear how the mechanical behavior of these particles affects collective transport through porous networks. To address this issue, we present a method to reconcile the pore-scale mechanics of the particles with the Darcy scale to understand the motion of a deformable particle through a porous network. We first show that particle transport is governed by a mechanical instability occurring at the pore scale, which leads to a binary permeation response on each pore. Then, using the principles of directed bond percolation, we are able to link this microscopic behavior to the probability of permeating through a random porous network. We show that this instability, together with network uniformity, are key to understanding the nonlinear permeation of particles at a given pressure gradient. The results are then summarized by a phase diagram that predicts three distinct permeation regimes based on particle properties and the randomness of the pore network.

Numerical Generation of Dense Plume Fingers in Unsaturated Homogeneous Porous Media

NASA Astrophysics Data System (ADS)

Cremer, C.; Graf, T.

2012-04-01

In nature, the migration of dense plumes typically results in the formation of vertical plume fingers. Flow direction in fingers is downwards, which is counterbalanced by upwards flow of less dense fluid between fingers. In heterogeneous media, heterogeneity itself is known to trigger the formation of fingers. In homogeneous media, however, fingers are also created even if all grains had the same diameter. The reason is that pore-scale heterogeneity leading to different flow velocities also exists in homogeneous media due to two effects: (i) Grains of identical size may randomly arrange differently, e.g. forming tetrahedrons, hexahedrons or octahedrons. Each arrangement creates pores of varying diameter, thus resulting in different average flow velocities. (ii) Random variations of solute concentration lead to varying buoyancy effects, thus also resulting in different velocities. As a continuation of previously made efforts to incorporate pore-scale heterogeneity into fully saturated soil such that dense fingers are realistically generated (Cremer and Graf, EGU Assembly, 2011), the current paper extends the research scope from saturated to unsaturated soil. Perturbation methods are evaluated by numerically re-simulating a laboratory-scale experiment of plume transport in homogeneous unsaturated sand (Simmons et al., Transp. Porous Media, 2002). The following 5 methods are being discussed: (i) homogeneous sand, (ii) initial perturbation of solute concentration, (iii) spatially random, time-constant perturbation of solute source, (iv) spatially and temporally random noise of simulated solute concentration, and (v) random K-field that introduces physically insignificant but numerically significant heterogeneity. Results demonstrate that, as opposed to saturated flow, perturbing the solute source will not result in plume fingering. This is because the location of the perturbed source (domain top) and the location of finger generation (groundwater surface) do not coincide. Alternatively, similar to saturated flow, applying either a random concentration noise (iv) or a random K-field (v) generates realistic plume fingering. Future work will focus on the generation mechanisms of plume finger splitting.

Heat of adsorption, adsorption stress, and optimal storage of methane in slit and cylindrical carbon pores predicted by classical density functional theory.

PubMed

Hlushak, Stepan

2018-01-03

Temperature, pressure and pore-size dependences of the heat of adsorption, adsorption stress, and adsorption capacity of methane in simple models of slit and cylindrical carbon pores are studied using classical density functional theory (CDFT) and grand-canonical Monte-Carlo (MC) simulation. Studied properties depend nontrivially on the bulk pressure and the size of the pores. Heat of adsorption increases with loading, but only for sufficiently narrow pores. While the increase is advantageous for gas storage applications, it is less significant for cylindrical pores than for slits. Adsorption stress and the average adsorbed fluid density show oscillatory dependence on the pore size and increase with bulk pressure. Slit pores exhibit larger amplitude of oscillations of the normal adsorption stress with pore size increase than cylindrical pores. However, the increase of the magnitude of the adsorption stress with bulk pressure increase is more significant for cylindrical than for slit pores. Adsorption stress appears to be negative for a wide range of pore sizes and external conditions. The pore size dependence of the average delivered density of the gas is analyzed and the optimal pore sizes for storage applications are estimated. The optimal width of slit pore appears to be almost independent of storage pressure at room temperature and pressures above 10 bar. Similarly to the case of slit pores, the optimal radius of cylindrical pores does not exhibit much dependence on the storage pressure above 15 bar. Both optimal width and optimal radii of slit and cylindrical pores increase as the temperature decreases. A comparison of the results of CDFT theory and MC simulations reveals subtle but important differences in the underlying fluid models employed by the approaches. The differences in the high-pressure behaviour between the hard-sphere 2-Yukawa and Lennard-Jones models of methane, employed by the CDFT and MC approaches, respectively, result in an overestimation of the heat of adsorption by the CDFT theory at higher loadings. However, both adsorption stress and adsorption capacity appear to be much less sensitive to the differences between the models and demonstrate excellent agreement between the theory and the computer experiment.

Molecular Dynamic Simulation of Water Vapor and Determination of Diffusion Characteristics in the Pore

NASA Astrophysics Data System (ADS)

Nikonov, Eduard G.; Pavluš, Miron; Popovičová, Mária

2018-02-01

One of the varieties of pores, often found in natural or artificial building materials, are the so-called blind pores of dead-end or saccate type. Three-dimensional model of such kind of pore has been developed in this work. This model has been used for simulation of water vapor interaction with individual pore by molecular dynamics in combination with the diffusion equation method. Special investigations have been done to find dependencies between thermostats implementations and conservation of thermodynamic and statistical values of water vapor - pore system. The two types of evolution of water - pore system have been investigated: drying and wetting of the pore. Full research of diffusion coefficient, diffusion velocity and other diffusion parameters has been made.

The beneficial effects of straight open large pores in the support on steam electrolysis performance of electrode-supported solid oxide electrolysis cell

NASA Astrophysics Data System (ADS)

Lin, Jie; Chen, Long; Liu, Tong; Xia, Changrong; Chen, Chusheng; Zhan, Zhongliang

2018-01-01

This study is aimed at improving the electrochemical performance of electrode-supported solid oxide electrolysis cells (SOECs) by optimizing the pore structure of the supports. Two planar NiO-8 mol% yttria-stabilized zirconia supports are prepared, one by the phase-inversion tape casting, and the other by conventional tape casting method using graphite as the pore former. The former contains finger-like straight open large pores, while the latter contains randomly distributed and tortuous pores. The steam electrolysis of the cells with different microstructure cathode supports is measured. The cell supported on the cathode with straight pores shows a high current density of 1.42 A cm-2 and a H2 production rate of 9.89 mL (STP) cm-2 min-1 at 1.3 V and 50 vol % humidity and 750 °C, while the cell supported on the cathode with tortuous pores shows a current density of only 0.91 A cm-2 and a H2 production rate of 6.34 mL cm-2min-1. It is concluded that the introduction of large straight open pores into the cathode support allows fast gas phase transport and thus minimizes the concentration polarization. Furthermore, the straight pores could provide better access to the reaction site (the electrode functional layer), thereby reducing the activation polarization as well.

Pore-Scale Determination of Gas Relative Permeability in Hydrate-Bearing Sediments Using X-Ray Computed Micro-Tomography and Lattice Boltzmann Method

NASA Astrophysics Data System (ADS)

Chen, Xiongyu; Verma, Rahul; Espinoza, D. Nicolas; Prodanović, Maša.

2018-01-01

This work uses X-ray computed micro-tomography (μCT) to monitor xenon hydrate growth in a sandpack under the excess gas condition. The μCT images give pore-scale hydrate distribution and pore habit in space and time. We use the lattice Boltzmann method to calculate gas relative permeability (krg) as a function of hydrate saturation (Shyd) in the pore structure of the experimental hydrate-bearing sand retrieved from μCT data. The results suggest the krg - Shyd data fit well a new model krg = (1-Shyd)·exp(-4.95·Shyd) rather than the simple Corey model. In addition, we calculate krg-Shyd curves using digital models of hydrate-bearing sand based on idealized grain-attaching, coarse pore-filling, and dispersed pore-filling hydrate habits. Our pore-scale measurements and modeling show that the krg-Shyd curves are similar regardless of whether hydrate crystals develop grain-attaching or coarse pore-filling habits. The dispersed pore filling habit exhibits much lower gas relative permeability than the other two, but it is not observed in the experiment and not compatible with Ostwald ripening mechanisms. We find that a single grain-shape factor can be used in the Carman-Kozeny equation to calculate krg-Shyd data with known porosity and average grain diameter, suggesting it is a useful model for hydrate-bearing sand.

The exocytotic fusion pore modeled as a lipidic pore.

PubMed Central

Nanavati, C; Markin, V S; Oberhauser, A F; Fernandez, J M

1992-01-01

Freeze-fracture electron micrographs from degranulating cells show that the lumen of the secretory granule is connected to the extracellular compartment via large (20 to 150 nm diameter) aqueous pores. These exocytotic fusion pores appear to be made up of a highly curved bilayer that spans the plasma and granule membranes. Conductance measurements, using the patch-clamp technique, have been used to study the fusion pore from the instant it conducts ions. These measurements reveal the presence of early fusion pores that are much smaller than those observed in electron micrographs. Early fusion pores open abruptly, fluctuate, and then either expand irreversibly or close. The molecular structure of these early fusion pores is unknown. In the simplest extremes, these early fusion pores could be either ion channel like protein pores or lipidic pores. Here, we explored the latter possibility, namely that of the early exocytotic fusion pore modeled as a lipid-lined pore whose free energy was composed of curvature elastic energy and work done by tension. Like early exocytotic fusion pores, we found that these lipidic pores could open abruptly, fluctuate, and expand irreversibly. Closure of these lipidic pores could be caused by slight changes in lipid composition. Conductance distributions for stable lipidic pores matched those of exocytotic fusion pores. These findings demonstrate that lipidic pores can exhibit the properties of exocytotic fusion pores, thus providing an alternate framework with which to understand and interpret exocytotic fusion pore data. PMID:1420930

An Image-based Micro-continuum Pore-scale Model for Gas Transport in Organic-rich Shale

NASA Astrophysics Data System (ADS)

Guo, B.; Tchelepi, H.

2017-12-01

Gas production from unconventional source rocks, such as ultra-tight shales, has increased significantly over the past decade. However, due to the extremely small pores ( 1-100 nm) and the strong material heterogeneity, gas flow in shale is still not well understood and poses challenges for predictive field-scale simulations. In recent years, digital rock analysis has been applied to understand shale gas transport at the pore-scale. An issue with rock images (e.g. FIB-SEM, nano-/micro-CT images) is the so-called "cutoff length", i.e., pores and heterogeneities below the resolution cannot be resolved, which leads to two length scales (resolved features and unresolved sub-resolution features) that are challenging for flow simulations. Here we develop a micro-continuum model, modified from the classic Darcy-Brinkman-Stokes framework, that can naturally couple the resolved pores and the unresolved nano-porous regions. In the resolved pores, gas flow is modeled with Stokes equation. In the unresolved regions where the pore sizes are below the image resolution, we develop an apparent permeability model considering non-Darcy flow at the nanoscale including slip flow, Knudsen diffusion, adsorption/desorption, surface diffusion, and real gas effect. The end result is a micro-continuum pore-scale model that can simulate gas transport in 3D reconstructed shale images. The model has been implemented in the open-source simulation platform OpenFOAM. In this paper, we present case studies to demonstrate the applicability of the model, where we use 3D segmented FIB-SEM and nano-CT shale images that include four material constituents: organic matter, clay, granular mineral, and pore. In addition to the pore structure and the distribution of the material constituents, we populate the model with experimental measurements (e.g. size distribution of the sub-resolution pores from nitrogen adsorption) and parameters from the literature and identify the relative importance of different physics on gas production. Overall, the micro-continuum model provides a novel tool for digital rock analysis of organic-rich shale.

Kinetic models of controllable pore growth of anodic aluminum oxide membrane

NASA Astrophysics Data System (ADS)

Huang, Yan; Zeng, Hong-yan; Zhao, Ce; Qu, Ye-qing; Zhang, Pin

2012-06-01

An anodized Al2O3 (AAO) membrane with apertures about 72 nm in diameter was prepared by two-step anodic oxidation. The appearance and pore arrangement of the AAO membrane were characterized by energy dispersive x-ray spectroscopy and scanning electron microscopy. It was confirmed that the pores with high pore aspect ratio were parallel, well-ordered, and uniform. The kinetics of pores growth in the AAO membrane was derived, and the kinetic models showed that pores stopped developing when the pressure ( σ) trended to equal the surface tension at the end of anodic oxidation. During pore expansion, the effects of the oxalic acid concentration and expansion time on the pore size were investigated, and the kinetic behaviors were explained with two kinetic models derived in this study. They showed that the pore size increased with extended time ( r= G· t+ G'), but decreased with increased concentration ( r = - K·ln c- K') through the derived mathematic formula. Also, the values of G, G', K, and K' were derived from our experimental data.

Coupling root architecture and pore network modeling - an attempt towards better understanding root-soil interactions

NASA Astrophysics Data System (ADS)

Leitner, Daniel; Bodner, Gernot; Raoof, Amir

2013-04-01

Understanding root-soil interactions is of high importance for environmental and agricultural management. Root uptake is an essential component in water and solute transport modeling. The amount of groundwater recharge and solute leaching significantly depends on the demand based plant extraction via its root system. Plant uptake however not only responds to the potential demand, but in most situations is limited by supply form the soil. The ability of the plant to access water and solutes in the soil is governed mainly by root distribution. Particularly under conditions of heterogeneous distribution of water and solutes in the soil, it is essential to capture the interaction between soil and roots. Root architecture models allow studying plant uptake from soil by describing growth and branching of root axes in the soil. Currently root architecture models are able to respond dynamically to water and nutrient distribution in the soil by directed growth (tropism), modified branching and enhanced exudation. The porous soil medium as rooting environment in these models is generally described by classical macroscopic water retention and sorption models, average over the pore scale. In our opinion this simplified description of the root growth medium implies several shortcomings for better understanding root-soil interactions: (i) It is well known that roots grow preferentially in preexisting pores, particularly in more rigid/dry soil. Thus the pore network contributes to the architectural form of the root system; (ii) roots themselves can influence the pore network by creating preferential flow paths (biopores) which are an essential element of structural porosity with strong impact on transport processes; (iii) plant uptake depend on both the spatial location of water/solutes in the pore network as well as the spatial distribution of roots. We therefore consider that for advancing our understanding in root-soil interactions, we need not only to extend our root models, but also improve the description of the rooting environment. Until now there have been no attempts to couple root architecture and pore network models. In our work we present a first attempt to join both types of models using the root architecture model of Leitner et al., (2010) and a pore network model presented by Raoof et al. (2010). The two main objectives of coupling both models are: (i) Representing the effect of root induced biopores on flow and transport processes: For this purpose a fixed root architecture created by the root model is superimposed as a secondary root induced pore network to the primary soil network, thus influencing the final pore topology in the network generation. (ii) Representing the influence of pre-existing pores on root branching: Using a given network of (rigid) pores, the root architecture model allocates its root axes into these preexisting pores as preferential growth paths with thereby shape the final root architecture. The main objective of our study is to reveal the potential of using a pore scale description of the plant growth medium for an improved representation of interaction processes at the interface of root and soil. References Raoof, A., Hassanizadeh, S.M. 2010. A New Method for Generating Pore-Network Models. Transp. Porous Med. 81, 391-407. Leitner, D, Klepsch, S., Bodner, G., Schnepf, S. 2010. A dynamic root system growth model based on L-Systems. Tropisms and coupling to nutrient uptake from soil. Plant Soil 332, 177-192.

Root induced changes of effective 1D hydraulic properties in a soil column.

PubMed

Scholl, P; Leitner, D; Kammerer, G; Loiskandl, W; Kaul, H-P; Bodner, G

Roots are essential drivers of soil structure and pore formation. This study aimed at quantifying root induced changes of the pore size distribution (PSD). The focus was on the extent of clogging vs. formation of pores during active root growth. Parameters of Kosugi's lognormal PSD model were determined by inverse estimation in a column experiment with two cover crops (mustard, rye) and an unplanted control. Pore dynamics were described using a convection-dispersion like pore evolution model. Rooted treatments showed a wider range of pore radii with increasing volumes of large macropores >500 μm and micropores <2.5 μm, while fine macropores, mesopores and larger micropores decreased. The non-rooted control showed narrowing of the PSD and reduced porosity over all radius classes. The pore evolution model accurately described root induced changes, while structure degradation in the non-rooted control was not captured properly. Our study demonstrated significant short term root effects with heterogenization of the pore system as dominant process of root induced structure formation. Pore clogging is suggested as a partial cause for reduced pore volume. The important change in micro- and large macropores however indicates that multiple mechanic and biochemical processes are involved in root-pore interactions.

The Effect of the Pore Entrance on Particle Motion in Slit Pores: Implications for Ultrathin Membranes

PubMed Central

Delavari, Armin; Baltus, Ruth

2017-01-01

Membrane rejection models generally neglect the effect of the pore entrance on intrapore particle transport. However, entrance effects are expected to be particularly important with ultrathin membranes, where membrane thickness is typically comparable to pore size. In this work, a 2D model was developed to simulate particle motion for spherical particles moving at small Re and infinite Pe from the reservoir outside the pore into a slit pore. Using a finite element method, particles were tracked as they accelerated across the pore entrance until they reached a steady velocity in the pore. The axial position in the pore where particle motion becomes steady is defined as the particle entrance length (PEL). PELs were found to be comparable to the fluid entrance length, larger than the pore size and larger than the thickness typical of many ultrathin membranes. Results also show that, in the absence of particle diffusion, hydrodynamic particle–membrane interactions at the pore mouth result in particle “funneling” in the pore, yielding cross-pore particle concentration profiles focused at the pore centerline. The implications of these phenomena on rejection from ultrathin membranes are examined. PMID:28796197

Three-Dimensional Modeling of the Brain's ECS by Minimum Configurational Energy Packing of Fluid Vesicles

PubMed Central

Nandigam, Ravi K.; Kroll, Daniel M.

2007-01-01

The extracellular space of the brain is the heterogeneous porous medium formed by the spaces between the brain cells. Diffusion in this interstitial space is the mechanism by which glucose and oxygen are delivered to the brain cells from the vascular system. It is also a medium for the transport of certain informational substances between the cells (called volume transmission), and for drug delivery. This work involves three-dimensional modeling of the extracellular space as void space in close-packed arrays of fluid membrane vesicles. These packings are generated by minimizing the configurational energy using a Monte Carlo procedure. Both regular and random packs of vesicles are considered. A random walk algorithm is then used to compute the geometric tortuosities, and the results are compared with published experimental data. For the random packings, it is found that although the absolute values for the tortuosities differ, the dependence of the tortuosity on pore volume fraction is very similar to that observed in experiment. The tortuosities we measure are larger than those computed in previous studies of packings of convex polytopes, and modeling improvements, which require higher resolution studies and an improved modeling of brain cell shapes and mechanical properties, could help resolve remaining discrepancies between model simulations and experiment. It is also shown that the specular reflection scheme is the appropriate technique for implementing zero-flux boundary conditions in random walk simulations commonly encountered in diffusion problems. PMID:17307830

Three-dimensional modeling of the brain's ECS by minimum configurational energy packing of fluid vesicles.

PubMed

Nandigam, Ravi K; Kroll, Daniel M

2007-05-15

The extracellular space of the brain is the heterogeneous porous medium formed by the spaces between the brain cells. Diffusion in this interstitial space is the mechanism by which glucose and oxygen are delivered to the brain cells from the vascular system. It is also a medium for the transport of certain informational substances between the cells (called volume transmission), and for drug delivery. This work involves three-dimensional modeling of the extracellular space as void space in close-packed arrays of fluid membrane vesicles. These packings are generated by minimizing the configurational energy using a Monte Carlo procedure. Both regular and random packs of vesicles are considered. A random walk algorithm is then used to compute the geometric tortuosities, and the results are compared with published experimental data. For the random packings, it is found that although the absolute values for the tortuosities differ, the dependence of the tortuosity on pore volume fraction is very similar to that observed in experiment. The tortuosities we measure are larger than those computed in previous studies of packings of convex polytopes, and modeling improvements, which require higher resolution studies and an improved modeling of brain cell shapes and mechanical properties, could help resolve remaining discrepancies between model simulations and experiment. It is also shown that the specular reflection scheme is the appropriate technique for implementing zero-flux boundary conditions in random walk simulations commonly encountered in diffusion problems.

Pore-scale and continuum simulations of solute transport micromodel benchmark experiments

DOE PAGES

Oostrom, M.; Mehmani, Y.; Romero-Gomez, P.; ...

2014-06-18

Four sets of nonreactive solute transport experiments were conducted with micromodels. Three experiments with one variable, i.e., flow velocity, grain diameter, pore-aspect ratio, and flow-focusing heterogeneity were in each set. The data sets were offered to pore-scale modeling groups to test their numerical simulators. Each set consisted of two learning experiments, for which our results were made available, and one challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the transverse dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing,more » and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice Boltzmann (LB) approach, and one used a computational fluid dynamics (CFD) technique. Furthermore, we used the learning experiments, by the PN models, to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used the learning experiments to appropriately discretize the spatial grid representations. For the continuum modeling, the required dispersivity input values were estimated based on published nonlinear relations between transverse dispersion coefficients and Peclet number. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values, resulting in reduced dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models, which account for the micromodel geometry and underlying flow and transport physics, needed up to several days on supercomputers to resolve the more complex problems.« less

Modelling in conventional electroporation for model cell with organelles using COMSOL Multiphysics

NASA Astrophysics Data System (ADS)

Sulaeman, M. Y.; Widita, R.

2016-03-01

Conventional electroporation is a formation of pores in the membrane cell due to the external electric field applied to the cell. The purpose of creating pores in the cell using conventional electroporation are to increase the effectiveness of chemotherapy (electrochemotherapy) and to kill cancer tissue using irreversible electroporation. Modeling of electroporation phenomenon on a model cell had been done by using software COMSOL Multiphysics 4.3b with the applied external electric field with intensity at 1.1 kV/cm to find transmembrane voltage and pore density. It can be concluded from the results of potential distribution and transmembrane voltage, it show that pores formation only occurs in the membrane cells and it could not penetrate into inside the model cell so there is not pores formation in its organells.

A spatial dynamic model to assess piospheric land degradation processes of SW Iberian rangelands

NASA Astrophysics Data System (ADS)

Herguido Sevillano, Estela; Ibáñez, Javier; Francisco Lavado Contador, Joaquín; Pulido-Fernández, Manuel; Schnabel, Susanne

2015-04-01

Iberian open wooded rangelands (known as dehesas or montados) constitute valuable agro-silvo-pastoral systems traditionally considered as highly sustainable. Nevertheless, in the recent decades, those systems are undergoing changes of land use and management practices that compromise its sustainability. Some of those changes, as the rising construction of watering points and the high spatial fragmentation and livestock movement restriction associated to fencing, show an aggregated effect with livestock, producing an impact gradient over soil and vegetation. Soil compaction related to livestock pressure is higher around watering points, with bare soil halos and patches of scarce vegetation or nude soil developing with higher frequency in areas close to them. Using the freeware Dinamica EGO as environmental modeling platform, we have developed a theoretic spatial dynamic model that represents some of the processes of land degradation associated to livestock grazing in dehesa fenced enclosures. Spatial resolution is high since every cell in the model is a square unit area of 1 m2. We paid particular attention to the relationships between soil degradation by compaction (porosity), livestock pressure, rainfall, pasture growth and shrub cover and bare soil generation. The model considers pasture growth as related to soil compaction, measured by the pore space in the top 10 cm soil layer. Annual precipitation is randomly generated following a normal distribution. When annual precipitation and pore space increase, also does pasture growth. Besides, there is a feedback between pasture growth and pore space, given that pasture roots increases soil porosity. The cell utility for livestock function has been defined as an exponential function of the distance of a cell to watering points and the amount of pasture present in it. The closer the cell to a pond and the higher the amount of pasture, the higher is cell utility. The latter is modulated by a normal random variable to capture accidental effects. This variable has zero mean and a standard deviation linearly related to the distance to the pond. Livestock utilization of a cell is a function of its relative utility, the stocking rate and the time that animals spend at the enclosure. Since livestock trampling promotes soil compaction, livestock utilization has a negative effect on pore space. The probability of transition from herbaceous to shrubs is also modulated by pore space, and thus livestock utilization, as shrub development needs a minimum porosity value for seeds to successfully germinate. In addition, it is influenced by the proportion of cells occupied by shrubs in a radius where seed dispersal or exclusion by competition may occur. The model contemplates the probability of transition from shrubs to herbaceous through shrub mortality, and the age of the shrubs, which influences seed production and shrub cover. Pasture consumption by livestock and pasture remaining at the end of summer were also modeled, so that it is possible to obtain maps of bare soil at that time. Likewise, the model generates maps of vegetation state (shrubs or herbaceous) and pasture growth. The values of the set of 31 parameters were obtained from field measurements and from publications. Those parameters lacking quantitative information were calibrated by comparing model performance with the dynamics of true enclosures analyzed between 1984 and 2009 in ortophotographs. Stocking rates were inferred from farmers' interviews performed in 2009 about present and past land use and management practices. The model developed is intended to analyze strategies of livestock management in dehesas. Particularly, soil conservation practices as related to livestock pressure can be simulated looking for optimized schemes. Moreover, the model provides the possibility of generating simulations for future climate scenarios, studying the effects of climate change on livestock carrying capacity on these systems. Thanks to the Spanish Ministerio de Economía y Competitividad for financially supporting this study through AMID (CGL2011-23361) project.

Hydraulic conductivity of variably saturated porous media: Film and corner flow in angular pore space

NASA Astrophysics Data System (ADS)

Tuller, Markus; Or, Dani

2001-05-01

Many models for hydraulic conductivity of partially saturated porous media rely on oversimplified representation of the pore space as a bundle of cylindrical capillaries and disregard flow in liquid films. Recent progress in modeling liquid behavior in angular pores of partially saturated porous media offers an alternative framework. We assume that equilibrium liquid-vapor interfaces provide well-defined and stable boundaries for slow laminar film and corner flow regimes in pore space comprised of angular pores connected to slit-shaped spaces. Knowledge of liquid configuration in the assumed geometry facilitates calculation of average liquid velocities in films and corners and enables derivation of pore-scale hydraulic conductivity as a function of matric potential. The pore-scale model is statistically upscaled to represent hydraulic conductivity for a sample of porous medium. Model parameters for the analytical sample-scale expressions are estimated from measured liquid retention data and other measurable medium properties. Model calculations illustrate the important role of film flow, whose contribution dominates capillary flow (in full pores and corners) at relatively high matric potentials (approximately -100 to -300 J kg-1, or -1 to 3 bars). The crossover region between film and capillary flow is marked by a significant change in the slope of the hydraulic conductivity function as often observed in measurements. Model predictions are compared with the widely applied van Genuchten-Mualem model and yield reasonable agreement with measured retention and hydraulic conductivity data over a wide range of soil textural classes.

Modeling the interaction of ultrasound with pores

NASA Technical Reports Server (NTRS)

Lu, Yichi; Wadley, Haydn N. G.; Parthasarathi, Sanjai

1991-01-01

Factors that affect ultrasonic velocity sensing of density during consolidation of metal powders are examined. A comparison is made between experimental results obtained during the final stage of densification and the predictions of models that assume either a spherical or a spheroidal pore shape. It is found that for measurements made at low frequencies during the final stage of densification, relative density (pore fraction) and pore shape are the two most important factors determining the ultrasonic velocity, the effect of pore size is negligible.

A universal model for nanoporous carbon supercapacitors applicable to diverse pore regimes, carbon materials, and electrolytes.

PubMed

Huang, Jingsong; Sumpter, Bobby G; Meunier, Vincent

2008-01-01

Supercapacitors, commonly called electric double-layer capacitors (EDLCs), are emerging as a novel type of energy-storage device with the potential to substitute batteries in applications that require high power densities. In response to the latest experimental breakthrough in nanoporous carbon supercapacitors, we propose a heuristic theoretical model that takes pore curvature into account as a replacement for the EDLC model, which is based on a traditional parallel-plate capacitor. When the pore size is in the mesopore regime (2-50 nm), counterions enter mesoporous carbon materials and approach the pore wall to form an electric double-cylinder capacitor (EDCC); in the micropore regime (<2 nm), solvated/desolvated counterions line up along the pore axis to form an electric wire-in-cylinder capacitor (EWCC). In the macropore regime (>50 nm) at which pores are large enough so that pore curvature is no longer significant, the EDCC model can be reduced naturally to the EDLC model. We present density functional theory calculations and detailed analyses of available experimental data in various pore regimes, which show the significant effects of pore curvature on the supercapacitor properties of nanoporous carbon materials. It is shown that the EDCC/EWCC model is universal for carbon supercapacitors with diverse carbon materials, including activated carbon materials, template carbon materials, and novel carbide-derived carbon materials, and with diverse electrolytes, including organic electrolytes, such as tetraethylammonium tetrafluoroborate (TEABF(4)) and tetraethylammonium methylsulfonate (TEAMS) in acetonitrile, aqueous H(2)SO(4) and KOH electrolytes, and even an ionic liquid electrolyte, such as 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide (EMI-TFSI). The EDCC/EWCC model allows the supercapacitor properties to be correlated with pore size, specific surface area, Debye length, electrolyte concentration and dielectric constant, and solute ion size It may lend support for the systematic optimization of the properties of carbon supercapacitors through experiments. On the basis of the insight obtained from the new model, we also discuss the effects of the kinetic solvation/desolvation process, multimodal (versus unimodal) pore size distribution, and exohedral (versus endohedral) capacitors on the electrochemical properties of supercapacitors.

Is trabecular bone permeability governed by molecular ordering-induced fluid viscosity gain? Arguments from re-evaluation of experimental data in the framework of homogenization theory.

PubMed

Abdalrahman, T; Scheiner, S; Hellmich, C

2015-01-21

It is generally agreed on that trabecular bone permeability, a physiologically important quantity, is governed by the material׳s (vascular or intertrabecular) porosity as well as by the viscosity of the pore-filling fluids. Still, there is less agreement on how these two key factors govern bone permeability. In order to shed more light onto this somewhat open issue, we here develop a random homogenization scheme for upscaling Poiseuille flow in the vascular porosity, up to Darcy-type permeability of the overall porous medium "trabecular bone". The underlying representative volume element of the macroscopic bone material contains two types of phases: a spherical, impermeable extracellular bone matrix phase interacts with interpenetrating cylindrical pore channel phases that are oriented in all different space directions. This type of interaction is modeled by means of a self-consistent homogenization scheme. While the permeability of the bone matrix equals to zero, the permeability of the pore phase is found through expressing the classical Hagen-Poiseuille law for laminar flow in the format of a "micro-Darcy law". The upscaling scheme contains pore size and porosity as geometrical input variables; however, they can be related to each other, based on well-known relations between porosity and specific bone surface. As two key results, validated through comprehensive experimental data, it appears (i) that the famous Kozeny-Carman constant (which relates bone permeability to the cube of the porosity, the square of the specific surface, as well as to the bone fluid viscosity) needs to be replaced by an again porosity-dependent rational function, and (ii) that the overall bone permeability is strongly affected by the pore fluid viscosity, which, in case of polarized fluids, is strongly increased due to the presence of electrically charged pore walls. Copyright © 2014 Elsevier Ltd. All rights reserved.

Clinical outcome and patient satisfaction using biodegradable (NasoPore) and non-biodegradable packing, a double-blind, prospective, randomized study.

PubMed

Burduk, Pawel Krzysztof; Wierzchowska, Malgorzata; Grześkowiak, Blazej; Kaźmierczak, Wojciech; Wawrzyniak, Katarzyna

Nasal packing after endoscopic sinus surgery is used as a standard procedure. The optimum solution to minimize or eliminate all disadvantages of this procedure may be accomplished using biodegradable packs. The aim of this study was to compare patient satisfaction and clinical outcome associated with absorbable and non-absorbable packing after FESS. In total, 50 patients were included in a prospective, double-blind, randomized trial. One side was packed with polyurethane foam, while the opposite side was packed with gauze packing. On the 2nd, 10th, and 30th postoperative day, the patients were questioned with the aid of a visual analog scale. The standardized questionnaires for bleeding, nasal breathing, feeling of pressure, and headache were used. The presence of synechiae, infection, or granulation was noted and recorded with the video-endoscopy. A significant difference according to lower pressure was found in the NasoPore group compared to the controls on day ten after surgery. The NasoPore packing had lower scores with respect to postoperative nose blockage on the 2nd and 10th days. Mucosal healing was better for the NasoPore group, both at day ten and 30 compared with the control group. The overall patient comfort is higher when using NasoPore compared to non-resorbable traditional impregnated gauze packing. Intensive saline douches applied three to four times per day are mandatory after the operation to prevent synechiae formation and fluid resorption by the packing. Published by Elsevier Editora Ltda.

Investigating the relative permeability behavior of microporosity-rich carbonates and tight sandstones with multiscale pore network models

NASA Astrophysics Data System (ADS)

Bultreys, Tom; Stappen, Jeroen Van; Kock, Tim De; Boever, Wesley De; Boone, Marijn A.; Hoorebeke, Luc Van; Cnudde, Veerle

2016-11-01

The relative permeability behavior of rocks with wide ranges of pore sizes is in many cases still poorly understood and is difficult to model at the pore scale. In this work, we investigate the capillary pressure and relative permeability behavior of three outcrop carbonates and two tight reservoir sandstones with wide, multimodal pore size distributions. To examine how the drainage and imbibition properties of these complex rock types are influenced by the connectivity of macropores to each other and to zones with unresolved small-scale porosity, we apply a previously presented microcomputed-tomography-based multiscale pore network model to these samples. The sensitivity to the properties of the small-scale porosity is studied by performing simulations with different artificial sphere-packing-based networks as a proxy for these pores. Finally, the mixed-wet water-flooding behavior of the samples is investigated, assuming different wettability distributions for the microporosity and macroporosity. While this work is not an attempt to perform predictive modeling, it seeks to qualitatively explain the behavior of the investigated samples and illustrates some of the most recent developments in multiscale pore network modeling.

Effect of Pore Pressure on Slip Failure of an Impermeable Fault: A Coupled Micro Hydro-Geomechanical Model

NASA Astrophysics Data System (ADS)

Yang, Z.; Juanes, R.

2015-12-01

The geomechanical processes associated with subsurface fluid injection/extraction is of central importance for many industrial operations related to energy and water resources. However, the mechanisms controlling the stability and slip motion of a preexisting geologic fault remain poorly understood and are critical for the assessment of seismic risk. In this work, we develop a coupled hydro-geomechanical model to investigate the effect of fluid injection induced pressure perturbation on the slip behavior of a sealing fault. The model couples single-phase flow in the pores and mechanics of the solid phase. Granular packs (see example in Fig. 1a) are numerically generated where the grains can be either bonded or not, depending on the degree of cementation. A pore network is extracted for each granular pack with pore body volumes and pore throat conductivities calculated rigorously based on geometry of the local pore space. The pore fluid pressure is solved via an explicit scheme, taking into account the effect of deformation of the solid matrix. The mechanics part of the model is solved using the discrete element method (DEM). We first test the validity of the model with regard to the classical one-dimensional consolidation problem where an analytical solution exists. We then demonstrate the ability of the coupled model to reproduce rock deformation behavior measured in triaxial laboratory tests under the influence of pore pressure. We proceed to study the fault stability in presence of a pressure discontinuity across the impermeable fault which is implemented as a plane with its intersected pore throats being deactivated and thus obstructing fluid flow (Fig. 1b, c). We focus on the onset of shear failure along preexisting faults. We discuss the fault stability criterion in light of the numerical results obtained from the DEM simulations coupled with pore fluid flow. The implication on how should faults be treated in a large-scale continuum model is also presented.

Investigating Hydrophilic Pores in Model Lipid Bilayers using Molecular Simulations: Correlating Bilayer Properties with Pore Formation Thermodynamics

PubMed Central

Hu, Yuan; Sinha, Sudipta Kumar

2015-01-01

Cell-penetrating and antimicrobial peptides show remarkable ability to translocate across physiological membranes. Along with factors such as electric potential induced-perturbations of membrane structure and surface tension effects, experiments invoke pore-like membrane configurations during the solute transfer process into vesicles and cells. The initiation and formation of pores are associated with a non-trivial free energy cost, thus necessitating consideration of the factors associated with pore formation and attendant free energetics. Due to experimental and modeling challenges related to the long timescales of the translocation process, we use umbrella-sampling molecular dynamics simulations with a lipid-density based order parameter to investigate membrane pore-formation free energy employing Martini coarse-grained models. We investigate structure and thermodynamic features of the pore in 18 lipids spanning a range of head-groups, charge states, acyl chain lengths and saturation. We probe the dependence of pore-formation barriers on area per lipid, lipid bilayer thickness, membrane bending rigidities in three different lipid classes. The pore formation free energy in pure bilayers and peptide translocating scenarios are significantly coupled with bilayer thickness. Thicker bilayers require more reversible work to create pores. Pore formation free energy is higher in peptide-lipid systems relative to the peptide-free lipid systems due to penalties to maintain solvation of charged hydrophilic solutes within the membrane environment. PMID:25614183

Generalized network modeling of capillary-dominated two-phase flow

NASA Astrophysics Data System (ADS)

Raeini, Ali Q.; Bijeljic, Branko; Blunt, Martin J.

2018-02-01

We present a generalized network model for simulating capillary-dominated two-phase flow through porous media at the pore scale. Three-dimensional images of the pore space are discretized using a generalized network—described in a companion paper [A. Q. Raeini, B. Bijeljic, and M. J. Blunt, Phys. Rev. E 96, 013312 (2017), 10.1103/PhysRevE.96.013312]—which comprises pores that are divided into smaller elements called half-throats and subsequently into corners. Half-throats define the connectivity of the network at the coarsest level, connecting each pore to half-throats of its neighboring pores from their narrower ends, while corners define the connectivity of pore crevices. The corners are discretized at different levels for accurate calculation of entry pressures, fluid volumes, and flow conductivities that are obtained using direct simulation of flow on the underlying image. This paper discusses the two-phase flow model that is used to compute the averaged flow properties of the generalized network, including relative permeability and capillary pressure. We validate the model using direct finite-volume two-phase flow simulations on synthetic geometries, and then present a comparison of the model predictions with a conventional pore-network model and experimental measurements of relative permeability in the literature.

Quantifying the impact of early calcite cementation on the reservoir quality of carbonate rocks: A 3D process-based model

NASA Astrophysics Data System (ADS)

Hosa, Aleksandra; Wood, Rachel

2017-06-01

The reservoir properties of carbonate rocks are controlled by both deposition and diagenesis. The latter includes the early precipitation of calcite cements, which can exert a strong control on the evolution of subsequent diagenetic pathways. We quantify the impact of early marine cement growth in grainstones on evolving pore space by examining trends in the relationship between cementation and permeability using a 3D process-based model (Calcite3D). The model assumes varying proportions of polycrystalline and monocrystalline grain types, upon which we grow isopachous and syntaxial calcite cement types, respectively. We model two syntaxial cement shapes, compact and elongated, that approximate the geometries of typical rhombohedral calcite forms. Results demonstrate the effect of cement competition: an increasing proportion of monocrystalline grains creates stronger competition and a reduction in the impact of individual grains on final calcite cement volume and porosity. Isopachous cement is effective in closing pore throats and limiting permeability. We also show that the impact of syntaxial cement on porosity occlusion and therefore flow is highly dependent on monocrystalline grain location and the orientation of crystal axes. This demonstrates the importance of diagenetic overprint in controlling the evolution of rock properties, but also that this process can be essentially random. We also show that diagenesis alone can create notable heterogeneity in the permeability of carbonates. While Calcite3D is successful in modelling realistic changes in cement volumes and pore space morphology, modelled permeabilities (0.01 -30D) are above the range reported in reservoir grainstones due to the very high permeability of the initial synthetic sediment deposit (58.9D). Poroperm data generated by Calcite3D, however, exhibits a linear relationship between the logarithms of porosity and permeability with a high coefficient of determination, as observed in natural media.

Application of Mortar Coupling in Multiscale Modelling of Coupled Flow, Transport, and Biofilm Growth in Porous Media

NASA Astrophysics Data System (ADS)

Laleian, A.; Valocchi, A. J.; Werth, C. J.

2017-12-01

Multiscale models of reactive transport in porous media are capable of capturing complex pore-scale processes while leveraging the efficiency of continuum-scale models. In particular, porosity changes caused by biofilm development yield complex feedbacks between transport and reaction that are difficult to quantify at the continuum scale. Pore-scale models, needed to accurately resolve these dynamics, are often impractical for applications due to their computational cost. To address this challenge, we are developing a multiscale model of biofilm growth in which non-overlapping regions at pore and continuum spatial scales are coupled with a mortar method providing continuity at interfaces. We explore two decompositions of coupled pore-scale and continuum-scale regions to study biofilm growth in a transverse mixing zone. In the first decomposition, all reaction is confined to a pore-scale region extending the transverse mixing zone length. Only solute transport occurs in the surrounding continuum-scale regions. Relative to a fully pore-scale result, we find the multiscale model with this decomposition has a reduced run time and consistent result in terms of biofilm growth and solute utilization. In the second decomposition, reaction occurs in both an up-gradient pore-scale region and a down-gradient continuum-scale region. To quantify clogging, the continuum-scale model implements empirical relations between porosity and continuum-scale parameters, such as permeability and the transverse dispersion coefficient. Solutes are sufficiently mixed at the end of the pore-scale region, such that the initial reaction rate is accurately computed using averaged concentrations in the continuum-scale region. Relative to a fully pore-scale result, we find accuracy of biomass growth in the multiscale model with this decomposition improves as the interface between pore-scale and continuum-scale regions moves downgradient where transverse mixing is more fully developed. Also, this decomposition poses additional challenges with respect to mortar coupling. We explore these challenges and potential solutions. While recent work has demonstrated growing interest in multiscale models, further development is needed for their application to field-scale subsurface contaminant transport and remediation.

Product layer development during sulfation and sulfidation of uncalcined limestone particles at elevated pressures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zevenhoven, C.A.P.; Yrjas, K.P.; Hupa, M.M.

1998-07-01

Fluidized bed combustion or gasification allows for in-bed sulfur capture with a calcium-based sorbent such as limestone or dolomite. Sorbent particle size, porosity, internal surface, and their variation during conversion have great influence on the conversion of the sorbent. The uptake of SO{sub 2} and H{sub 2}S by five physically different limestones is discussed, for typical pressurized fluidized bed combustor or gasifier conditions: 850/950 C, 15/20 bar. Tests were done in a pressurized thermogravimetric apparatus (P-TGA), the size of the limestone particles was 250--300 {micro}m. It is stressed that the limestones remain uncalcined. A changing internal structure (CIS) model ismore » presented in which reaction kinetics and product layer diffusion are related to the intraparticle surface of reaction, instead of the outer particle surface as in unreacted shrinking core (USC)-type models. The random pore model was used for describing the changing internal pore and reaction surfaces. Rate parameters were extracted for all five limestones using the CIS model and a USC model with variable effective diffusivity. Differences in the sulfur capture performance of the limestones were evaluated. Plots of the CaSO{sub 4} or CaS product layer thickness as a function of conversion are given, and the relative importance of limestone porosity and internal surface is discussed.« less

Analyzing and modeling the kinetics of amyloid beta pores associated with Alzheimer’s disease pathology

DOE PAGES

Ullah, Ghanim; Demuro, Angelo; Parker, Ian; ...

2015-09-08

Amyloid beta (Aβ) oligomers associated with Alzheimer’s disease (AD) form Ca 2+-permeable plasma membrane pores, leading to a disruption of the otherwise well-controlled intracellular calcium (Ca 2+) homeostasis. The resultant up-regulation of intracellular Ca 2+ concentration has detrimental implications for memory formation and cell survival. The gating kinetics and Ca 2+ permeability of Aβ pores are not well understood. We have used computational modeling in conjunction with the ability of optical patch-clamping for massively parallel imaging of Ca 2+ flux through thousands of pores in the cell membrane of Xenopus oocytes to elucidate the kinetic properties of Aβ pores. Themore » fluorescence time-series data from individual pores were idealized and used to develop data-driven Markov chain models for the kinetics of the Aβ pore at different stages of its evolution. Our study provides the first demonstration of developing Markov chain models for ion channel gating that are driven by optical-patch clamp data with the advantage of experiments being performed under close to physiological conditions. As a result, we demonstrate the up-regulation of gating of various Ca 2+ release channels due to Aβ pores and show that the extent and spatial range of such up-regulation increases as Aβ pores with low open probability and Ca 2+ permeability transition into those with high open probability and Ca 2+ permeability.« less

Pore Formation During Solidification of Aluminum: Reconciliation of Experimental Observations, Modeling Assumptions, and Classical Nucleation Theory

NASA Astrophysics Data System (ADS)

Yousefian, Pedram; Tiryakioğlu, Murat

2018-02-01

An in-depth discussion of pore formation is presented in this paper by first reinterpreting in situ observations reported in the literature as well as assumptions commonly made to model pore formation in aluminum castings. The physics of pore formation is reviewed through theoretical fracture pressure calculations based on classical nucleation theory for homogeneous and heterogeneous nucleation, with and without dissolved gas, i.e., hydrogen. Based on the fracture pressure for aluminum, critical pore size and the corresponding probability of vacancies clustering to form that size have been calculated using thermodynamic data reported in the literature. Calculations show that it is impossible for a pore to nucleate either homogeneously or heterogeneously in aluminum, even with dissolved hydrogen. The formation of pores in aluminum castings can only be explained by inflation of entrained surface oxide films (bifilms) under reduced pressure and/or with dissolved gas, which involves only growth, avoiding any nucleation problem. This mechanism is consistent with the reinterpretations of in situ observations as well as the assumptions made in the literature to model pore formation.

Pore growth in U-Mo/Al dispersion fuel

NASA Astrophysics Data System (ADS)

Kim, Yeon Soo; Jeong, G. Y.; Sohn, D.-S.; Jamison, L. M.

2016-09-01

U-Mo/Al dispersion fuel is currently under development in the DOE's Material Management and Minimization program to convert HEU-fueled research reactors to LEU-fueled reactors. In some demanding conditions in high-power and high-performance reactors, large pores form in the interaction layers between the U-Mo fuel particles and the Al matrix, which pose a potential to cause fuel failure. In this study, comprehension of the formation and growth of these pores was explored. As a product, a model to predict pore growth and porosity increase was developed. The model includes three major topics: fission gas release from the U-Mo and the IL to the pores, stress evolution in the fuel meat, and the effect of amorphous IL growth. Well-characterized in-pile data from reduced-size plates were used to fit the model parameters. A data set from full-sized plates, independent and distinctively different from those used to fit the model parameters, was used to examine the accuracy of the model. The model showed fair agreement with the measured data. The model suggested that the growth of the IL has a critical effect on pore growth, as both its material properties and energetics are favorable to pore formation. Therefore, one area of the current effort, focused on suppressing IL growth, appears to be on the right track to improve the performance of this fuel.

Prediction of porosity of food materials during drying: Current challenges and directions.

PubMed

Joardder, Mohammad U H; Kumar, C; Karim, M A

2017-07-18

Pore formation in food samples is a common physical phenomenon observed during dehydration processes. The pore evolution during drying significantly affects the physical properties and quality of dried foods. Therefore, it should be taken into consideration when predicting transport processes in the drying sample. Characteristics of pore formation depend on the drying process parameters, product properties and processing time. Understanding the physics of pore formation and evolution during drying will assist in accurately predicting the drying kinetics and quality of food materials. Researchers have been trying to develop mathematical models to describe the pore formation and evolution during drying. In this study, existing porosity models are critically analysed and limitations are identified. Better insight into the factors affecting porosity is provided, and suggestions are proposed to overcome the limitations. These include considerations of process parameters such as glass transition temperature, sample temperature, and variable material properties in the porosity models. Several researchers have proposed models for porosity prediction of food materials during drying. However, these models are either very simplistic or empirical in nature and failed to consider relevant significant factors that influence porosity. In-depth understanding of characteristics of the pore is required for developing a generic model of porosity. A micro-level analysis of pore formation is presented for better understanding, which will help in developing an accurate and generic porosity model.

Investigation of Coupled model of Pore network and Continuum in shale gas

NASA Astrophysics Data System (ADS)

Cao, G.; Lin, M.

2016-12-01

Flow in shale spanning over many scales, makes the majority of conventional treatment methods disabled. For effectively simulating, a coupled model of pore-scale and continuum-scale was proposed in this paper. Based on the SEM image, we decompose organic-rich-shale into two subdomains: kerogen and inorganic matrix. In kerogen, the nanoscale pore-network is the main storage space and migration pathway so that the molecular phenomena (slip and diffusive transport) is significant. Whereas, inorganic matrix, with relatively large pores and micro fractures, the flow is approximate to Darcy. We use pore-scale network models (PNM) to represent kerogen and continuum-scale models (FVM or FEM) to represent matrix. Finite element mortars are employed to couple pore- and continuum-scale models by enforcing continuity of pressures and fluxes at shared boundary interfaces. In our method, the process in the coupled model is described by pressure square equation, and uses Dirichlet boundary conditions. We discuss several problems: the optimal element number of mortar faces, two categories boundary faces of pore network, the difference between 2D and 3D models, and the difference between continuum models FVM and FEM in mortars. We conclude that: (1) too coarse mesh in mortars will decrease the accuracy, while too fine mesh will lead to an ill-condition even singular system, the optimal element number is depended on boundary pores and nodes number. (2) pore network models are adjacent to two different mortar faces (PNM to PNM, PNM to continuum model), incidental repeated mortar nodes must be deleted. (3) 3D models can be replaced by 2D models under certain condition. (4) FVM is more convenient than FEM, for its simplicity in assigning interface nodes pressure and calculating interface fluxes. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the 973 Program (2014CB239004), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the National Natural Science Foundation of China (41574129).

Estimation and modeling of coal pore accessibility using small angle neutron scattering

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Rui; Liu, Shimin; Bahadur, Jitendra

Gas diffusion in coal is controlled by nano-structure of the pores. The interconnectivity of pores not only determines the dynamics of gas transport in the coal matrix but also influences the mechanical strength. In this study, small angle neutron scattering (SANS) was employed to quantify pore accessibility for two coal samples, one of sub-bituminous rank and the other of anthracite rank. Moreover, a theoretical pore accessibility model was proposed based on scattering intensities under both vacuum and zero average contrast (ZAC) conditions. Our results show that scattering intensity decreases with increasing gas pressure using deuterated methane (CD 4) at lowmore » Q values for both coals. Pores smaller than 40 nm in radius are less accessible for anthracite than sub-bituminous coal. On the contrary, when the pore radius is larger than 40 nm, the pore accessibility of anthracite becomes larger than that of sub-bituminous coal. Only 20% of pores are accessible to CD 4 for anthracite and 37% for sub-bituminous coal, where the pore radius is 16 nm. For these two coals, pore accessibility and pore radius follows a power-law relationship.« less

Estimation and modeling of coal pore accessibility using small angle neutron scattering

DOE PAGES

Zhang, Rui; Liu, Shimin; Bahadur, Jitendra; ...

2015-09-04

Gas diffusion in coal is controlled by nano-structure of the pores. The interconnectivity of pores not only determines the dynamics of gas transport in the coal matrix but also influences the mechanical strength. In this study, small angle neutron scattering (SANS) was employed to quantify pore accessibility for two coal samples, one of sub-bituminous rank and the other of anthracite rank. Moreover, a theoretical pore accessibility model was proposed based on scattering intensities under both vacuum and zero average contrast (ZAC) conditions. Our results show that scattering intensity decreases with increasing gas pressure using deuterated methane (CD 4) at lowmore » Q values for both coals. Pores smaller than 40 nm in radius are less accessible for anthracite than sub-bituminous coal. On the contrary, when the pore radius is larger than 40 nm, the pore accessibility of anthracite becomes larger than that of sub-bituminous coal. Only 20% of pores are accessible to CD 4 for anthracite and 37% for sub-bituminous coal, where the pore radius is 16 nm. For these two coals, pore accessibility and pore radius follows a power-law relationship.« less

Pore-throat radius and tortuosity estimation from formation resistivity data for tight-gas sandstone reservoirs

NASA Astrophysics Data System (ADS)

Ziarani, Ali S.; Aguilera, Roberto

2012-08-01

A new model is proposed for estimation of pore-throat aperture size from formation resistivity factor and permeability data. The model is validated with data from the Mesaverde sandstone using brine salinities ranging from 20,000 to 200,000 ppm. The data analyzed includes various basins such as Green River, Piceance, Sand Wash, Powder River, Uinta, Washakie and Wind River, available in the literature. For pore-throat radii analysis the methodology involves the use of log-log plots of pore-throat radius versus the product of formation resistivity factor and permeability (rT = a(FK)b + c). The model fits over 280 samples from the Mesaverde formation with coefficients of determination varying between 0.95 and 0.99 depending primarily on the type of model used for pore throat radius calculation. The brine salinity has some minor effects on the results. The model can provide better estimates of pore-throat radii if it is calibrated with experimental techniques such as mercury porosimetry. The results show pore-throat radii varying between 0.001 and 5 μm for the Mesaverde tight sandstone; however, most of the samples fall in a range between 0.01 and 1 μm. For tortuosity analysis, the calculation involves the use of product of formation factor and porosity data. Results indicate that the estimated tortuosity values range mainly between 1 and 5. For samples with lower porosities (< 5%), tortuosity values show a wider scatter (between 1 and 8); whereas for samples with larger porosities (> 15%), the scattering in tortuosity decreases significantly. In general, for tortuosity calculation in tight gas sandstone formations, a square root model with a parameter (bf) representing various types of connecting pores, i.e., sheet-like and tubular pores, is recommended.

Multi-scale Pore Imaging Techniques to Characterise Heterogeneity Effects on Flow in Carbonate Rock

NASA Astrophysics Data System (ADS)

Shah, S. M.

2017-12-01

Digital rock analysis and pore-scale studies have become an essential tool in the oil and gas industry to understand and predict the petrophysical and multiphase flow properties for the assessment and exploitation of hydrocarbon reserves. Carbonate reservoirs, accounting for majority of the world's hydrocarbon reserves, are well known for their heterogeneity and multiscale pore characteristics. The pore sizes in carbonate rock can vary over orders of magnitudes, the geometry and topology parameters of pores at different scales have a great impact on flow properties. A pore-scale study is often comprised of two key procedures: 3D pore-scale imaging and numerical modelling techniques. The fundamental problem in pore-scale imaging and modelling is how to represent and model the different range of scales encountered in porous media, from the pore-scale to macroscopic petrophysical and multiphase flow properties. However, due to the restrictions of image size vs. resolution, the desired detail is rarely captured at the relevant length scales using any single imaging technique. Similarly, direct simulations of transport properties in heterogeneous rocks with broad pore size distributions are prohibitively expensive computationally. In this study, we present the advances and review the practical limitation of different imaging techniques varying from core-scale (1mm) using Medical Computed Tomography (CT) to pore-scale (10nm - 50µm) using Micro-CT, Confocal Laser Scanning Microscopy (CLSM) and Focussed Ion Beam (FIB) to characterise the complex pore structure in Ketton carbonate rock. The effect of pore structure and connectivity on the flow properties is investigated using the obtained pore scale images of Ketton carbonate using Pore Network and Lattice-Boltzmann simulation methods in comparison with experimental data. We also shed new light on the existence and size of the Representative Element of Volume (REV) capturing the different scales of heterogeneity from the pore-scale imaging.

Indirect Reconstruction of Pore Morphology for Parametric Computational Characterization of Unidirectional Porous Iron.

PubMed

Kovačič, Aljaž; Borovinšek, Matej; Vesenjak, Matej; Ren, Zoran

2018-01-26

This paper addresses the problem of reconstructing realistic, irregular pore geometries of lotus-type porous iron for computer models that allow for simple porosity and pore size variation in computational characterization of their mechanical properties. The presented methodology uses image-recognition algorithms for the statistical analysis of pore morphology in real material specimens, from which a unique fingerprint of pore morphology at a certain porosity level is derived. The representative morphology parameter is introduced and used for the indirect reconstruction of realistic and statistically representative pore morphologies, which can be used for the generation of computational models with an arbitrary porosity. Such models were subjected to parametric computer simulations to characterize the dependence of engineering elastic modulus on the porosity of lotus-type porous iron. The computational results are in excellent agreement with experimental observations, which confirms the suitability of the presented methodology of indirect pore geometry reconstruction for computational simulations of similar porous materials.

The Effect of Porosity on Fatigue of Die Cast AM60

NASA Astrophysics Data System (ADS)

Yang, Zhuofei; Kang, Jidong; Wilkinson, David S.

2016-07-01

AM60 high-pressure die castings are known to contain significant porosity which can affect fatigue life. We have studied this using samples drawn from prototype AM60 shock towers by conducting strain-controlled fatigue tests accompanied by X-ray computed tomography analysis. The results show that the machined surface is the preferential location for fatigue crack development, with pores close to these surfaces serving as initiation sites. Fatigue life shows a strong inverse correlation with the size of the fatigue-crack-initiating pore. Pore shape and pore orientation also influence the response. A supplemental study on surface roughness shows that porosity is the dominant factor in fatigue. Tomography enables the link between porosity and fatigue crack initiation to be clearly identified. These data are complemented by SEM observations of the fracture surfaces which are generally flat and full of randomly oriented serration patterns but without long-range fatigue striations.

Measurement and modelling of reactive transport in geological barriers for nuclear waste containment

DOE PAGES

Xiong, Qingrong; Joseph, Claudia; Schmeide, Katja; ...

2015-10-26

Compacted clays are considered as excellent candidates for barriers to radionuclide transport in future repositories for nuclear waste due to their very low hydraulic permeability. Diffusion is the dominant transport mechanism, controlled by a nano-scale pore system. Assessment of the clays' long-term containment function requires adequate modelling of such pore systems and their evolution. Existing characterisation techniques do not provide complete pore space information for effective modelling, such as pore and throat size distributions and connectivity. Special network models for reactive transport are proposed here using the complimentary character of the pore space and the solid phase. Here, this balancesmore » the insufficient characterisation information and provides the means for future mechanical–physical–chemical coupling. The anisotropy and heterogeneity of clays is represented using different length parameters and percentage of pores in different directions. Resulting networks are described as mathematical graphs with efficient discrete calculus formulation of transport. Opalinus Clay (OPA) is chosen as an example. Experimental data for the tritiated water (HTO) and U(VI) diffusion through OPA are presented. Calculated diffusion coefficients of HTO and uranium species are within the ranges of the experimentally determined data in different clay directions. This verifies the proposed pore network model and validates that uranium complexes are diffusing as neutral species in OPA. In the case of U(VI) diffusion the method is extended to account for sorption and convection. Finally, rather than changing pore radii by coarse grained mathematical formula, physical sorption is simulated in each pore, which is more accurate and realistic.« less

Measurement and modelling of reactive transport in geological barriers for nuclear waste containment.

PubMed

Xiong, Qingrong; Joseph, Claudia; Schmeide, Katja; Jivkov, Andrey P

2015-11-11

Compacted clays are considered as excellent candidates for barriers to radionuclide transport in future repositories for nuclear waste due to their very low hydraulic permeability. Diffusion is the dominant transport mechanism, controlled by a nano-scale pore system. Assessment of the clays' long-term containment function requires adequate modelling of such pore systems and their evolution. Existing characterisation techniques do not provide complete pore space information for effective modelling, such as pore and throat size distributions and connectivity. Special network models for reactive transport are proposed here using the complimentary character of the pore space and the solid phase. This balances the insufficient characterisation information and provides the means for future mechanical-physical-chemical coupling. The anisotropy and heterogeneity of clays is represented using different length parameters and percentage of pores in different directions. Resulting networks are described as mathematical graphs with efficient discrete calculus formulation of transport. Opalinus Clay (OPA) is chosen as an example. Experimental data for the tritiated water (HTO) and U(vi) diffusion through OPA are presented. Calculated diffusion coefficients of HTO and uranium species are within the ranges of the experimentally determined data in different clay directions. This verifies the proposed pore network model and validates that uranium complexes are diffusing as neutral species in OPA. In the case of U(vi) diffusion the method is extended to account for sorption and convection. Rather than changing pore radii by coarse grained mathematical formula, physical sorption is simulated in each pore, which is more accurate and realistic.

Energetics and Self-Assembly of Amphipathic Peptide Pores in Lipid Membranes

PubMed Central

Zemel, Assaf; Fattal, Deborah R.; Ben-Shaul, Avinoam

2003-01-01

We present a theoretical study of the energetics, equilibrium size, and size distribution of membrane pores composed of electrically charged amphipathic peptides. The peptides are modeled as cylinders (mimicking α-helices) carrying different amounts of charge, with the charge being uniformly distributed over a hydrophilic face, defined by the angle subtended by polar amino acid residues. The free energy of a pore of a given radius, R, and a given number of peptides, s, is expressed as a sum of the peptides' electrostatic charging energy (calculated using Poisson-Boltzmann theory), and the lipid-perturbation energy associated with the formation of a membrane rim (which we model as being semitoroidal) in the gap between neighboring peptides. A simple phenomenological model is used to calculate the membrane perturbation energy. The balance between the opposing forces (namely, the radial free energy derivatives) associated with the electrostatic free energy that favors large R, and the membrane perturbation term that favors small R, dictates the equilibrium properties of the pore. Systematic calculations are reported for circular pores composed of various numbers of peptides, carrying different amounts of charge (1–6 elementary, positive charges) and characterized by different polar angles. We find that the optimal R's, for all (except, possibly, very weakly) charged peptides conform to the “toroidal” pore model, whereby a membrane rim larger than ∼1 nm intervenes between neighboring peptides. Only weakly charged peptides are likely to form “barrel-stave” pores where the peptides essentially touch one another. Treating pore formation as a two-dimensional self-assembly phenomenon, a simple statistical thermodynamic model is formulated and used to calculate pore size distributions. We find that the average pore size and size polydispersity increase with peptide charge and with the amphipathic polar angle. We also argue that the transition of peptides from the adsorbed to the inserted (membrane pore) state is cooperative and thus occurs rather abruptly upon a change in ambient conditions. PMID:12668433

Poring over two-pore channel pore mutants

PubMed Central

Penny, Christopher J.; Patel, Sandip

2016-01-01

Two-pore channels are members of the voltage-gated ion channel superfamily. They localise to the endolysosomal system and are likely targets for the Ca2+ mobilising messenger NAADP. In this brief review, we relate mutagenesis of the TPC pore to a recently published homology model and discuss how pore mutants are informing us of TPC function. Molecular physiology of these ubiquitous proteins is thus emerging. PMID:27226934

Resolving structural influences on water-retention properties of alluvial deposits

USGS Publications Warehouse

Winfield, K.A.; Nimmo, J.R.; Izbicki, J.A.; Martin, P.M.

2006-01-01

With the goal of improving property-transfer model (PTM) predictions of unsaturated hydraulic properties, we investigated the influence of sedimentary structure, defined as particle arrangement during deposition, on laboratory-measured water retention (water content vs. potential [??(??)]) of 10 undisturbed core samples from alluvial deposits in the western Mojave Desert, California. The samples were classified as having fluvial or debris-flow structure based on observed stratification and measured spread of particle-size distribution. The ??(??) data were fit with the Rossi-Nimmo junction model, representing water retention with three parameters: the maximum water content (??max), the ??-scaling parameter (??o), and the shape parameter (??). We examined trends between these hydraulic parameters and bulk physical properties, both textural - geometric mean, Mg, and geometric standard deviation, ??g, of particle diameter - and structural - bulk density, ??b, the fraction of unfilled pore space at natural saturation, Ae, and porosity-based randomness index, ??s, defined as the excess of total porosity over 0.3. Structural parameters ??s and Ae were greater for fluvial samples, indicating greater structural pore space and a possibly broader pore-size distribution associated with a more systematic arrangement of particles. Multiple linear regression analysis and Mallow's Cp statistic identified combinations of textural and structural parameters for the most useful predictive models: for ??max, including Ae, ??s, and ??g, and for both ??o and ??, including only textural parameters, although use of Ae can somewhat improve ??o predictions. Textural properties can explain most of the sample-to-sample variation in ??(??) independent of deposit type, but inclusion of the simple structural indicators Ae and ??s can improve PTM predictions, especially for the wettest part of the ??(??) curve. ?? Soil Science Society of America.

Pore‐Scale Hydrodynamics in a Progressively Bioclogged Three‐Dimensional Porous Medium: 3‐D Particle Tracking Experiments and Stochastic Transport Modeling

PubMed Central

Carrel, M.; Dentz, M.; Derlon, N.; Morgenroth, E.

2018-01-01

Abstract Biofilms are ubiquitous bacterial communities that grow in various porous media including soils, trickling, and sand filters. In these environments, they play a central role in services ranging from degradation of pollutants to water purification. Biofilms dynamically change the pore structure of the medium through selective clogging of pores, a process known as bioclogging. This affects how solutes are transported and spread through the porous matrix, but the temporal changes to transport behavior during bioclogging are not well understood. To address this uncertainty, we experimentally study the hydrodynamic changes of a transparent 3‐D porous medium as it experiences progressive bioclogging. Statistical analyses of the system's hydrodynamics at four time points of bioclogging (0, 24, 36, and 48 h in the exponential growth phase) reveal exponential increases in both average and variance of the flow velocity, as well as its correlation length. Measurements for spreading, as mean‐squared displacements, are found to be non‐Fickian and more intensely superdiffusive with progressive bioclogging, indicating the formation of preferential flow pathways and stagnation zones. A gamma distribution describes well the Lagrangian velocity distributions and provides parameters that quantify changes to the flow, which evolves from a parallel pore arrangement under unclogged conditions, toward a more serial arrangement with increasing clogging. Exponentially evolving hydrodynamic metrics agree with an exponential bacterial growth phase and are used to parameterize a correlated continuous time random walk model with a stochastic velocity relaxation. The model accurately reproduces transport observations and can be used to resolve transport behavior at intermediate time points within the exponential growth phase considered. PMID:29780184

Saturation-dependent solute dispersivity in porous media: Pore-scale processes

NASA Astrophysics Data System (ADS)

Raoof, A.; Hassanizadeh, S. M.

2013-04-01

It is known that in variably saturated porous media, dispersion coefficient depends on Darcy velocity and water saturation. In one-dimensional flow, it is commonly assumed that the dispersion coefficient is a linear function of velocity. The coefficient of proportionality, called the dispersivity, is considered to depend on saturation. However, there is not much known about its dependence on saturation. In this study, we investigate, using a pore network model, how the longitudinal dispersivity varies nonlinearly with saturation. We schematize the porous medium as a network of pore bodies and pore throats with finite volumes. The pore space is modeled using the multidirectional pore-network concept, which allows for a distribution of pore coordination numbers. This topological property together with the distribution of pore sizes are used to mimic the microstructure of real porous media. The dispersivity is calculated by solving the mass balance equations for solute concentration in all network elements and averaging the concentrations over a large number of pores. We have introduced a new formulation of solute transport within pore space, where we account for different compartments of residual water within drained pores. This formulation makes it possible to capture the effect of limited mixing due to partial filling of the pores under variably saturated conditions. We found that dispersivity increases with the decrease in saturation, it reaches a maximum value, and then decreases with further decrease in saturation. To show the capability of our formulation to properly capture the effect of saturation on solute dispersion, we applied it to model the results of a reported experimental study.

Pore-network model of evaporation-induced salt precipitation in porous media: The effect of correlations and heterogeneity

NASA Astrophysics Data System (ADS)

Dashtian, Hassan; Shokri, Nima; Sahimi, Muhammad

2018-02-01

Salt transport and precipitation in porous media constitute a set of complex and fascinating phenomena that are of considerable interest to several important problems, ranging from storage of CO2 in geological formations, to soil fertility, and protection of pavements and roads, as well as historical monuments. The phenomena occur at the pore scale and are greatly influenced by the heterogeneity of the pore space morphology. We present a pore-network (PN) model to study the phenomena. Vapor diffusion, capillary effect at the brine-vapor interface, flow of brine, and transport of salt and its precipitation in the pores that plug the pores partially or completely are all accounted for. The drying process is modeled by the invasion percolation, while transport of salt in brine is accounted for by the convective-diffusion equation. We demonstrate that the drying patterns, the clustering and connectivity of the pore throats in which salt precipitation occurs, the saturation distribution, and the drying rate are all strongly dependent upon the pore-size distribution, the correlations among the pore sizes, and the anisotropy of the pore space caused by stratification that most natural porous media contain. In particular, if the strata are more or less parallel to the direction of injection of the gas that dries out the pore space (air, for example) and/or causes salt precipitation (CO2, for example), the drying rate increases significantly. Moreover, salt tends to precipitate in clusters of neighboring pores that are parallel to the open surface of the porous medium.

Investigating Hydrophilic Pores in Model Lipid Bilayers Using Molecular Simulations: Correlating Bilayer Properties with Pore-Formation Thermodynamics.

PubMed

Hu, Yuan; Sinha, Sudipta Kumar; Patel, Sandeep

2015-06-23

Cell-penetrating and antimicrobial peptides show a remarkable ability to translocate across physiological membranes. Along with factors such as electric-potential-induced perturbations of membrane structure and surface tension effects, experiments invoke porelike membrane configurations during the solute transfer process into vesicles and cells. The initiation and formation of pores are associated with a nontrivial free-energy cost, thus necessitating a consideration of the factors associated with pore formation and the attendant free energies. Because of experimental and modeling challenges related to the long time scales of the translocation process, we use umbrella sampling molecular dynamics simulations with a lipid-density-based order parameter to investigate membrane-pore-formation free energy employing Martini coarse-grained models. We investigate structure and thermodynamic features of the pore in 18 lipids spanning a range of headgroups, charge states, acyl chain lengths, and saturation. We probe the dependence of pore-formation barriers on the area per lipid, lipid bilayer thickness, and membrane bending rigidities in three different lipid classes. The pore-formation free energy in pure bilayers and peptide translocating scenarios are significantly coupled with bilayer thickness. Thicker bilayers require more reversible work to create pores. The pore-formation free energy is higher in peptide-lipid systems than in peptide-free lipid systems due to penalties to maintain the solvation of charged hydrophilic solutes within the membrane environment.

Multi-Scale Multi-Physics Modeling of Matrix Transport Properties in Fractured Shale Reservoirs

NASA Astrophysics Data System (ADS)

Mehmani, A.; Prodanovic, M.

2014-12-01

Understanding the shale matrix flow behavior is imperative in successful reservoir development for hydrocarbon production and carbon storage. Without a predictive model, significant uncertainties in flowback from the formation, the communication between the fracture and matrix as well as proper fracturing practice will ensue. Informed by SEM images, we develop deterministic network models that couple pores from multiple scales and their respective fluid physics. The models are used to investigate sorption hysteresis as an affordable way of inferring the nanoscale pore structure in core scale. In addition, restricted diffusion as a function of pore shape, pore-throat size ratios and network connectivity is computed to make correct interpretation of the 2D NMR maps possible. Our novel pore network models have the ability to match sorption hysteresis measurements without any tuning parameters. The results clarify a common misconception of linking type 3 nitrogen hysteresis curves to only the shale pore shape and show promising sensitivty for nanopore structre inference in core scale. The results on restricted diffusion shed light on the importance of including shape factors in 2D NMR interpretations. A priori "weighting factors" as a function of pore-throat and throat-length ratio are presented and the effect of network connectivity on diffusion is quantitatively assessed. We are currently working on verifying our models with experimental data gathered from the Eagleford formation.

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.; Gao, Shengyan

2015-01-01

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir. PMID:26310236

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N; Gao, Shengyan

2015-08-27

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir.

Effects of Pore Distributions on Ductility of Thin-Walled High Pressure Die-Cast Magnesium

DOE Office of Scientific and Technical Information (OSTI.GOV)

Choi, Kyoo Sil; Li, Dongsheng; Sun, Xin

2013-06-01

In this paper, a microstructure-based three-dimensional (3D) finite element modeling method is adopted to investigate the effects of porosity in thin-walled high pressure die-cast (HPDC) Magnesium alloys on their ductility. For this purpose, the cross-sections of AM60 casting samples are first examined using optical microscope and X-ray tomography to obtain the general information on the pore distribution features. The experimentally observed pore distribution features are then used to generate a series of synthetic microstructure-based 3D finite element models with different pore volume fractions and pore distribution features. Shear and ductile damage models are adopted in the finite element analyses tomore » induce the fracture by element removal, leading to the prediction of ductility. The results in this study show that the ductility monotonically decreases as the pore volume fraction increases and that the effect of ‘skin region’ on the ductility is noticeable under the condition of same local pore volume fraction in the center region of the sample and its existence can be beneficial for the improvement of ductility. The further synthetic microstructure-based 3D finite element analyses are planned to investigate the effects of pore size and pore size distribution.« less

Computational fluid dynamics models of conifer bordered pits show how pit structure affects flow.

PubMed

Schulte, Paul J

2012-02-01

• The flow of xylem sap through conifer bordered pits, particularly through the pores in the pit membrane, is not well understood, but is critical for an understanding of water transport through trees. • Models solving the Navier-Stokes equation governing fluid flow were based on the geometry of bordered pits in black spruce (Picea mariana) and scanning electron microscopy images showing details of the pores in the margo of the pit membrane. • Solutions showed that the pit canals contributed a relatively small fraction of resistance to flow, whereas the torus and margo pores formed a large fraction, which depended on the structure of the individual pit. The flow through individual pores in the margo was strongly dependent on pore area, but also on the radial location of the pore with respect to the edge of the torus. • Model results suggest that only a few per cent of the pores in the margo account for nearly half of the flow and these pores tend to be located in the inner region of the margo where their contribution will be maximized. A high density of strands in outer portions of the margo (hence narrower pores) may be more significant for mechanical support of the torus. © 2011 The Author. New Phytologist © 2011 New Phytologist Trust.

Micro-CT Pore Scale Study Of Flow In Porous Media: Effect Of Voxel Resolution

NASA Astrophysics Data System (ADS)

Shah, S.; Gray, F.; Crawshaw, J.; Boek, E.

2014-12-01

In the last few years, pore scale studies have become the key to understanding the complex fluid flow processes in the fields of groundwater remediation, hydrocarbon recovery and environmental issues related to carbon storage and capture. A pore scale study is often comprised of two key procedures: 3D pore scale imaging and numerical modelling techniques. The essence of a pore scale study is to test the physics implemented in a model of complicated fluid flow processes at one scale (microscopic) and then apply the model to solve the problems associated with water resources and oil recovery at other scales (macroscopic and field). However, the process of up-scaling from the pore scale to the macroscopic scale has encountered many challenges due to both pore scale imaging and modelling techniques. Due to the technical limitations in the imaging method, there is always a compromise between the spatial (voxel) resolution and the physical volume of the sample (field of view, FOV) to be scanned by the imaging methods, specifically X-ray micro-CT (XMT) in our case In this study, a careful analysis was done to understand the effect of voxel size, using XMT to image the 3D pore space of a variety of porous media from sandstones to carbonates scanned at different voxel resolution (4.5 μm, 6.2 μm, 8.3 μm and 10.2 μm) but keeping the scanned FOV constant for all the samples. We systematically segment the micro-CT images into three phases, the macro-pore phase, an intermediate phase (unresolved micro-pores + grains) and the grain phase and then study the effect of voxel size on the structure of the macro-pore and the intermediate phases and the fluid flow properties using lattice-Boltzmann (LB) and pore network (PN) modelling methods. We have also applied a numerical coarsening algorithm (up-scale method) to reduce the computational power and time required to accurately predict the flow properties using the LB and PN method.

Passive advection-dispersion in networks of pipes: Effect of connectivity and relationship to permeability

NASA Astrophysics Data System (ADS)

Bernabé, Y.; Wang, Y.; Qi, T.; Li, M.

2016-02-01

The main purpose of this work is to investigate the relationship between passive advection-dispersion and permeability in porous materials presumed to be statistically homogeneous at scales larger than the pore scale but smaller than the reservoir scale. We simulated fluid flow through pipe network realizations with different pipe radius distributions and different levels of connectivity. The flow simulations used periodic boundary conditions, allowing monitoring of the advective motion of solute particles in a large periodic array of identical network realizations. In order to simulate dispersion, we assumed that the solute particles obeyed Taylor dispersion in individual pipes. When a particle entered a pipe, a residence time consistent with local Taylor dispersion was randomly assigned to it. When exiting the pipe, the particle randomly proceeded into one of the pipes connected to the original one according to probabilities proportional to the outgoing volumetric flow in each pipe. For each simulation we tracked the motion of at least 6000 solute particles. The mean fluid velocity was 10-3 ms-1, and the distance traveled was on the order of 10 m. Macroscopic dispersion was quantified using the method of moments. Despite differences arising from using different types of lattices (simple cubic, body-centered cubic, and face-centered cubic), a number of general observations were made. Longitudinal dispersion was at least 1 order of magnitude greater than transverse dispersion, and both strongly increased with decreasing pore connectivity and/or pore size variability. In conditions of variable hydraulic radius and fixed pore connectivity and pore size variability, the simulated dispersivities increased as power laws of the hydraulic radius and, consequently, of permeability, in agreement with previously published experimental results. Based on these observations, we were able to resolve some of the complexity of the relationship between dispersivity and permeability.

An improved method for permeability estimation of the bioclastic limestone reservoir based on NMR data

NASA Astrophysics Data System (ADS)

Ge, Xinmin; Fan, Yiren; Liu, Jianyu; Zhang, Li; Han, Yujiao; Xing, Donghui

2017-10-01

Permeability is an important parameter in formation evaluation since it controls the fluid transportation of porous rocks. However, it is challengeable to compute the permeability of bioclastic limestone reservoirs by conventional methods linking petrophysical and geophysical data, due to the complex pore distributions. A new method is presented to estimate the permeability based on laboratory and downhole nuclear magnetic resonance (NMR) measurements. We divide the pore space into four intervals by the inflection points between the pore radius and the transversal relaxation time. Relationships between permeability and percentages of different pore intervals are investigated to investigate influential factors on the fluid transportation. Furthermore, an empirical model, which takes into account of the pore size distributions, is presented to compute the permeability. 212 core samples in our case show that the accuracy of permeability calculation is improved from 0.542 (SDR model), 0.507 (TIM model), 0.455 (conventional porosity-permeability regressions) to 0.803. To enhance the precision of downhole application of the new model, we developed a fluid correction algorithm to construct the water spectrum of in-situ NMR data, aiming to eliminate the influence of oil on the magnetization. The result reveals that permeability is positively correlated with percentages of mega-pores and macro-pores, but negatively correlated with the percentage of micro-pores. Poor correlation is observed between permeability and the percentage of meso-pores. NMR magnetizations and T2 spectrums after the fluid correction agree well with laboratory results for samples saturated with water. Field application indicates that the improved method provides better performance than conventional models such as Schlumberger-Doll Research equation, Timur-Coates equation, and porosity-permeability regressions.

A universal model for nanoporous carbon supercapacitors

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Jingsong; Sumpter, Bobby G; Meunier, Vincent

2009-01-01

Supercapacitors based on nanoporous carbon materials, commonly called electric double-layer capacitors (EDLCs), are emerging as a novel type of energy-storage device with the potential to substitute batteries in applications that require high power densities. Nanoporous carbon supercapacitors are generally viewed as a parallel-plate capacitor since supercapacitors store energy by charge separation in an electric double layer formed at the electrode/electrolyte interface. The EDLC model has been used to characterize the energy storage of supercapacitors for decades. We comment in this chapter on the shortcomings of the EDLC model when applied to nanoporous carbon supercapacitors. In response to the latest experimentalmore » breakthrough in nanoporous carbon supercapacitors, we have proposed a heuristic model that takes pore curvature into account as a replacement for the EDLC model. When the pore size is in the mesopore regime (2 50 nm), electrolyte counterions enter mesoporous carbons and approach the pore wall to form an electric double-cylinder capacitor (EDCC); in the micropore regime (< 2 nm), solvated/desolvated counterions line up along the pore axis to form an electric wire-in-cylinder capacitor (EWCC). In the macropore regime (> 50 nm), where pores are large enough so that pore curvature is no longer significant, the EDCC model can be reduced to the EDLC model. With the backing of experimental data and quantum density functional theory calculations, we have shown that the EDCC/EWCC model is universal for carbon supercapacitors with diverse carbon materials and electrolytes. The strengths and limitations of this new model are discussed. The new model allows the supercapacitor properties to be correlated with pore size, specific surface area, Debye length, electrolyte concentration, dielectric constant, and solute ion size, and may lend support to the systematic optimization of the properties of carbon supercapacitors through experiments.« less

An Engineered Membrane to Measure Electroporation: Effect of Tethers and Bioelectronic Interface

PubMed Central

Hoiles, William; Krishnamurthy, Vikram; Cranfield, Charles G.; Cornell, Bruce

2014-01-01

This article reports on the construction and predictive models for a platform comprised of an engineered tethered membrane. The platform provides a controllable and physiologically relevant environment for the study of the electroporation process. The mixed self-assembled membrane is formed via a rapid solvent exchange technique. The membrane is tethered to the gold electrode and includes an ionic reservoir separating the membrane and gold surface. Above the membrane, there is an electrolyte solution, and a gold counterelectrode. A voltage is applied between the gold electrodes and the current measured. The current is dependent on the energy required to form aqueous pores and the conductance of each pore. A two-level predictive model, consisting of a macroscopic and a continuum model, is developed to relate the pore dynamics to the measured current. The macroscopic model consists of an equivalent circuit model of the tethered membrane, and asymptotic approximations to the Smoluchowski-Einstein equation of electroporation that is dependent on the pore conductance and the energy required to form aqueous pores. The continuum model is a generalized Poisson-Nernst-Planck (GPNP) system where an activity coefficient to account for steric effects of ions is added to the standard PNP system. The GPNP is used to evaluate the conductance of aqueous pores, and the electrical energy required to form the pores. As an outcome of the setup of the device and the two-level model, biologically important variables can be estimated from experimental measurements. To validate the accuracy of the two-level model, the predicted current is compared with experimentally measured current for different tethering densities. PMID:25229142

Fabrication of channeled scaffolds with ordered array of micro-pores through microsphere leaching and indirect Rapid Prototyping technique.

PubMed

Tan, J Y; Chua, C K; Leong, K F

2013-02-01

Advanced scaffold fabrication techniques such as Rapid Prototyping (RP) are generally recognized to be advantageous over conventional fabrication methods in terms architectural control and reproducibility. Yet, most RP techniques tend to suffer from resolution limitations which result in scaffolds with uncontrollable, random-size pores and low porosity, albeit having interconnected channels which is characteristically present in most RP scaffolds. With the increasing number of studies demonstrating the profound influences of scaffold pore architecture on cell behavior and overall tissue growth, a scaffold fabrication method with sufficient architectural control becomes imperative. The present study demonstrates the use of RP fabrication techniques to create scaffolds having interconnected channels as well as controllable micro-size pores. Adopted from the concepts of porogen leaching and indirect RP techniques, the proposed fabrication method uses monodisperse microspheres to create an ordered, hexagonal closed packed (HCP) array of micro-pores that surrounds the existing channels of the RP scaffold. The pore structure of the scaffold is shaped using a single sacrificial construct which comprises the microspheres and a dissolvable RP mold that were sintered together. As such, the size of pores as well as the channel configuration of the scaffold can be tailored based on the design of the RP mold and the size of microspheres used. The fabrication method developed in this work can be a promising alternative way of preparing scaffolds with customized pore structures that may be required for specific studies concerning cell-scaffold interactions.

Influence of pore structure on compressive strength of cement mortar.

PubMed

Zhao, Haitao; Xiao, Qi; Huang, Donghui; Zhang, Shiping

2014-01-01

This paper describes an experimental investigation into the pore structure of cement mortar using mercury porosimeter. Ordinary Portland cement, manufactured sand, and natural sand were used. The porosity of the manufactured sand mortar is higher than that of natural sand at the same mix proportion; on the contrary, the probable pore size and threshold radius of manufactured sand mortar are finer. Besides, the probable pore size and threshold radius increased with increasing water to cement ratio and sand to cement ratio. In addition, the existing models of pore size distribution of cement-based materials have been reviewed and compared with test results in this paper. Finally, the extended Bhattacharjee model was built to examine the relationship between compressive strength and pore structure.

Influence of Pore Structure on Compressive Strength of Cement Mortar

PubMed Central

Zhao, Haitao; Xiao, Qi; Huang, Donghui

2014-01-01

This paper describes an experimental investigation into the pore structure of cement mortar using mercury porosimeter. Ordinary Portland cement, manufactured sand, and natural sand were used. The porosity of the manufactured sand mortar is higher than that of natural sand at the same mix proportion; on the contrary, the probable pore size and threshold radius of manufactured sand mortar are finer. Besides, the probable pore size and threshold radius increased with increasing water to cement ratio and sand to cement ratio. In addition, the existing models of pore size distribution of cement-based materials have been reviewed and compared with test results in this paper. Finally, the extended Bhattacharjee model was built to examine the relationship between compressive strength and pore structure. PMID:24757414

Generation of dense plume fingers in saturated-unsaturated homogeneous porous media

NASA Astrophysics Data System (ADS)

Cremer, Clemens J. M.; Graf, Thomas

2015-02-01

Flow under variable-density conditions is widespread, occurring in geothermal reservoirs, at waste disposal sites or due to saltwater intrusion. The migration of dense plumes typically results in the formation of vertical plume fingers which are known to be triggered by material heterogeneity or by variations in source concentration that causes the density variation. Using a numerical groundwater model, six perturbation methods are tested under saturated and unsaturated flow conditions to mimic heterogeneity and concentration variations on the pore scale in order to realistically generate dense fingers. A laboratory-scale sand tank experiment is numerically simulated, and the perturbation methods are evaluated by comparing plume fingers obtained from the laboratory experiment with numerically simulated fingers. Dense plume fingering for saturated flow can best be reproduced with a spatially random, time-constant perturbation of the solute source. For unsaturated flow, a spatially and temporally random noise of solute concentration or a random conductivity field adequately simulate plume fingering.

Evaluating the hydraulic and transport properties of peat soil using pore network modeling and X-ray micro computed tomography

NASA Astrophysics Data System (ADS)

Gharedaghloo, Behrad; Price, Jonathan S.; Rezanezhad, Fereidoun; Quinton, William L.

2018-06-01

Micro-scale properties of peat pore space and their influence on hydraulic and transport properties of peat soils have been given little attention so far. Characterizing the variation of these properties in a peat profile can increase our knowledge on the processes controlling contaminant transport through peatlands. As opposed to the common macro-scale (or bulk) representation of groundwater flow and transport processes, a pore network model (PNM) simulates flow and transport processes within individual pores. Here, a pore network modeling code capable of simulating advective and diffusive transport processes through a 3D unstructured pore network was developed; its predictive performance was evaluated by comparing its results to empirical values and to the results of computational fluid dynamics (CFD) simulations. This is the first time that peat pore networks have been extracted from X-ray micro-computed tomography (μCT) images of peat deposits and peat pore characteristics evaluated in a 3D approach. Water flow and solute transport were modeled in the unstructured pore networks mapped directly from μCT images. The modeling results were processed to determine the bulk properties of peat deposits. Results portray the commonly observed decrease in hydraulic conductivity with depth, which was attributed to the reduction of pore radius and increase in pore tortuosity. The increase in pore tortuosity with depth was associated with more decomposed peat soil and decreasing pore coordination number with depth, which extended the flow path of fluid particles. Results also revealed that hydraulic conductivity is isotropic locally, but becomes anisotropic after upscaling to core-scale; this suggests the anisotropy of peat hydraulic conductivity observed in core-scale and field-scale is due to the strong heterogeneity in the vertical dimension that is imposed by the layered structure of peat soils. Transport simulations revealed that for a given solute, the effective diffusion coefficient decreases with depth due to the corresponding increase of diffusional tortuosity. Longitudinal dispersivity of peat also was computed by analyzing advective-dominant transport simulations that showed peat dispersivity is similar to the empirical values reported in the same peat soil; it is not sensitive to soil depth and does not vary much along the soil profile.

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.