Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Analysis”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-dimensional control co-design of a wave energy converter with a novel pitch resonator power takeoff system

Researchers are exploring adding wave energy converters to existing oceanographic buoys to provide a predictable source of renewable power. A ”pitch resonator” power take-off system has been developed that generates power using a geared flywheel system designed to match resonance with the pitching motion of the buoy. However, the novelty of the concept leaves researchers uncertain about various design aspects of the system. This work presents a novel design study of a pitch resonator to inform design decisions for an upcoming deployment of the system. The assessment uses control co-design via WecOptTool to optimize control trajectories for maximal electrical power production while varying five design parameters of the pitch resonator. Given the large search space of the problem, the control trajectories are optimized within a Monte Carlo analysis to identify optimal designs, followed by parameter sweeps around the optimum to identify trends between the design parameters. The gear ratio between the pitch resonator spring and flywheel are found to be the most sensitive design variables to power performance. Finally, the assessment also finds similar power generation for various sizes of resonator components, suggesting that correctly designing for optimal control trajectories at resonance is more critical to the design than component sizing.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Planar XY magnetic glass state in the Gd 2 ScNbO 7 pyrochlore

Here a spin glass system with emergent planar ordered spin clusters is investigated. The mixed B-site pyrochlore Gd 2 ScNbO 7 has been synthesized and characterized through a variety of techniques, including x-ray diffraction, magnetic susceptibility, muon spin relaxation, heat capacity and neutron scattering. Despite a Curie–Weiss temperature of -3.93(3) K, indicating net antiferromagnetic interactions, no signs of long ranged magnetic ordering are found down to T = 0.3 K. Instead, a disordered magnetic state emerges with a small correlation length of 2.1(1) Å of single tetrahedra. A Reverse Monte Carlo analysis of the polarized neutron scattering data reveals short-range antiferromagnetic order with emergent XY spin ordering similar to the parent pyrochlore compounds. Muon spin relaxation, and AC susceptibility measurements confirm that the magnetization condenses into a glass, with 10 % of the potential entropy missing in the specific heat. This magnetic ground state is similar to what is observed in Gd 2 Sn 2 O 7 just above the ordering temperature, without the eventual long-range ordering at low temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Improved modeling of β electronic recoils in liquid xenon using LUX calibration data

We report here methods and techniques for creating and improving a model that reproduces the scintillation and ionization response of a dual-phase liquid and gaseous xenon time-projection chamber. Starting with the recent release of the Noble Element Simulation Technique (NEST v2.0), electronic recoil data from the $\beta$ decays of ${}^3$H and ${}^{14}$C in the Large Underground Xenon (LUX) detector were used to tune the model, in addition to external data sets that allow for extrapolation beyond the LUX data-taking conditions. This paper also presents techniques used for modeling complicated temporal and spatial detector pathologies that can adversely affect data using a simplified model framework. The methods outlined in this report show an example of the robust applications possible with NEST v2.0, while also providing the final electronic recoil model and detector parameters that will used in the new analysis package, the LUX Legacy Analysis Monte Carlo Application (LLAMA), for accurate reproduction of the LUX data. As accurate background reproduction is crucial for the success of rare-event searches, such as dark matter direct detection experiments, the techniques outlined here can be used in other single-phase and dual-phase xenon detectors to assist with accurate ER background reproduction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Quantifying airborne fraction trends and the destination of anthropogenic CO 2 by tracking carbon flows in a simple climate model

Abstract Atmospheric carbon dioxide (CO 2 ) concentrations have increased as a direct result of human activity and are at their highest level over the last 2 million years, with profound impacts on the Earth system. However, the magnitude and future dynamics of land and ocean carbon sinks are not well understood; therefore, the amount of anthropogenic fossil fuel emissions that remain in the atmosphere (the airborne fraction) is poorly constrained. This work aims to quantify the sources and controls of atmospheric CO 2 , the fate of anthropogenic CO 2 over time, and the likelihood of a trend in the airborne fraction. We use Hector v3.0, a coupled simple climate and carbon cycle model with the novel ability to explicitly track carbon as it flows through the Earth system. We use key model parameters in a Monte Carlo analysis of 15 000 model runs from 1750 to 2300. Results are filtered for physical realism against historical observations and CMIP6 projection data, and we calculate the relative importance of parameters controlling how much anthropogenic carbon ends up in the atmosphere. Modeled airborne fraction was roughly 52%, consistent with observational studies. The overwhelming majority of model runs exhibited a negative trend in the airborne fraction from 1960–2020, implying that current-day land and ocean sinks are proportionally taking up more carbon than the atmosphere. However, the percentage of atmospheric CO 2 derived from anthropogenic origins can be much higher because of Earth system feedbacks. We find it peaks at over 90% between 2010–2050. Moreover, when looking at the destination of anthropogenic fossil fuel emissions, only a quarter ends up in the atmosphere while more than half of emissions are taken up by the land sink on centennial timescales. This study evaluates the likelihood of airborne fraction trends and provides insights into the dynamics of anthropogenic CO 2 in the Earth system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Oscillating scalar fields and the Hubble tension: A resolution with novel signatures

Here, we present a detailed investigation of a subdominant oscillating scalar field [“early dark energy” (EDE)] in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until log 10 ($z_c$) ~ 3.5, reaching $ρ_{\text{EDE}}(z_c)/ρ_{\text{tot}}$ ~ 10 % and diluting faster than matter afterwards, can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES Collaboration into agreement. A scalar field potential that scales as $V (\phi) \propto \phi ^{2n}$ with 2 ≲ $\textit{n}$ ≲ 3.4 around the minimum is preferred at the 68% confidence level, and the Planck polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefer a potential that flattens at large field displacements. A Markov-chain Monte Carlo analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at a very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the E-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $\textit{n}$ ≈ 2. This instability and ensuing potentially nonlinear, spatially inhomogeneous, dynamics may provide unique signatures of this scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Five percent measurement of the gravitational constant in the Large Magellanic Cloud

We perform a novel test of General Relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids’ periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler’s third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be G LMC /G SS = $0.93{^{+0.05}_{-0.04}}$. We discuss the implications of this 5% measurement of Newton’s constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers G LMC < G SS at 2.6σ. We also obtain new bounds on critical parameters that appear in semi-analytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $α = 0.90{^{+0.36}_{-0.26}}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Planck constraints and gravitational wave forecasts for primordial black hole dark matter seeded by multifield inflation

Here, we perform a Markov Chain Monte Carlo analysis of a simple yet generic multifield inflation model characterized by two scalar fields coupled to each other and nonminimally coupled to gravity, fit to Planck 2018 cosmic microwave background (CMB) data. In particular, model parameters are constrained by data on the amplitude of the primordial power spectrum of scalar curvature perturbations on CMB scales A s , the spectral index n s , and the ratio of power in tensor to scalar modes r, with a prior that the primordial power spectrum should also lead to primordial black hole (PBH) production sufficient to account for the observed dark matter abundance. We find that n s in particular controls the constraints on our model. Whereas previous studies of PBH formation from an ultraslow-roll phase of inflation have highlighted the need for at least one model parameter to be highly fine-tuned, we identify a degeneracy direction in parameter space such that shifts by ~10% of one parameter can be compensated by comparable shifts in other parameters while preserving a close fit between model predictions and observations. Furthermore, we find this allowed parameter region produces observable gravitational wave signals in the frequency ranges to which upcoming experiments are projected to be sensitive, including Advanced LIGO and Virgo, the Einstein Telescope, Cosmic Explorer, DECIGO, and LISA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Considering the Variability of Soiling in Long-Term PV Performance Forecasting

This study presents the development of a methodology for evaluating the variability associated with soiling on long-term photovoltaic (PV) forecasting. Independent engineering firms typically build forecasts for large PV plants through the use of the PVsyst software, where monthly soiling losses are one of many inputs to the P50 model. Subsequently, long-term performance distributions are constructed through a Monte Carlo analysis that includes various factors, such as satellite irradiance modeling uncertainty, uncertainty in the PVsyst model, and long-term irradiance variability. Often the PVsyst model uncertainty is increased to account for sites with significant soiling concerns but no systematic method has been presented in the literature to specifically include soiling variability within long-term performance uncertainty. In this work soiling information from 16 sites in the U.S. Southwest are combined with 24 years of rainfall data to generate 24 years of energy production with soiling losses and then subsequently generate probability of exceedance values (e.g., P50, P90, P95…). The results show that the size of the 90% confidence interval (P5–P95) can increase from –0.7% to 10.1% when interannual soiling variability and soiling rate uncertainty is included.

14 SOLAR ENERGY↗
DOE OSTI2025

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Monte-Carlo Method for Evaluating the Economic Performance of Plastics Recycling Systems Using Historical Pricing

The management of plastic waste is a considerable emerging global concern. Conventional plastics recycling rates range from 5–10% in the US, which results in considerable material being sent to landfill. New sorting technologies using artificial intelligence have the potential to dramatically increase the ability of waste management companies to sort out valuable fractions of plastic waste and create high-purity streams for secondary markets. However, several fractions of the plastics recycling stream, i.e., low-density polyethylene, polypropylene, and polystyrene, do not have well-developed markets. These materials may be sent to pyrolysis facilities for thermal recycling, but the economics of the pyrolysis industry are currently uncertain. This work aims to identify the breakeven price of the plastic fractions that would be sent to pyrolysis in order to determine the economic viability of the sorting facility. The work will use Monte Carlo analysis to evaluate several scenarios and understand how compositional variation, price variability, and facility attributes affect the breakeven price. Here, the presorting conditions in Materials Recovery Facilities strongly affects the profitability of the sorting facility. High residual high density polyethylene (HDPE) and polyethylene terephthalate (PET) content in presorted bales are important revenue drivers for the facility, even in presorted 3-7 bales. Key variables that lead to profitability include a high sale price of polyethylene terephthalate, a high sale price of high-density polyethylene, and if polypropylene can be sold to recycling markets (as opposed to pyrolysis facilities).

Plastics↗
DOE OSTI2020

HEPfit: a code for the combination of indirect and direct constraints on high energy physics models

HEPfit is a flexible open-source tool which, given the Standard Model or any of its extensions, allows to (i) fit the model parameters to a given set of experimental observables; (ii) obtain predictions for observables. HEPfit can be used either in Monte Carlo mode, to perform a Bayesian Markov Chain Monte Carlo analysis of a given model, or as a library, to obtain predictions of observables for a given point in the parameter space of the model, allowing HEPfit to be used in any statistical framework. In the present version, around a thousand observables have been implemented in the Standard Model and in several new physics scenarios. In this paper, we describe the general structure of the code as well as models and observables implemented in the current release.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

scANN

SAND2025-00656O scANN, also known as sampling by coinflips artificial neural networks, is a software tool that estimates uncertainty in artificial intelligence by performing Monte Carlo analysis on the weight matrices of feedforward neural networks. This computationally intensive process aims to explore the potential value added by future probabilistic hardware. The program’s output helps researchers gain insights into how probabilistic neural networks work. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2020

Flaw Tolerance of Octet Truss Lattices with Random Flaws

We investigate the behavior of octet truss lattices in which some proportion of the bars are randomly removed from a perfect lattice. Using a lattice specimen of 9 × 9 × 9 octet truss unit cells as the nominal design, we simulate the response of prefect and flawed lattices with a finite element method implementation of rod theory. For each of two simple load cases, we perform a Monte Carlo analysis where we consider many instantiations of flawed lattices where between 1% and 20% of the bars are missing, and compare the axial loads in the remaining bars to those in the perfect lattice. We also compare the overall stiffness of the flawed lattices to the perfect lattice. Our results show progressive deterioration of lattice strength and stiffness with increasing proportion of missing bars, as expected, with the location and magnitude of the highest-loaded (most failure-prone) bars depending more on the load case and bar orientation than on the location and connectivity of missing bars. As such, the results reported herein may not be directly transferable to different lattice geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Boiling-Water Reactor Testing Capability in the Advanced Test Reactor

I-Loop is an irradiation facility that is currently being installed at the Advanced Test Reactor. It is a two-loop test facility capable of performing Light Water Reactor (LWR) irradiations in prototypic coolant conditions. The two loops are being installed to be capable of both Boiling Water Reactor (BWR) and Pressurized Water Reactor (PWR) pressure, temperature, and chemistry environments. Each loop is nominally dedicated as a BWR or PWR for simplicity of operations. In-reactor water loop testing that an I-Loop provides is key to the deployment of new accident tolerant fuel technologies and other advanced LWR fuel concepts. Currently, pressurized water loops are the only testing facilities available to test BWR fuel concepts. Their test environments are non-prototypic at higher pressure/temperature and at single-phase fluid flow conditions. This void in the LWR test bed capabilities is one that the I-Loop is uniquely situated to provide. This report discusses the mechanical design, thermal hydraulic calculations, and neutronic calculations of a proposed standard experiment of accident tolerant BWR fuel concepts. Mechanical design examines the geometry and features of the main components. Thermal hydraulic calculations examine the modeling and results of the two-phase flow options available. Lastly, neutronic calculations examine the Monte Carlo analysis of enrichment, heat rates, and flux spectrum.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Considering the Variability of Soiling in Long-Term PV Performance Forecasting: Preprint

This study presents the development of a methodology for evaluating the variability associated with soiling on long-term PV forecasting. Independent engineering firms typically build P50 forecasts for large PV plants through the use of the PVsyst software, where monthly soiling losses are one of many inputs to the P50 model. Subsequently, long-term performance distributions, or Pvalues, are constructed through a Monte Carlo analysis that includes various factors such as: satellite irradiance modeling uncertainty, uncertainty in the PVsyst model, and long-term irradiance variability. Often the PVsyst model uncertainty is increased to account for sites with significant soiling concerns but no systematic method has been presented in the literature to specifically include soiling variability within Pvalues. In this work soiling information from 16 sites in the U.S. Southwest are combined with 20 years of rainfall data to generate 20 years of energy production with soiling losses and then subsequently generate Pvalues. The results show that the spread of Pvalues (P1-P99) can increase from 0-13% when interannual soiling variability is included.

interannual variability↗
DOE OSTI2023

High-Performance Accelerator Modeling: Toward Improving Controls and Diagnostics for High-Brightness Beams in Experiment

One of the main charges of beam physics is to improve the overall beam quality delivered to experiments and applications, which is generally quantified by the beam brightness. High brightness beams have numerous applications, including x-ray free electron lasers and ultrafast electron diffraction (UED). To these ends, this thesis details efforts to use high-performance models --- high-fidelity models that execute quickly --- for controls and diagnostics. To demonstrate the generality of these techniques, this thesis focuses on three of the most successful photoinjector designs currently used worldwide: the UCLA/SLAC/BNL-type high-gradient S-band gun, the continuous-wave high-repetition rate VHF APEX gun and the L-band DESY-PITZ-type gun. Work will be shown from Pegasus (UCLA), HiRES (LBNL) and FAST (FNAL).Specifically, data-driven models for online virtual diagnostics are presented, in this case, in the context of UED at HiRES, leading to a temporal resolution improvement. Methods for improving the fidelity of physics-based models are discussed, with examples at HiRES, Pegasus and FAST. Markov-chain Monte Carlo analysis is applied to match simulations, in the context of photocathode studies at HiRES and Pegasus. Lastly, the augmentation of online model-based predictions with model-independent optimization is explored in a fluctuating environment at Pegasus and HiRES. A central theme of this dissertation is working with parameter fluctuations when modeling an accelerator beamline. Long-term drifts and shot-to-shot jitter exist in every accelerator to a varying degree and therefore play an important role in every chapter of this dissertation. This thesis attempts to address the issues associated with these fluctuations when trying to develop faithful model representations of the system.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

In-situ relative calibration of Bragg crystals with Monte Carlo line ratio analysis

X-ray line emission spectra can thoroughly characterize hot plasmas, especially when line shapes and ratios convey distinct aspects of plasma conditions. However, the high spectral resolution required for observing line shapes is often at odds with the large bandwidth required to observe many line ratios across a wide spectral range. One strategy to obtain high spectral resolution over a wide bandwidth is to use multiple crystals with calibrated reflectivity so that line intensities across different crystals can be compared. Here, we explore the use of a low-resolution, wide-bandwidth mica survey spectrometer to infer relative reflectivity of two high-resolution, narrow-bandwidth quartz crystals. Finally, a Monte Carlo error analysis determines comparable x-ray line ratios measured from both spectrometers, resulting in an in situ calibration factor and associated uncertainty for the relative reflectivity of the high-resolution crystals.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Proton Spin Structure from Simultaneous Monte Carlo Global QCD Analysis

Despite the great effort and achievements made towards understanding proton spin structure in the past few decades, a complete picture is still elusive. Parton distribution functions (PDFs), which in quantum chromodynamics (QCD) encode the momentum and helicity distributions of quarks and gluons inside a proton, provide the means by which to quantify the proton structure information. Being inherently nonperturbative, PDFs have to be extracted from unpolarized and polarized lepton-hadron and hadron-hadron scattering data. In particular, experiments that measure unpolarized and polarized jet observables can provide insight into the momentum and helicity distributions of gluons, which have generally been more difficult to determine reliably than those of quarks. In the past, extraction of the spin-averaged and spin-dependent (or helicity) PDFs has been performed in separate analyses. In this thesis, we perform the first simultaneous extraction of both types of quantities from deep-inelastic scattering (DIS), Drell-Yan and single jet observables, within the Monte Carlo global QCD analysis framework developed by the Jefferson Lab Angular Momentum (JAM) Collaboration. The results from this work indicate that the gluon helicity distributions depend rather strongly on the theory assumptions on which the global analysis is based, which calls for the need of measurements with higher precision. As an application of the new simultaneous JAM analysis, we perform an impact study for future Electron-Ion Collider (EIC) data with parity-conserving and parity-violating polarization asymmetries on quark and gluon helicity distributions in the proton. The extrapolation of structure functions from the current data is studied for the first time in the context of the impact study. Theory assumptions, such as SU(2) and SU(3) flavor symmetries, are also studied to give a more thorough understanding of the impact of EIC pseudodata on proton spin structure.

Zhou, Yiyu↗