Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deformation estimation”

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 55 records · Page 3

Search for Infrasound Signals in InSight Data Using Coupled Pressure/Ground Deformation Methods

The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.

58 GEOSCIENCES↗

Overcoming Anode Instability in Solid‐State Batteries through Control of the Lithium Metal Microstructure

Abstract Enabling the lithium metal anode (LMA) in solid‐state batteries (SSBs) is the key to developing high energy density battery technologies. However, maintaining a stable electrode–electrolyte interface presents a critical challenge to high cycling rate and prolonged cycle life. One such issue is the interfacial pore formation in LMA during stripping. To overcome this, either higher stack pressure or binary lithium alloy anodes are used. Herein, it is shown that fine‐grained ( d = 20 µm) polycrystalline LMA can avoid pore formation by exploiting the microstructural dependence of the creep rates. In a symmetric cell set‐up, i.e., LiǀLi 6.25 Al 0.25 La 3 Zr 2 O 12 (LLZO)ǀLi, fine‐grained LMA achieves > 11.0 mAh cm −2 compared to ≈ 3.6 mAh cm −2 for coarse‐grained LMA ( d = 295 µm) at 0.1 mA cm −2 and at moderate stress of 2.0 MPa. Smaller diffusion lengths (≈ 20 µm) and higher diffusivity pathway along dislocations ( D d ≈ 10 −7 cm 2 s −1 ), generated during cell fabrication, result in enhanced viscoplastic deformation in fine‐grained polycrystalline LMA. The electrochemical performances corroborate well with estimated creep rates. Thus, microstructural control of LMA can significantly reduce the required stack pressure during stripping. These results are particularly relevant for “anode‐free” SSBs wherein both the microstructure and the mechanical state of the lithium are critical parameters.

25 ENERGY STORAGE↗

Comprehensive Review of Multi-arm Caliper Data for the Big Hill SPR Site

The Big Hill SPR site has a rich data set consisting of multi-arm caliper (MAC) logs collected from the cavern wells. This data set provides insight into the on-going casing deformation at the Big Hill site. This report summarizes the MAC surveys for each well and presents well longevity estimates where possible. Included in the report is an examination of the well twins for each cavern and a discussion on what may or may not be responsible for the different levels of deformation between some of the well twins. The report also takes a systematic view of the MAC data presenting spatial patterns of casing deformation and deformation orientation in an effort to better understand the underlying causes. The conclusions present a hypothesis suggesting the small-scale variations in casing deformation are attributable to similar scale variations in the character of the salt-caprock interface. These variations do not appear directly related to shear zones or faults.

58 GEOSCIENCES↗

2024 Update of Comprehensive Review of Multi-arm Caliper Data for the Big Hill SPR Site

The Big Hill SPR site has a rich data set consisting of multi-arm caliper (MAC) logs collected from the cavern wells. This data set provides insight into the on-going casing deformation at the Big Hill site. This report summarizes the MAC surveys for each well and presents well longevity estimates where possible. Included in the report is an examination of the well twins for each cavern and a discussion on what may or may not be responsible for the different levels of deformation between some of the well twins. The report also takes a systematic view of the MAC data presenting spatial patterns of casing deformation and deformation orientation in an effort to better understand the underlying causes. The conclusions present a hypothesis suggesting the small-scale variations in casing deformation are attributable to similar scale variations in the character of the salt-caprock interface. These variations do not appear directly related to shear zones or faults. In addition, the deformation orientation shows no preferred directionality. This 2024 edition of this report represents an update to the original, 2023 edition. The updates primarily focus on the inclusion of MAC log data run since the December 2021 threshold date for the original report, but some new analyses are also included.

58 GEOSCIENCES↗

Tensile and compressive plastic deformation behavior of medium-entropy Cr-Co-Ni single crystals from cryogenic to elevated temperatures

The equiatomic Cr-Co-Ni medium-entropy alloy has the face-centered cubic (FCC) structure. Bulk single crystals of this alloy were grown and tested in tension and compression between 14 K and 1373 K with the loading axis parallel to [23]. At room temperature, the critical resolved shear stress (CRSS) for {111}<10> slip is 65 ± 5 MPa and does not exhibit a tension-compression asymmetry. It does, however, increase significantly as the test temperature decreases. A dulling of this temperature dependence occurs below 50 K, which may be due to the inertia effect. When the measured values above 50 K are extrapolated to lower temperatures, a value of 225 MPa is estimated for the CRSS at 0 K. This is larger than that (168 MPa) previously determined for the equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy using a similar procedure. The stacking fault energy of the present Cr-Co-Ni is estimated to be about 14 mJm-2, which is sufficiently low to account for deformation twinning both at 77 K and room temperature. Twinning at 77 K occurs on conjugate (1) planes at an onset shear stress of 482 MPa after primary slip and propagates in the form of Lüders deformation. At room temperature, twinning occurs uniformly throughout the gauge section on primary (111) planes at an onset shear stress of 381 MPa after primary and subsequent conjugate slip. Thin layers with the hexagonal close-packed stacking are observed in association with twinning both at 77 K and room temperature.

36 MATERIALS SCIENCE↗

Development of a Stress Sensor for In-Situ High-Pressure Deformation Experiments Using Radial X-Ray Diffraction

We developed a stress sensor for in-situ deformation experiments using synchrotron radial X-ray diffraction. This stress sensor provided nearly diffraction-plane-independent stress that, when used in series with a sample, reduced the uncertainty of the average stress estimation acting on a sample. Here, we present the results of a study where pyrope was used as a stress sensor. Using a Deformation-DIA (D-DIA) high-pressure deformation apparatus, pyrope, olivine and alumina were deformed in the same run/cell assembly placed in series along the compression direction. Deformation experiments were conducted at pressures between 4 and 5 GPa and temperatures between 730 and 1273 K with strain-rates between 10 –5 and 10 –6 s –1 . Stresses estimated from various (hkl) planes in pyrope were nearly the same; i.e., pyrope is plastically isotropic with ≤10 % variation with (hkl). However, stresses from various (hkl) planes in olivine and alumina varied by approximately a factor of 3. Comparisons between average stresses inferred from pyrope and those from different diffraction planes in olivine and alumina showed that the average stress in these materials evolved from low-end stress, estimated from various (hkl) planes at small strain, to high-end stress at a large strain. This suggests that the rate-controlling slip system in these materials changes from the soft to the hard slip system with strain.

58 GEOSCIENCES↗

Microstructures and deformation mechanisms of the medium-entropy alloy (NiCoCr) 76 (Ni 6 AlTi) 3

Here the microstructures and deformation mechanisms of the f.c.c. medium-entropy alloy (MEA) (NiCoCr) 76 (Ni 6 AlTi) 3 has been analyzed after various thermo-mechanical treatments. The solutionized, single-phase MEA, which had a grain size of 93.4 ± 31.9 μm, was cold-rolled (CR) to a 80% thickness reduction after which it showed both a high yield strength (YS) of 1539 MPa and a high ultimate tensile strength (UTS) of 1602 MPa, but an elongation to failure, ε, of only 16%. The CR MEA was then subjected to one of two heat treatments: (1) the CR MEA was recrystallized at 1100°C for 24 h, which produced a single-phase material with 85 ± 61.7 μm grain size that exhibited a much lower YS and UTS of 364 MPa and 747 MPa, respectively, but a much greater ε of 73%; and (2) the MEA was aged at 700°C for 24 h, which produced a fine-grained (1.1 ± 0.9 μm) material and a high volume fraction (0.35) of 12 nm dia. L1 2 nanoparticles that exhibited an excellent combination of strength and ductility, viz., YS~1501 MPa, UTS~1651 MPa, elongation ~26%. The grain boundary strengthening and precipitation strengthening were together estimated to account for ~58% of the YS in the latter material. After deformation, the recrystallized MEA contained numerous stacking faults and a Taylor lattice structure containing domain boundaries and microbands, while the aged MEA exhibited numerous stacking faults and sheared particles. The density of low-angle grain boundaries (LAGBs) in the recrystallized MEA increased by 100 times to 2.2 × 10 –1 μm/μm 2 while the density of CSL Σ3 n boundaries decreased by 71% to 7.8 × 10 –3 μm/μm 2 after deformation. In contrast, the density of LAGBs in the aged MEA increased by 5 times to 2.6 μm/μm 2 while the density of CSL Σ3 n boundaries decreased by 57% to 0.6 μm/μm 2 after deformation.

36 MATERIALS SCIENCE↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗

Geophysical Signatures of Crack Network Coalescence in Rocks at Multiple Length Scales

The main goal of the research project was to identify the geophysical signatures of fracture growth in natural rocks by utilizing novel geophysical techniques. The research objectives were to (a) investigate the potential for geophysical methods to determine when cracks initiate, the types and locations of propagated cracks, and the coalescence of networks of cracks in natural rocks at multiple scales, (b) determine how damage at the microscale evolved into damage at the macroscale and then link the microscopic and macroscopic observations, (c) quantify crack coalescence in rocks under realistic stress conditions using coupled mechanical-geophysical-optical visualization, and (d) identify the precursors in geophysical signals to crack coalescence. The following research thrusts were explored to achieve the research objectives: (1) uniaxial compression testing of rock specimens with and without a set of pre‐existing flaws and (2) triaxial compression testing of natural rock specimens. These thrusts allowed for exploring fracturing in rocks under realistic in situ environments and at multiple scales. This project provided educational opportunities for nine graduate and undergraduate students and resulted in 27 peer-reviewed publications. This first research thrust focused on investigating the micromechanics of fractures in rocks through uniaxial compression testing combined with advanced geophysical and imaging techniques, specifically acoustic emission (AE) monitoring, ultrasonic imaging, and 2-dimensional Digital Image Correlation (2D-DIC). By examining damage processes under time-independent and time-dependent loading conditions, insights into damage localization, crack initiation, and fracturing mechanisms were gained. It was observed that the AE signals and the strain-based measurements directly reflect the state of damage in the rock specimen and could be used to identify the cracking levels, such as the crack initiation (CI) and crack damage (CD), and the mode of deformation. A novel calibration apparatus was developed to enhance the accuracy of AE sensors, allowing for the estimation of key parameters such as magnitude, source dimension, stress drop, and radiated seismic energy associated with the fractures. The findings highlighted significant variations in the temporal evolution of AE source parameters during the primary, secondary, and tertiary stages of creep, identifying tensile cracking as the primary deformation mode. The second research thrust focused on enhancing the understanding of fracturing processes in natural rocks through triaxial compression testing, real-time AE monitoring, and ultrasonic monitoring. We investigated the impact of various factors such as fracture propagation regimes, injection parameters, rock types, and pre-existing conditions on the hydraulic fracture (HF) behavior using scaled true-triaxially loaded specimens of Barre granite and Lyons sandstone. Custom sensor housing facilitated concurrent active and passive monitoring to analyze hydro-mechanical responses and microseismicity associated with different HF scenarios. A coupled investigation of passive microseismicity and active signal attributes permitted a detailed comprehension of the various HF processes (aseismic deformation, fracture initiation and propagation, fluid permeation, and leak-off) and their dependence on the specific rock type. The findings of this research demonstrated the effectiveness of AE monitoring techniques in providing valuable insights into the impact of various factors on the behavior and dynamics of HF processes. The advancements in monitoring techniques, offering a more thorough and precise approach, represent a significant step towards optimizing HF practices and ensuring sustainable resource extraction.

58 GEOSCIENCES↗

Ac losses in field-cooled type I superconducting cavities

As superconductors are cooled below their critical temperature, stray magnetic flux can become trapped in regions that remain normal. The presence of trapped flux facilitates dissipation of ac current in a superconductor, leading to losses in superconducting elements of microwave devices. In type II superconductors, dissipation is well-understood in terms of the dynamics of vortices hosting a single flux quantum. In contrast, the ac response of type I superconductors with trapped flux has not received much attention. Building on Andreev's early work (Andreev 1967 Sov. Phys. JETP 24 1019), here we show theoretically that the dominant dissipation mechanism is the absorption of the ac field at the exposed surfaces of the normal regions, while the deformation of the superconducting/normal interfaces is unimportant. We use the developed theory to estimate the degradation of the quality factors in field-cooled cavities, and we satisfactorily compare these theoretical estimates to the measured field dependence of the quality factors of two aluminum cavities. We also identify a regime in which the dissipated power depends weakly on the Ginzburg-Landau parameter; this makes it possible to apply our findings to cavities made of other materials, such as niobium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗

Bayesian analysis for estimating statistical parameter distributions of elasto-viscoplastic material models

High temperature design methods rely on constitutive models for inelastic deformation and failure typically calibrated against the mean of experimental data without considering the associated scatter. Variability may arise from the experimental data acquisition process, from heat-to-heat material property variations, or both and need to be accurately captured to predict parameter bounds leading to efficient component design. Applying the Bayesian Markov Chain Monte Carlo (MCMC) method to produce statistical models capturing the underlying uncertainty in the experimental data is an area of ongoing research interest. This work varies aspects of the Bayesian MCMC method and explores their effect on the posterior parameter distributions for a uniaxial elasto-viscoplastic damage model using synthetically generated reference data. From our analysis with the uniaxial inelastic model we determine that an informed prior distribution including different types of test conditions results in more accurate posterior parameter distributions. The parameter posterior distributions, however, do not improve when increasing the number of similar experimental data. Additionally, changing the amount of scatter in the data affects the quality of the posterior distributions, especially for the less sensitive model parameters. Moreover, we perform a sensitivity study of the model parameters against the likelihood function prior to the Bayesian analysis. The results of the sensitivity analysis help to determine the reliability of the posterior distributions and reduce the dimensionality of the problem by fixing the insensitive parameters. The comprehensive study described in this work demonstrates how to efficiently apply the Bayesian MCMC methodology to capture parameter uncertainties in high temperature inelastic material models. Quantifying these uncertainties in inelastic models will improve high temperature engineering design practices and lead to safer, more effective component designs.

42 ENGINEERING↗

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE↗

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

Experimental Evaluation of Deformation and Fracture Mechanisms in Highly Irradiated Austenitic Steels

The present report documents recent experimental results of analysis using scanning electron microscopy/electron backscatter diffraction (SEM-EBSD) of plastic deformation mechanisms and strain localization phenomena in austenitic steels irradiated by neutrons. Experiments were performed with specimens irradiated to 125 dpa and, additionally, with specimens that experienced radiation-induced swelling up to 3%. Section 1 briefly analyzes the deformation localization in irradiated steels and its consequences on the material performance. The section describes the advantages and importance of the SEM-EBSD approach combined with in situ mechanical testing capability. Section 2 briefly introduces the experimental tools and methods (i.e., SEM/EBSD in situ tensile frame, electric discharge machine to manufacture irradiated specimens) and describes the investigated materials (i.e., element composition, irradiation conditions, and general microstructure). Section 3 describes the key experimental results and provides a brief analysis and comparison with the datasets obtained earlier within the same task (i.e., low-dose specimens). The discussion focuses on EBSD microstructure maps with strain localization features, misorientation evolution as a function of strain, and observed deformation mechanisms. Section 4 evaluates data collected in recent years on highly irradiated steel and estimates the possible misorientation evolution under irradiation. The section introduces and discusses the concept of in-service-induced damage as an irradiation-assisted stress-corrosion cracking precursor. Section 5 summarizes the work performed. As expected, the present work results are beneficial for exploring and understanding degradation mechanisms in highly irradiated in-core materials found in light water reactors after long-term in-service life.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Isotropic seismic radiation from rock damage and dilatancy

SUMMARY We present semi-analytical and numerical results on isotropic radiation expected to be generated by seismic sources with damaged dilatant materials that sustain during failure an overall shear deformation. The results are derived for inclusions with different stiffness components, including those responsible for shear-volumetric coupling in dilatant material, and different geometries, in a surrounding elastic solid. The obtained relations indicate that the isotropic seismic potency and seismic moment components are proportional to the shear potency components multiplied by the aspect ratio between the thickness and length of the source zone. For a typical rupture zone with a thin tabular geometry, the expected small aspect ratio explains the general seismological observations of small isotropic radiation. However, isotropic radiation induced by dynamic dilation for short and thick sources is more pronounced and may reach for spherical sources sustaining shear deformation 18 per cent of the shear components. Sources that include also tensile components and explosions can produce larger isotropic radiation. The results provide rigorous relations that can be used to estimate different source properties in numerical simulations, experiments and field studies.

Geochemistry & Geophysics↗