SEARCH · Engineering Papers
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.
Measured cation, anion, and silica are presented as well as calculations evaluating mass balance performed for 4 fracture sustainability tests under EGS conditions: Stripa granite, Metamudstone, Rhyolitic Ashflow Tuff, and Silicified Rhyolitic Tuff. Calculations of charge balance are included for water samples in the 4 tests.
Using laboratory slide-hold-slide experiments, at temperatures from 22 to 200 degrees C, to examine effects of fracture reactivation and quasi-static loading on the evolution of fluid transport properties of simulated fractures in Westerly granite. At all temperatures, the in-plane hydraulic transmissivity consistently decays during hold periods resulting in an overall reduction in transmissivity. During the first three to fifteen hours of an experiment, transmissivity decreases rapidly due to the generation of wear products, development of a sliding surface, and compaction of the resulting gouge. Once the sliding surface has developed, the long-term transmissivity decay rate at 22 and 100 degrees C is significantly lower than the transmissivity decay rate during the initial 3-15 hours of the experiment. However, at 200 degrees C, the decay of hydraulic transmissivity remains high throughout the experiment. The long-term decay of hydraulic transmissivity can be fitted with a power law model with more rapid reduction of hydraulic transmissivity at higher temperature. Periods of sliding on the fracture surface result in transient increases in the transmissivity, due to shear dilation, as is expected for Coulomb materials. These transients are superimposed on the long-term decay. When sliding ceases and a new hold period commences, there is a rapid reduction in transmissivity and return to the long-term rate of transmissivity decay. The rate of decay of the transmissivity transients is inversely correlated with temperature, in contrast to the long-term decay and the expected behavior for processes like subcritical crack growth and indentation creep. The higher decay rates that are observed during the initial 3-15 hours of the tests and following sliding, are associated with times that the porosity of the gouge is expected to be high. The difference in decay rates suggests that when the gouge is driven far from equilibrium by active shearing, densification may be dominated by a different mechanism from long-term compaction.
Included are results from shear reactivation experiments on laboratory faults pre-loaded close to failure and reactivated by the injection of fluid into the fault. The sample comprises a single-inclined-fracture (SIF) transecting a cylindrical sample of Westerly granite. All experiments are conducted at ambient temperature and follow a similar protocol: (i) application of confining stresses (3MPa) on the fault fully saturated with DI water, (ii) shear-mobilization through the increase of axial loading at a constant displacement rate until a post-peak steady-state condition is reached, (iii) reduction of axial loading and related shear stress to a prescribed fraction of the peak steady-state frictional strength (typically 60% to 90%, representing intermediate to high magnitudes) and (iv), fault reactivation triggered by a stepwise increase of pore pressure on the fault in 0.1 MPa increments held constant for 1-5 minutes. Mechanical data from three ISCO pumps connected to a Temco pressure vessel measure axial, confining, and fault-related parameters, including fluid pressure (kPa), fluid flow rate (mL/min), and axial displacement (mm). See included code for initial data analysis and visualization for select experiments. Resource names represent experiment numbers found in the "Read Me" file, which describes each experimental setup and parameters.
Included are a report and data from triaxial shear tests on simulated Sierra White fault gouge and a borehole stimulation experiment. This serves as a progress report for the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University (Utah FORGE R&D project 5-2557). Tests were done with Sierra White granite (SWG) to investigate slip mechanisms and to evaluate the fracture initiation and propagation from a geothermal borehole. Raw data is provided alongside a report on all experimental testing protocols, experimental design, and results.
This dataset contains experimental data from fluid injection experiments conducted to investigate the influence of temperature on fracture seismicity. The experiments were performed on granite samples from Utah FORGE. The samples were prepared with a 30-degree inclined fracture and subjected to controlled stress and temperature conditions. Data were collected under three distinct temperature settings: 24 C, 78 C, and 137 C. During the experiments, a constant confining pressure of 10 MPa and a constant shear stress at 80% of the shear strength of the sample were maintained. Pore pressure was incrementally increased at a rate of 300 kPa every three minutes to simulate fluid injection. Temperature was raised rapidly and then stabilized for the duration of each test. The dataset includes shear stress and displacement measurements under each temperature condition, along with supplementary figures illustrating the experimental setup and time-series plots of pressures and temperature.
In the present study, we extensively explored the high-pressure behaviors and vibrational properties of amblygonite LiAlPO 4 F with elevated pressures up to 34.3 GPa based on single-crystal X-ray diffraction measurements, Raman spectroscopy, and DFT calculations. The compressibility and elastic properties of amblygonite were determined first. Specifically, the obtained isothermal bulk modulus of LiAlPO 4 F is 128(4) GPa and the triclinic phase exhibited anisotropic compression with axial compressibility β c > β a > β b with a ratio of 1.11:1.00:1.20. The Raman spectra showed no indication of phase transformation and were used to obtained mode Grüneisen parameters. The average Grüneisen parameter for PO 4 tetrahedral sites was smaller than for the LiO 4 F sites. Our results provide new insights into the phase stability and elastic properties of lithium-fluorite granites at extreme conditions.
In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)
Around the repository for HLW underground, cementitious materials used for the repository may complicate the migration of radionuclides in the host rock by the alteration of groundwater to a high alkaline condition. Specially, the alkaline groundwater dissolves silicate minerals of host rocks but the plenty of silicic acid dissolved from host rocks in groundwater will become supersaturated and deposit on the surface of flow-paths after decreasing the pH value to about 8 with dispersion and dilution. The deposition of silicic acid is presumed to recover and improve the barrier performance of host rocks for the radionuclide migration because of the narrowing of flow paths and the decrease of the permeability. In this study, the deposition behavior of supersaturated silicic acid was examined by the deposition experiments and the deposition rate constants, k [m/s], were evaluated from the experimental data. As the minerals of the surface of flow-paths, this study focused on chlorite and sericite, hydrothermally-altered minerals distributed to fill fractures of granite. In the experiments, deposition behavior of silicic acid on the powder of chlorite or sericite was observed after the pH of the solution was adjusted to 8 in order to make the silicic acid supersaturated. The solution of CaCl{sub 2} and NaCl were added into a silicic acid to simulate Ca{sup 2+} from the cementitious materials of the repository, Na{sup +} of the saline groundwater in Japan. From the results, the value of k was estimated to be about the order of 10{sup -10}∼10{sup -9} m/s for both chlorite and sericite and the value was almost the same as those for amorphous silica in the previous study. Therefore, the deposition of silicic acid may be assessed by similar apparent deposition rate constants even in actual underground containing various silicate minerals. Moreover, the Damkoehler number which is the dimensionless number expressed by the ratio of characteristic time of mass transfer and chemical reaction was calculated to be 10{sup 4} or more from the value of k. This value shows that the deposition reaction of silicic acid was estimated as much more predominant than the groundwater flow under the situation simulated in this study. These results suggest that the deposition of supersaturated silicic acid on the flow paths composed with various silicate minerals will be able to retard the migration of radionuclides because of narrowing micro flow-paths. (authors)
Abstract The 2085 km 2 Jordan River Basin, and its seven sub‐catchments draining the Central Wasatch Range immediately east of Salt Lake City, UT, are home to an array of hydrologic, atmospheric, climatic and chemical research infrastructure that collectively forms the Wasatch Environmental Observatory (WEO). WEO is geographically nested within a wildland to urban land‐use gradient and built upon a strong foundation of over a century of discharge and climate records. A 2200 m gradient in elevation results in variable precipitation, temperature and vegetation patterns. Soil and subsurface structure reflect systematic variation in geology from granitic, intrusive to mixed sedimentary clastic across headwater catchments, all draining to the alluvial or colluvial sediments of the former Lake Bonneville. Winter snowfall and spring snowmelt control annual hydroclimate, rapid population growth dominates geographic change in lower elevations and urban gas and particle emissions contribute to episodes of severe air pollution in this closed‐basin. Long‐term hydroclimate observations across this diverse landscape provide the foundation for an expanding network of infrastructure in both montane and urban landscapes. Current infrastructure supports both basic and applied research in atmospheric chemistry, biogeochemistry, climate, ecology, hydrology, meteorology, resource management and urban redesign that is augmented through strong partnerships with cooperating agencies. These features allow WEO to serve as a unique natural laboratory for addressing research questions facing seasonally snow‐covered, semi‐arid regions in a rapidly changing world and an excellent facility for providing student education and research training.
Abstract Passage of downstream-migrating fishes through hydropower dams poses a management issue in rivers throughout the world. As such, regulations have been enacted at many locales to understand and limit fish losses by requiring dam passage survival estimation. However, limiting estimates of survival to the dam and immediate tailrace presents a challenge to researchers. In this study, two alternative release–recapture methods were used to estimate dam passage survival of juvenile salmonids through Lower Granite Dam on the Snake River, Washington, in 2018. One approach, the virtual/paired-release (ViPRe) model, was extensively used in the Columbia–Snake River basin at federally operated hydroelectric dams during 2010–2014. This existing approach uses three releases of tagged smolts to isolate dam passage survival, defined as survival from the upstream dam face to the tailrace mixing zone 1–2 km downstream of the dam. An alternative approach, the virtual release/dead-fish correction (ViRDCt) model, uses one live-release group paired with a release of dead tagged fish at the dam to estimate the same survival parameter. The alternative estimation approaches were tested on two spring out-migrating stocks, yearling Chinook Salmon Oncorhynchus tshawytscha and steelhead Oncorhynchus mykiss, and on a summer out-migrating stock of subyearling Chinook Salmon. The alternative estimates for these stocks were all within 1 SE. However, the ViRDCt model produced survival estimates with SEs that were 59% smaller than those of the existing ViPRe model and did so using 42% fewer tagged fish and one less downstream acoustic detection array. Because of the reduced sample size, nearly US$540,000 less in tag costs and 4,000 fewer labor-hours were required to implement the ViRDCt model compared to the ViPRe model. As such, the ViRDCt model represents a cost-effective and precise approach to characterizing dam passage survival.
Abstract Current noble gas detection systems for nuclear explosion monitoring are based on the detection of four radioxenon isotopes—Xe-131m, -133, -133m and -135. The data provided by radioxenon detection could be enhanced by other radionuclide signatures such as Ar-37. Activation of Ca-40 in rock by neutrons produces Ar-37, and monitoring for this additional nuclide could help distinguish detections of nuclear explosions from background sources, such as medical isotope production. This work studies the capabilities of a hypothetical argon detection network. A 10 kt explosion was modeled using MCNP and SCALE to determine the inventory of Ar-37 created in a representative granite rock layer, assuming either 0.1, 1 or 10% of the total inventory was released. The Ar-37 inventory was combined with atmospheric transport data from HYSPLIT compiled in a previous study, along with the detection limits of standard Ar-37 detection systems, to determine how many hypothetical monitoring stations would detect Ar-37 from an explosion. This method was repeated for 365 HYSPLIT data sets to create a year’s worth of hypothetical explosions, releases, and detections. The study quantified the average number of detections per release, the number of stations detecting Ar-37, and the possibility of detecting Ar-37 in coincidence with xenon.
In this work, we aim to verify the predictions of the numerical simulators, which are used for designing field-scale hydraulic stimulation experiments. Although a strong theoretical understanding of this process has been gained over the past few decades, numerical predictions of fracture propagation in low-permeability rocks still remains a challenge. Against this background, we performed controlled laboratory-scale hydraulic fracturing experiments in granite samples, which not only provides high-quality experimental data but also a well-characterized experimental set-up. Using the experimental pressure responses and the final fracture sizes as benchmark, we compared the numerical predictions of two coupled hydraulic fracturing simulators—CSMP and GEOS. Both the simulators reproduced the experimental pressure behavior by implementing the physics of Linear Elastic Fracture Mechanics (LEFM) and lubrication theory within a reasonable degree of accuracy. The simulation results indicate that even in the very low-porosity (1–2 %) and low-permeability (10 -18 m 2 - 10 -19 m 2 ) crystalline rocks, which are usually the target of EGS, fluid-loss into the matrix and unsaturated flow impacts the formation breakdown pressure and the post-breakdown pressure trends. Therefore, underestimation of such parameters in numerical modeling can lead to significant underestimation of breakdown pressure. The simulation results also indicate the importance of implementing wellbore solvers for considering the effect of system compressibility and pressure drop due to friction in the injection line. The varying injection rate as a result of decompression at the instant of fracture initiation affects the fracture size, while the entry friction at the connection between the well and the initial notch may cause an increase in the measured breakdown pressure.
Abstract In 2017 and 2019, injection testing was carried out in three zones in a vertical well in granite at the Frontier Observatory for Research in Geothermal Energy site near Milford, Utah, USA. In several injection cycles, flowback was implemented rather than shut-in. The goal was to explore an alternative to prolonged shut-in periods for inferring closure stress, formation compressibility, and formation permeability (permeability thickness product). The flowback procedures involved a cyclic flowback/shut-in, while pressure decreased. The flowback data are presented, and analyses are shown. The inferred closure stress(es) from flowback analyses are lower than for equivalent injection cycles that were strictly shut-in. Relatively high formation compressibility obtained from the flowback analyses indicates an extensive, fractured system. This study also includes numerical simulation of the flowback events. The numerical model shows that the rebound pressure is not necessarily the lower bound of the minimum principal stress. The signature of stiffness change can be identified as the process when the depletion mainly transitions from hydraulic fracture to natural fractures from numerical analysis. Overall, flowback potentially has advantages over shut-in because of the reduced time to closure.
Rock fracturing (cracking) is a universal process that drives and limits chemical degradation, sediment production and erosion, and deterioration of infrastructure. Despite extensive research gains in rock mechanics on one hand and geochronology on the other, there remains a glaring gap in our ability to understand the long term evolution of natural, in situ fractures. Here we develop a novel fracture exposure dating technique, grounded in modern advances in luminescence geochronology. We apply our new dating method to a granitic boulder from a glacial outwash terrace in California, US. We conclude that the longest, clast-splitting E-W fracture appeared shortly after the boulder's deposit, whereas the secondary N-S fracture appeared 5 ka after the deposition, approximately correlating with the Last Glacial Maximum and Younger Dryas periods of the region, respectively. However, dating of the third fracture (<< 50 µm width) which does not fully split the rock, is ambiguous due to negligible daylight penetration and poor determination of fracture width. The fracture dating method presented herein brings with it the potential to decipher relationships that are crucial for the interpretation and modeling of, for example, long-term landscape and atmospheric evolution relating rock weathering to climate change and erosion.
Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.
Allanite is commonly encountered as an accessory rare-earth silicate mineral in association with minerals such as garnet, biotite, and feldspar. It is distributed globally and occurs in igneous formations such as granites, pegmatites, and syenites, as well as in various metamorphic rocks such as schist, gneiss, and amphibolite. Moreover, it can be found in mineral veins formed through hydrothermal activity. While allanite has not yet been extensively utilized for the production of rare-earth elements, recent discoveries of high-grade rare-earth-rich allanite deposits in Wyoming, USA, highlight its economic potential. However, despite ongoing research on the mineralogy and processing of rare-earth minerals, allanite has not received widespread attention in mineral processing. To achieve economical extraction of rare-earth elements from allanite in the future, systematic studies on processing techniques (e.g., density separation, magnetic separation, flotation, leaching) are imperative to fully unlock the potential of allanite as a rare-earth element source. To pave the way for future investigation of allanite and address the unique processing challenges, this review article aims to comprehensively summarize previous studies, encompassing properties, occurrences, and processing technologies of allanite.
Nonlinearity of Mohr envelopes confirmed by past laboratory test data can be of considerable importance in theoretical interpretation of rock strength behaviors and many practical applications. Current nonlinear failure criteria are introduced with parameters to empirically fit test data, and the physical and mechanical meanings of these fitting parameters are often not known. Predicted strengths by theoretical criteria (e.g., Griffith criterion) generally do not agree well with test data. Considering the nonlinear properties of a Mohr envelope, a nonlinear failure criterion in a generalized power-law form is derived in this paper based on strength magnitude coefficient (a) and failure envelope curvature coefficient (n) which are expressed in terms of uniaxial tensile strength, uniaxial compressive strength, and angle of internal friction. Test data on granite, marble, dolomite, and andesite from published literature were chosen to test the performance of the proposed criterion along with Griffith and Hoek-Brown criteria. Statistical analysis was performed to evaluate the application of different criteria in predicting the strength of rocks. In addition, some Griffith-typed criteria are mathematically shown to be the special cases of the proposed criterion. The proposed and Hoek-Brown criteria have similar power-law forms with a significant difference in the coefficient functions that govern the shear strength magnitude of rock. The different tensile and compressive failure mechanisms are also investigated from the perspective of the relationship between the uniaxial tension Mohr circle and the proposed failure envelope.