Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ROCK”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

ISRM Suggested Method: Determining Deformation and Failure Characteristics of Rocks Subjected to True Triaxial Compression

The purpose of this ISRM Suggested Method is to introduce a guideline on determining deformation and failure characteristics of rocks subjected to true triaxial compression on different stress path. Here, the true triaxial testing apparatus was reviewed by means of the function and engineering application. Some key techniques, such as stress and strain measurements, and reduction of end effect between specimen and metal platens, preventing metal platens interference, were stated and suggested in detail. Methodology of specimen processing, specimen shape, and testing procedure are characterized. There is an explanation of the experimental data processing on stress–strain curves, strength, and fracture mode.

42 ENGINEERING↗

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗

Evaluation of the predictive capability of coupled thermo-hydro-mechanical models for a heated bentonite/clay system (HE-E) in the Mont Terri Rock Laboratory

Process understanding and parameter identification using numerical methods based on experimental findings are a key aspect of the international cooperative project DECOVALEX. Comparing the predictions from numerical models against experimental results increases confidence in the site selection and site evaluation process for a radioactive waste repository in deep geological formations. In the present phase of the project, DECOVALEX-2015, eight research teams have developed and applied models for simulating an in-situ heater experiment HE-E in the Opalinus Clay in the Mont Terri Rock Laboratory in Switzerland. The modelling task was divided into two study stages, related to prediction and interpretation of the experiment. A blind prediction of the HE-E experiment was performed based on calibrated parameter values for both the Opalinus Clay, that were based on the modelling of another in-situ experiment (HE-D), and modelling of laboratory column experiments on MX80 granular bentonite and a sand/bentonite mixture .. After publication of the experimental data, additional coupling functions were analysed and considered in the different models. Moreover, parameter values were varied to interpret the measured temperature, relative humidity and pore pressure evolution. The analysis of the predictive and interpretative results reveals the current state of understanding and predictability of coupled THM behaviours associated with geologic nuclear waste disposal in clay formations.

54 ENVIRONMENTAL SCIENCES↗

A New Workflow of X-ray CT Image Processing and Data Analysis of Structural Features in Rock Using Open-Source Software

X-ray computed tomography (CT) images of rock specimens often contain artifacts which must be corrected before scientific analyses are performed. Here, we present a new workflow of automated image processing to utilize poor-quality X-ray CT scan images. The workflow runs on the open-source image analysis software and efficiently separates desired features from low-contrast scanned images. The new workflow is a two-step technique using contrast enhancement and automated feature segmentation to generate noise-free binary images. The results of binary images using the proposed workflow and using a conventional thresholding technique are analyzed to show the quality of the proposed method. The paper also presents a workflow of estimating the structural geometries of features in two and three dimensions. The results of the structural feature analyses and computational time were compared between the open-source (ImageJ) and commercial image analysis software (Bruker Computed Tomography Analyzer). The commercial software was more computationally efficient, but the task-specific macros in open-source software enabled the user-desired automation in image processing and data extraction of desired structural features of comparable quality.

47 OTHER INSTRUMENTATION↗

Rapid imbibition of water in fractures within unsaturated sedimentary rock

The spontaneous imbibition of water and other liquids into gas-filled fractures in variably-saturated porous media is important in a variety of engineering and geological contexts. However, surprisingly few studies have investigated this phenomenon. In this paper, we present a theoretical framework for predicting the 1-dimensional movement of water into air-filled fractures within a porous medium based on early-time capillary dynamics and spreading over the rough surfaces of fracture faces. The theory permits estimation of sorptivity values for the matrix and fracture zone, as well as a dispersion parameter which quantifies the extent of spreading of the wetting front. Quantitative data on spontaneous imbibition of water in unsaturated Berea sandstone cores were acquired to evaluate the proposed model. The cores with different permeability classes ranging from 50 to 500 mD and were fractured using the Brazilian method. Spontaneous imbibition in the fractured cores was measured by dynamic neutron radiography at the Neutron Imaging Prototype Facility (beam line CG-1D, HFIR), Oak Ridge National Laboratory. Water uptake into both the matrix and the fracture zone exhibited square-root-of-time behavior. The matrix sorptivities ranged from 2.9 to 4.6 mm s -0.5 , and increased linearly as the permeability class increased. The sorptivities of the fracture zones ranged from 17.9 to 27.1 mm s -0.5 , and increased linearly with increasing fracture aperture width. The dispersion coefficients ranged from 23.7 to 66.7 mm 2 s -1 and increased linearly with increasing fracture aperture width and damage zone width. Both theory and observations indicate that fractures can significantly increase spontaneous imbibition in unsaturated sedimentary rock by capillary action and surface spreading on rough fracture faces. Fractures also increase the dispersion of the wetting front. In conclusion, further research is needed to investigate this phenomenon in other natural and engineered porous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mineral-specific heterogeneous neptunium sorption onto geological repository rocks in oxic and anoxic conditions and various temperatures

Neptunium poses challenges for radioactive wastes management due to its long half-life and ingrowth from americium decay, as well as its high mobility in water in its oxidized forms. Np(V)-bearing aqueous solutions were contacted with polished and crushed gneiss and dolerite from the Yeniseisky site, Russia, which is under consideration as a location for an underground repository, both in oxic and anoxic conditions at 20 and 90 °C. Neptunium sorption on crushed gneiss and dolerite in oxic conditions at 20 °C was quite low, but increased substantially in anoxic conditions. Autoradiography on samples treated at 20 °C revealed that sorbed Np was associated mostly with dark colored iron-bearing minerals, and solvent extraction indicated partial reduction of Np(V) to Np(IV). Sorption of Np from aqueous solutions contacted with dolerite and gneiss at 90 °C was in all cases higher than at the lower temperature. Autoradiography of polished rock samples treated in anoxic conditions at 90 °C revealed that Fe-bearing minerals had 3–4 orders of magnitude more Np sorbed on them than the quartz and feldspar grains. X-ray photoelectron spectroscopy indicated this sorbed Np was present as NpO 2 . This study supports the concept that the geologic conditions present at the planned repository site will retard Np migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fate and transport of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and its degradation products in sedimentary and volcanic rocks, Los Alamos, New Mexico

High-explosive compounds including hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) were used extensively in weapons research and testing at Los Alamos National Laboratory (LANL). Liquid effluents containing RDX were released to an outfall pond that flowed to Cañon de Valle at LANL's Technical Area 16 (TA-16), resulting in the contamination of the alluvial, intermediate and regional groundwater bodies. Monitoring of groundwater within Cañon de Valle has shown persistent RDX in the intermediate perched zone located between 225 and 311 m below ground surface. Monitoring data also show detectable levels of RDX putative anaerobic degradation products. Batch and column experiments were conducted to determine the extent of adsorption-desorption and transport of RDX and its degradation products (MNX, DNX, and TNX) in major rock types that are within the RDX plume. All experiments in this paper were performed in the dark using water obtained from a well located at the center of the plume, which is fairly oxic and has a neutral pH of 7.5. Retardation factors and partitioning coefficient (K d ) values for RDX were calculated from batch experiments. Additionally, retardation factors and K d values for RDX and its degradation products were calibrated from column experiments using a one-dimensional transport model with equilibrium sorption (linear isotherm). Results from the column and batch experiments showed little to no sorption of RDX to the aquifer materials tested, with retardation factors ranging from 1.0 to 1.8 and K d values varying from 0 to 0.70 L/kg. Finally, results also showed no measurable differences between the transport properties of RDX and its degradation products.

54 ENVIRONMENTAL SCIENCES↗

A mixed-mode phase field fracture model in anisotropic rocks with consistent kinematics

Under a pure tensile loading, cracks in brittle, isotropic, and homogeneous materials often propagate such that pure mode I kinematics are maintained at the crack tip. However, experiments performed on geo-materials, such as sedimentary rock, shale, mudstone, concrete and gypsum, often lead to the conclusion that the mode I and mode II critical fracture energies/surface energy release rates are distinctive. This distinction has great influence on the formation and propagation of wing cracks and secondary cracks from pre-existing flaws under a combination of shear and tensile or shear and compressive loadings. To capture the mixed-mode fracture propagation, a mixed-mode I/II fracture model that employs multiple critical energy release rates based on Shen and Stephansson, IJRMMS, 1993 is reformulated in a regularized phase field fracture framework. Here we obtain the mixed-mode driving force of the damage phase field by balancing the microforce. Meanwhile, the crack propagation direction and the corresponding kinematics modes are determined via a local fracture dissipation maximization problem. Several numerical examples that demonstrate mode II and mixed-mode crack propagation in brittle materials are presented. Possible extensions of the model capturing degradation related to shear/compressive damage, as commonly observed in sub-surface applications and triaxial compression tests, are also discussed.

42 ENGINEERING↗

Assessing the carbon sequestration potential of basalt using X-ray micro-CT and rock mechanics

Mineral carbonation in basaltic rock provides a permanent storage solution for the mitigation of anthropogenic CO 2 emissions in the atmosphere. 3D X-ray micro-CT (XCT) image analysis is applied to a core sample from the main basaltic reservoir of the CarbFix site in Iceland, which obtained a connected porosity of 2.05–8.76%, a reactive surface area of 0.10–0.33 mm -1 and a larger vertical permeability (2.07 × 10 -10 m 2 ) compared to horizontal permeability (5.10 × 10 -11 m 2 ). The calculations suggest a CO 2 storage capacity of 0.33 Gigatonnes at the CarbFix pilot site. The XCT results were compared to those obtained from a hydromechanical test applied to the same sample, during which permeability, electrical resistivity and volumetric deformation were measured under realistic reservoir pressure conditions. It was found that permeability is highly stress sensitive, dropping by two orders of magnitude for a -0.02% volumetric deformation change, equivalent to a mean pore diameter reduction of 5 μm. This pore contraction was insufficient to explain such a permeability reduction according to the XCT analysis, unless combined with the effects of clay swelling and secondary mineral pore clogging. Here the findings provide important benchmark data for the future upscaling and optimisation of CO 2 storage in basalt formations.

58 GEOSCIENCES↗

Noble and Major Gases Released from Rock Core Materials as Intrinsic Tracers for Detecting Carbon Dioxide Leakage – Laboratory Evaluation

The elemental and isotopic compositions of noble gases constitute a set of non-reactive geochemical tracers. Noble gases released from geological formations in contact with pore water acidified by leaked carbon dioxide (CO2) can potentially be employed as tracers for detection of CO2 leakage from storage reservoirs in carbon sequestration sites. In a set of laboratory tests, core materials from a stratigraphic borehole at the proposed FutureGen 2.0 carbon storage site (Jacksonville, IL, USA) were reacted with an aqueous acid to simulate contact with aqueous CO2 and the gas released from the reaction was collected. The major gas components, the noble gas concentrations and isotopic ratios were determined to evaluate the feasibility of using the released noble gases as tracers for crustal water-rock interactions. Experiments yielded several promising tracer candidates for intrinsic monitoring of CO2 leakage at the proposed FutureGen 2.0 site. The most robust signals stem from CH4 and 4He released by carbonates associated with dolomite and shale formations from the primary containment zone. No useful signal was obtained from the sandstone formation due to its very low yield of gas when acidified. Xenon released from the primary-containment zone also showed promise as a tracer, but a more robust analysis requires a different experimental set up. A general enrichment of CH4, N2 and light noble gas isotopes with respect to Ar abundance was observed and is consistent with the addition of major gases (CO2, CH4, and to a lesser extent N2) released by a combination of the decomposition of carbonate mineral phases and/or organic matter in the dolomite and shale formations. The enrichment of these major gases leads to enhanced signals from heterogeneous molecular and isotopic ratios including 3He/4He, 4He/22Ne, 20Ne/22Ne, 21Ne*/22Ne, 20Ne/36Ar, 84Kr/36Ar, 132Xe/36Ar, and 132Xe/84Kr, amongst others.

noble gases, tracers, carbon sequestration, CO2, i↗

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

Reactive transport of uranium in fractured crystalline rock: Upscaling in time and distance

In this study, batch adsorption and breakthrough column experiments were conducted to evaluate uranium transport through altered material that fills fractures in a granite rock system at the Grimsel Test Site in Switzerland at pH 6.9 and 7.9. The role of adsorption and desorption kinetics was evaluated with reactive transport modeling by comparing one-, two-, and three-site models. Emphasis was placed on describing long desorption tails that are important for upscaling in time and distance. The effect of increasing pH in injection solutions was also evaluated. For pH 6.9, a three-site model with forward rate constants between 0.07 and 0.8 ml g –1 h –1 , reverse rate constants between 0.001 and 0.06 h –1 , and site densities of 1.3, 0.104, and 0.026 μmol g –1 for ‘weak/fast’, ‘strong/slow’, and ‘very strong/very slow’ sites provided the best fits. For pH 7.9, a three-site model with forward rate constants between 0.05 and 0.8 mL g –1 h –1 , reverse rate constants between 0.001 and 0.6 h –1 , and site densities of 1.3, 0.039, and 0.013 μmol g –1 for a ‘weak/fast’, ‘strong/slow’, and ‘very strong/very slow’ sites provided the best fits. Column retardation coefficients (R d ) were 80 for pH 6.9 and 10.3 for pH 7.9. Model parameters determined from the batch and column experiments were used in 50 year large-scale simulations for continuous and pulse injections and indicated that a three-site model is necessary at pH 6.9, although a K d -type equilibrium partition model with one-site was adequate for large scale predictions at pH 7.9. Batch experiments were useful for predicting early breakthrough times in the columns while column experiments helped differentiate the relative importance of sorption sites and desorption rate constants on transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗

Machine learning-assisted upscaling analysis of reservoir rock core properties based on micro-computed tomography imagery

Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.

58 GEOSCIENCES↗

Optimization of Eisenia fetida stocking density for the bioconversion of rock phosphate enriched cow dung–waste paper mixtures

Highlights: • Vermidegradation of RP-enriched waste mixtures is dependent on E. fetida stocking density. • A stocking density of 12.5 g-worms kg{sup -1} resulted in highly humified vermicomposts. • P release from RP-enriched waste vermicomposts increases with E. fetida stocking density. • RP-enriched waste vermicomposts had no inhibitory effect on seed germination. - Abstract: Vermitechnology is gaining recognition as an environmental friendly waste management strategy. Its successful implementation requires that the key operational parameters like earthworm stocking density be established for each target waste/waste mixture. One target waste mixture in South Africa is waste paper mixed with cow dung and rock phosphate (RP) for P enrichment. This study sought to establish optimal Eisenia fetida stocking density for maximum P release and rapid bioconversion of RP enriched cow dung–paper waste mixtures. E. fetida stocking densities of 0, 7.5, 12.5, 17.5 and 22.5 g-worms kg{sup −1} dry weight of cow dung–waste paper mixtures were evaluated. The stocking density of 12.5 g-worms kg{sup −1} resulted in the highest earthworm growth rate and humification of the RP enriched waste mixture as reflected by a C:N ratio of <12 and a humic acid/fulvic acid ratio of >1.9 in final vermicomposts. A germination test revealed that the resultant vermicompost had no inhibitory effect on the germination of tomato, carrot, and radish. Extractable P increased with stocking density up to 22.5 g-worm kg{sup −1} feedstock suggesting that for maximum P release from RP enriched wastes a high stocking density should be considered.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reaction-Induced Fracturing of Porous Carbonate Rocks during Volume-Increasing Replacement by Witherite

Volume-increasing replacement reactions can lead either to fracturing of the parent phase or to the formation of cohesive layers that passivate further reaction, but the factors that drive one outcome or the other are not understood. In this experimental study, we investigated the volume-increasing replacement of carbonate rocks (consisting of CaCO 3 and CaMg(CO 3 ) 2 ) of different porosity by witherite (BaCO 3 ). Samples were characterized using scanning electron microscopy, Raman spectroscopy, small-angle neutron/X-ray scattering, X-ray tomography, and scanning transmission electron microscopy. We observed the formation of a witherite reaction layer and witherite formation within pores and along grain boundaries. Despite this being a volume-increasing replacement reaction, newly formed witherite was porous, potentially allowing further replacement. Filled fractures were observed in the low-porosity carbonates, whereas witherite formed within pores in high-porosity carbonates. We conclude that fracturing of the parent phase versus passivation is contingent on the initial microstructure of the parent with an optimal degree of porosity required for fracturing.

Calcite↗

Eu 5 Al 3 Sb 6 : Al 4 Tetrahedra Embedded in a Rock-Salt-Like Structure

The new Eu 5 Al 3 Sb 6 phase has been successfully synthesized as a pure phase through Sn flux methods yielding large, high-quality crystals. This structure type features disordered Al clusters that appear in the form of dual tetrahedra. It crystallizes in the monoclinic C2/m space group exhibiting a rock-salt-like Eu–Sb framework with [Al 4 ] tetrahedra replacing some of the cationic Eu atoms (space group: C2/m, a = 8.151(1) Å, b = 14.181(2) Å, c = 8.145(1) Å, β = 109.577(2)°). The structure models the [Al 4 ] as dual tetrahedra with the Al atom sites 37.5% occupied along with Eu present on the central site at 8% occupancy and the remainder of the site being vacant. The presence of the [Al 4 ] cluster is further supported by HRTEM. Electronic structure calculations show that this material is a semimetal with observed band crossings close to the Fermi level. Strong Al–Sb antibonding interactions were found from COHP calculations close to the Fermi level and provide the rationale for the deficiency of the Al cluster. Mössbauer spectroscopy on Eu-151 and Sb-121 provides oxidation states of 2+ and 3– along with the local environment. Magnetic susceptibility measurements can be described well with a Curie–Weiss law where an effective moment of 7.80 μB/mol Eu is obtained, consistent with Eu 2+ , and show canted antiferromagnetic behavior below 10 K. Temperature dependent resistivity shows a Kondo-like low-temperature upturn caused by enhanced scattering of the itinerant electrons with the 4f orbitals of Eu.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reaction Mechanism of the Synthesis of a Disordered Rock Salt Cathode Material

Li-rich disordered rock salt (DRX) materials have recently emerged as a class of promising cathode materials for cobalt- and nickel-free Li-ion batteries that can achieve a high energy density and high electrochemical storage capacity. Here, the synthesis conditions of DRX were found to have a significant impact on the electrochemical performance of the final product formed and hence need to be carefully controlled. To date, the formation pathway of DRX materials has remained largely unexplored. In this work, the phase evolution of a DRX material, Li 1.1 Mn 0.8 Ti 0.1 O 1.9 F 0.1 (LMTF1811), during solid-state synthesis, was monitored using in situ synchrotron X-ray diffraction (sXRD). We have observed several crystalline intermediates formed during the synthesis process before the full formation of DRX. Ex situ X- ray absorption spectroscopy (XAS) was used to assess the local structure and charge states of the transition metals. F K-edge XAS collected from the final product LMTF1811 shows that the degree of F incorporation into the DRX phase is low.

Diffraction↗