Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity 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 1,621 records · Page 90

Sensitivity and Uncertainty Analysis of Plutonium and Cesium Isotopes in Modeling of BR3 Reactor Spent Fuel

Nuclear forensics has a broad task to characterize recovered nuclear or radiological material and interpret the results of investigation. One approach to isotopic characterization of nuclear material obtained from a reactor is to chemically separate and perform isotopic measurements on the sample and verify the results with modeling of the sample history, for example, operation of a nuclear reactor. The major actinide plutonium and fission product cesium are commonly measured signatures of the fuel history in a reactor core. This study investigates the uncertainty of the plutonium and cesium isotope ratios of a fuel rod discharged from a research pressurized water reactor when the location of the sample is not known a priori. A sensitivity analysis showed overpredicted values for the 240 Pu/ 239 Pu ratio toward the axial center of the rod and revealed a lower probability of the rod of interest (ROI) being on the periphery of the assembly. The uncertainty analysis found the relative errors due to only the rod position and boron concentration to be 17% to 36% and 7% to 15% for the 240 Pu/ 239 Pu and 137 Cs/ 135 Cs ratios, respectively. Lastly, this study provides a method for uncertainty quantification of isotope concentrations due to the location of the ROI. Similar analyses can be performed to verify future chemical and isotopic analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Adding stable carbon isotopes improves model representation of the role of microbial communities in peatland methane cycling

Abstract Climate change is expected to have significant and uncertain impacts on methane (CH 4 ) emissions from northern peatlands. Biogeochemical models can extrapolate site‐specificCH 4 measurements to larger scales and predict responses of CH 4 emissions to environmental changes. However, these models include considerable uncertainties and limitations in representing CH 4 production, consumption, and transport processes. To improve predictions of CH 4 transformations, we incorporated acetate and stable carbon (C) isotopic dynamics associated with CH 4 cycling into a biogeochemistry model, DNDC. By including these new features, DNDC explicitly simulates acetate dynamics and the relative contribution of acetotrophic and hydrogenotrophic methanogenesis (AM and HM) to CH 4 production, and predicts the C isotopic signature (δ 13 C) in soil C pools and emitted gases. When tested against biogeochemical and microbial community observations at two sites in a zone of thawing permafrost in a subarctic peatland in Sweden, the new formulation substantially improved agreement with CH 4 production pathways and δ 13 C in emitted CH 4 (δ 13 C‐CH 4 ), a measure of the integrated effects of microbial production and consumption, and of physical transport. We also investigated the sensitivity of simulated δ 13 C‐CH 4 to C isotopic composition of substrates and, to fractionation factors for CH 4 production (α AM and α HM ), CH 4 oxidation (α MO ), and plant‐mediated CH 4 transport (α TP ). The sensitivity analysis indicated that the δ 13 C‐CH 4 is highly sensitive to the factors associated with microbial metabolism (α AM , α HM , and α MO ). The model framework simulating stable C isotopic dynamics provides a robust basis for better constraining and testing microbial mechanisms in predicting CH 4 cycling in peatlands.

54 ENVIRONMENTAL SCIENCES↗

Simulation process and data flow for a large system dynamics model

This paper documents the workflow and supporting technologies that a large system dynamics model, the biomass scenario model, employs to streamline the data preparation, simulation, quality control, and analysis process at the National Renewable Energy Laboratory. The workflow centers on automation of routine aspects of the flow of data between data stores, simulations, and visualizations. It enforces quality checks on data, reproducibility of computations, and traceability of results, while maintaining complete archives of modeling and analysis artifacts. The resulting frictionless simulation/analysis environment supports large-scale sensitivity analysis, interactive creation of ensembles of simulations, and rapid visualization-based exploration of simulation results.

09 BIOMASS FUELS↗

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗

DME-Propane Ignition Delay Time Measurements at Mixing Controlled Compression Ignition Engine-Relevant Conditions

The blend of dimethyl ether (DME, CH 3 OCH 3 ) and propane (C 3 H 8 ) is a potentially renewable fuel mixture that has the potential to replace diesel in compression ignition engines. The combination can potentially reduce particulate and greenhouse gas emissions compared to a conventional diesel engine operating under similar conditions. However, detailed conceptual and simulation studies must be conducted before adopting a new fuel on a compression ignition engine. For these simulations, accurate chemical kinetic models are necessary. However, the validity of chemical kinetic mechanisms in the literature is unknown for mixing controlled compression ignition (MCCI) engine operating conditions. Hence, in this work, we studied the ignition of dimethyl ether (DME) and propane blends in a shock tube at MCCI engine conditions. Ignition delay time (IDT) data was collected behind the reflected shock for DME-propane mixtures for heavy-duty compression ignition (CI) engine parameters. Undiluted experiments spanning temperatures of 700 to 1100 K and pressures of 55 to 84 bar for various blends (100% CH 3 OCH 3 / 0% C 3 H 8, 100% C 3 H 8 / 0% CH 3 OCH 3 , 60% CH 3 OCH 3 / 40% C 3 H 8 ) of DME and propane were combusted in synthetic air (21% O 2 / 79% N 2 ). Some experiments were conducted at higher pressures (90-120 bar) to understand model performance at these conditions. Comparisons of IDT were made with the predictions of recent chemical kinetic mechanisms for DME-propane mixture, including the Aramco3.0, NUIG, and Dames et al. mechanisms. All mechanisms overpredicted IDT compared to experimental values. Sensitivity analysis was conducted with Dames et al. model, and critical reactions sensitive to IDT of DME-propane mixture near 100 bar are outlined.

Mohammed, Zuhayr Pasha↗

Using a surrogate-assisted Bayesian framework to calibrate the runoff-generation scheme in the Energy Exascale Earth System Model (E3SM) v1

Abstract. Runoff is a critical component of the terrestrial water cycle, and Earth system models (ESMs) are essential tools to study its spatiotemporal variability. Runoff schemes in ESMs typically include many parameters so that model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at a global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using a surrogate-assisted Bayesian framework. First, the polynomial chaos expansion machinery with Bayesian compressed sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5∘ × 0.5∘ for 1991–2010. The error metric between the ELM simulations and the benchmark data is selected to construct the surrogates, which facilitates efficient calibration and avoids the more conventional, but challenging, construction of high-dimensional surrogates for the ELM simulated runoff. Second, the Sobol' index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that, in most regions, ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multimodel ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.

58 GEOSCIENCES↗

Solar Electric Propulsion Stage design

This paper summarizes the analyses and design studies performed for the Solar Electric Propulsion Stage (SEPS). Mission analysis results were used to establish general design requirements. A performance sensitivity analysis was then used to define an operational power level of 21 kW for accomplishing a broad range of interplanetary and earth-orbital missions. With this base, propulsion system and stage subsystem design analyses were performed to support concept and configuration studies. Design study areas include solar array, thrust subsystem, thermal control, structures, and astrionics. Results are presented primarily for areas not covered by more specific technical papers (i.e., thrust subsystem, solar array, design requirements and configuration studies).

Masek, T. D.↗

Characterization of iron-doped lithium niobate for holographic storage applications

A comprehensive characterization of chemical and holographic properties of eight systematically chosen Fe:LiNbO3 crystals is performed in order to determine optimum performance of the crystals in holographic storage and display applications. The discussion covers determination of Fe(2+) and Fe(3+) ion concentrations in Fe:LiNbO3 system from optical absorption and EPR measurements; establishment of the relation between the photorefractive sensitivity of Fe(2+) and Fe(3+) concentrations; study of the spectral dependence, the effect of oxygen annealing, and of other impurities on the photorefractive sensitivity; analysis of the diffraction efficiency curves for different crystals and corresponding sensitivities with the dynamic theory of hologram formation; and determination of the bulk photovoltaic fields as a function of Fe(2+) concentrations. In addition to the absolute Fe(2+) concentration, the relative concentrations of Fe(2+) and Fe(3+) ions are also important in determining the photorefractive sensitivity. There exists an optimal set of crystal characteristics for which the photorefractive sensitivity is most favorable.

Shah, R. R.↗

Shape optimization of three-dimensional stamped and solid automotive components

The shape optimization of realistic, 3-D automotive components is discussed. The integration of the major parts of the total process: modeling, mesh generation, finite element and sensitivity analysis, and optimization are stressed. Stamped components and solid components are treated separately. For stamped parts a highly automated capability was developed. The problem description is based upon a parameterized boundary design element concept for the definition of the geometry. Automatic triangulation and adaptive mesh refinement are used to provide an automated analysis capability which requires only boundary data and takes into account sensitivity of the solution accuracy to boundary shape. For solid components a general extension of the 2-D boundary design element concept has not been achieved. In this case, the parameterized surface shape is provided using a generic modeling concept based upon isoparametric mapping patches which also serves as the mesh generator. Emphasis is placed upon the coupling of optimization with a commercially available finite element program. To do this it is necessary to modularize the program architecture and obtain shape design sensitivities using the material derivative approach so that only boundary solution data is needed.

Botkin, M. E.↗

Dual Extended Kalman Filter for the Identification of Time-Varying Human Manual Control Behavior

A Dual Extended Kalman Filter was implemented for the identification of time-varying human manual control behavior. Two filters that run concurrently were used, a state filter that estimates the equalization dynamics, and a parameter filter that estimates the neuromuscular parameters and time delay. Time-varying parameters were modeled as a random walk. The filter successfully estimated time-varying human control behavior in both simulated and experimental data. Simple guidelines are proposed for the tuning of the process and measurement covariance matrices and the initial parameter estimates. The tuning was performed on simulation data, and when applied on experimental data, only an increase in measurement process noise power was required in order for the filter to converge and estimate all parameters. A sensitivity analysis to initial parameter estimates showed that the filter is more sensitive to poor initial choices of neuromuscular parameters than equalization parameters, and bad choices for initial parameters can result in divergence, slow convergence, or parameter estimates that do not have a real physical interpretation. The promising results when applied to experimental data, together with its simple tuning and low dimension of the state-space, make the use of the Dual Extended Kalman Filter a viable option for identifying time-varying human control parameters in manual tracking tasks, which could be used in real-time human state monitoring and adaptive human-vehicle haptic interfaces.

manual control↗

Analysis of the uncertainties in velocity measurements with triple hot-wire probes

A detailed computerized sensitivity analysis of the triple hot-wire equations has been performed in order to delineate the uncertainties associated with measurements of the velocity components. Absolute and relative uncertainties for the instantaneous hot-wire outputs are calculated as functions of roll and pitch angles, based on a constant probability combination of the uncertainties in the inputs. From the results, it is concluded that the small inherent difficulties associated with the triple hot-wire data do not reflect artifacts introduced by the data processing. Fixed errors present in the V and W channels of the output are due to the nonzero measuring volume of the triple wire probe, and are entirely predictable.

Frota, M. N.↗

Aerodynamic Sensitivities over Separable Shape Tensors

Here, we present a comprehensive aerodynamic sensitivity analysis of airfoil parameterization informed by separable shape tensors. This parameterization approach uniquely benefits the design process by isolating various well-studied shape characteristics, such as airfoil thickness, and providing a well-regulated low-dimensional parameter domain for aerodynamic designs. Exploring the aerodynamic sensitivities of this novel parameterization can provide valuable insights for more robust designs and future manufacturing efforts. We construct a data-driven parameter space of airfoils using principal geodesic analysis of separable shape tensors informed by a curated database containing almost 20,000 suitable engineering airfoils. Analyzing the shape reconstruction error and the maximum mean discrepancy between joint distributions of aerodynamic quantities, we study the dimensionality of the learned parameter space. This simple numerical experiment demonstrates a dramatic dimension reduction that retains design effectiveness and promotes regularity of the shape representations. Finally, we generate new airfoils and use the HAM2D Reynolds-averaged Navier–Stokes solver to predict lift, drag, and moment coefficients. We compute multiple sensitivity metrics to quantify and assert the consistency of parameter influence on the aerodynamic quantities. We also explore low-dimensional polynomial ridge approximations to motivate physical intuitions and offer explanations of the approximated sensitivities.

17 WIND ENERGY↗

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than "building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of MCNP, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers. In this paradigm, reduced complexity computational and analytical models still have an important purpose. The explicit form of an analytic solution is arguably the best way to understand the qualitative properties of simple models. In contrast to "building up" a complex problem through understanding simpler problems, results from detailed computational scenarios can be better explained by "building down" the complex model through simple models rooted in the fundamental or essential phenomenology. Simplified analytic and computational models can be used to 1) increase a user's confidence in the computational solution of a complex model, 2) confirm there are no user input errors, and 3) ensure essential assumptions of the simulation tool are preserved. This process of using analytic models to develop a more valuable analysis of simulation results is named the results analysis methodology. The utility of the results assessment methodology and a complimentary sensitivity analysis is exemplified through the analysis of the neutron flux in a dry used fuel storage cask. This application was chosen due to current scientific interest in used nuclear fuel storage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sensitivity and uncertainty of the IFR-1 BISON benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered for adoption by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This paper summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment.For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. Here, the resulting variations in benchmark predictions were normalized to calculate sensitivities. These sensitivities were then used to inform input selections for the UQ.The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of the benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because the confidence intervals associated with benchmark predictions were wide. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A simple, mass balance model of carbon flow in a controlled ecological life support system

Internal cycling of chemical elements is a fundamental aspect of a Controlled Ecological Life Support System (CELSS). Mathematical models are useful tools for evaluating fluxes and reservoirs of elements associated with potential CELSS configurations. A simple mass balance model of carbon flow in CELSS was developed based on data from the CELSS Breadboard project at Kennedy Space Center. All carbon reservoirs and fluxes were calculated based on steady state conditions and modelled using linear, donor-controlled transfer coefficients. The linear expression of photosynthetic flux was replaced with Michaelis-Menten kinetics based on dynamical analysis of the model which found that the latter produced more adequate model output. Sensitivity analysis of the model indicated that accurate determination of the maximum rate of gross primary production is critical to the development of an accurate model of carbon flow. Atmospheric carbon dioxide was particularly sensitive to changes in photosynthetic rate. The small reservoir of CO2 relative to large CO2 fluxes increases the potential for volatility in CO2 concentration. Feedback control mechanisms regulating CO2 concentration will probably be necessary in a CELSS to reduce this system instability.

Garland, Jay L.↗

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗

mvBayesR

SAND2025-11559O The mvBayesR tool performs multivariate Bayesian analysis on generic data. It includes tools for regression modeling, diagnosis, basis decomposition, sensitivity analysis, and visualization. The tool compiles state-of-the-art methodology into one easy-to-use package. 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.

Tucker, James [Sandia National Lab. (SNL-CA), Live↗