SEARCH · Search NASA

Results for “error bars”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Nuclear data uncertainty and sensitivity analysis of the VHTRC benchmark using SCALE

The Very High Temperature Critical Assembly (VHTRC) experiment represents one of the few data sets available for the validation of HTGR lattice physics. Within the framework of the IAEA Coordinated Research Project on HTGR Uncertainty Analysis in Modeling, uncertainty and sensitivity analyses of this graphite-moderated facility are performed as the validation reference to the prismatic MHTGR-350 fuel block calculations. Nominal multi-group and continuous-energy criticality calculations with the KENO-VI Monte Carlo code of the SCALE code package are compared with the continuous-energy Monte Carlo Code Serpent. Good agreement with a maximum difference of 250 pcm is obtained. The experimental data set, however, differs by several hundred pcm with the results of both codes when using the ENDF-VII.0 library. When using ENDF/B-VII.1 data, this difference is reduced to a few hundred pcm and the calculations lie within the experimental error bars. Uncertainties of the VHTRC multiplication factors due to uncertainties in nuclear data are determined with SAMPLER/KENO-VI of SCALE 6.2 and with continuous-energy TSUNAMI of SCALE 6.2. For all experimental configurations, the obtained uncertainty is found to be 0.58% when using ENDF/B-VII.0 data. The top contributor to this uncertainty is the average number of neutrons per fission event of U-235. With ENDF/B-VII.1 data, uncertainties of about 0.66% are obtained in particular due to an increased uncertainty of U-235 nu-bar in the latest ENDF release. In conclusion, when considering these nuclear data uncertainties, the obtained 1s uncertainty intervals are overlapping with the experimental error bars.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Validation of SOLPS-ITER simulations with kinetic, fluid, and hybrid neutral models for JET-ILW low-confinement mode plasmas

For JET L-mode plasmas in low-recycling conditions (electron temperature at the outer strike point, T e,ot ≳ 30 eV ), SOLPS-ITER simulations agree within the error bars for the experimental profiles at the low-field side (LFS) divertor target. The peak Balmer-α (Dα ) emission in the LFS divertor agrees within the error bars of the KS3 filterscope diagnostic, but is approximately 30% lower than the peak value of the KT1 spectrometer. Simulations have been performed with fluid, kinetic, and hybrid models for the neutrals. The large fluid-kinetic discrepancies of more than a factor 2 are successfully corrected by using a hybrid fluid-kinetic approach, for which kinetic atoms are transferred to the fluid population when the local Knudsen number of the atom becomes smaller than a user-defined transition Knudsen number Kn t . The hybrid-kinetic discrepancies are limited to a few % for Kn t ≤ 100 . When increasing the upstream density to high-recycling conditions, at the onset of detachment ( T e,ot ≈ 5 eV ), the simulations predict more than a factor 2 lower peak ion saturation current to the LFS divertor than the experiments. Also the Dα emission is underpredicted with approximately a factor 2. For these high-recycling conditions, the fluid-kinetic discrepancies are limited to maximum 50%, which are again corrected by using the hybrid approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗

DOE OSTI2024

The BOSS bispectrum analysis at one loop from the Effective Field Theory of Large-Scale Structure

Abstract We analyze the BOSS power spectrum monopole and quadrupole, and the bispectrum monopole and quadrupole data, using the predictions from the Effective Field Theory of Large-Scale Structure (EFTofLSS). Specifically, we use the one loop prediction for the power spectrum and the bispectrum monopole, and the tree level for the bispectrum quadrupole. After validating our pipeline against numerical simulations as well as checking for several internal consistencies, we apply it to the observational data. We find that analyzing the bispectrum monopole to higher wavenumbers thanks to the one-loop prediction, as well as the addition of the tree-level quadrupole, significantly reduces the error bars with respect to our original analysis of the power spectrum at one loop and bispectrum monopole at tree level. After fixing the spectral tilt to Planck preferred value and using a Big Bang Nucleosynthesis prior, we measureσ 8 = 0.794 ± 0.037,h = 0.692 ± 0.011, and Ω m = 0.311 ± 0.010 to about 4.7%, 1.6%, and 3.2%, at 68% CL, respectively.This represents an error bar reduction with respect to the power spectrum-only analysis of about 30%, 18%, and 13% respectively. Remarkably, the results are compatible with the ones obtained with a power-spectrum-only analysis, showing the power of the EFTofLSS in simultaneously predicting several observables. We find no tension with Planck.

Astronomy & Astrophysics↗

DOE OSTI2025

Measuring and unbiasing the BAO shift in the Ly α forest with AbacusSummit

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) places sub- per cent constraints on measurements of the Baryon Acoustic Oscillation (BAO) scaling parameters from the Ly $\alpha$ forest. However, no systematic error budget stemming from non-linearities in the three-dimensional clustering of the Ly $\alpha$ forest is included in the DESI-Ly $\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly $\alpha$ forest mocks produced on the N-body simulation suite AbacusSummit, which adopt the Fluctuating–Gunn–Peterson Approximation (FGPA). Specifically, we measure the Ly $\alpha$ autocorrelation and the Ly $\alpha$-quasar cross-correlation functions. To mitigate the noise, we adopt a linear control variates technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the autocorrelation, we detect a small positive shift in radial direction of $\Delta \alpha _{\parallel }= 0.35~{{\ \rm per\ cent}}$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha _\perp$. From the cross-correlation, we see a similar shift to $\Delta \alpha _\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha _{\perp }=\sim$0.25 per cent, at the 2$\sigma$ level. We also make a connection with the Ly $\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly $\alpha$ auto- and the Ly $\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly $\alpha$ forest analyses based on EFT. This work paves the way for novel full-shape analyses of the currently observing DESI and future surveys such as the PFS, WEAVE-QSO, and 4MOST.

Hadzhiyska, Boryana↗

DOE OSTI2017

Superiority in Robustness of Multifield Optimization Over Single-Field Optimization for Pencil-Beam Proton Therapy for Oropharynx Carcinoma: An Enhanced Robustness Analysis

Purpose: To compare the difference in robustness of single-field optimized (SFO) and robust multifield optimized (rMFO) proton plans for oropharynx carcinoma patients by an improved robustness analysis. Methods and Materials: We generated rMFO proton plans for 11 patients with oropharynx carcinoma treated with SFO intensity modulated proton therapy with simultaneous integrated boost prescription. Doses from both planning approaches were compared for the initial plans and the worst cases from 20 optimization scenarios of setup errors and range uncertainties. Expected average dose distributions per range uncertainty were obtained by weighting the contributions from the respective scenarios with their expected setup error probability, and the spread of dose parameters for different range uncertainties were quantified. Using boundary dose distributions created from 56 combined setup error and range uncertainty scenarios and considering the vanishing influence of setup errors after 30 fractions, we approximated realistic worst-case values for the total treatment course. Error bar metrics derived from these boundary doses are reported for the clinical target volumes (CTVs) and organs at risk (OARs). Results: The rMFO plans showed improved CTV coverage and homogeneity while simultaneously reducing the average mean dose to the constrictor muscles, larynx, and ipsilateral middle ear by 5.6 Gy, 2.0 Gy, and 3.9 Gy, respectively. We observed slightly larger differences during robustness evaluation, as well as a significantly higher average brainstem maximum and ipsilateral parotid mean dose for SFO plans. For rMFO plans, the range uncertainty-related spread in OAR dose parameters and many error bar metrics were found to be superior. The SFO plans showed a lower global maximum dose for single-scenario worst cases and a slightly lower mean oral cavity dose throughout. Conclusions: An enhanced robustness analysis has been proposed and implemented into clinical systems. The benefit of better CTV coverage and OAR dose sparing in oropharynx carcinoma patients by rMFO compared with SFO proton plans is preserved in a robustness analysis with consideration of setup error and range uncertainty.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2023

Current data are consistent with flat spatial hypersurfaces in the Λ CDM cosmological model but favor more lensing than the model predicts

Here, we study the performance of three pairs of tilted, and a pair of untilited, ΛCDM cosmological models, with three of these four pairs allowing for non-flat spatial hypersurfaces, against cosmic microwave background (CMB) temperature and polarization power spectrum data (P18), measurements of the Planck 2018 lensing potential power spectrum (lensing), and a large compilation of non-CMB data (non-CMB). For the eight models, we measure cosmological parameters and study whether or not pairs of the data sets (as well as subsets of them) are mutually consistent in these models. Half of these models allow the lensing consistency parameter A L , which re-scales the gravitational potential power spectrum, to be an additional free parameter to be determined from data, while the other three have A L = 1 which is the theoretically expected value. The pair of untilted non-flat ΛCDM models are incompatible with P18 data. The tilted spatially-flat models assume the usual primordial spatial inhomogeneity power spectrum that is a power law in wave number. The tilted non-flat models assume either the primordial power spectrum used in the Planck group anal yses [Planck P(q)], that has recently been numerically shown to be a good approximation to what is quantum-mechanically generated from a particular choice of closed inflation model initial conditions, or a recently computed power spectrum [new P(q)] that quantum-mechanically follows from a different set of non-flat inflation model initial conditions. In the tilted non-flat models with A L = 1 we find differences between P18 data and non-CMB data cosmological parameter constraints, which are large enough to rule out the Planck P(q) model at 3σ but not the new P(q) model. No significant differences are found when cosmological parameter constraints obtained with two different data sets are compared within the standard tilted flat ΛCDM model. While both P18 data and non-CMB data separately favor a closed geometry, with spatial curvature density parameter Ω k < 0, when P18+non-CMB data are jointly analyzed the evidence in favor of non-flat hypersurfaces subsides. Differences between P18 data and non-CMB data cosmological constraints subside when A L is allowed to vary. From the most restrictive P18+lensing+non-CMB data combination we get almost model-independent constraints on the cosmological parameters and find that the A L > 1 option is preferred over the Ω k < 0 one, with the A L parameter, for all models, being larger than unity by ~ 2.5σ. According to the deviance information criterion, in the P18+lensing+non-CMB analysis, the varying A L option is on the verge of being strongly favored over the A L = 1 one, which could indicate a problem for the standard tilted flat ΛCDM model. These data are consistent with flat spatial hypersurfaces but more and better data could improve the constraints on Ω k and might alter this conclusion. Error bars on some cosmological parameters are significantly reduced when non-CMB data are used jointly with P18+lensing data. For example, in the tilted flat ΛCDM model for P18+lensing+non-CMB data the Hubble constant H 0 = 68.09 ± 0.38 km s -1 Mpc -1 , which is consistent with that from a median statistics analysis of a large compilation of H 0 measurements as well as with a number of local measurements of the cosmological expansion rate. This H 0 error bar is 31% smaller than that from P18+lensing data alone.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Exploring experimental conditions to reduce uncertainties in the optical potential

Uncertainty quantification for nuclear theories has gained a more prominent role in the field, with more groups attempting to understand the uncertainties in their calculations. However, recent studies have shown that the uncertainties on the optical potentials are too large for the theory to be useful. The goal of this work is to explore possible experimental conditions that may reduce the uncertainties on elastic scattering and single-nucleon transfer cross sections that come from the fitting of the optical model parameters to experimental data. Using Bayesian methods, we explore the effect of the uncertainties of optical model parameters on the angular grid of the differential cross section, including cross-section data at nearby energies, and changes in the experimental error bars. We also study the effect on the resulting uncertainty when other observables are included in the fitting procedure, particularly the total (reaction) cross sections. We study proton and neutron elastic scattering on 48 Ca and 208 Pb. We explore the parameter space with the Markov-chain Monte Carlo method, produce posterior distributions for the optical model parameters, and construct the corresponding 95% confidence intervals on the elastic-scattering cross sections. We also propagate the uncertainties on the optical potentials to the 48 Ca(d,p) 49 Ca(g.s.) and 208 Pb(d,p) 209 Pb(g.s.) cross sections. We find little sensitivity to the angular grid and an improvement of up to a factor of 2 on the uncertainties by including data at a nearby energy. Although reducing the error bars on the data does reduce the uncertainty, the gain is often considerably smaller than one would naively expect. We also find that the inclusion of total reaction cross section can improve the uncertainty although the magnitude of the effect depends strongly on the cases considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

S -factor and scattering-parameter extractions from ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}{ \rightarrow }^{7}\mathrm{Be}+\gamma $

Previous studies of the reaction ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ have mainly focused on providing the best central value and error bar for the S factor at solar energies. Experimental measurements of this capture reaction at higher energies, the ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering phase shifts, as well as properties of ${}^{7}\mathrm{Be}$ and its excited state, have been used to constrain the theoretical models employed for this purpose. Here we show that much more information than was previously appreciated can be extracted from angle-integrated capture data alone. We use the next-to-leading-order (NLO) amplitude in an effective field theory (EFT) for ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ to perform the extrapolation. At this order the EFT describes the capture process using an s-wave scattering length and effective range, the asymptotic properties of ${}^{7}\mathrm{Be}$ and its excited state, and short-distance contributions to the E1 capture amplitude. We extract the multi-dimensional posterior of all these parameters via a Bayesian analysis that uses capture data below 2 MeV. We find that properties of the ${}^{7}\mathrm{Be}$ ground and excited states are well constrained. The total S factor $S(0)\,={0.577}_{-0.016}^{+0.015}$ keV b, while the branching ratio for excited- to ground-state capture at zero energy, ${Br}(0)={0.406}_{-0.011}^{+0.013}$, both at 68% degree of belief. This S(0) is broadly consistent with other recent evaluations, and agrees with the previously recommended value $S(0)=0.56\pm 0.03\,\,\mathrm{keV}$ b, but has a smaller error bar. We also find significant constraints on ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering parameters, and we obtain constraints on the angular distribution of capture gamma rays, which is important for interpreting experiments. The path forward for this reaction seems to lie with better measurements of the scattering phase shifts and S(E)'s angular dependence away from zero energy, together with better understanding of the asymptotic normalization coefficients of the ${}^{7}$Be bound states' wave functions. Lastly, data on these could further reduce the uncertainty on S(0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

The development of explosions in axisymmetric ab initio core-collapse supernova simulations of 12–25 M{sub ⊙} STARS

We present four ab initio axisymmetric core-collapse supernova simulations initiated from 12, 15, 20, and 25 M{sub ⊙} zero-age main sequence progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 s after core bounce and 1 s after material first becomes unbound. These simulations were computed with our Chimera code employing RbR spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 s after core bounce allows the explosions to develop more fully and the processes involved in powering the explosions to become more clearly evident. We compute explosion energy estimates, including the negative gravitational binding energy of the stellar envelope outside the expanding shock, of 0.34, 0.88, 0.38, and 0.70 Bethe (B ≡ 10{sup 51} erg) and increasing at 0.03, 0.15, 0.19, and 0.52 B s{sup −1}, respectively, for the 12, 15, 20, and 25 M{sub ⊙} models at the endpoint of this report. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We discuss how the explosion energies may be subject to stochastic variations as exemplfied by the effect of the explosion geometry of the 20 M{sub ⊙} model in reducing its explosion energy. We compute the proto-neutron star masses and kick velocities. We compare our results for the explosion energies and ejected {sup 56}Ni masses against some observational standards despite the large error bars in both models and observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

( Hubble Space Telescope Hubble Space Telescope) proper motion (HSTPROMO) catalogs of galactic globular clusters. IV. Kinematic profiles and average masses of blue straggler stars

We make use of the Hubble Space Telescope proper-motion catalogs derived by Bellini et al. to produce the first radial velocity dispersion profiles σ(R) for blue straggler stars (BSSs) in Galactic globular clusters (GCs), as well as the first dynamical estimates for the average mass of the entire BSS population. We show that BSSs typically have lower velocity dispersions than stars with mass equal to the main-sequence turnoff mass, as one would expect for a more massive population of stars. Since GCs are expected to experience some degree of energy equipartition, we use the relation σ∝M{sup −η}, where η is related to the degree of energy equipartition, along with our velocity dispersion profiles to estimate BSS masses. We estimate η as a function of cluster relaxation from recent Monte Carlo cluster simulations by Bianchini et al. and then derive an average mass ratio M{sub BSS}/M{sub MSTO}=1.50±0.14 and an average mass M{sub BSS}=1.22±0.12 M {sub ⊙} from 598 BSSs across 19 GCs. The final error bars include any systematic errors that are random between different clusters, but not any potential biases inherent to our methodology. Our results are in good agreement with the average mass of M{sub BSS}=1.22±0.06 M {sub ⊙} for the 35 BSSs in Galactic GCs in the literature with properties that have allowed individual mass determination.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Shock compression response of highly reactive Ni + Al multilayered thin foils

The shock-compression response of Ni + Al multilayered thin foils is investigated using laser-accelerated thin-foil plate-impact experiments over the pressure range of 2 to 11 GPa. The foils contain alternating Ni and Al layers (parallel but not flat) of nominally 50 nm bilayer spacing. The goal is to determine the equation of state and shock-induced reactivity of these highly reactive fully dense thin-foil materials. The laser-accelerated thin-foil impact set-up involved combined use of photon-doppler-velocimetry to monitor the acceleration and impact velocity of an aluminum flyer, and VISAR interferometry was used to monitor the back free-surface velocity of the impacted Ni + Al multilayered target. The shock-compression response of the Ni + Al target foils was determined using experimentally measured parameters and impedance matching approach, with error bars identified considering systematic and experimental errors. Meso-scale CTH shock simulations were performed using real imported microstructures of the cross-sections of the multilayered Ni + Al foils to compute the Hugoniot response (assuming no reaction) for correlation with their experimentally determined equation of state. It was observed that at particle velocities below ∼150 m/s, the experimentally determined equation of state trend matches the CTH-predicted inert response and is consistent with the observed unreacted state of the recovered Ni + Al target foils from this velocity regime. At higher particle velocities, the experimentally determined equation of state deviates from the CTH-predicted inert response. A complete and self-sustained reaction is also seen in targets recovered from experiments performed at these higher particle velocities. The deviation in the measured equation of state, to higher shock speeds and expanded volumes, combined with the observation of complete reaction in the recovered multilayered foils, confirmed via microstructure characterization, is indicative of the occurrence of shock-induced chemical reaction occurring in the time-scale of the high-pressure state. TEM characterization of recovered shock-compressed (unreacted) Ni + Al multilayered foils exhibits distinct features of constituent mixing revealing jetted layers and inter-mixed regions. These features were primarily observed in the proximity of the undulations present in the alternating layers of the Ni + Al starting foils, suggesting the important role of such instabilities in promoting shock-induced intermetallic-forming reactions in the fully dense highly exothermic multilayered thin foils.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Applying the noiseless extrapolation error mitigation protocol to calculate real-time quantum field theory scattering phase shifts

Not Available

Parks, Zachary↗

DOE OSTI2022

Numerical error analysis of SOLPS-ITER simulations of EAST

Abstract Plasma edge simulations with codes like SOLPS-ITER are widely employed to interpret fusion experiments. However, numerical errors appearing in such simulations are rarely investigated, despite their potential large impact on simulation results. These errors consist of the statistical error and the bias, both resulting from the finite number of employed EIRENE Monte Carlo particles and incomplete convergence, and the discretization error due to the finite resolution of the computational grids. In this contribution, the resulting numerical errors on simulations of pure deuterium and neon seeded H-mode EAST discharges are examined. The statistical error can be kept small compared to other numerical error contributions by averaging the plasma profiles. This allows investigating the bias and discretization errors using Richardson extrapolation. It is shown that grid refinement and the number of employed Monte Carlo particles have the largest influence on the result, in agreement with similar studies of an ITER deuterium case. For the first time, numerical error bars on the entire simulated target profiles are determined showing that the largest numerical error is 17.9%, mainly due to the plasma grid discretization. On top, also numerical errors on simulated neutral pressures are investigated in detail, for which the statistical error is dominant. The analysis demonstrates which setup is needed to keep numerical errors limited: the SOLPS-ITER averaging procedure should be employed including enough EIRENE particles, and the involved grids should be sufficiently refined to reduce discretization errors.

Boeyaert, Dieter (ORCID:0000000309208660)↗

DOE OSTI2023

Count Every Trip: Finding the Uncertainty in Energy Estimates Made from Inferred Travel Modes

To properly inform transport policy and infrastructure changes, transportation related metrics need both measured values and uncertainties of those values. Travel monitoring smartphone apps can record people's travel behavior, but trip data quality is limited by sensor errors, user labeling rates and the accuracy of inference algorithms used for travel diary creation. We discuss the use of phone app recorded travel diary data to estimate energy consumption, and propose the use of propagation of variance to find error bars for such estimates. We define energy consumption for one trip as trip length times the energy intensity per distance unit of the travel mode used. We characterize trip length errors with relative error and inferred trip mode errors with confusion matrix columns. The resulting variances of each measurement are then propagated to the final calculated energy consumption. We tested our uncertainty methods on a dataset that used phone app data combined with prompted recall, consisting of 92,234 labeled trips for over 500,000 miles. Accounting for uncertainty using expected energy intensities and variance propagation gives a dataset-wide aggregate energy consumption percent error of about 9%, within one standard deviation from the truth. Future work could involve applying similar methods to other travel diary based metrics.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗

DOE OSTI2023

Count Every Trip: Finding the Uncertainty in Energy Estimates Made from Inferred Travel Modes

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Estimating Travel Energy Consumption Uncertainty Based on Inferred Travel Mode and Sensed Travel Length

To properly inform transport policy and infrastructure changes, transportation related metrics need both measured values and uncertainties of those values. Travel monitoring smartphone apps can record people's travel behavior, but trip data quality is limited by sensor errors, user labeling rates and the accuracy of inference algorithms used for travel diary creation. We discuss the use of phone app recorded travel diary data to estimate energy consumption, and propose the use of propagation of variance to find error bars for such estimates. We define energy consumption for one trip as trip length times the energy intensity per distance unit of the travel mode used. We characterize trip length errors with relative error and inferred trip mode errors with confusion matrix columns. The resulting variances of each measurement are then propagated to the final calculated energy consumption. We tested our uncertainty methods on a dataset that used phone app data combined with prompted recall, consisting of 92,234 labeled trips for over 500,000 miles. Accounting for uncertainty using expected energy intensities and variance propagation gives a dataset-wide aggregate energy consumption percent error of about 8%, within one standard deviation from the truth. Future work could involve applying similar methods to other travel diary based metrics.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗