Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress Modeling”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An Investigation of Chemo-Mechanical Phenomena and Li Metal Penetration in All-Solid-State Lithium Metal Batteries Using In Situ Optical Curvature Measurements

Solid-electrolytes (SEs) can provide a pathway to increase energy-density in lithium metal batteries. However, lithium metal penetration through garnet based LLZO solid electrolytes has been identified as a critical failure process. This phenomenon is related to chemo-mechanical processes which are difficult to probe. In particular, characterizing the dynamic mechanical deformations that occur in electrode-SE structures is very challenging. Here study reports in situ curvature measurements that are thus designed to probe chemo-mechanical phenomena that occur during lithium plating. The novel experimental cell configuration created for this work shows that pressure builds up in the Li metal during plating, up until the point where short circuits occur. The resulting data are analyzed with a detailed finite element model (FEM) to quantitatively evaluate stress evolution. The results show that Li metal plating within a surface flaw can produce stress build-up prior to short-circuiting. The combined results from both the experiments and the FEM suggest that it is critical to minimize surface defects and flaws during the manufacturing processes.

25 ENERGY STORAGE↗
DOE OSTI2020

Effect of heat treatment on the tensile behavior of selective laser melted Ti-6Al-4V by in situ X-ray characterization

Selective Laser Melted Ti-6Al-4V (as-SLMed) exhibits decreased yield strength, increased work hardening, and increased ductility after heat treatment at 730 degrees C (HT-730) or 900 degrees C (HT-900) for 2 h. To understand the change of mechanical properties, in situ high energy X-ray diffraction (HEXRD) is used to examine the phase composition, load partitioning, slip system activity, and dislocation density evolution in all three specimens. The as-SLMed specimen is dominated by martensitic α'. After heat treatment, α' partly or fully decomposes into α+β, reducing the yield strength. In HT-730, beta precipitates with confined size show much higher lattice strain than the α'/α matrix during deformation; in HT-900, the lattice strain difference is mostly eliminated. This is a key reason for the increased ductility in HT-900. From the anisotropic lattice strain development, basal slip is identified as the easiest slip system in α'/α. Using an elasto-plastic self-consistent (EPSC) model, the critical resolved shear stress ratio between prismatic slip and basal slip (CRSS prism /CRSS basal ) is estimated to be 1.31 and 1.16 in the as -SLMed and the HT-900 specimens, respectively. Further, the alpha phase in HT-900 is able to activate multiple slip systems and accumulate more dislocations during plastic deformation. This explains why HT-900 has better ductility and higher work hardening rate than the other two specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ultrastructure of the green alga Dunaliella salina strain CCAP19/18 (Chlorophyta) as investigated by quick-freeze deep-etch electron microscopy

The single-celled green alga Dunaliella salina is a model system for studies on stress biology, in particular regarding secondary carotenoid accumulation. Under non-stress conditions the cells are green, but under abiotic stress the cells turn orange, because they switch their metabolism and accumulate β-Carotene in globules in the chloroplast. For the first time, Quick-freeze deep-etch electron microscopy was used to visualize cellular structures in green and orange cells of D. salina strain CCAP19/18. This allowed us to present an in-depth analysis of the cellular ultrastructure describing and comparing the features of the two cell types. Our images illustrate the presence of a pericellular matrix for this strain of D. salina. The pericellular matrix was spongy and strands of unknown material anchored it into the plasma membrane. The cytoplasm contained a variety of vesicles, vacuoles, and acidocalcisomes. We could show for strain CCAP19/18 that cytoplasmic lipid bodies were often in close proximity to and sometimes in contact with the outer chloroplast envelope membrane and with the endoplasmic reticulum. Major visible differences between green and orange cells were in the chloroplast: the orange cells have greatly reduced amounts of thylakoid membranes and greatly increased numbers of β-Carotene globules. We showed that the β-Carotene globules often made point contacts with thylakoid membranes, and frequently laid side-by-side along the thylakoid membrane surface, providing support to studies that indicated exchange of molecules between β-Carotene globules and thylakoid membranes. A unique finding was the β-Carotene globule duplets, suggesting intermediate stages in β-Carotene globule morphogenesis. Overall, the β-Carotene globules appear to be similar to plant plastoglobuli. We provide this description of the cellular ultrastructure features as a resource in the context of the recent publication of the genome of D. salina strain CCAP19/18 to expand on the knowledge regarding this novel reference strain.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Micromechanical response of SiC-OPyC layers in TRISO fuel particles

Tristructural isotropic (TRISO)–coated particle fuel is a proposed fuel for multiple advanced reactor concepts. The performance of the particle depends on whether the silicon carbide (SiC) layer remains intact to prevent the release of metallic and gaseous fission products. Mechanical fracture of the SiC layer is a potential failure mode under various fuel configurations and operating environments, including the potential transmission of matrix-originating cracks through TRISO particles. Furthermore, this study uses instrumented indentation techniques on cross-sectioned surrogate particles to examine the mechanical stability of the critical interface between SiC and the outer pyrolytic carbon (OPyC) layer. The observed behavior at the interface is rationalized by examining the radially dependent fracture behavior of the SiC layer and performing a numerical analysis to quantify the residual stresses that develop during the processing and cross-sectioning of the as-fabricated particle. Characterizing the SiC-OPyC interface of surrogate TRISO particles using nanoindentation provides unique insight into the interface's room-temperature residual stress and mechanical stability. The modeling efforts were used to investigate the experimental procedure further, and the results are presented herein to validate this fuel form's potential mechanical failure modes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The formation and evolution of Ni 2 Cr precipitates in Ni–Cr model alloys as a function of stoichiometry characterized by synchrotron x-ray diffraction

Ni-based alloys containing significant amounts of Cr, may precipitate new phases during long-term service at elevated temperatures. In Ni–Cr model alloys, the formation of Ni 2 Cr has been found to impact the material properties, including the strength and ductility after isothermal aging below the critical temperature of 590°C. In this work, we quantify the formation and evolution of long-range ordering in four Ni–Cr binary model alloys (Ni/Cr = 1.8, 2.0, 2.2, 2.4) after isothermal aging up to 10,000h at temperatures between 373 and 475°C. The alloys were characterized by hardness testing and synchrotron-based x-ray diffraction to quantify the impact of Ni 2 Cr phase fraction and precipitate size on mechanical properties. After 500h of isothermal aging, the formation of Ni 2 Cr was detected in all four alloys at 475°C. In the stochiometric alloy samples (Ni/Cr = 2.0), the formation of Ni 2 Cr was found after 500 and 3000h aging at 418 and 373°C, respectively. We found that the matrix lattice contraction and Ni 2 Cr phase fraction both saturate after early aging times. This is in stark contrast to the hardness and Ni 2 Cr precipitate size that both continue to increase with increasing aging time. Our results highlight that changes in hardness correlate linearly with Ni 2 Cr precipitate size rather than phase fraction. This important structure-property relationship can potentially help define Ni–Cr-based component lifetimes directly through an understanding of how Ni 2 Cr formation impacts strength and ductility. Here, we find that a precipitation hardening model for critical resolved shear stress with weakly coupled dislocations shows good agreement with the material property changes quantified from experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio calculations of third-order elastic coefficients

Third-order elasticity (TOE) theory predicts strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an “effective” third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. Additionally, we extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Predicting Rupture Pressure and Leak Rate of Alloy 600 Tubes with ODSCC and Comparing with Test Results (Technical Letter Report)

This report summarizes structural integrity work performed for the U.S. NRC that predicts the rupture pressure and leak rate of flawed steam generator (SG) tubes. The work includes pressure testing of straight tubes with lab-induced outer diameter stress corrosion cracks (ODSCC). Existing models for predicting rupture pressure and leak rate of straight tubes with well-defined, rectangular, electro-discharge-machined (EDM) flaws were modified and used to predict the measured rupture pressures and leak rates from straight tubes with ODSCC. Models previously developed for EDM notches in straight tubes were modified for ODSCC by using the equivalent rectangular crack (ERC) method.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

MRCI Task 2: Addressing Key Technical Challenges - Executive Summary Report

This report provides a high-level overview of Task 2 products that were developed through defining carbon storage systems, Precambrian basement structure and stress, developing regional technical collaboration, modeling and risk assessment. Subtask 2.1: Defining Sub-Regional CS/CCUS Systems Subtask 2.2: Defining Precambrian basement faulting/stress Subtask 2.3: Developing industrial partnerships and regional technical collaborations Subtask 2.4: Conducting regional/subregional analysis Subtask 2.5: Assessing and managing risk for potential commercial-scale storage complexes

MRCI,Technical Challenges↗
DOE OSTI2025

Virtual Reality for Shoot/No-Shoot Decision Training in Law Enforcement: A Literature Review and Research Agenda

Virtual reality (VR) can materially improve “shoot / no-shoot” (SNS) training by giving officers realistic, repeatable practice making high-stakes decisions under pressure. Traditional tools—live-fire ranges and video simulators—build basics, but they cannot adapt to each officer in real time or fully mirror the complexity of the field. VR closes that gap by creating immersive scenarios that are safer, more flexible, easier to scale across units, and able to capture objective performance data. SNS decisions are not just about marksmanship; they rely on perception, judgment, memory, and the ability to hold fire when a threat is uncertain. Effective training therefore needs realism, decision complexity, and branching outcomes that reflect the true consequences of choices. These elements strengthen recognition of hostile intent while reducing false positives and building the self-control required in ambiguous situations. VR brings specific advantages: dynamic environments, full-body interaction, and the ability to measure performance with precision—enabling targeted feedback and better transfer of learning to the street. At the same time, responsible deployment must address scenario quality (credible environments and behaviors), lawful decision models, and user wellbeing (appropriate stress levels, comfort, and safety). Sandia’s VIPER Lab is positioned to lead this work. The team combines human-performance science, AI/ML, and VR/AR development with a deep equipment bench (e.g., omnidirectional treadmill, eye-tracking, haptics, multiple HMDs). This ecosystem supports building and validating next-generation SNS training that is immersive, measurable, and trustworthy. Bottom line: Investment in VR-enabled SNS training that blends evidence-based design with careful validation and legal safeguards is expected to pay off in safer, more consistent decision-making and improved community trust, delivered through training that is practical to deploy at scale.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
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 OSTI2022

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Burst pressure models and validations for thick-walled pipelines containing corrosion defects

Corrosion is one major threat to pipeline integrity. Over the past decades, many corrosion models have been developed for determining the remaining strength of corroded pipelines, including ASME B31.G, Modified B31.G, LPC, PCORRC and their modified models. All these corrosion models are applicable only to large diameter, thin-walled pipelines with a diameter to wall thickness ratio D/t ≥ 20. In practice, many pipelines have a small diameter and thick wall with a D/t ratio < 20, and thus an adequate corrosion model is needed for assessing remaining strength for corroded thick-walled pipelines. This paper briefly reviews the theoretical burst pressure models for defect-free thin and thick-walled pipelines and four representative corrosion assessment models for thin-walled corroded pipelines. On this basis, two modified corrosion models are proposed to thick-walled pipelines in terms of the average shear stress yield theory. To verify the proposed corrosion models, comprehensive validations are performed. Numerical validations include the elastic-plastic finite element analysis to determine burst pressure for pipelines without and with corrosion defects and the model evaluation using a large dataset of available FEA results of burst pressure for machined defects. Experimental validations include a set of burst pressure tests for defect-free thick-walled pipes with different thicknesses and the model evaluation using one large burst dataset for machined defects with flat bottoms and another large dataset for real corrosion defects with curved river bottom profiles. Both numerical and experimental validations show that the proposed corrosion models can more accurately predict the remaining strength for corroded thin and thick-walled pipelines.

Pipeline↗
DOE OSTI2023

Crack-Free Joining of High-Strength AA7055 Sheets by Friction Based Self-Piercing Riveting with the Aid of Numerical Design

Unique friction-based self-piercing riveting (F-SPR) was employed to join high-strength, low-ductility aluminum alloy 7055 for lightweight vehicle applications. This study aimed to maximize the joint strength of the AA7055 F-SPR joint while avoiding cracking issues due to low ductility at room temperature. A fully coupled Eulerian–Lagrangian (CEL) model was employed to predict the process temperature during F-SPR, and the temperature field was then mapped onto a 2D axisymmetric equivalent model for accelerated numerical analysis. The geometry, dimensions, and material strength of the rivet, as well as the depth of the die cavity and plunging depth, were investigated to enhance joint formation. Also, a static finite-element analysis model was developed to predict and analyze the stress distribution in the rivet under different mechanical testing loading conditions. Overall, the numerical model showed good agreement with the experiment results, such as joint formation and mechanical joint strength. With the aid of virtual fabrication through numerical modeling, the joint design iterations and process development time of F-SPR were greatly reduced regarding the goal of lightweight, high-strength aluminum joining.

36 MATERIALS SCIENCE↗
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 OSTI2022

Modeling Validation of Mechanical and Thermal Borehole Breakouts From Polyaxial Laboratory Tests to Determine Maximum Horizontal in Situ Stress

To accurately assess rock behavior in the deep subsurface, it is necessary to measure the in situ stress directions and magnitudes. The current methods for measuring the in situ stress state in the deep subsurface primarily include hydraulic fracturing tests (i.e., the minimum horizontal stress) or occasionally observing existing compressive borehole breakouts (i.e., the maximum horizontal stress) that have occurred naturally from drilling. If there are no existing compressive breakouts in a borehole, the maximum horizontal in situ stress cannot be estimated with much confidence. In response to this data gap, a new thermal breakout technology is being developed that will provide a method for thermally inducing borehole breakouts and obtaining consistent measurements of the maximum horizontal stress magnitude. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which is directly correlated to the maximum horizontal stress magnitude. In support of developing the thermal breakout technology, polyaxial laboratory tests have been performed on small-scale boreholes within rock blocks where mechanically- and thermally-induced borehole breakouts have been created. Numerical models along with the principal of superposition were created and used to analyze the polyaxial laboratory tests to predict the maximum horizontal stress, given the same data that would be obtained in an actual subsurface borehole thermal breakout test. Multiple failure criteria were used to evaluate the best prediction of the breakout onset and maximum horizontal stress. The maximum horizontal stress predictions were compared to the actual maximum horizontal stress applied at breakout using acoustic events recorded from emission sensors in the polyaxial tests. The results showed consistent results that could be used to refine the modeling approach and failure criterion that are used to make the maximum horizontal stress predictions. This study provided insight and validation for the thermal breakout stress measurement concept.

Jones, Matt↗
DOE OSTI2025

Importance of Considering Near-Surface Attenuation in Earthquake Source Parameter Estimation: Insights from Kappa at a Dense Array in Oklahoma

ABSTRACT Separating earthquake source spectra from propagation effects is challenging. The propagation effect contains a site-dependent term related to the high attenuation of shallow sediments. Neglecting the site-dependent attenuation can cause large biases and scattering in the corner-frequency (fc) estimates, resulting in significant stress-drop deviations. In this study, we investigate shallow attenuation at the LArge-n Seismic Survey in Oklahoma (LASSO) and site-related biases and scattering in source parameter measurements due to simplified attenuation models. We measure the high-frequency spectral decay parameter kappa on the vertical acceleration spectra of regional earthquakes (125 km away). The site-dependent kappa (κ0,acc) suggests that attenuation increases rapidly at shallow depth and is highly site-dependent. 10%–75% of the attenuation is site-dependent for S waves and even larger for P waves. The quality factor for S waves (QS) ranges from 10 to 100 in the upper 400 m. QP for P waves is mostly below 10 within the same depth. The Quaternary sediments tend to be more attenuating (QS<30), but the Permian rocks also can have high attenuation. We demonstrate that using a non-site-dependent attenuation model in single-spectra fitting leads to large scattering in fc estimates among stations with apparent good fits. The apparent fc can significantly deviate when the range of site-dependent kappa is large or with a higher assumed source spectral fall-off rate n. The biases in apparent fc depend on site condition and distance; however, the correlation between fc and these factors might not be obvious, depending on model assumptions. An apparent increase of stress drop with magnitude in a previous study for local microearthquakes (1.3

Chang, Hilary↗