SEARCH · Engineering Papers

Results for “spatial constraints”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗

NASA NTRS2006

On Ambiguities in SAR Design

Ambiguities are an aliasing effect caused by the periodic sampling of the scene backscatter inherent to pulsed radar systems such as Synthetic Aperture radar (SAR). In this paper we take a fresh look at the relationship between SAR range and azimuth ambiguity constraints on the allowable pulse repetition frequency (PRF) and the antenna length. We show that for high squint angles smaller antennas may be feasible in some cases. For some applications, the ability to form a synthetic aperture at high squint angles is desirable, but the size of the antenna causes problems in the design of systems capable of such operation. This is because the SAR system design is optimized for a side-looking geometry. In two examples design examples we take a suboptimum antenna size and examine the performance in terms of azimuth resolution and swath width as a function of squint angle. We show that for stripmap SARs, the swath width is usually worse for off-boresight squint angles, because it is severely limited by range walk, except in cases where we relax the spatial resolution. We consider the implications for the design of modest-resolution, narrow swath, scanning SAR scatterometers .

Synthetic Aperture radar (SAR)↗

NASA NTRS2012

Dynamic Range and Sensitivity Requirements of Satellite Ocean Color Sensors: Learning from the Past

Sensor design and mission planning for satellite ocean color measurements requires careful consideration of the signal dynamic range and sensitivity (specifically here signal-to-noise ratio or SNR) so that small changes of ocean properties (e.g., surface chlorophyll-a concentrations or Chl) can be quantified while most measurements are not saturated. Past and current sensors used different signal levels, formats, and conventions to specify these critical parameters, making it difficult to make cross-sensor comparisons or to establish standards for future sensor design. The goal of this study is to quantify these parameters under uniform conditions for widely used past and current sensors in order to provide a reference for the design of future ocean color radiometers. Using measurements from the Moderate Resolution Imaging Spectroradiometer onboard the Aqua satellite (MODISA) under various solar zenith angles (SZAs), typical (L(sub typical)) and maximum (L(sub max)) at-sensor radiances from the visible to the shortwave IR were determined. The Ltypical values at an SZA of 45 deg were used as constraints to calculate SNRs of 10 multiband sensors at the same L(sub typical) radiance input and 2 hyperspectral sensors at a similar radiance input. The calculations were based on clear-water scenes with an objective method of selecting pixels with minimal cross-pixel variations to assure target homogeneity. Among the widely used ocean color sensors that have routine global coverage, MODISA ocean bands (1 km) showed 2-4 times higher SNRs than the Sea-viewing Wide Field-of-view Sensor (Sea-WiFS) (1 km) and comparable SNRs to the Medium Resolution Imaging Spectrometer (MERIS)-RR (reduced resolution, 1.2 km), leading to different levels of precision in the retrieved Chl data product. MERIS-FR (full resolution, 300 m) showed SNRs lower than MODISA and MERIS-RR with the gain in spatial resolution. SNRs of all MODISA ocean bands and SeaWiFS bands (except the SeaWiFS near-IR bands) exceeded those from prelaunch sensor specifications after adjusting the input radiance to L(sub typical). The tabulated L(sub typical), L(sub max), and SNRs of the various multiband and hyperspectral sensors under the same or similar radiance input provide references to compare sensor performance in product precision and to help design future missions such as the Geostationary Coastal and Air Pollution Events (GEO-CAPE) mission and the Pre-Aerosol-Clouds-Ecosystems (PACE) mission currently being planned by the U.S. National Aeronautics and Space Administration (NASA).

chlorophyll fluorescence↗

DOE OSTI2023

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Observational Study of the Thermodynamics and Morphological Characteristics of a Midlatitude Continental Cold Pool Event

Abstract Convectively generated cold pools are important to the Earth system as they exert strong controls on deep convective-storm initiation, intensity, and life cycle. Despite their importance, efforts to introduce such cold pool controls into weather and climate models lack guidance and/or physical constraints from cold pool observations. This work presents a detailed, purely observational analysis of a cold pool event that took place on 23–24 May 2011 in north-central Oklahoma. The characteristics of the cold pool, and the spatiotemporal evolution of the hydrometeors and dynamics in the proximity of the cold pool, are studied with high-resolution observations. The unprecedented dataset used in this work to study cold pool characteristics includes an enhanced network of surface weather stations, a high-temporal-frequency sounding array, and the NEXRAD and Atmospheric Radiation Measurement (ARM) Southern Great Plains radar networks. The potential use of NEXRAD surveillance scans to estimate height and propagation speed of the leading edge of the cold pool (LECP) is presented in this work. Manual identification and tracking of the LECP from NEXRAD imagery shows a spatial and temporal heterogeneity of the LECP properties. Surprisingly, over its detected life cycle, the LECP speed remains almost constant, even though the strength of the cold pool diminishes in time and its height varies. Radar analysis shows that pulses of graupel and hail within downdrafts in the convective system generating the cold pool appeared to be related to temporary increases in the LECP height.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Interseismic Strain Accumulation between the Colorado Plateau and the Eastern California Shear Zone: Implications for the Seismic Hazard near Las Vegas, Nevada

ABSTRACT Crustal deformation in the central Basin and Range between the Colorado plateau and the Eastern California Shear Zone is active but slow, making it a challenge to assess how strain is distributed and crustal motion transferred. However, knowledge of strain rates is very important, particularly for addressing the seismic hazard for both the Las Vegas urban area and the site of the proposed Yucca Mountain nuclear waste repository, in southern Nevada. Global Positioning System (GPS) data provide important constraints, particularly now that the GPS network in the area has substantially expanded in recent years. However, because deformation is slow, it is important to mitigate any transient tectonic and nontectonic signals to obtain the most accurate long-term interseismic motion and robust estimation of strain rates. We use data from all GPS stations in the region including both long-running continuous and semicontinuous stations. We model and remove postseismic displacements at these stations using source parameters for 41 events, dating back to the 1700 Cascadia megathrust earthquake, which contribute significantly to the deformation field within the central Basin and Range. We also remove correlated noise from the time series with the common-mode component imaging technique. We find that removal of both the postseismic transients and common-mode noise substantially reduces the uncertainties and spatial variation in the velocities. We find east–west extension across the Las Vegas Valley of 0.5–0.6 mm/yr. The interseismic strain rate field, calculated with the final velocities, reveals higher strain rates through southern Nevada than in previous studies, with rates within Las Vegas Valley of 8.5±2.4×10−9 yr−1. Our results also confirm shear along the Pahranagat shear zone, but the estimated amplitude is strongly affected by postseismic relaxation.

Geochemistry & Geophysics↗

NASA NTRS2020

Global-scale distribution of ozone in the remote troposphere from ATom and HIPPO airborne field missions

Ozone is a key constituent of the troposphere, where it drives photochemical processes, impacts air quality, and acts as a climate forcer. Large-scale in situ observations of ozone commensurate with the grid resolution of current Earth system models are necessary to validate model outputs and satellite retrievals. In this paper, we examine measurements from the Atmospheric Tomography (ATom; four deployments in 2016–2018) and the HIAPER Pole-to-Pole Observations (HIPPO; five deployments in 2009–2011) experiments, two global-scale airborne campaigns covering the Pacific and Atlantic basins. ATom and HIPPO represent the first global-scale, vertically resolved measurements of O3 distributions throughout the troposphere, with HIPPO sampling the atmosphere over the Pacific and ATom sampling both the Pacific and Atlantic. Given the relatively limited temporal resolution of these two campaigns, we first compare ATom and HIPPO ozone data to longer-term observational records to establish the representativeness of our dataset. We show that these two airborne campaigns captured on average 53 %, 54 %, and 38 % of the ozone variability in the marine boundary layer, free troposphere, and upper troposphere–lower stratosphere (UTLS), respectively, at nine well-established ozonesonde sites. Additionally, ATom captured the most frequent ozone concentrations measured by regular commercial aircraft flights in the northern Atlantic UTLS. We then use the repeated vertical profiles from these two campaigns to confirm and extend the existing knowledge of tropospheric ozone spatial and vertical distributions throughout the remote troposphere. We highlight a clear hemispheric gradient, with greater ozone in the Northern Hemisphere, consistent with greater precursor emissions and consistent with previous modeling and satellite studies. We also show that the ozone distribution below 8 km was similar in the extra-tropics of the Atlantic and Pacific basins, likely due to zonal circulation patterns. However, twice as much ozone was found in the tropical Atlantic as in the tropical Pacific, due to well-documented dynamical patterns transporting continental air masses over the Atlantic. Finally, we show that the seasonal variability of tropospheric ozone over the Pacific and the Atlantic basins is driven year-round by transported continental plumes and photochemistry, and the vertical distribution is driven by photochemistry and mixing with stratospheric air. This new dataset provides additional constraints for global climate and chemistry models to improve our understanding of both ozone production and loss processes in remote regions, as well as the influence of anthropogenic emissions on baseline ozone.

ATom, aircraft profiles, Atlantic pollution, ozone↗

NASA NTRS2024

Testing of Advanced Capabilities to Enable In-time Safety Management and Assurance for Future Flight Operations

In order to refine an initial Concept of Operations, explore Concepts of Use, and expose/validate requirements for future In-Time Aviation Safety Management Systems (IASMS), testing architectures were created, along with a set of capabilities and underlying information exchange protocols. These systems were conceived and developed based on hazards associated with two envisioned urban area flight domains: (1) highly autonomous small uncrewed aerial systems (sUAS) operating at low altitudes, and (2) highly autonomous air taxis. The initial scope of this development is described in [1]; this report provides an update, focusing on the subsequent developments and test activities. As stated in [1], it is important to note that there are many capabilities already in use by the industry (or soon to be in use) that will play critical roles in future IASMS designs. Those reported here were developed to address a gap in the current state-of-the-art regarding specific hazards/risks, and/or to allow for investigation of the interplay between and across hazard types — particularly regarding how overall safety risk can be reduced or managed effectively. Results of testing and development activities are organized by the operational phase wherein a particular capability would be employed (i.e., preflight, in-flight, and post-flight/off-line). Pre-flight: A set of capabilities were developed to help mitigate safety risk prior to flight (e.g., during flight and mission planning). Results of testing summarize (1) validation activities to raise the Technology Readiness Level (TRL) and (2) evaluation activities where the capabilities were applied to flight/mission planning procedures and used by operators/pilots. For the latter, flight plans were automatically assessed, and operators/pilots were notified of hazardous flight segments so as to enable adjustment of the flight plan and re-evaluation, and/or to better inform go/no-go decisions. Capabilities addressed hazards associated with power consumption, third-party risk, wind, navigation system performance, radiofrequency interference, and proximity to geo-spatial threats (e.g., buildings, trees, and no-fly zones). In-flight: Flight experiments tested capabilities that detect and respond to hazards encountered during flight. In the first series, safety hazards were monitored and assessed onboard, and system-generated mitigation maneuvers were recorded (but not acted upon by the vehicle). In the second series, mitigation maneuver commands directed the aircraft in response to safety hazards (i.e., auto-mitigation). The sUAS used for testing is described in full, as is the test architecture, which included commercial avionics, research avionics, and onboard software designed to detect, assess, and respond to hazards. The onboard system was designed as a run-time assurance framework, consistent with [2] and supportive of both supervisory and automated modes. The primary functions included: real-time risk assessment (RTRA), auto-pilot monitoring, constraint monitoring, and contingency select/triggering. RTRA performs integrated risk assessment considering data from several hazard-related monitors (e.g., battery, motors, navigation, communications, population density, and loss-of-control). Post-flight/off-line: Data monitored and recorded during flights can enable IASMS capabilities that execute after flights have completed (or “off-line”). These include: (1) the ability to identify anomalies and trends that may only be observable when comparing data spanning a number of similar flights; (2) the ability to update and validate pre-flight and in-flight capabilities and any underlying models to improve their performance; (3) the ability to report anomalies/off-nominals that may indicate design changes or maintenance actions are needed; and (4) the ability for humans involved in operations to report safety-relevant observations to help in understanding the flight data and/or the operational context of a flight. Progress on three such capabilities is summarized; the first investigates anomaly detection given a limited set of flight logs and applies an approach previously used for space operations. The second explores what could be identified using a larger set of flight logs, including from web-based forums where flight logs are posted by sUAS autopilot users. The third creates a new means of collecting information on UAS incidents and accidents via the Aviation Safety Reporting System (ASRS).

sUAS↗

DOE OSTI2023

Development of a Geothermal Module in reV: Quantifying the Geothermal Potential while Accounting for the Geospatial Intersection of the Grid Infrastructure and Land Use Characteristics

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. This paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results should be considered with care due to the highly uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for five scenarios: four hydrothermal (3.5km depth) and one EGS (4.5km depth). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

exclusions↗

DOE OSTI2024

Development of a Geothermal Module in reV: Quantifying the Geothermal Potential While Accounting for the Geospatial Intersection of the Grid Infrastructure and Land Use Characteristics: Preprint

exclusions↗

DOE OSTI2023

Renewable Energy Potential Model: Geothermal Supply Curves

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. The included paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results provided here should be considered with care due to the high uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for three scenarios: two hydrothermal (3.5km depth, USGS heat flow & SMU temperatures respectively) and one EGS (4.5km depth, SMU temperatures). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Confronting Dark Matter with the Multiwavelength Sky (Final Technical Report)

The existence of dark matter is a key clue to new fundamental physics, yet its nature and non-gravitational interactions remain mysterious. In the course of this five-year project (spanning July 2015 - July 2020), we have developed novel analysis techniques to separate dark matter signals from astrophysical backgrounds in the inner Galaxy, and measured hitherto-unknown properties of the Fermi Bubbles and the Galactic Center excess in gamma-rays from the inner Milky Way. We have improved predictions for the spatial distribution and energy spectra of photon signals from dark matter, and set new limits from astrophysical and cosmological datasets on dark matter inhabiting complex dark sectors, multi-particle dark matter annihilation, and tiny primordial black holes. We have successfully adapted effective field theory techniques from collider physics to make high-precision predictions for signals from heavy weakly-interacting DM, allowing us to essentially exclude the possibility that DM is a thermal wino, and have discovered novel theoretical and observational implications for dark matter bound states and self-interactions. We have developed a comprehensive public code package for modeling the effects of exotic energy injections, such as from dark matter annihilation and decay, on the ionization and temperature history of the early cosmos, and have used these tools to set robust, powerful and broadly-applicable constraints on models of light dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗

DOE OSTI2025

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Directivity Modes of Earthquake Populations with Unsupervised Learning

We present a novel approach for resolving modes of rupture directivity in large populations of earthquakes. A seismic spectral decomposition technique is used to first produce relative measurements of radiated energy for earthquakes in a spatially compact cluster. The azimuthal distribution of energy for each earthquake is then assumed to result from one of several distinct modes of rupture propagation. Rather than fitting a kinematic rupture model to determine the most likely mode of rupture propagation, we instead treat the modes as latent variables and learn them with a Gaussian mixture model. The mixture model simultaneously determines the number of events that best identify with each mode. The technique is demonstrated on four datasets in California, each with compact clusters of several thousand earthquakes with comparable slip mechanisms. We show that the datasets naturally decompose into distinct rupture propagation modes that correspond to different rupture directions, and the fault plane is unambiguously identified for all cases. We find that these small earthquakes exhibit unilateral ruptures 63–73% of the time on average. Here, the results provide important observational constraints on the physics of earthquakes and faults.

58 GEOSCIENCES↗

DOE OSTI2026

Stable and tunable MeV $$\gamma$$-ray generation via dual-laser inverse Thomson scattering from a laser-plasma accelerator

Abstract Inverse Thomson scattering from laser-plasma accelerators offers a pathway to compact, tunable MeV $$\gamma$$ -ray sources for reduced-dose radiography and enhanced performance in nuclear resonance fluorescence (NRF)-based isotope identification. However, photon yield and spectral quality are often limited by constraints on interaction geometry and scatter-laser tunability. Here we demonstrate a MeV $$\gamma$$ -ray source based on a dual-laser inverse Thomson scattering configuration driven by a 100-TW laser-plasma accelerator. Electron beams tunable from 122 to 204 MeV with $$<5$$ mrad divergence and $$<1$$ mrad pointing stability generate $$\gamma$$ rays with peak energies from 276 keV to 1.2 MeV and yields up to $$2\times 10^{7}$$ photons per shot. By independently controlling the interaction position and the scatter-pulse duration, we experimentally match the scatter pulse to the walk-off-limited interaction length. Extending the scatter pulse to 200 fs increases photon production by approximately $$15\%$$ while maintaining operation in the linear Thomson regime, thereby preserving narrow spectral bandwidth and controlled radiation divergence. Radiographic characterization demonstrates MeV-level penetration and $$\approx 0.1$$ mm spatial resolution, while stable operation is sustained over multi-hour timescales across multiple days. These results show that interaction-length optimization provides a scalable strategy for improving photon yield, spectral control, and operational stability in compact laser-plasma-accelerator-driven $$\gamma$$ -ray sources.

Tsai, Hai-En↗

NASA NTRS1981

The Gemini-Monoceros X-ray enhancement - A giant X-ray ring

A 1.5 deg spatial resolution map of the 1/4 KeV diffuse X-ray background enhancement in the Gemini and Monoceros constellations is found to show a circular ring-shaped emission feature with a diameter of 20 deg. The feature and possible X-ray contributions from the Mon OB1 association and neighboring supernova remnants are discussed. From 300 pc, the region has a radius of 50 pc, with an emitting electron density of 0.01 per cu cm. A shell of expanding neutral hydrogen and nonthermal radio spur is observed outside the ring with the X-ray emitting pulsar PSR 0656 + 14 lying close to the center of the ring. Origins of the ring are discussed, ruling out formation by the association Mon OB1. The ring is considered to be a field supernova remnant formed by the progenitor of the central pulsar, and providing constraints on theories of remnant evolution. This conclusion is found to agree with estimated supernova rates, and the absence of additional examples of this stage of evolution is an observational selection effect.

Nousek, J. A.↗

NASA NTRS1999

Non-Intrusive Laser-Induced Imaging for Speciation and Patternation in High Pressure Gas Turbine Combustors

The next generation of was turbine combustors for aerospace applications will be required to meet increasingly stringent constraints on fuel efficiency, noise abatement, and emissions. The power plants being designed to meet these constraints will operate at extreme conditions of temperature and pressure, thereby generating unique challenges to the previously employed diagnostic methodologies. Current efforts at NASA Glenn Research Center (GRC) utilize optically accessible, high pressure flametubes and sector combustor rigs to probe, via advanced nonintrusive laser techniques, the complex flowfields encountered in advanced combustor designs. The fuel-air mixing process is of particular concern for lowering NO(x) emissions generated in lean, premixed engine concepts. Using planar laser-induced fluorescence (PLIF) we have obtained real-time, detailed imaging of the fuel spray distribution for a number of fuel injector over a wide range of operational conditions that closely match those expected in the proposed propulsion systems. Using a novel combination of planar imaging, of fuel fluorescence and computational analysis that allows an examination of the flowfield from any perspective, we have produced spatially and temporally resolved fuel-air distribution maps. These maps provide detailed insight into the fuel injection at actual conditions never before possible, thereby greatly enhancing the evaluation of fuel injector performance and combustion phenomena.

Locke, Randy J.↗