Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity and Uncertainty 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 289 records · Page 16

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗

10° off-axis testing of CFRP using DIC: A study on strength, strain and modulus

For the determination of the shear properties such as strength, strain and modulus of unidirectionally reinforced composites, nearly a dozen measurement techniques exist, several of which have found their way into standards. One non-standardized method for sheet materials is the 10° off-axis tensile test of unidirectional coupons. The current work focuses on the experimental investigation of tensile testing of CFRP-UD coupon specimens using the 10° off-axis test and DIC strain mapping to measure the full-field deformation response. A sensitivity analysis on the measured stresses, strains and moduli is performed by varying the aspect ratio of the specimens and the location of strain measurement. Test related uncertainties such as the loading angle and specimen preparation details are investigated and quantified with the help of these full-field DIC measurements.

42 ENGINEERING↗

Modeling Riverine Dissolved and Particulate Organic Carbon Fluxes from Two Small Watersheds in the Northeastern United States

The coupled carbon (C) cycle across terrestrial and aquatic environments at the watershed scale has been identified as an important, but poorly constrained component of the global carbon budget. Here, we extended Soil and Water Assessment Tool (SWAT) with coupled riverine particulate organic carbon (POC) and dissolved organic carbon (DOC) modules (referred to as SWAT-C hereafter). Results show that SWAT-C reproduced daily POC and DOC fluxes well in two watersheds in the Northeastern United States. We found that SWAT-C tended to underestimate high flows and peak DOC and POC fluxes. Uncertainty analysis indicated flux uncertainties associated with POC and DOC simulation were larger than those for flow simulation. Sensitive parameters controlling POC and DOC biogeochemical processes were identified along with how these parameters influence mechanisms underlying C cycling. We anticipate that the tool developed and applied here will inform C related ecosystem services in watershed assessment and planning.

Dissolved organic carbon, Particulate organic carb↗

Seasonal variability of wake impacts on US mid-Atlantic offshore wind plant power production

Abstract. The mid-Atlantic will experience rapid wind plant development due to its promising wind resource located near large population centers. Wind turbines and wind plants create wakes, or regions of reduced wind speed, that may negatively affect downwind turbines and plants. We evaluate wake variability and annual energy production with the first yearlong modeling assessment using the Weather Research and Forecasting model, deploying 12 MW turbines across the domain at a density of 3.14 MW km−2, matching the planned density of 3 MW km−2. Using a series of simulations with no wind plants, one wind plant, and complete build-out of lease areas, we calculate wake effects and distinguish the effect of wakes generated internally within one plant from those generated externally between plants. We also provide a first step towards uncertainty quantification by testing the amount of added turbulence kinetic energy (TKE) by 0 % and 100 %. We provide a sensitivity analysis by additionally comparing 25 % and 50 % for a short case study period. The strongest wakes, propagating 55 km, occur in summertime stable stratification, just when New England's grid demand peaks in summer. The seasonal variability of wakes in this offshore region is much stronger than the diurnal variability of wakes. Overall, yearlong simulated wake impacts reduce power output by a range between 38.2 % and 34.1 % (for 0 %–100 % added TKE). Internal wakes cause greater yearlong power losses, from 29.2 % to 25.7 %, compared to external wakes, from 14.7 % to 13.4 %. The overall impact is different from the linear sum of internal wakes and external wakes due to non-linear processes. Additional simulations quantify wake uncertainty by modifying the added amount of turbulent kinetic energy from wind turbines, introducing power output variability of 3.8 %. Finally, we compare annual energy production to New England grid demand and find that the lease areas can supply 58.8 % to 61.2 % of annual load. We note that the results of this assessment are not intended to make nor are they suitable to make commercial judgments about specific wind projects.

17 WIND ENERGY↗

BioSTEAMDevelopmentGroup/thermosteam

BioSTEAM is a fast and flexible package for the design, simulation, and techno-economic analysis of biorefineries under uncertainty. BioSTEAM’s framework is built to streamline and automate early-stage technology evaluations and to enable rigorous sensitivity and uncertainty analyses. Complete biorefinery configurations are available at the Bioindustrial-Park GitHub repository, BioSTEAM’s premier repository for biorefinery models and results. The long-term growth and maintenance of BioSTEAM is supported through both community-led development and the research institutions invested in BioSTEAM. Through the open-source and community-lead platform, BioSTEAM aims to foster communication and transparency within the biorefinery research community for an integrated effort to expedite the evaluation of candidate biofuels and bioproducts. Additionally, an agile life cycle assessment (LCA) platform has been designed to interface with BioSTEAM, BioSTEAM-LCA. This open-source, installable package allows users to perform streamlined LCAs of biorefineries. The focus of BioSTEAM-LCA is to streamline and automate early-stage environmental impact analyses of processes and technologies, and to enable rigorous sensitivity and uncertainty analyses linking process design, performance, economics, and environmental impacts. ThermoSTEAM is a standalone thermodynamic engine capable of estimating mixture properties, solving thermodynamic phase equilibria, and modeling stoichiometric reactions. ThermoSTEAM builds upon chemicals, the chemical properties component of the Chemical Engineering Design Library, with a robust and flexible framework that facilitates the creation of property packages. The Biorefinery Simulation and Techno-Economic Analysis Modules (BioSTEAM) is dependent on ThermoSTEAM for the simulation of unit operations.

Cortes-Peña, Yoel↗

Uniformly Ordered Binary Decision Algorithm for Benchmark Experiment Correlations in Whisper Validation

When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sensitivity Analysis of Salt Storage Cavern Mechanical Integrity Test Parameters

This report examines the sensitivity of salt cavern Mechanical Integrity Tests (MIT) to uncertainties in key test parameters. MIT's are used by cavern operators to detect and quantify leak rates in access wells in underground salt storage caverns and involve a suite of measured and assumed parameters that have a direct impact on the sensitivity of the testing to detect actual leaks from the cavern storage system. Determining the sensitivity of the testing to these different parameters provides a basis for understanding the results from, and informing the design criteria for, this type of testing. Without fully understanding the sensitivity of the test to the testing parameters, it is possible that the test results may not accurately reflect the integrity of the cavern system; an actual leak may be missed, or an intact system may be interpreted as leaking. This report reviews the main parameters included in MITs and examines how selected changes in their values can impact test results. The deviations used in the sensitivity analyses were designed to be within the ranges believed to be similar to those which may be encountered during testing. The results show that small, plausible fluctuations in some of the parameters measured values can have a significant impact on the testing results. Of the parameters studied here, the sensitivity analyses showed the order of importance to be (from highest to lowest): nitrogen-oil interface depth measurement, well bore temperature, well head pressure, and finally the internal geometry of the testing interval.

42 ENGINEERING↗

Making coal relevant for small scale applications: Modular gasification for syngas/engine CHP applications in challenging environments

Small-scale coal gasification technology, coupled to a reciprocating engine generator, has the potential for making coal a cost competitive resource for meeting the flexible energy needs and resiliency requirements of utilities across the United States. To maintain grid stability and reliability, electrical generation must be regulated to match the load at the proper voltage and frequency. With the expansion of intermittent sources into the grid, such as wind and solar, frequency and voltage regulation become increasingly important and challenging. This work presents the results of a Front End Engineering Design (FEED) effort to detail the engineering and preliminary economics of a small-scale, air blown, fixed-bed gasification process, operating at near-atmospheric pressure, with gas cleanup to provide syngas and pyrolysis liquid fuels for use in reciprocating engine generators for combined heat and power at the University of Alaska Fairbanks. A very detailed assessment of capital and operating costs allows the evaluation of a levelized cost of electricity of 208.06 $/MWh, which can be reduced through the selling of by-products (steam for heating purposes and pyrolysis liquids to be used for power generation in an existing diesel engine). Here, a combined sensitivity analysis, based on the Monte Carlo approach, has been carried out to evaluate the effects of the uncertainties (capital and operating costs and plant annual availability) on the LCOE.

01 COAL, LIGNITE, AND PEAT↗

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geo- and reactor antineutrino sensitivity at THEIA

We present the sensitivity of the THEIA experiment to low-energy geo- and reactor antineutrinos. For this study, we consider one of the possible proposed designs, a 17.8-ktonne fiducial volume THEIA-25 detector filled with water-based liquid scintillator placed at Sanford Underground Research Facility (SURF). We demonstrate THEIA’s sensitivity to measure the geo- and reactor antineutrinos via inverse-beta decay interactions after one year of data taking with $11.9 × 10^{32}$ free target protons. Considering all uncertainties on input throughout the whole analysis chain, the expected number of geo- and reactor antineutrinos is $220^{+30}_{–24}$ (stats+syst) and $168^{+26}_{–24}$ (stat+sys), respectively, after one year of data taking. The corresponding expected fit precision of a sole experiment is evaluated at 8.7% and 10.1%, respectively. We also demonstrate the sensitivity towards fitting individual Th and U contributions, with best fit values of $N_{\text{Th}} = 40^{+26}_{–22}$ (stat+sys) and $N_\text{U}=180^{+30}_{–24}$ (stat+sys). Finally, from the fit results of individual Th and U contributions, we evaluate the mantle signal to be $S_{\text{mantle}} = 9.3±[5.2,5.4]$ NIU (stat+sys). This was obtained assuming a full-range positive correlation $(ρ_c ϵ[0,1])$ between Th and U, and the projected uncertainties on the crust contributions of 8.3% (Th) and 7.0% (U).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Probabilistic R&D Portfolio Management Applied to Energy Storage

Here this paper proposes a methodology to analyze the potential of several alternative technologies in advancing a research and development portfolio. The methodology considers technology maturation, use case requirements, performance goals, cost projections, and existing research and development (R&D) efforts to predict the impact of various R&D investment tracks across several energy storage technologies. Technology maturation and cost projections are derived from a probabilistic Monte Carlo simulation, yielding results that incorporate uncertainty in modeling the impact of future R&D. This framework allows the researcher to perform a robust sensitivity analysis to uncover critical pathways and areas where investment will have a high impact.

25 ENERGY STORAGE↗

High-Pressure Laminar Flame Speeds and Markstein Lengths of Syngas Flames Diluted in Carbon Dioxide and Helium

Abstract New laminar flame speed and burned-gas Markstein length data for H2–CO–O2–CO2–He mixtures have been measured from spherically expanding flames. Experiments were conducted at 10 atm and room temperature for H2:CO ratios ranging from 2:1 to 1:4 and for overall CO2 mole fractions from 0% to 30%. CO2 dilution had little effect on Markstein length, but CO2 dilutions of 10%, 20%, and 30% caused average reductions in flame speed of 47%, 73%, and 89%, respectively, regardless of H2:CO ratio. The study was designed to isolate the dilution effect of CO2 on flame speed, and a detailed analysis using the FCO2 method was used to show that the chemical-kinetic participation of CO2 was responsible for up to 20% of the reduction in flame speed. Hence, the majority (80% or more) of the reduction in flame speed due to CO2 is from the thermal effect. Accurate flame speed predictions were produced by five different chemical kinetics mechanisms for most conditions, with the slight exception of high-CO, high-CO2 mixtures. A thorough sensitivity analysis highlighted the larger effect of CO2 dilution on the important kinetics reactions than the effect of changing H2:CO. Sensitivity analysis also showed that the chain branching reaction H2O + O ⇌ OH + OH could be modified (albeit beyond its uncertainty) to achieve more accurate flame speed predictions, but also indicated that further improvement of flame speed modeling would require changes to many lesser reactions.

Engineering↗

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗

Accounting for uncertainty sources in pyrometric measurements: Variables affecting Surface Temperature measurements using the LLNL Streak Spectrometers

Here we present a detailed analysis of the error sensitivities of the LLNL temperature analysis package used in analyzing streak spectrometer data. These sensitivities include noise, temperature, wavelength, emissivity, missing bands, and missing data points. The analysis package was tested for its ability to accurately reconstruct randomly generated noisy data and found to perform within specifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗

SCALE inventory and reactivity analysis as part of the Hermes 2021 PSAR review

The readiness of SCALE for comprehensive studies of pebble-bed reactors has been demonstrated through detailed analysis of a fluoride salt–cooled, high-temperature pebble-bed reactor (PB-FHR). The methods developed for pebble-bed reactor modeling in SCALE, particularly for inventory generation, have proven effective in gaining insights into the reactor physics of this advanced reactor. Excellent agreement with another code package has been observed, further highlighting SCALE’s strong performance. The SCALE results supported the US Nuclear Regulatory Commission’s construction permit application review of the Hermes low-power PB-FHR demonstration reactor. A SCALE model of the Hermes reactor was developed at Oak Ridge National Laboratory using information from the Preliminary Safety Analysis Report (PSAR) and supplemented with publicly available data. SCALE reactivity coefficient simulations reproduced PSAR results within 1σ statistical uncertainties. Sensitivity studies emphasized the importance of graphite specifications for accurate keff predictions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Uncertainty Quantification in High-Low Dynamic System Coupling using RAVEN and TRANSFORM

This work demonstrates new functionality and applications stemming from the development of high-fidelity to low-fidelity (high-low) coupling for system simulations and to further explore the capabilities of the Risk Analysis Virtual Environment (RAVEN) in the performance of uncertainty quantification in this kind of high-low coupled system models. The work builds from previous work on high-low coupling that utilized COBRA-TF (CTF), the high-fidelity subchannel analysis code, with a low fidelity model built in ORNL’s TRANSFORM, the system analysis code, utilizing the Functional Mock-Up Interface (FMI). Steady-state and transient analysis examples using the high/low coupled models generated from CTF and TRANSFORM/FMI are investigated. The workflows for both steady-state and transient coupled simulations are described. A steady-state parameter sweep and uncertainty analysis of the primary flow rates and reactor power are demonstrated. Likewise, a transient pump trip and power ramp sensitivity studies are also demonstrated. This work elucidates some of the potential benefits and future needs of using RAVEN for high/low system coupling analysis of energy systems. It also shows some of the difficulties that can be encountered in coupling system simulations.

Williams, Wesley↗

A non-intrusive optical (NIO) approach to characterize heliostats in utility-scale power tower plants: Sensitivity study

Power towers, a type of concentrating power tower technology, use a large number of heliostats to concentrate sunlight and produce renewable energy. Optical errors of heliostats can cause drastic losses in the efficiency of power-tower plants. Accurately measuring optical errors is crucial to assessing and improving plant performance. This analysis discusses an innovate non-intrusive optical (NIO) approach to measure slope, canting, and tracking errors by detecting distortions in the reflected tower structure in heliostat images. Further in this work a sensitivity study is carefully conducted to determine the uncertainty requirements to allow the method to calculate slope errors with an accuracy of 0.25 mrad. Measurement uncertainty sources include camera resolution and position uncertainty, tower position uncertainty, and number of collected images. Each uncertainty source is investigated to determine its impact on the accuracy of the slope-error calculation. A combination of theoretical results and experimental results obtained from data collected on a heliostat at Sandia National Laboratories is used to determine and validate uncertainty requirements. The analysis shows that a measurement uncertainty of 0.25 mrad can be realized by realistically controlling uncertainty sources when implementing the NIO method. It demonstrates the superior performance of NIO in performing in-situ optical characterization.

14 SOLAR ENERGY↗