Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initial state”

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 847 records · Page 47

Domain Compilation for Embedded Real-Time Planning

A recently conceived approach to automated real-time control of the actions of a robotic system enables an embedded real-time planning algorithm to develop plans that are more robust than they would otherwise be, without imposing an excessive computational burden. This approach occupies a middle ground between two prior approaches known in the art as the universal-plan and hybrid approaches. Ever since discovering the performance limitations of taking a sense-plan-act approach to controlling robots, the robotics community has endeavored to follow a behavior-based approach in which a behavior includes a rapid feedback loop between state estimation and motor control. Heretofore, system architectures following this approach have been based, variously, on algorithms that implement universal plans or algorithms that function as hybrids of planners and executives. In a typical universal-plan case, a set of behaviors is merged into the plan, but the system must be restricted to relatively small problem domains to avoid having to reason about too many states and represent them in the plan. In the hybrid approach, one implements actions as small sets of behaviors, each applicable to a limited set of circumstances. Each action is intended to bring the system to a subgoal state. A planning algorithm is used to string these actions together into a sequence to traverse the state space from an initial or current state to a goal state. The hybrid approach works well in a static environment, but it is inherently brittle in a dynamic environment because a failure can occur when the environment strays beyond the region of applicability of the current activity. In the present approach, a system can vary from the hybrid approach to the universal-plan approach, depending on a single integer parameter, denoted n, which can range from 1 to a maximum domain-dependent value of M. As illustrated in the figure, n = 1 represents the hybrid approach, in which each linked action covers a small part of the state space of the system. As n increases, the portion of state space associated with each action and its subgoal grows. When n reaches M, coverage extends over the full state space, so that the system contains a universal plan.

Barrett, Anthony↗

Solute Nucleation and Growth in Supercritical Fluid Mixtures

This research effort is directed toward two primary scientific objectives: (1) to determine the gravitational effect on the measurement of nucleation and growth rates near a critical point and (2) to investigate the nucleation process in supercritical fluids to aid in the evaluation and development of existing theoretical models and practical applications. A nucleation pulse method will be employed for this investigation using a rapid expansion to a supersaturated state that is maintained for approximately 1 ms followed by a rapid recompression to a less supersaturated state that effectively terminates nucleation while permitting growth to continue. Nucleation, which occurs during the initial supersaturated state, is decoupled from growth by producing rapid pressure changes. Thermodynamic analysis, condensation modeling, apparatus design, and optical diagnostic design necessary for the initiation of a theoretical and experimental investigation of naphthalene nucleation from supercritical CO2 have been completed.

Smedley, Gregory T.↗

High Energy Density Shape Memory Polymers Using Strain-Induced Supramolecular Nanostructures

Shape memory polymers are promising materials in many emerging applications due to their large extensibility and excellent shape recovery. However, practical application of these polymers is limited by their poor energy densities (up to ~1 MJ/m 3 ). Here, we report an approach to achieve a high energy density, one-way shape memory polymer based on the formation of strain-induced supramolecular nanostructures. As polymer chains align during strain, strong directional dynamic bonds form, creating stable supramolecular nanostructures and trapping stretched chains in a highly elongated state. Upon heating, the dynamic bonds break, and stretched chains contract to their initial disordered state. This mechanism stores large amounts of entropic energy (as high as 19.6 MJ/m 3 or 17.9 J/g), almost six times higher than the best previously reported shape memory polymers while maintaining near 100% shape recovery and fixity. The reported phenomenon of strain-induced supramolecular structures offers a new approach toward achieving high energy density shape memory polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Motif magnetism and quantum many-body scars

In this work, we generally expect quantum systems to thermalize and satisfy the eigenstate thermalization hypothesis (ETH), which states that finite-energy-density eigenstates are thermal. However, some systems, such as many-body localized systems and systems with quantum many-body scars, violate ETH and have high-energy athermal eigenstates. In systems with scars, most eigenstates thermalize, but a few atypical scar states do not. Scar states can give rise to a periodic revival when time-evolving particular initial product states, which can be detected experimentally. Recently, a family of spin Hamiltonians was found with magnetically ordered three-colored eigenstates that are quantum many-body scars [Lee et al., Phys. Rev. B 101, 241111(R) (2020)]. These models can be realized in any lattice that can be tiled by triangles, such as the triangular or kagome lattices, and have been shown to have close connections to the physics of quantum spin liquids in the Heisenberg kagome antiferromagnet. In this paper, we introduce a generalized family of n-colored Hamiltonians with “spiral colored” eigenstates made from n-spin motifs such as polygons or polyhedra. We show how these models can be realized in many different lattice geometries and provide numerical evidence that they can exhibit quantum many-body scars with periodic revivals that can be observed by time-evolving simple product states. The simple structure of these Hamiltonians makes them promising candidates for future experimental studies of quantum many-body scars.

36 MATERIALS SCIENCE↗

Phase-Shifting Liquid Crystal Point-Diffraction Interferometry

Microgravity fluid physics experiments frequently measure concentration and temperature. Interferometers such as the Twyman Green illustrated have performed full-field measurement of these quantities. As with most such devices, this interferometer uses a reference path that is not common with the path through the test section. Recombination of the test and reference wavefronts produces interference fringes. Unfortunately, in order to obtain stable fringes, the alignment of both the test and reference paths must be maintained to within a fraction of the wavelength of the light being used for the measurement. Otherwise, the fringes will shift and may disappear. Because these interferometers are extremely sensitive to bumping, jarring and transmitted vibration, they are typically mounted on optical isolation tables. Schlieren deflectometers or the more recent Shack-Hartmann wavefront sensors also measure concentration and temperature in laboratory fluid flows. Ray optics describe the operation of both devices. In a schlieren system, an expanded, collimated beam passes through a test section where refractive index gradients deflect rays. A lens focuses the beam to a filter placed in the rear focal plane of the decollimating lens. In a quantitative color schlieren system, gradients in the index of refraction appear as colors in the field of view due to the action of the color filter. Since sensitivity is a function of the focal length of the decollimating lens, these systems are rather long and filter fabrication and calibration is rather difficult. A Shack-Hartmann wavefront sensor is an array of small lenslets. Typical diameters are on the order of a few hundred microns. Since these lenslets divide the test section into resolution elements, the spatial resolution can be no smaller than an individual lenslet. Such a device was recently used to perform high-speed tomography of heated air exiting a 1.27 cm diameter nozzle. While these wavefront sensors are very compact, the limited spatial resolution and the methods required for data reduction suggest that a more useful instrument needs to be developed. The category of interferometers known as common path interferometers can eliminate much of the vibration sensitivity associated with traditional interferometry as described above. In these devices, division of the amplitude of the wavefront following the test section produces the reference beam. Examples of these instruments include shearing and point diffraction interferometers. In the latter case, shown schematically, a lens focuses light passing through the test section onto a small diffracting object. Such objects are typically either a circle of material on a high quality glass plate or a small sphere in a glass cell. The size of the focused spot is several times larger than the object so that the light not intercepted by the diffracting object forms the test beam while the diffracted light generates a spherical reference beam. While this configuration is mechanically stable, phase shifting one beam with respect to the other is difficult due to the common path. Phase shifting enables extremely accurate measurements of the phase of the interferogram using only gray scale intensity measurements and is the de facto standard of industry. Mercer and Creath 2 demonstrated phase shifting in a point diffraction interferometer using a spherical spacer in a liquid crystal cell as the diffracting object. By changing the voltage across the cell, they were able to shift the phase of the undiffracted beam relative to the reference beam generated by diffraction from the sphere. While they applied this technology to fluid measurements, the device shifted phase so slowly that it was not useful for studying transient phenomena. We have identified several technical problems that precluded operation of the device at video frame rates and intend to solve them to produce a phase-shifting liquid crystal point-diffraction interferometer operating at video frame rates. The first task is to produce high contrast fringes. Since the diffracted beam is much weaker than the transmitted beam, interferograms have poor contrast unless a dye is added to the liquid crystal to reduce the intensity of the undiffracted light. Dyes previously used were not rigorously characterized and suffered from hysteresis in both the initial alignment state of the device and the electro-optic switching characteristics. Hence, our initial effort will identify and characterize dyes that do not suffer from these difficulties and are readily soluble in the liquid crystal host. Since the ultimate goal of this research is to produce interferometers capable of phase shifting at video frame rates, we will quantify the difference in switching times between ferroelectric and nematic liquid crystals. While we have more experience with nematic crystals, they typically switch more slowly than ferroelectric cells. As part of that effort, we will investigate the difference in the modulation of the interferograms as a function of the type of liquid crystal in the cell. Because the temporal switching response of a liquid crystal cell is directly related its thickness, we intend to explore techniques required to produce cells that are as thin as possible. However, the cells must still produce a total phase shift of two pi radians.

Griffin, DeVon W.↗

Advanced Diagnostic System on Earth Observing One

In this infusion experiment, the Livingstone 2 (L2) model-based diagnosis engine, developed by the Computational Sciences division at NASA Ames Research Center, has been uploaded to the Earth Observing One (EO-1) satellite. L2 is integrated with the Autonomous Sciencecraft Experiment (ASE) which provides an on-board planning capability and a software bridge to the spacecraft's 1773 data bus. Using a model of the spacecraft subsystems, L2 predicts nominal state transitions initiated by control commands, monitors the spacecraft sensors, and, in the case of failure, isolates the fault based on the discrepant observations. Fault detection and isolation is done by determining a set of component modes, including most likely failures, which satisfy the current observations. All mode transitions and diagnoses are telemetered to the ground for analysis. The initial L2 model is scoped to EO-1's imaging instruments and solid state recorder. Diagnostic scenarios for EO-1's nominal imaging timeline are demonstrated by injecting simulated faults on-board the spacecraft. The solid state recorder stores the science images and also hosts: the experiment software. The main objective of the experiment is to mature the L2 technology to Technology Readiness Level (TRL) 7. Experiment results are presented, as well as a discussion of the challenging technical issues encountered. Future extensions may explore coordination with the planner, and model-based ground operations.

Hayden, Sandra C.↗

Toward Realistic Dynamics of Rotating Orbital Debris, and Implications for Lightcurve Interpretation

Optical observations of rotating space debris near GEO contain important information on size, shape, composition, and rotational states, but these aspects are difficult to extract due to data limitations and the high number of degrees of freedom in the modeling process. For tri-axial rigid debris objects created by satellite fragmentations, the most likely initial rotation state has a large component of initial angular velocity directed along the intermediate axis of inertia, leading to large angular reorientations of the body on the timescale of the rotation period. This lends some support to the simplest possible interpretation of light curves -- that they represent sets of random orientations of the objects of study. However, effects of internal friction and solar radiation are likely to cause significant modification of rotation states within a time as short as a few orbital periods. In order to examine the rotational dynamics of debris objects under the influences of these effects, a set of seven first-order coupled equations of motion were assembled in state form: three are Euler equations describing the rates of change of the components of angular velocity in the body frame, and four describe the rates of change of the components of the unit quaternion. Quaternions are a four-dimensional extension of complex numbers that form a seamless, singularity-free representation of body orientation on S3. The Euler equations contain explicit terms describing torque from solar radiation in terms of spherical harmonics, and terms representing effects of a prescribed rate of internal friction. Numerical integrations of these equations of motion are being performed, and results will be presented. Initial tests show that internal friction without solar radiation torque leads to rotation about the maximum principal axis of inertia, as required, and solar radiation torque is expected to lead to spin-up of objects. Because the axis of maximum rotational inertia tends to be roughly coincident with the normal to the largest projected cross-sectional area, internal friction is expected to lead to reduced variation of light curve amplitudes at a given phase angle, but a large dependence of the same on phase angle. At a given phase angle, databases are generated which contain reflected intensities for comprehensive sets of equally-likely orientations, represented as unit quaternions. When projected onto three dimensions (S2) and color-coded by intensity, the set is depicted as points within a solid, semi-transparent unit sphere, within which all possible reflected intensities for an object at a given phase angle may be inspected simultaneously. Rotational sequences are represented by trajectories through the sphere. Databases are generated for each of a set of phase angles separately, forming a comprehensive dataset of reflected intensities spanning all object orientations and solar phase angles. Symmetries in the problem suggest that preferred rotation states are likely, defined relative to the object-sun direction in inertial space and relative to the maximum principal axis of inertia in the body coordinate system. Such rotation states may greatly simplify the problem of light curve interpretation by reducing the number of degrees of freedom in the problem.

Ojakangas, Gregory W.↗

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic strain evolution in an optically excited Pt thin film

The structural evolution of a Pt thin film following photo-thermal excitation by 1 ps optical laser pulses was studied with a time resolution of 100 ps over a total time period of 1 ms. Laser pulse fluences below 50 mJ/cm 2 were insufficient to relax the residual stress state of the as-prepared film even after 10 000 pulses. In this fluence regime, a rapid initial lattice expansion and a decrease in the lattice coherence length due to ultrafast photo-thermal heating were observed. The lattice expansion reached a maximum, and the coherence length reached a minimum, 100–200 ps after excitation before monotonically decaying back to their initial values in about 1 µs. Laser pulse fluences greater than 50 mJ/cm 2 produced irreversible stress relaxation within the first 10 optical pulses. In this regime, the lattice expansion was qualitatively similar to that in the low fluence regime, except that the initial structural state was not recovered. The evolution in the coherence length, however, was more complex. Following an initial decrease similar to that observed at low fluence, the coherence length then increased to a broad maximum greater than the initial value, before recovery.

36 MATERIALS SCIENCE↗

Ultrafast Exciton Dynamics of CH 3 NH 3 PbBr 3 Perovskite Nanoclusters

Exciton dynamics of perovskite nanoclusters has been investigated for the first time using femtosecond transient absorption (TA) and time-resolved photoluminescence (TRPL) spectroscopy. The TA results show two photoinduced absorption signals at 420 and 461 nm and a photoinduced bleach (PB) signal at 448 nm. The analysis of the PB recovery kinetic decay and kinetic model uncovered multiple processes contributing to electron-hole recombination. The fast component (~8 ps) is attributed to vibrational relaxation within the initial excited state, and the medium component (~60 ps) is attributed to shallow carrier trapping. The slow component is attributed to deep carrier trapping from the initial conduction band edge (~666 ps) and the shallow trap state (~40 ps). The TRPL reveals longer time dynamics, with modeled lifetimes of 6.6 and 93 ns attributed to recombination through the deep trap state and direct band edge recombination, respectively. In conclusion, the significant role of exciton trapping processes in the dynamics indicates that these highly confined nanoclusters have defect-rich surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Opportunities and Challenges from Artificial Intelligence and Machine Learning for the Advancement of Science, Technology, and the Office of Science Missions

In February 2019, the President signed Executive Order 13859, Maintaining American Leadership in Artificial Intelligence. This order launched the American Artificial Intelligence Initiative, a concerted effort to promote and protect AI technology and innovation in the United States. The Initiative implements a government-wide strategy in collaboration and engagement with the private sector, academia, the public, and like-minded international partners. Among other actions, key directives in the Initiative called for Federal agencies to: Prioritize AI research and development investments, Enhance access to high-quality cyberinfrastructure and data, Ensure that the US maintains an international leadership role in the development of technical standards for AI, and Provide education and training opportunities to prepare the American workforce for the new era of AI. The mission of the Department of Energy (DOE) is to ensure America’s security and prosperity by addressing its energy, environmental, and nuclear challenges through transformative science and technology solutions. In terms of Science and Innovation, the DOE’s mission is to maintain a vibrant US effort in science and engineering as a cornerstone of our economic prosperity with clear leadership in strategic areas. From July to October in 2019, the Argonne, Oak Ridge, and Berkeley National Laboratories hosted a series of four AI for Science Town Hall meetings in Chicago, Oak Ridge, Berkeley, and Washington DC. The four meetings were attended by over 1300 scientists from the 17 DOE Labs, 39 companies, and over 90 universities. The goal of the Town Hall series was ‘to examine scientific opportunities in the areas of artificial intelligence, Big Data, and high-performance computing (HPC) in the next decade, and to capture the big ideas, grand challenges, and next steps to realizing these.’ The discussions at the meetings were captured in the final report of the AI for Science Town Hall meetings.

42 ENGINEERING↗

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗

Quantum magic and computational complexity in the neutrino sector

We consider the quantum magic in systems of dense neutrinos undergoing coherent flavor transformations, relevant for supernova and neutron-star binary mergers. Mapping the three-flavor-neutrino system to qutrits, the evolution of quantum magic is explored in the single scattering angle limit for a selection of initial tensor-product pure states for 𝑁 𝜈 ≤ 8 neutrinos. For |𝜈𝑒⟩ ⊗𝑁𝜈 initial states, the magic, as measured by the 𝛼 = 2 stabilizer Renyi entropy ℳ 2 , is found to decrease with radial distance from the neutrino sphere, reaching a value that lies below the maximum for tensor-product qutrit states. Further, the asymptotic magic per neutrino, ℳ 2 /𝑁 𝜈 , decreases with increasing 𝑁 𝜈 . In contrast, the magic evolving from states containing all three flavors reaches values only possible with entanglement, with the asymptotic ℳ 2 /𝑁 𝜈 increasing with 𝑁 𝜈 . These results highlight the connection between the complexity in simulating quantum physical systems and the parameters of the Standard Model.

computational complexity↗

Optical Response of Soda-Lime Glass Shocked to 14 GPa

To investigate the potential use of soda-lime glass (SLG) as an optical window in shock wave experiments, laser interferometry (VISAR) was used to determine the window correction and the refractive index in SLG samples subjected to plane shock wave compression. Symmetric plate impact experiments were conducted at peak stresses varying between 4 and 14 GPa and laser interferometry ($λ_0$ = 532$\textit {nm}$ wavelength) measurements were used to monitor the particle velocity histories at the impact surface. The higher stress experiments displayed a constant state prior to stress unloading. In contrast, the lower stress experiments (4.3 and 7.8 GPa) displayed an initial constant state followed by weakly time-dependent response prior to stress unloading. Because the actual particle velocities were independently known from the projectile velocity measurements, the window correction factor and the refractive index changes were determined over the stress range examined. Above 7.5 GPa, our window correction and refractive index results show significant differences with previously published data for the same material. Potential reasons for these differences are indicated.

36 MATERIALS SCIENCE↗

Covalency in Fe 2 O 3 and FeO: Consequences for XPS satellite intensity

The covalent character of the interaction between the metal cation and the oxygen ligands has been examined for two Fe oxides with different nominal oxidation states, Fe(II)O, and Fe(III) 2 O 3 . The covalent character is examined for the initial, ground state configuration and for the ionic states involving the removal of a shallow core, Fe 3p, and a deep core, Fe 2p, electron. The covalency is assessed based on novel theoretical analyses of wave functions for the various cases. It is found that the covalency is considerably different for different oxidation states and for different ionized and non-ionized configurations. The changes in covalency for the ions are shown to be responsible for important changes in relaxation energies for X-Ray Photoelectron Spectroscopy (XPS) spectra and in the intensity lost from main XPS peaks to shake satellites. While these consequences are not observables themselves, they are important for the interpretation of the XPS spectra, in particular, for efforts to extract stoichiometries of these iron oxides from XPS data. This is a finding likely applicable across various 3d transition metal oxide materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Current understanding of the physics of type III solar radio bursts

One of the most exciting plasma physics investigations of recent years has been connected with the understanding of a new strong turbulent plasma state excited by propagating electron beams. This new state is initiated on the linear level by parametric instabilities (OTS, modulational, etc.) and results in a very dynamic state composed of collective clusters of modes called solitons, cavitons, spikons, etc. Introduction of these concepts into the classic beam-plasma interaction problem has rendered quasi-linear and weak turbulence theories inapplicable over most of the interesting parameter range, and helped explain many paradoxes connected with the propagation of beams in the laboratory and space. Following a brief review of these nonlinear notions, the means by which their application to type III solar radiobursts has revolutionized understanding of their propagation, radioemission and scaling properties and has guided the in situ observations towards a more complete understanding are demonstrated. A particular burst (May 16, 1971) is analyzed in detail and compared with numerical predictions.

Papadopoulos, K.↗

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗