Statistical Synthesis of Imaging and Porosimetry Data for the Characterization of Microstructure and Transport Properties of Sandstones

The microstructure of a suite of sandstone samples is quantitatively analyzed using a method which combines information from thin section micrographs of the pore space with mercury injection porosimetry in a statistical framework. This method enables the determination of a continuous distribution of pore sizes ranging from few nanometre to several hundred micrometre. The data obtained unify fractal and Euclidean aspects of the void space geometry, yield estimates of the pore-to-throat aspect ratio and challenge the ability of commonly used network models to describe fluid percolation in multiscale porous media. Application of critical path analysis to the prediction of flow permeability and electrical conductivity of sandstone core samples using the new information produces results comparable to those obtained by the classical approach—a fact attributed to the presence of macroscopic heterogeneity at the scale of several millimetres.     

Characterization of the Mixing Layer Resulting from the Detonation of Heterogeneous Explosive Charges

A dense, two-phase numerical methodology is used to study the mixing layer developing behind the detonation of a heterogeneous explosive charge, i.e., a charge comprising of a high explosive with metal particles. The filtered Navier–Stokes equations are solved in addition to a sub-grid kinetic energy equation, along with a recently developed Eulerian–Lagrangian formulation to handle dense flow-fields. The mixing layer resulting from the post-detonation phase of the explosion of a nitromethane charge consisting of inert steel particles is of interest in this study. Significant mixing and turbulence effects are observed in the mixing layer, and the rms of the radial velocity component is found to be about 25% higher than that of the azimuthal and zenith velocity components due to the flow being primarily radial. The mean concentration profiles are self-similar in shape at different times, based on a scaling procedure used in the past for a homogeneous explosive charge. The peak rms of concentration profiles are 23–30% in intensity and decrease in magnitude with time. The behavior of concentration gradients in the mixing layer is investigated, and stretching along the radial direction is observed to decrease the concentration gradients along the azimuth and zenith directions faster than the radial direction. The mixing and turbulence effects in the mixing layer subsequent to the detonation of the heterogeneous explosive charge are superior to that of a homogeneous explosive charge containing the same amount of the high explosive, exemplifying the role played by the particles in perturbing the flow-field. The non-linear growth of the mixing layer width starts early for the heterogeneous explosive charge, and the rate is reduced during the implosion phase in comparison with the homogeneous charge. The turbulence intensities in the mixing layer for the heterogeneous explosive charge are found to be nearly independent of the particle size for two different sizes considered in the initial charge. Overall, this study has provided some useful insights on the mixing layer characteristics subsequent to the detonation of heterogeneous explosives, and has also demonstrated the efficacy of the dense, multiphase formulation for such applications.     

Experimental Characterization of Porosity Structure and Transport Property Changes in Limestone Undergoing Different Dissolution Regimes

Limestone dissolution by $\hbox {CO}_2$ -rich brine induces critical changes of the pore network geometrical parameters such as the pore size distribution, the connectivity, and the tortuosity which govern the macroscopic transport properties (permeability and dispersivity) that are required to parameterize the models, simulating the injection and the fate of $\hbox {CO}_2$ . A set of four reactive core-flood experiments reproducing underground conditions ( $T = 100\,^{\circ }\hbox {C}$ and $P = 12$ MPa) has been conducted for different $\hbox {CO}_2$ partial pressures $(0.034 < P_{\mathrm{CO}_2}< 3.4\; \hbox {MPa})$ in order to study the different dissolution regimes. X-ray microtomographic images have been used to characterize the changes in the structural properties from pore scale to Darcy scale, while time-resolved pressure loss and chemical fluxes enabled the determination of the sample-scale change in porosity and permeability. The results show the growth of localized dissolution features associated with high permeability increase for the highest $P_{\mathrm{CO}_2}$ , whereas dissolution tends to be more homogeneously distributed for lower values of $P_{\mathrm{CO}_2}$ . For the latter, the higher the $P_{\mathrm{CO}_2}$ , the more the dissolution patterns display ramified structures and permeability increase. For the lowest value of $P_{\mathrm{CO}_2}$ , the preferential dissolution of the calcite cement associated with the low dissolution kinetics triggers the transport that may locally accumulate and form a microporous material that alters permeability and produces an anti-correlated porosity–permeability relationship. The combined analysis of the pore network geometry and the macroscopic measurements shows that $P_{\mathrm{CO}_2}$ regulates the tortuosity change during dissolution. Conversely, the increase of the exponent value of the observed power law permeability–porosity trend while $P_{\mathrm{CO}_2}$ increases, which appears to be strongly linked to the increase of the effective hydraulic diameter, depends on the initial rock structure.     

Nonlinear Inversion of an Unconfined Aquifer: Simultaneous Estimation of Heterogeneous Hydraulic Conductivities,Recharge Rates,and Boundary Conditions

A new inverse method is developed to simultaneously estimate heterogeneous hydraulic conductivities, source/sink rates, and unknown boundary conditions for steady-state flow in an unconfined aquifer. Unlike objective function-based techniques, the new method does not optimize any data-model misfits. Instead, its formulation is developed by honoring physical flow principles as well as observation data at sampled locations. Under the Dupuit–Forchheimer assumption of negligible vertical flow, accuracy and stability of the new method are demonstrated using synthetic heterogeneous aquifer problems with increasingly complex flow: (1) aquifer domains without source/sink effects; (2) aquifer domains with a point sink (a pumping well operating under a constant discharge rate); (3) aquifer domains with constant or spatially variable recharge; (4) aquifer domains with constant or spatially variable recharge undergoing single-well pumping. For all problems, inversion yields stable solutions under increasing head measurement errors (up to $\pm $ 10 % of the total head variation in a problem), although accuracy of the estimated parameters degrades with the increasing errors. The inverse method is successfully tested on problems with high hydraulic conductivity contrasts—up to 10,000 times between the maximum and minimum values. In inverting several heterogeneous problems, if the aquifer is assumed homogeneous with a constant recharge rate, physically meaningful parameter estimates (i.e., equivalent conductivities and mean recharge rates) can be determined. Alternatively, if the inverse parameterization contains spurious parameters, inversion can identify such parameters, while the simultaneous estimation of non-spurious parameters is not affected. The method obviates the well-known issues associated with model “structure errors”, when inverse parameterization either simplifies or complexifies the true parameter field.     

Characterization and Monitoring of Porous Media with Electrical Imaging: A Review

Transport in Porous Media - Electrical geophysical imaging is a widely used noninvasive technology for visualizing porous media at scales larger than individual pores. Originally developed for...     

In situ observations of unusual drop deformation and wobbling in simple shear flow

Deformation and wobbling of a liquid drop immersed in a liquid matrix were studied under mild shear conditions for various viscosity ratios. In situ visualization experiments were conducted on a homemade transparent Couette cell incorporated to the Paar Physica MCR500 shear rheometer. The effect of drop or matrix elasticity was examined and was found to play a major role in both deformation and wobbling processes. Experimental results were compared to Jackson and Tucker (J Rheol 47:659–682, 2003), Maffettone and Minale (J Non-Newton Fluid Mech 78:227–241, 1998) and Yu and Bousmina (J Rheol 47:1011–1039, 2003) ellipsoidal models. It was found that the agreement between the Newtonian models and the experimental results required an increase in the drop viscosity. Such increment in viscosity was found to scale with the first normal stress difference.     

In situ synthesis of magnetic mesoporous silica via sol-gel process coupled with precipitation and oxidation

A novel method was proposed for the in situ synthesis of magnetite-containing mesoporous silica SBA- 16 via a sol-gel process coupled with precipitation and oxidation. The effect of the added amounts of reactants on the mesostructural and magnetic properties of the magnetic mesoporous silica was investigated. It was determined that the synthesized magnetic mesoporous silica with a total pore volume of 0.64-0.96 cm 3 /g and an average pore diameter of 4.0-14.9 nm had a relatively high saturation magnetizatio...     

Impact of Capillary Pressure, Salinity and In situ Conditions on CO2 Injection into Saline Aquifers

This article presents a numerical investigation of the combined effects of capillary pressure, salinity and in situ thermodynamic conditions on CO₂-brine-rock interactions in a saline aquifer. We demonstrate that the interrelations between capillary pressure, salinity, dissolution and drying-out affect CO₂ injectivity and storage capacity of a saline aquifer. High capillary forces require a high injection pressure for a given injection rate. Depending on salinity, the increase in injection pressure due to capillary forces can be offset by the dissolution of CO₂ in formation water and its compressibility. Higher capillary forces also reduce gravity segregation, and this gives a more homogeneous CO₂ plume which improves the dissolution of CO₂. The solubility of CO₂ in formation water decreases with increasing salinity which requires an increased injection pressure. Higher salinity and capillary pressure can even block the pores, causing an increased salt precipitation. Simulations with various pressure-temperature conditions and modified salinity and capillary pressure curves demonstrate that, with the injection pressures similar for both cold and warm basins at a given injection rate, CO₂ dissolves about 10% more in the warm basin water than in the cold basin. The increase in dissolution lowers the injection pressure compensating the disadvantage of low CO₂ density and compressibility for storage in warm basins.     

In situ 3D Synchrotron Laminography Assessment of Edge Fracture in Dual-Phase Steels: Quantitative and Numerical Analysis

The mechanical performance of automotive structures made of advanced high strength steels (AHSS) is often seen reduced by the presence of cut edges. An attempt is made to assess and quantify the initial damage state and the damage evolution during mechanical testing of a punched edge and a machined edge via a recently developed 3D imaging technique called synchrotron radiation computed laminography. This technique allows us to observe damage in regions of interest in thin sheet-like objects at micrometer resolution. In terms of new experimental mechanics, steel sheets having sizes and mechanical boundary conditions of engineering relevance can be tested for the first time with in situ 3D damage observation and quantification. It is found for the investigated DP600 steel that the fracture zone of the punched edge is rough and that needle-shape voids at the surface and in the bulk follow ferrite-martensite flow lines. During mechanical in situ testing the needle voids grow from the fracture zone surface and coalesce with the sheared zone. In contrast, during in situ mechanical testing of a machined edge the damage starts away from the edge (～800μm) where substantial necking has occurred. Three-dimensional image analysis was performed to quantify the initial damage and its evolution. These data can be used as input and validation data for micromechanical damage models. To interpret the experimental findings in terms of mechanical fields, combined surface digital image correlation and 3D finite element analysis were carried out using an elasto-plastic constitutive law of the investigated DP steel. The stress triaxiality and the accumulated plastic strain were calculated in order to understand the influence of the edge profile and the hardening of the cutting-affected zone on the mechanical fields.     

Effect of Sudden Expansion on Entrainment and Spreading Rates of a Jet Issuing from Asymmetric Nozzles

Turbulent free jets issuing from five different nozzle geometries; smooth pipe, contracted circular, rectangular, triangular, and square, are experimentally investigated by using TSI 2-D laser Doppler velocimetry (LDV) to assess the effect of nozzle geometry and quarl (i.e. a cylindrical sudden expansion) on jet entrainment and spreading. The centerline mean velocity decay and the jet half-velocity width, which are indicators of jet entrainment and spreading rates, are determined for each nozzle’s flow configuration, i.e. with and without sudden expansion. Furthermore, turbulence quantities, such as the flow mean velocities and their mean fluctuating components, as well as Reynolds shear stresses, are all measured along the centerline plane of the jet to facilitate understanding the extent of the effect of nozzle’s geometry (i.e. nozzle’s orifice shape and sudden expansion) on jet’s entrainment and spreading. The main results show that the jet flow with the presence of sudden expansion exhibits higher rates of entrainment and spreading than without. In addition, these results reveal that sudden expansion exercises a greater effect on the asymmetric jet characteristics, especially for the triangular and rectangular nozzles compared to their axisymmetric counterparts (i.e. circular contracted nozzle).     

Statistical modeling and optimization of a multistage gas-solid fluidized bed for removing pollutants from flue gases

The present paper describes the statistical modeling and optimization of a multistage gas-solid fluidized bed reactor for the control of hazardous pollutants in flue gas.In this work,we study the hydrodynamics of the pressure drop and minimum fluidization velocity.The hydrodynamics of a three-stage fluidized bed are then compared with those for a single-stage unit.It is observed that the total pressure drop over all stages of the three-stage fluidized bed is less than that of an identical single-stage system.However,the minimum fluidization velocity is higher in the single-stage unit.Under identical conditions,the minimum fluidization velocity is highest in the top bed,and lowest in the bottom bed.This signifies that the behavior of solids changes from a well-mixed flow to a plug-flow,with intermediate behavior in the middle bed.     

In situ synthesis of CNTs/Fe-Ni/TiO2 nanocomposite by fluidized bed chemical vapor deposition and the synergistic effect in photocatalysis

Carbon nanotube(CNTs)/Fe-Ni/TiO_2 nanocomposite photocatalysts have been synthesized by an in situ fluidized bed chemical vapor deposition(FBCVD) method.The composite photocatalysts were characterized by XRD,Raman spectroscopy,BET,FESEM,TEM,UV-vis spectroscopy,and XPS.The results showed that the CNTs were grown in situ on the surface of TiO_2.Fe(lll) in TiO_2 showed no chemical changes in the growth of CNTs.Ni(Ⅱ) was partly reduced to metal Ni in the FBCVD process,and the metal Ni acted as a catalyst for the growth of CNTs.The photocatalytic activities of CNTs/Fe-Ni/TiO_2 decreased with the rise of the FBCVD reaction temperature.For the sample synthesized at low FBCVD temperature(500℃),more than 90% and nearly 50% of methylene blue were removed under UV irradiation in 180 min and under visible light irradiation in 300 min,respectively.The probable mechanism of synergistic enhancement of photocatalysis on the CNTs/Fe-Ni/TiO_2 nanocomposite is proposed.     

A comparison of ground vehicle mobility analysis based on soil moisture time series datasets from WindSat,LIS, and in situ sensors

Soil moisture is a key terrain variable in ground vehicle off-road mobility. Historically, models of the land water balance have been used to estimate soil moisture. Recently, satellites have provided another source of soil moisture estimates that can be used to estimate soil-limited vehicle mobility. In this study, we compared the off-road vehicle mobility estimates based on three soil moisture sources: WindSat (a satellite source), LIS (a computer model source), and in situ ground sensors (to represent ground truth). Mobility of six vehicles, each with different ranges of sensitivity to soil moisture, was examined in three test sites. The results demonstrated that the effect of the soil moisture error on mobility predictions is complex and may produce very significant errors in off-road mobility analysis for certain combinations of vehicles, seasons, and climates. This is because soil moisture biases vary in both direction and magnitude with season and location. Furthermore, vehicles are sensitive to different ranges of soil moistures. Modeled vehicle speeds in the dry time periods were limited by the interaction between soil traction and the vehicles’ powertrain characteristics. In the wet season, differences in soil strength resulted in more significant differences in mobility predictions.