Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model error”

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 199 records · Page 11

Assimilation of Satellite-Derived Soil Moisture for Improved Forecasts of the Great Plains Low-Level Jet

In the context of forecasting societally impactful Great Plains low-level jets (GPLLJs), the potential added value of satellite soil moisture (SM) data assimilation (DA) is high. GPLLJs are both sensitive to regional soil moisture gradients and frequent drivers of severe weather, including mesoscale convective systems. An untested hypothesis is that SM DA is more effective in forecasts of weakly synoptically forced, or uncoupled GPLLJs, than in forecasts of cyclone-induced coupled GPLLJs. Using the NASA Unified Weather Research and Forecasting (NU-WRF) Model, 75 GPLLJs are simulated at 9-km resolution both with and without NASA Soil Moisture Active Passive SM DA. Differences in modeled SM, surface sensible (SH) and latent heat (LH) fluxes, 2-m temperature (T2), 2-m humidity (Q2), PBL height (PBLH), and 850-hPa wind speed (W850) are quantified for individual jets and jet-type event subsets over the south-central Great Plains, as well as separately for each GPLLJ sector (entrance, core, and exit). At the GPLLJ core, DA-related changes of up to 5.4 kg m–2 in SM can result in T2, Q2, LH, SH, PBLH, and W850 differences of 0.68°C, 0.71 g kg –2 , 59.9 W m –2 , 52.4 W m –2 , 240 m, and 4 m s –1 , respectively. W850 differences focus along the jet axis and tend to increase from south to north. Further, jet-type differences are most evident at the GPLLJ exit where DA increases and decreases W850 in uncoupled and coupled GPLLJs, respectively. Data assimilation marginally reduces negative wind speed bias for all jets, but the correction is greater for uncoupled GPLLJs, as hypothesized.

54 ENVIRONMENTAL SCIENCES↗

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗

2020 natural gas LCA appendices Rev1

This collection is the data-centric appendices for the report Life Cycle Analysis of Natural Gas Extraction and Power Generation: U.S. 2020 Emissions Profile. It consists of Appendix A: Additional Modeling Parameters [spreadsheet]; Appendix B: Water Burdens [spreadsheet]; Appendix D: Simulation of Liquids Unloading [python script and spreadsheet]; Appendix E: Detailed GHG Results for All Scenarios [spreadsheet]; Appendix F: Full Inventory Results [spreadsheet]; and Appendix I: Stage-Level Natural Gas Loss and Consumption Rates [spreadsheet]. These results have been updated from the previous version (https://edx.netl.doe.gov/dataset/2020-natural-gas-lca-data-appendices) to correct a modeling error where the same post-processing natural gas composition was used instead of the intended regional compositions.

Appendices↗

NETL Natural Gas Lifecycle Model 2020 Rev1

This is the excel-based life cycle model that contains all the parameters and Monte Carlo simulation capabilities to model the techno-regions contained in the report: Life Cycle Analysis of Natural Gas Extraction and Power Generation: U.S. 2020 Emissions Profile. This is an update to a previous version (https://edx.netl.doe.gov/dataset/netl-natural-gas-lifecycle-model) and corrects a modeling error in the post-processing natural gas composition.

Calcium Chloride↗

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗

Life Cycle Analysis of Natural Gas Extraction and Power Generation: U.S. 2020 Emissions Profile

This analysis expands upon previous life cycle analyses (LCAs) of natural gas systems performed by the National Energy Technology Laboratory (NETL). It provides a complete inventory of emissions to air and water, water consumption, and land use change. These environmental burdens are detailed for all supply chain steps from natural gas production through natural gas distribution. This package includes the report, the NETL Natural Gas Model, and appendices that include several Excel workbooks and a python script to provide transparent access to the calculations and resulting data. This is revision 1 of the 2024 study (published December 17, 2024 and updated on January 24, 2025) and corrects a modeling error in natural gas composition. See the errata on page 2 for more information.

03 NATURAL GAS↗

Emerging low-cloud feedback and adjustment in global satellite observations

From mid-2003 to mid-2024, a global decrease in low-cloud amount enhanced the absorption of solar radiation by 0.22±0.07 W m −2 per decade (±1σ range), accelerating the energy imbalance trend during that period (0.44 W m −2 per decade). Through controlling factor analysis, here we show that the low-cloud trend is due to a combination of cloud feedback and adjustments to greenhouse gases and aerosols (respectively 0.09±0.02, 0.05±0.03, and 0.03±0.03 W m −2 per decade), which jointly account for 74 % of the trend. The contribution of natural climate variability is weak but uncertain (0.01±0.08 W m −2 per decade), owing to a poorly constrained trend in boundary-layer inversion strength. Importantly, the observed low-cloud radiative trend lies well within the range of values simulated by contemporary global climate models under conditions close to present day. Any systematic model error in the representation of present-day global energy imbalance trends is thus likely to originate in processes unrelated to low clouds.

Geosciences↗

\texttt{qec\_code\_sim}: An open-source Python framework for estimating the effectiveness of quantum-error correcting codes on superconducting qubits

Quantum computers are highly susceptible to errors due to unintended interactions with their environment. It is crucial to correct these errors without gaining information about the quantum state, which would result in its destruction through back-action. Quantum Error Correction (QEC) provides information about occurred errors without compromising the quantum state of the system. However, the implementation of QEC has proven to be challenging due to the current performance levels of qubits -- break-even requires fabrication and operation quality that is beyond the state-of-the-art. Understanding how qubit performance factors into the success of a QEC code is a valuable exercise for tracking progress towards fault-tolerant quantum computing. Here we present \texttt{qec\_code\_sim}, an open-source, lightweight Python framework for studying the performance of small quantum error correcting codes under the influence of a realistic error model appropriate for superconducting transmon qubits, with the goal of enabling useful hardware studies and experiments. \texttt{qec\_code\_sim} requires minimal software dependencies and prioritizes ease of use, ease of change, and pedagogy over execution speed. As such, it is a tool well-suited to small teams studying systems on the order of one dozen qubits.

Lopez, Santiago↗

DUNE-PRISM – A New Method to Measure Neutrino Oscillations

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long baseline neutrino oscillation experiment designed to make precision measurements in a 1.2 2.4 MW neutrino beam, which is directed 1285 km from the Fermi National Accelerator Laboratory (Fermilab) to the Sandford Underground Research Facility (SURF) in South Dakota. Neutrinos are measured at two detector facilities: a near detector located at Fermilab close where the beam is produced and a far detector at SURF. The neutrino beam can be configured to be composed primarily of either muon or anti-muon neutrinos. DUNE measures the disappearance of muon and anti-muon neutrinos and appearance electron and anti-electron neutrinos in the neutrino beam. Measuring these neutrino flavour transitions provides DUNE with sensitivity to the neutrino mass ordering, $\delta_{CP}$, $\theta_{13}$, $\theta_{23}$ and the magnitude of $\Delta m^2_{23}$. The DUNE Precision Reaction Independent Spectrum Measurement (DUNE-PRISM) concept presents a novel way to perform a neutrino oscillation analysis, which has the potential to significantly reduce the impact of large systematic uncertainties in the neutrino interaction model. The PRISM method linearly combines measurements of off-axis neutrino interactions at the DUNE near detector to produce data-driven predictions of the oscillated neutrino event rate spectrum at the far detector. By building an oscillated far detector prediction directly from data, any unknown or poorly modelled neutrino interaction effects will be naturally incorporated into the measurement of the parameters of the neutrino oscillation model. This thesis presents the first complete neutrino oscillation analysis for DUNE using the PRISM method. Details of the methodology are fully explained and the prospects for further improvements to the techniques described are highlighted. The expected impact and relative importance of the neutrino flux, cross section and detector systematic uncertainties are described in detail. Finally, this thesis demonstrates that the PRISM method is capable of performing a measurement of the oscillation parameters that is robust against neutrino interaction modelling errors.

Hasnip, Ciaran↗

diffReplication - An Energy-Aware Fault Tolerance Model for Silent Error Detection and Mitigation in Heterogeneous Extreme-scale Computing Environment

At extreme scale, the frequency of silent errors – a class of errors that remain undetected by low-level error detection mechanisms – increases significantly with the computational complexity of the application and the scale of the computing infrastructure. As hardware and software advances are made to usher in the next scientific era of computing, developing new approaches to mitigate the impact of silent errors remains a challenging problem. In this work, we propose an energy-aware fault-tolerance model, referred to diffReplication to overcome silent errors. In the proposed model, the main process is associated with one replica that executes at the same rate as the main process, and one diffReplica that is executed at a fraction of the main process' execution rate. If the main and its replica reach consensus at the end of a computation phase, the state of the diffReplica is updated and computation is resumed. If the synchronization attempt results in a disagreement, however, the diffReplica increases its execution speed to complete the computation and quickly reach the synchronization barrier. Assuming a single error over any given synchronization interval, a majority voting is used to reach consensus and tolerate silent errors. To further enhance its performance, diffReplication is augmented with speculative execution, whereby the main or its fast replica is selected to continue execution without waiting for the diffReplica. The selection process is based on the previous behaviour of the main and its replica. A performance analysis study is carried out to assess the performance of diffReplication, in terms of the energy saving and time-to-completion reduction achieved by the diffReplication scheme. The experiment shows that speculative execution reduces the time to completion with additional energy, and dynamic decision-making balances the energy consumption and time to completion.

97 MATHEMATICS AND COMPUTING↗

Assessment of errors in analytic modeling of permanent magnet electron spectrometers for laser-plasma accelerators

Here, we assess the error in three treatments of a critical electron beam diagnostic for laser-plasma acceleration (LPA) experiments: a permanent magnet electron spectrometer. Since LPA electron spectrometers are often difficult to calibrate due to the scarcity of well-characterized, tunable electron beam sources in the appropriate energy range and the mechanical complexity of electron spectrometers, the standard of calibration is Hall probe measurements. We first compare the electron spectrometer performance between SIMION calculations and the Hall probe measurements. We find up to a 60% (<8%) error in determining the absolute energy for electrons below (above) 80 MeV when using the SIMION-modeled magnetic field vs the measured magnetic field due to SIMION overpredicting fringe field strength. The difference in spectrometer energy resolution is within ∼10% between the modeled and measured fields. We then assess a simple block model, which is commonly used. We find that the simple block model presented here sufficiently captures the predicted as-built spectrometer performance for any application provided that the user can tolerate uncertainties on absolute energy determination up to 5% and errors on energy resolution up to 1%.

Electromagnetic optics↗

Multi-phenomenology Yield Characterization

This report serves as the first delivery of a four-year applied science effort to transform and advance the error bounds for the yield estimate of an explosion. Each year’s delivery will be in this form, culminating in the submission of this work for peer review to a scientific journal. Importantly, the yearly progress reports can then also be viewed as expanding drafts working towards a formal journal article submission. For the first tranche of funding, we collaborated with Air Force Technical Applications Center (AFTAC) scientists to identify unclassified real-world data that demonstrate and validate our advanced error propagation methods. Collaboration includes visits to AFTAC and telecons. For this development, we illustrate the fusion of seismic, acoustic, optical, and surface effect signatures from an explosion. The mathematics and code being adapted to this specific application (Williams et al., 2021) involves physics models of multiple sensor signatures. We have also identified related physics models and have integrated them into code. Current methods of underground explosion yield estimation for the Threshold Test Ban Treaty (TTBT) have served the US treaty monitoring mission well for decades. A research objective of the Defense Nuclear Nonproliferation Research and Development (DNN R&D) office of the National Nuclear Security Administration (NNSA) has always been to provide new technical capabilities for monitoring lower thresholds. The general model and error propagation code to be developed in this project is based on significant advances in error modeling and propagation needed to analyze data at lower detection thresholds. The second tranche of funding for this project began on May 1, 2022, and planned work for the second tranche includes: i) completing the integration of physical model code into the general error model framework; this code accommodates a wide range of linear/nonlinear source models, fixed/ random effects, and frequentist/Bayesian analyses (the purpose of which is not to dictate to users how to analyze data, but instead to allow users the maximum flexibility in their work); ii) illustrative application of code to identified data, and; iii) initial planning with AFTAC researchers on delivery of code to the Common Development Environment at AFTAC, and continued writing of the planned final journal article submission (year two of this progress report), with particular emphasis on descriptions of data identified for this effort.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗

Model validation and error attribution for a drifting qubit

Qubit performance is often reported in terms of a variety of single-value metrics, each providing a facet of the underlying noise mechanism limiting performance. However, the value of these metrics may drift over long timescales, and reporting a single number for qubit performance fails to account for the low-frequency noise processes that give rise to this drift. Here, in this work, we demonstrate how we can use the distribution of these values to validate or invalidate candidate noise models. We focus on the case of randomized benchmarking (RB), where typically a single error rate is reported but this error rate can drift over time when multiple passes of RB are performed. We show that using a statistical test as simple as the Kolmogorov-Smirnov statistic on the distribution of RB error rates can be used to rule out noise models, assuming the experiment is performed over a long enough time interval to capture relevant low frequency noise. With confidence in a noise model, we show how care must be exercised when performing error attribution using the distribution of drifting RB error rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗