An introduction to flow and transport in fractal models of porous media: Part I

This special issue gathers together a number of recent papers on fractal geometry and its applications to the modeling of flow and transport in porous media. The aim is to provide a systematic approach for analyzing the statics and dynamics of fluids in fractal porous media by means of theory, modeling and experimentation. The topics covered include lacunarity analyses of multifractal and natural grayscale patterns, random packing's of self-similar pore/particle size distributions, Darcian and non-Darcian hydraulic flows, diffusion within fractals, models for the permeability and thermal conductivity of fractal porous media and hydrophobicity and surface erosion properties of fractal structures.

Comparison Between Dcrit Considering the Abrupt Variation and Inflexion in the Concrete Mercury Intrusion Porosimetry Curve

Mercury intrusion porosimetry (MIP) has been widely used to evaluate the quality of concrete through the pore size distribution parameters. Two of these parameters are the critical pore diameter (Dcrit) and the percentage of the most interconnected net of pores compared to the total volume of pores. Some researchers consider Dcrit as the diameter obtained from the inflexion point of the cumulative mercury intrusion curve while others consider Dcrit as the diameter obtained from the point of abrupt variation in the same curve. This study aims to analyze two groups of concretes of varying w/c ratios, one cast with pozzolanic cement and another with high initial strength cement, in order to determine which of these diameters feature a better correlation with the quality parameters of the concretes. The concrete quality parameters used for the evaluations were (1) the w/c ratios and (2) chloride diffusion coefficients measured at approximately 90 days. MIP cumulative distributions of the same concretes were also measured at about 90 days, and Dcrit values were determined (1) from the point of abrupt variation and alternatively, (2) from the inflexion point of each of these plots. It was found that Dcrit values measured from the point of abrupt variation were useful indicators of the quality of the concrete, but the Dcrit values based on the inflexion points were not. Hence, it is recommended that Dcrit and the percentage of the most interconnected volume of pores should be obtained considering the point of abrupt variation of the cumulative curve of pore size distribution.

Petrophysical measurement of ODP Leg 122 and 123 samples (Tables 1, 2)

The capillary-pressure characteristics of 22 samples of lithified post-Paleozoic Indian-Ocean carbonates were compared to published data from older carbonate rocks (lower Paleozoic Hunton Group of Texas and Oklahoma). The Indian-Ocean samples are considerably more porous than are the Paleozoic samples, yet all of the Indian-Ocean samples fit readily into a descriptive petrofacies scheme previously established for the Hunton Group. The Indian-Ocean samples may be assigned to four petrophysical facies (petrofacies) based on the shapes of their capillary-pressure curves, their pore-throat-size distributions, their estimated recovery efficiency values (for nonwetting fluids), and the visual characteristics of their pore systems, as observed with a scanning-electron microscope. Petrofacies assignments for the Indian-Ocean samples are as follows. Petrofacies I includes six samples collected from the coarse basal portions of event deposits (primarily turbidites). These samples have large throats, leptokurtic throat-size distributions, low- to moderate recovery efficiency values, concave cumulative-intrusion capillary-pressure curves, and high porosity values. Petrofacies II includes two sedimentologically dissimilar samples that have medium-size throats, platykurtic throat-size distributions, moderate- to-high recovery efficiency values, gently sloping cumulative-intrusion capillary-pressure curves, and high porosity values. Petrofacies III includes two polymictic sandstones and a skeletal packstone that have small throats, polymodal throat-size distributions, moderate recovery efficiency values, gently sloping cumulative-intrusion capillary-pressure curves, and high porosity values. Petrofacies IV includes 11 samples, mostly recrystallized neritic carbonates, that have small throats, leptokurtic throat-size distributions, high recovery efficiency values, convex cumulative-intrusion capillary-pressure curves, and low porosity values. Comparison of petrofacies assignment to core-, thin-section-, and smear-slide data, and to inferred depositional setting, suggests that pore systems in most samples from Holes 765C and 766A result from primary depositional features, whereas pore systems in samples from Hole 761C and one sample from Hole 765C have been strongly influenced by diagenetic processes. For Hole 761C, prediction of petrophysical parameters should be most successful if based on diagenetic facies patterns. By contrast, the distribution of favorable reservoir facies and of permeability barriers in less highly altered rocks collected from Holes 765C and 766A is related to depositional patterns. Recovery efficiency is inversely related to both porosity and median throat size for the present data set. This relationship is similar to that observed for carbonates of the lower Paleozoic Hunton Group and the Ordovician Ellenburger dolomite, but opposite of that observed for some other ancient carbonates. The coarse deposits of the massive basal units of turbidites are petrophysically distinct and form a coherent petrophysical group (Petrofacies I) with substantial reservoir potential. Two samples assigned to Petrofacies I have extremely large throats (median throat size at least 4 ?m, and at least six times that of any other sample) and therefore high permeability values. These two samples come from thin, coarse turbidites that lack or have poorly developed fine divisions and are interpreted to have been deposited on channeled suprafan lobes in a proximal mid-fan setting. The restriction of extremely high permeability values to a single depositional facies suggests that careful facies mapping of deep-sea fans in a deliberate search for such coarse turbidites could dramatically enhance the success of exploration for aquifers or hydrocarbon reservoirs. Such reservoirs should have substantial vertical heterogeneity. They should have high lateral permeability values but low vertical permeability values, and reservoir sections should include numerous thin units having widely differing petrophysical characteristics.

Physical Characteristics and Estimated Hydraulic Conductivity from Grain-Size Distribution for Vashon Till, South-Central Puget Lowland, Washington

Population growth, urban development, and increased commercial and industrial activity in the south-central Puget Lowlands of Washington State has led to an increased demand for groundwater. The Vashon till is a glacially consolidated, low-permeability unit comprising unstratified clay, silt, cobbles and boulders with ubiquitous coarse-grained lenses and is an extensive surficial unit throughout the south-central Puget Lowland. Thus, understanding the physical and hydrological characteristics – specifically, the hydraulic conductivity – of this unit is a necessary component of a groundwater model. This study provides (1) a record of the physical characteristics of Vashon till deposits within the study area; and (2) an estimate of the highest, lowest, and average value of saturated hydraulic conductivity based on the grain-size distribution of Vashon till samples collected from six field sites in the Puyallup River Watershed. Analysis shows that the average moisture content ranges between about 1 and 6%, average dry bulk density is about 2.20 g/cm3, and average porosity is about 17%. Grain-size distributions show that half of the samples analyzed are well graded, while the other half is poorly graded. Grain-size distributions also show an average d10 value of about 0.20 mm, and average ff values ≤ 16%, which are key values in estimating the saturated hydraulic conductivity of over-consolidated glacial deposits. Based on these observed values, the estimates of hydraulic conductivity range from a minimum of 0.02 m/d to a maximum of 1.38 m/d in within the general Vashon till.

SMPS particle size data from the Mt. Bachelor Observatory for summer 2015

Aerosol size distributions from the Mt. Bachelor Observatory for 7/30/2015 to 9/23/2015.

Pore characterization of carbonaceous materials by DFT and GCMC simulations: A review

A review is given of the pore characterization of carbonaceous materials, including activated carbon, carbon fibres, carbon nanotubes, etc., using adsorption techniques. Since the pores of carbon media are mostly of molecular dimensions, the appropriate modem tools for the analysis of adsorption isotherms are grand canonical Monte Carlo (GCMC) simulations and density functional theory (DFT). These techniques are presented and applications of such tools in the derivation of pore-size distribution highlighted.

Study of hexane adsorption in nanoporous MCM-41 silica

We study here the adsorption of hexane on nanoporous MCM-41 silica at 303, 313, and 323 K, for various pore diameters between 2.40 and 4.24 nm. Adsorption equilibria, measured thermogravimetrically, show that all the isotherms, that are somewhat akin to those of type V, exhibit remarkably sharp capillary adsorption phase transition steps and are reversible. The position of the phase transition step gradually shifts from low to high relative pressure with an increase in the temperature as well as the pore sizes. The isosteric heats of adsorption derived from the equilibrium information using the Clapeyron equation reveal a gradual decrease with increasing adsorbed amount because of the surface heterogeneity but approach a constant value near the phase transition. A decrease in the pore size results in an increase in the isosteric heat of adsorption because of the increased dispersion forces. A simple strategy, based on the Broekhoff and De Boer adsorption theory, successfully interprets the hexane adsorption isotherms for the different pore size MCM-41 samples. The parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting the monolayer region prior to capillary condensation and the experimental phase transition simultaneously, for some pore sizes. Subsequently, the parameters are used to predict the adsorption isotherm on other pore size samples, which showed good agreement with experimental data.

Flow regime, hydrodynamics, floc size distribution and sludge properties in activated sludge bubble column, air-lift and aerated stirred reactors

This paper presents a comparative study how reactor configuration, sludge loading and air flowrate affect flow regimes, hydrodynamics, floc size distribution and sludge solids-liquid separation properties. Three reactor configurations were studied in bench scale activated sludge bubble column reactor (BCR), air-lift reactor (ALR) and aerated stirred reactor (ASR). The ASR demonstrated the highest capacity of gas holdup and resistance, and homogeneity in flow regimes and shearing forces, resulting in producing large numbers of small and compact floes. The fluid dynamics in the ALR created regularly directed recirculation forces to enhance the gas holdup and sludge flocculation. The BCR distributed a high turbulent flow regime and non-homogeneity in gas holdup and mixing, and generated large numbers of larger and looser floes. The sludge size distributions, compressibility and settleability were significantly influenced by the reactor configurations associated with the flow regimes and hydrodynamics.

Influence of charge size distribution on net-power draw of tumbling mill based on DEM modelling

Modelling and optimization of the power draw of large SAG/AG mills is important due to the large power draw which modern mills require (5-10 MW). The cost of grinding is the single biggest cost within the entire process of mineral extraction. Traditionally, modelling of the mill power draw has been done using empirical models. Although these models are reliable, they cannot model mills and operating conditions which are not within the model database boundaries. Also, due to its static nature, the impact of the changing conditions within the mill on the power draw cannot be determined using such models. Despite advances in computing power, discrete element method (DEM) modelling of large mills with many thousands of particles could be a time consuming task. The speed of computation is determined principally by two parameters: number of particles involved and material properties. The computational time step is determined by the size of the smallest particle present in the model and material properties (stiffness). In the case of small particles, the computational time step will be short, whilst in the case of large particles; the computation time step will be larger. Hence, from the point of view of time required for modelling (which usually corresponds to time required for 3-4 mill revolutions), it will be advantageous that the smallest particles in the model are not unnecessarily too small. The objective of this work is to compare the net power draw of the mill whose charge is characterised by different size distributions, while preserving the constant mass of the charge and mill speed. (C) 2004 Elsevier Ltd. All rights reserved.

Characterization of activated carbon fibers using argon adsorption

In this paper, we present results of the internal structure (pore size and pore wall thickness distributions) of a series of activated carbon fibers with different degrees of burn-off, determined from interpretation of argon adsorption data at 87 K using infinite and finite wall thickness models. The latter approach has recently been developed in our laboratory. The results show that while the low bun-off samples have nearly uniform pore size (

Application of density functional theory to equilibrium adsorption of argon and nitrogen on amorphous silica surface

We present a new version of non-local density functional theory (NL-DFT) adapted to description of vapor adsorption isotherms on amorphous materials like non-porous silica. The novel feature of this approach is that it accounts for the roughness of adsorbent surface. The solid–fluid interaction is described in the same framework as in the case of fluid–fluid interactions, using the Weeks–Chandler–Andersen (WCA) scheme and the Carnahan–Starling (CS) equation for attractive and repulsive parts of the Helmholtz free energy, respectively. Application to nitrogen and argon adsorption isotherms on non-porous silica LiChrospher Si-1000 at their boiling points, recently published by Jaroniec and co-workers, has shown an excellent correlative ability of our approach over the complete range of pressures, which suggests that the surface roughness is mostly the reason for the observed behavior of adsorption isotherms. From the analysis of these data, we found that in the case of nitrogen adsorption short-range interactions between oxygen atoms on the silica surface and quadrupole of nitrogen molecules play an important role. The approach presented in this paper may be further used in quantitative analysis of adsorption and desorption isotherms in cylindrical pores such as MCM-41 and carbon nanotubes.