Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock deformation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Metasediments of the deep crustal section of Southern Karnataka

A review of the rocks of supracrustal origin in the amphibolite and granulite facies terrane in southern Karnataka was presented. In addition to introducing the metasediments in the field area of the workshop, a review was presented of the common occurrence of metasediments in amphibolite and granulite facies rocks worldwide. Models of granulite metamorphism must include a mechanism for the burial of these sediments to the depths recorded by the geobarometers in granulite metamorphism in addition to their reexposure at the surface. Unfortunately, the common occurrence of supracrustals in granulite facies rocks, sometimes with remarkably little deformation was deemed significant.

Devaraju, T. C.↗
NASA NTRS1986

Polyphase thrust tectonic in the Barberton greenstone belt

In the circa 3.5 by-old Barberton greenstone belt, the supracrustal rocks form a thick and strongly deformed thrust complex. Structural studies in the southern part of the belt have shown that 2 separate phases of over-thrusting (D sub 1 and D sub 2) successively dismembered the original stratigraphy. Thrust nappes were subsequently refolded during later deformations (D sub 3 and D sub 4). This report deals with the second thrusting event which, in the study region appears to be dominant, and (unlike the earlier thrusting), affects the entire supracrustal pile. The supracrustal rocks form a predominantly NE/SW oriented, SE dipping tectonic fan (the D sub 2 fan) in which tectonic slices of ophiolitic-like rocks are interleaved with younger sedimentary sequences of the Diepgezet and malalotcha groups. Structural and sedimentological data indicate that the D sub 2 tectonic fan was formed during a prolonged, multi-stage regional horizontal shortening event during which several types of internal deformation mechanisms were successively and/or simultaneously active. Movement appears to have been predominantly to the NW and to the N. During D sub 2, periods of quiescence and sedimentation followed periods of thrust propagation. Although the exact kinematics which led to the formation of this fan is not yet known, paleoenvironmental interpretations together with structural data suggest that D sub 2 was probably related to (an) Archean collision(s).

Paris, I. A.↗
DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗
NASA NTRS1986

Transpression as the main deformational event in an Archaean greenstone belt, northeastern Minnesota

Deformed and metamorphosed sedimentary and volcanic rocks of the Vermilion district constitute an Archean greenstone belt trending east-west between higher grade rocks of the Vermilion Granitic Complex to the north and the Giants Range batholith to the south. Metamorphic grade is low throughout, being lowest in the center of the belt (chlorite zone of the greenschist facies). All the measured strain, a cleavage or schistosity, and a mineral lineation in this belt are attributed to the main phase of deformation D sub 2 that followed an earlier nappe-forming event D sub 1, which left little evidence of penetrative fabric. Previous work assumed that the D sub 2 deformation resulted from north-south compression across the district. It is now believed that a significant component of this deformation resulted from dextral shear across the whole region. Thus the Vermilion fault, a late-state largely strike-slip structure that bounds the Vermilion district to the north, may simply be the latest, most brittle expression of a shear regime that was much more widespread in space and time. Features that are indicative of shear include ductile shear zones with sigmoidal foliation patterns, highly schistose zones with the development of shear bands, feldspar clasts or pyrite cubes with asymmetric pressure shadows, and the fact that the asymmetry of the F sub 2 folds is predominantly Z for at least 15 km south of the Vermilion fault.

Hudleston, P. J.↗
NASA NTRS1989

River Valley pluton, Ontario - A late-Archean/early-Proterozoic anorthositic intrusion in the Grenville Province

This paper presents Nd, Sr, and Pb isotopic data indicating a late-Archean/early-Proterozoic age for the River Valley anorthositic pluton of the southwestern Grenville Province of Sudbury, Ontario. Pb-Pb isotopic data on 10 whole-rock samples ranging in composition from anorthosite to gabbro yield an age of 2560 + or - 155 Ma. The River Valley pluton is thus the oldest anorthositic intrusive yet recognized within the Grenville Province. The Sm-Nd isotopic system records an age of 2377 + or - 68 Ma. High Pb-208/Pb-204 of deformed samples relative to igneous-textured rocks implies Th introduction and/or U loss during metamorphism in the River Valley area. Rb-Sr data from igneous-textured and deformed samples and from mineral separates give an age of 2185 + or - 105 Ma, indicating substantial disturbance of the Rb-Sr isotopic system.

Ashwal, Lewis D.↗
DOE OSTI2022

Creep fronts and complexity in laboratory earthquake sequences illuminate delayed earthquake triggering

Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters.

58 GEOSCIENCES↗
NASA NTRS1972

Shock metamorphic effects in the Luna 16 soil sample from Mare Fecunditatis.

Diverse lunar soil fragments returned from the Luna 16 mission were examined for evidence in shock metamorphism in order to evaluate the role of meteorite impact in forming the lunar regolith at a new site relatively distant from the Apollo landing sites. It was found that shock metamorphic effects, characteristic of meteorite impact and virtually identical to those observed in Apollo samples, are common in the Luna 16 fragments. Two types of shock effects are present - i.e., a deformation and partial melting features in rock and mineral fragments, and heterogeneous glasses and glassy breccias produced by shock melting. It is concluded that these shock metamorphic effects indicate that regolith formation by meteorite impact has occurred in Mare Fecunditatis and is a general process over the entire moon.

French, B. M.↗
NASA NTRS1988

Deformational models of rifting and folding on Venus

Features of presumed tectonic origin on Venus are reviewed, and lithospheric strength envelopes are derived based on laboratory measurements of the deformational properties of crustal and subcrustal rocks, extrapolated to conditions appropriate to Venus. Models for rifting and folding are developed that use this lithospheric structure and take into account both brittle and ductile yielding as well as finite elastic strength. For both rifting and folding, structures with characteristic widths and spacings are predicted whose size depends on the thickness of the lithosphere, density contrast, and elastic properties of the layer. Finally, the model predictions are compared with the widths and spacings of observed tectonic features, and it is concluded that they are consistent with a relatively strong mantle layer separated from a thin brittle surface layer by a ductile lower crust. These results allow constraints to be placed on the crustal thickness and thermal gradient on Venus.

Banerdt, W. Bruce↗
DOE OSTI2025

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗
DOE OSTI2023

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2020

Lattice modeling of excavation damage in argillaceous clay formations: Influence of deformation and strength anisotropy

This paper presents modeling of mechanical anisotropy in argillaceous rocks using an irregular lattice modeling approach, namely the rigid-body-spring network. To represent the mechanical anisotropy, new schemes are implemented in the modeling framework. The directionality of elastic deformation is resolved by modifying the element formulation with anisotropic elastic properties. The anisotropy of strength and failure characteristics is facilitated by adopting orientation-dependent failure criteria into the failure model. The verification of the improved modeling procedures is performed against theoretical model predictions for unconfined compression tests with various bedding orientations. Furthermore, excavation damage and fracturing processes in rock formations are simulated for different geomechanical configurations, such as rock anisotropy and tectonic heterogeneity. The simulated excavation damage characteristics are realistic and comparable with the actual field observation at a tunnel located in an argillaceous clay formation. The simulation results provide insights into the excavation damage zone phenomena with an explicit representation of fracturing processes.

58 GEOSCIENCES↗
DOE OSTI2022

Polytypism in semi-disordered lizardite and amesite by low-dose HAADF-STEM

Serpentine minerals exert important controls on the physical properties of ultramafic rocks and have the potential to influence deformation phenomena in fault zones and to control the release of water in subducted slabs. Sheet serpentine generally, and lizardite and amesite specifically, can adopt alternative crystallographic stacking arrangements called polytypes. Polytypism has been extensively studied in fully ordered crystals, but it remains largely enigmatic in the more common semi-disordered crystals that in long-range analyses such as X-ray diffraction only exhibit random combinations of 0b and ±1/3b interlayer shifts. To date, atomic-resolution imaging to identify locally ordered polytypes has been precluded by the beam-sensitive nature of this hydrous magnesium silicate mineral. Here, in this study, we employed low-dose high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) to study the polytypic structure of semi-disordered lizardite and amesite. Because the electron dose was as low as ~6000 e–/Å 2 , it was possible to directly resolve oxygen atomic columns and all the cations with a resolution of ~1 Å and reveal the short-range order. For lizardite, we identified long-period non-standard polytypes, including examples with 3, 4, 8, and 9 layers stemming from the ordering of the octahedral tilt along the a-axis. For amesite, we found short-range ordered polytypes with periodicities of up to 42 Å stemming from the ordering of interlayer shifts along the b-axis. The resolution was sufficient to determine the relative abundance of 6R 2 , 6R 1 , 2H 1 , and 2H 2 polytypes in amesite to be 46.1, 29.6, 7.7, and 1.9%, respectively. This is contrary to the expectation that the most common form of amesite is the 2H 2 polytype, which may be more likely to form macroscopic crystals suitable for conventional X-ray diffraction-based studies. We conclude that HAADF-STEM methods open the way for the characterization of beam-sensitive minerals and to resolve the structural details of less well-ordered (but possibly more abundant) minerals at a unit-cell scale.

58 GEOSCIENCES↗
NASA NTRS1978

Bedrock geology of the northern Columbia Plateau and adjacent areas

The Columbia Plateau is surrounded by a complex assemblage of highly deformed Precambrian to lower Tertiary continental and oceanic rocks that reflects numerous episodes of continental accretion. The plateau itself is comprised of the Columbia River basalt group formed between about 16.5 x 1 million years B.P. and 6 x 1 million years B.P. Eruptions were infrequent between about 14 and 6 x 1 million years B.P., allowing time for erosion and deformation between successive outpourings. The present-day courses of much of the Snake River, and parts of the Columbia River, across the plateau date from this time. Basalt produced during this waning activity is more heterogeneous chemically and isotopically than older flows, reflecting its prolonged period of volcanism.

Swanson, D. A.↗
DOE OSTI2021

Experimental workflow to estimate model parameters for evaluating long term viscoelastic response of CO2 storage caprocks

Understanding the time-dependent behavior of reservoir and sealing formations is critical to assessing risks associated with geological carbon storage since time-dependent deformation strongly influences mechanical responses of some rock types. Many studies have evaluated the risk of CO2 leakage and induced seismicity by assuming poroelastic rheology in sealing formations. Few have considered viscoelastic or other time-dependent responses, where the existing literature adopts 1D models to represent long-term time-dependent responses. This is primarily because to date, the general form of a reasonable 3D time-dependent model for rocks remains unclear. In this paper, we address this unclear issue by proposing a new workflow to select constitutive modeling parameters to evaluate if a 3D viscoelastic model is reasonable using several-hour-long experimental data and a power-law response to extrapolate to the decades-long time frames of interest in geologic carbon storage. To provide experimental data, we conducted multi-level loading/unloading triaxial relaxation tests with four rock types. The experimental results showed that the maximum load relaxation observed is approximately 49%, with some rock types showing as little as 1.4%. Using a simple linear viscoelastic model, parameters were chosen such that a maximum deviation of 1.5 MPa in axial stress and 7 MPa in radial stress was attained with the extrapolated 30-year data. We found that a reasonable parameter range for the normalized elastic modulus is 0.1~2 for rocks with significant time-dependent responses and 0.01~0.06 for those with small time-dependent responses. No matter how significant time-dependent responses are for rocks considered, our results showed that the relaxation time has a general range of 1~10^10 s, whose time scale can be one or two orders higher than a time frame typically envisioned for CO2 injection projects.

Stress relaxation, 3D time-dependent model, viscoe↗
DOE OSTI2026

Role of Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes for Enhanced Geothermal Systems

The ability to sustain high energy extraction efficiency from Enhanced Geothermal Systems (EGS) is affected by the ability to measure, characterize, and predict the effect of different perturbations that arise from fluid injection rates, fluid temperature, and shut-in conditions on existing and stimulated fracture systems throughout a subsurface reservoir’s lifecycle. The difficulty arises from limited knowledge of the interaction of fluid and temperature driven fractures with frictional interfaces that govern local deformation and frictional behavior in the surrounding rock. Thermo-poro-mechanical coupled processes tend to dominate this interaction and strongly affect the permeability and local stress distributions that impact efficiency and prevent optimal stimulation conditions over multi-decadal time frames.

Pyrak-Nolte, Laura [Purdue Univ., West Lafayette, ↗