Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “State 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 1,477 records · Page 82

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗

srlife : a software tool for estimating the life of high temperature concentrating solar receivers. Part I – metallic receivers

Here, this paper introduces srlife, a tool for estimating the structural service life of concentrating solar power (CSP) receivers operating at high temperatures. Supporting both metallic and ceramic receiver designs, srlife is available as open-source software at https://github.com/applied-material-modeling/srlife and can be installed via the PyPi package manager (https://pypi.org). Given basic receiver geometry and incident heat flux, the tool performs thermohydraulic and structural analysis and estimates the life of a receiver. Designed for easy integration into a software stack, including solar field and levelized cost analysis, the tool can be utilized for optimizing receiver designs to meet service life and economic targets. This paper is Part I in a two-part series. Part I discusses the analysis process used to estimate the life of metallic receivers, along with a description of the required input data. Additionally, several heuristics applied within srlife can reduce analysis time significantly while maintaining accurate life estimations for metallic receivers when compared to full analyses. Several examples demonstrating the utility of srlife in receiver design are also discussed. Part II focuses on the life estimation of ceramic receivers, using time-dependent reliability analysis and various ceramic failure models implemented in srlife.

Creep-fatigue analysis↗

Facile alignment estimation in carbon nanotube films using image processing

Whether a macroscopic assembly of carbon nanotubes can exhibit the one-dimensional properties expected from individual nanotubes critically depends on how well the nanotubes are aligned inside the assembly. Therefore, a simple and accurate method for assessing the degree of alignment is desired for the rapid characterization of carbon nanotube films and fibers. Here, we present an end-to-end solution for determining the global and local spatial orientation of carbon nanotubes in films within a short amount of time using a fast, precise, and economical approach based on an image processing method applied to scanning electron microscopy images. Further, we first use Laplacian edge enhancement filtering for improving the appearance of edge regions, which is followed by image partitioning into multiple blocks to capture the nanoscale orientation characteristics and total variation-based image decomposition of these image blocks. We then perform a 2D-fast Fourier transform on the image decomposed textural components of these edge-enhanced image blocks to determine the orientation distribution, which is utilized to estimate the 2D nematic order parameter. To show the effectiveness of our method, we corroborated our results against results obtained with other state-of-the-art image processing and experimental techniques.

2D-FFT↗

Fault-tolerant grid frequency measurement algorithm during transients

Many critical electric grid operations rely on accurate grid frequency measurements. Unfortunately, the measurement accuracy can be easily undermined by power system transient faults. During a power system transient fault, the power grid voltages and currents are usually highly distorted by high-frequency components. What is worse, the power grid signals could have discontinuity during some system transient faults such as phase angle jump, and the discontinuity could result in large measurement errors to state-of-the-art grid measurement algorithms. In this study, a fault-tolerant grid frequency measurement algorithm during transients is proposed. The new algorithm consists of two stages. The first stage is a transient detector, and it can detect the occurrence of system transient faults instantaneously. The second stage is the intelligent frequency estimator, and it will adapt its measurements according to the transient detector. The performance of the algorithm is evaluated under different steady-state and transient conditions. Both dependability and security of the fault-tolerant algorithm are assessed by using PSCAD simulation data and IEEE Standard test data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying Value and Representing Competitiveness of Electricity System Technologies in Economic Models

Evaluating competition between electricity technologies is challenging because it depends on both their costs and their values. While technology costs can typically be estimated from projections of the cost components - capital, fuel, and O&M - estimating a technology's value is more complex due to its dependence on its contributions to multiple different grid services, each with prices that can vary substantially over space and time. In this work, using an electricity model of the contiguous United States, we develop relationships between relative value and share of total generation for major electricity generation technologies which, when paired with projections of technology costs, can be used to estimate technology competitiveness. We identify significant differences in the relationship between relative value and generation share for variable renewable energy (VRE) and non-VRE sources, but we demonstrate that all technologies require consideration of their dynamic values (in addition to cost) when evaluating competitiveness. In addition, we demonstrate that relative value of a technology is substantially impacted by not only its own generation share but also other aspects of the system state, in particular the mix of other technologies present in the system. Finally, we use the developed relative value relationships in combination with projections of future technology costs in a coarse resolution model that competes technologies based on a comprehensive competitiveness metric: profitability-adjusted LCOE (PLCOE). We show that this simple representation of technology competition approximately recovers the generation mix from a detailed model, which is not possible using LCOE alone. Such an approach can be used to improve the representation of technology competition in coarse-resolution models such as integrated assessment models, for which simplified metrics are often needed.

electricity model↗

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

A semi-analytic estimate for the effective sound speed counterterm in the EFTofLSS

The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. Here, we apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value.

Cosmological perturbation theory in GR and beyond↗

The development of natural gas and hydrogen pipeline capital cost estimating equations

Natural gas pipeline cost data collected by the Oil and Gas Journal (O&GJ) for interstate pipelines constructed from 1980 through 2017 were used to develop capital cost estimating equations that are a function of pipeline diameter, length, and U.S. region. Equations were developed for material, labor, miscellaneous, and right-of-way costs, the four cost components in the O&GJ data, for six different regions of the United States (U.S.). Each equation is a function of pipeline diameter and length. Adjustment mechanisms were then developed for converting the natural gas pipeline equations into equations for estimating the costs of hydrogen pipelines. These adjustments were based in part on an analysis completed by the National Institute for Standards and Technology (NIST). Here the results of this work were used to update cost models in the Hydrogen Delivery Scenario Analysis Model (HDSAM), developed by Argonne National Laboratory for the U.S. Department of Energy’s Hydrogen Program. Our analysis shows a wide range of pipeline cost across different U.S. regions, especially with respect to labor and right-of-way costs. The developed cost formulas for hydrogen pipelines are both important and timely as hydrogen is being considered as a zero-carbon energy carrier with the potential to decarbonize all energy sectors, and the cost of hydrogen transportation is essential for techno-economic analysis of its potential use in these sectors.

03 NATURAL GAS↗

Diffuse boosted cosmic neutrino background

Energetic cosmic rays scatter off the cosmic neutrino background throughout the history of the Universe, yielding a diffuse flux of cosmic relic neutrinos boosted to high energies. We calculate this flux under different assumptions of the cosmic-ray flux spectral slope and redshift evolution. The nonobservation of the diffuse flux of boosted relic neutrinos with current high-energy neutrino experiments already excludes an average cosmic neutrino background overdensity larger than ∼10 4 over cosmological distances. We discuss the future detectability of the diffuse flux of boosted relic neutrinos in light of neutrino overdensity estimates and cosmogenic neutrino backgrounds.

Astronomy & Astrophysics↗

Combining Eddy Covariance Towers, Field Measurements, and the MEMS 2 Ecosystem Model Improves Confidence in the Climate Impacts of Bioenergy With Carbon Capture and Storage

ABSTRACT Carbon dioxide removal technologies such as bioenergy with carbon capture and storage (BECCS) are required if the effects of climate change are to be reversed over the next century. However, BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty. We place these net ecosystem carbon balance forcings in the context of associated radiative forcings from changes to soil N 2 O and CH 4 fluxes, land surface albedo, embedded fossil fuel use, and geologically stored carbon. Results from long term eddy covariance measurements, a soil and plant carbon inventory, and the MEMS 2 process‐based ecosystem model all agree that establishing perennials such as switchgrass or mixed prairie on former cropland resulted in net negative radiative forcing (i.e., global cooling) of −26.5 to −39.6 fW m −2 over 100 years. Establishing these perennials on former grassland sites had similar climate mitigation impacts of −19.3 to −42.5 fW m −2 . However, the largest climate mitigation came from establishing corn for BECCS on former cropland or grassland, with radiative forcings from −38.4 to −50.5 fW m −2 , due to its higher plant productivity and therefore more geologically stored carbon. Our results highlight the strengths and limitations of each method for quantifying the field scale climate impacts of BECCS and show that utilizing multiple methods can increase confidence in the final radiative forcing estimates.

Falvo, Grant [Department of Plant, Soil and Microb↗

Full‐Waveform Simulation of Infrasound Propagation in the Atmosphere: A Case Study of the 2023 April 20 SpaceX Starship Explosion

Infrasound, low‐frequency sound below 20 Hz, has been a key technology to monitor explosion events in the atmosphere. The International Monitoring System (IMS) of the Comprehensive Nuclear‐Test‐Ban Treaty Organization provides the means for continuous monitoring of infrasonic events worldwide. Infrasonic techniques for event location and size estimation can also complement other observational techniques for the detection and characterization of the entry of asteroids or large meteoroids. In this study, we describe the detection capability of IMS infrasound stations for an explosive event in the middle of the atmosphere. Full‐waveform simulations are performed with the specification of atmospheric conditions and incorporated into the event location and explosion yield estimation. We applied it to the 2023 April 20 SpaceX Starship explosion at 29 km altitude. Starship is a super heavy‐lift space vehicle constructed by SpaceX and known as the largest and most powerful rocket ever built. The Starship explosion created huge pressure disturbances in the atmosphere, and its infrasound was detected by the IMS arrays in North America between 2000 and 4000 km. Independent observational data and available ground‐truth information provide a rare opportunity to evaluate the monitoring capability of the IMS network for elevated sources in the atmosphere. We also demonstrate the capability of full‐waveform simulation for infrasound wavefield characterization and prediction to improve event location and yield estimation.

Geosciences↗

CosmoLensNRE

Cosmology Inference from Strong Gravitational Lensing using Neural Ratio Estimation

Jarugula, Sreevani [Fermi National Accelerator Lab↗

Quantification of regional net CO 2 flux errors in the Orbiting Carbon Observatory-2 (OCO-2) v10 model intercomparison project (MIP) ensemble using airborne measurements

Inverse model intercomparison projects (MIPs) provide a chance to assess the uncertainties in inversion estimates arising from various sources. However, accurately quantifying ensemble CO 2 flux errors remains challenging and often relies on the ensemble spread. This study proposes a method for quantifying the errors in regional net surface–atmosphere CO 2 flux estimates from models taken from the Orbiting Carbon Observatory-2 (OCO-2) v10 MIP by using independent airborne CO 2 measurements for the period 2015–2017. We first calculate the root mean square error (RMSE) between the ensemble mean of posterior CO 2 concentrations and airborne observations and then isolate the CO 2 concentration errors caused solely by the ensemble mean of posterior net fluxes by subtracting the observation, representation, and transport errors from seven regions. Our analysis reveals that the flux errors projected onto CO 2 space account for 55 %–85 % of the regional average RMSE over the 3 years, ranging from 0.88 to 1.91 ppm. In five regions, the error estimates based on observations exceed those computed from the ensemble spread of posterior fluxes by a factor of 1.33–1.93, implying an underestimation of the actual flux errors, while their magnitudes are comparable in two regions. The adjoint sensitivity analysis identifies that the underestimation of flux errors is prominent where the magnitudes of fossil fuel emissions exceed those of terrestrial-biosphere fluxes by a factor of 3–31 over the 3 years. This suggests the presence of systematic biases in the inversion estimates associated with errors in the prescribed fossil fuel emissions common to all models. Our study emphasizes the value of airborne measurements for quantifying regional errors in ensemble net CO 2 flux estimates.

54 ENVIRONMENTAL SCIENCES↗

Regional Oil and gas Aerial Methane Synthesis model (Analytica) (ROAMS Analytica) v1.5.2

The Regional Oil and gas Aerial Methane Synthesis model (Analytica) is a tool to convert the results of wide-area, source-resolved aerial methane remote sensing surveys of oil and natural gas infrastructure in a given region into methane emissions inventories (estimates of the magnitude and breakdown of methane emissions from the surveyed infrastructure). This version is written in the Analytica programming language, and this version accompanies a correction in preparation for submission to Sherwin et al. 2024 (Nature). The tool leverages databases of source-resolved methane emissions detected in aerial surveys, aerial survey coverage information (which areas were measured and when), data summarizing surveyed oil and natural gas infrastructure and production (derived from third-party databases), as well as state-of-the-art mechanistic emissions simulation tools to characterize emissions too small for the aerial system to see. The regional methane emissions estimates produced by this tool are much more granular in both space and asset type than common satellite- or flux tower-based regional estimates. Unlike other tools for converting site-level measurements into regional emissions estimates, our unique geostatistical approach integrates aerially measured emissions with limited need for statistical extrapolation, which can be highly sensitive to modeler assumptions. As a result, ROAMS-based estimates of regional methane emissions from oil and gas activity are widely viewed as highly credible, as evidenced by the success of Dr. Sherwin's recent paper in Nature.

Sherwin, Evan [Lawrence Berkeley National Laborato↗

Autonomous monitoring of algal biomass: Success stories and lessons learned from long-term field deployment

Autonomous, high-frequency monitoring of outdoor algal ponds is needed to quantify biomass productivity and detect culture decline in environments prone to contamination, grazers, and variable operating conditions. We report successes and lessons learned in translating a laboratory spectroradiometric monitoring approach to a multi-year autonomous field deployment at the Arizona Center for Algae Technology and Innovation (AzCATI). The system measures spectrally resolved pond reflectance by ratioing upwelling radiance from each raceway to simultaneous downwelling sky irradiance using fiber-coupled spectrometers. A physics-based reflectance model (ASHARP) is fit to each spectrum pair to estimate optical parameters, including a biomass-proxy coefficient (C a ) which enables near-real-time tracking of biomass accumulation and culture state at 2–5 min intervals. From May 2022 through September 2025 the platform operated continuously while scaling from two to six raceway ponds. Several strains of algae were monitored successfully, including the high productivity Tetraselmis striata and Picochlorum celeri. Transitioning data acquisition from a Windows laptop to a Raspberry Pi improved uptime from 57% (2022) to ~89% (2024–2025) and enabled routine real-time analysis. Further, we converted relative biomass estimates to absolute ash-free dry weight (AFDW) using experimentally-derived calibrations, providing field-relevant biomass predictions with conservative confidence bounds. These results demonstrate the feasibility of long-term, autonomous optical monitoring for well-mixed open-raceway algal cultivation and provide practical guidance for reliable field operation and scaling.

Katinas, Christopher Michael [Sandia National Labo↗

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of Regional Interconnections to Meet Southeast Asia's Renewable Energy Targets: Phase I of the Singapore – U.S. Collaboration

This report summarizes the results of a study of long-distance subsea electrical interconnections in Southeast Asia, conducted under the first phase of the Singapore-U.S. collaboration under the Net Zero World Initiative. In Phase 1, the Net Zero World team evaluated the current state of the renewable energy landscape in Southeast Asia; developed an open-source capacity expansion model using GridPath and evaluated several future scenarios including changes in interregional transmission; explored the technical feasibility of regional subsea interconnectors; and estimated socio-economic benefits of expanded regional transmission.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗