Search NASA⌕ Search

SEARCH · Search NASA

Results for “Input Sensitivity”

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

iMatTOUGH: An open-source Matlab-based graphical user interface for pre- and post-processing of TOUGH2 and iTOUGH2 models

TOUGH2 and iTOUGH2 are powerful models that simulate the heat and fluid flows in porous and fracture media, and perform parameter estimation, sensitivity analysis and uncertainty propagation analysis. However, setting up the input files is not only tedious, but error prone, and processing output files is time consuming. Here, we present an open source Matlab-based tool (iMatTOUGH) that supports the generation of all necessary inputs for both TOUGH2 and iTOUGH2 and visualize their outputs. The tool links the inputs of TOUGH2 and iTOUGH2, making sure the two input files are consistent. It supports the generation of rectangular computational mesh, i.e., it automatically generates the elements and connections as well as their properties as required by TOUGH2. The tool also allows the specification of initial and time-dependent boundary conditions for better subsurface heat and water flow simulations. The effectiveness of the tool is illustrated by an example that uses TOUGH2 and iTOUGH2 to estimate soil hydrological and thermal properties from soil temperature data and simulate the heat and water flows at the Rifle site in Colorado.

58 GEOSCIENCES↗

Effects of different vegetation drag parameterizations on the tidal propagation in coastal marshlands

Vegetation drag is a fundamental quantity directly affecting results for both long- and short-term coastal marsh and geomorphological studies. The vegetation drag in coastal marshland has been modeled by various two-dimensional (2D) and threedimensional (3D) numerical parameterizations. 2D parameterizations treat coastal marshes as bottom roughness elements, while 3D parameterizations resolve the vertically-variable vegetation drag through the water column. However, differences in tidal propagation arising from different drag parameterizations within a single model are largely unknown, and clear guidance on parameterization selection is still missing. In this study, we implemented four vegetation drag parameterizations into the Model for Prediction Across Scales-Ocean (MPAS-O), which include 1) a 2D parameterization using land-cover type-determined Manning’s n (2DLM); 2) a 2D parameterization using vegetation-determined Manning’s n (2DVM); 3) a 3D parameterization for stiff vegetation (3DSV); and 4) a 3D parameterization for flexible vegetation (3DFV). Estimates of the flow resistance effects from these parameterizations were compared using a series of idealized tidal propagation simulations. Given the same tidal condition, flooding depth and flooding distance are the largest in the 2DLM simulations and the smallest in the 3DSV simulations. 2DVM results are the closest to the 2DLM results. 3DFV results are the closest to the average of 2DVM, 3DSV, and 3DFV results. 2DVM and 3DSV results are the least and most sensitive to the vegetation aboveground biomass, respectively. Based on the input data requirement and computational efficiency of each parameterization, a comparison summary is provided to help inform parameterization selection for specific applications. Here, the effects of these parameterizations on coastal geomorphology are further discussed, and the results demonstrate that estimates of the long-term evolution of coastal marshes and coastal morphology depend upon the selection of the vegetation drag parameterization

54 ENVIRONMENTAL SCIENCES↗

Transient uncertainty quantification and Global Sensitivity Analysis of the open-source Molten Chloride Reactor Experiment (MCRE) using GP-PCA surrogate models

Uncertainties in the thermophysical properties of molten salts impact both the steady-state and transient behavior of Molten Salt Reactors (MSRs). In this work, we aim to quantify the influence of such uncertainties on the transient operation of the Molten Chloride Reactor Experiment (MCRE), utilizing the open-source specifications provided for this reactor. Seven representative transient scenarios are considered. For each scenario, we evaluate the impact of thermophysical property uncertainties on four key multiphysics model output variables of interest (VoIs): maximum power density, maximum fuel temperature, maximum reflector temperature, and average fuel velocity magnitude. In addition, we perform a Global Sensitivity Analysis (GSA) by computing Sobol’ indices for the uncertain input parameters to determine their contribution to the variability of each VoI. Conducting GSA is computationally intensive due to the large number of required evaluations of the high-fidelity multiphysics model. To mitigate this cost, we develop a surrogate modeling framework that combines Gaussian Process (GP) regression with Principal Component Analysis (PCA), enabling efficient sample generation for the GSA. Our results show that for energy-related VoIs, thermal conductivity is the dominant contributor to uncertainty. In contrast, for flow-related VoIs, density and dynamic viscosity are the primary sources of uncertainty. The specific heat of the fuel salt was found to play a secondary role in the transient analyses.

42 - ENGINEERING↗

Uncertainty evaluation for twist drilling stability model

This paper describes the first uncertainty analysis for drilling stability using a frequency-domain drilling stability model. The stability model inputs include: the modal parameters for the torsional-axial vibration mode from the twist drill-holder-spindle axial frequency response function; and the mechanistic coefficients that relate the torque and thrust force to chip area for the selected drill-workpiece material combination. Furthermore, Monte Carlo simulation is applied to propagate the input uncertainties to output uncertainty in the predicted stability map, which separates stable from unstable (chatter) zones in the spindle speed-chip width parameter space. Additionally, the mean stability boundary and its 95% confidence intervals are determined for five cases: varying all four inputs simultaneously and varying them individually. This enables the individual sensitivities to be compared. Experimental results from drilling tests are included for comparison to the prediction. Additionally, Matlab code is provided to implement the stability model and Monte Carlo uncertainty analysis.

42 ENGINEERING↗

Life-cycle greenhouse gas emissions analysis of battery-grade lithium production in Finland

Various countries are undertaking initiatives to domestically produce battery-related critical materials. Within Finland, Keliber Technology Oy is developing capabilities for battery-grade lithium hydroxide monohydrate (LHM) production from spodumene ores. A detailed life-cycle assessment (LCA) of this pathway is conducted to determine its life-cycle GHG impacts using Argonne's R&D GREET (Research and Development Greenhouse gases, Regulated Emissions, and Energy use in Technologies) model. The analysis shows life-cycle GHG emissions of similar to 9.2 kg CO 2-eq /kg LHM, dominated by contributions from three energy sources - diesel, natural gas, and electricity - and two material inputs - lime (CaO) and soda ash (Na 2 CO 3 ). Sensitivity analyses highlight the potential to reduce these impacts using low-carbon electricity, sequestration of process CO 2 emissions generated during CaO and Na 2 CO 3 production, and bio-based energy for LHM production (by similar to 15 % each). A comparative analysis shows lower impacts for Keliber's LHM than for existing LHM production from Australian spodumene ores processed in China (by similar to 40 %).

CCS↗

The effects of electron thermal radiation on laser ablative shock waves from aluminum plasma into ambient air

The effect of electron thermal radiation on 7 ns laser ablative shock waves from aluminum (Al) plasma into an ambient atmospheric air has been numerically investigated using a one-dimensional, three-temperature (electron, ion, and radiation) radiation hydrodynamic code MULTI. The governing equations in Lagrangian form are solved using an implicit scheme for planar, cylindrical, and spherical geometries. The shockwave velocities (V{sub sw}) obtained numerically are compared with our experimental values obtained over the intensity range of 2.0 × 10{sup 10} to 1.4 × 10{sup 11 }W/cm{sup 2}. It is observed that the numerically obtained V{sub sw} is significantly influenced by the thermal radiation effects which are found to be dominant in the initial stage up to 2 μs depending on the input laser energy. Also, the results are found to be sensitive to the co-ordinate geometry used in the simulation (planar, cylindrical, and spherical). Moreover, it is revealed that shock wave undergoes geometrical transitions from planar to cylindrical nature and from cylindrical to spherical nature with time during its propagation into an ambient atmospheric air. It is also observed that the spatio-temporal evolution of plasma electron and ion parameters such as temperature, specific energy, pressure, electron number density, and mass density were found to be modified significantly due to the effects of electron thermal radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward Verifying Fossil Fuel CO2 Emissions with the CMAQ Model: Motivation, Model Description and Initial Simulation

Motivated by the urgent need for emission verification of CO2 and other greenhouse gases, we have developed regional CO2 simulation with CMAQ over the contiguous U.S. Model sensitivity experiments have been performed using three different sets of inputs for net ecosystem exchange (NEE) and two fossil fuel emission inventories, to understand the roles of fossil fuel emissions, atmosphere-biosphere exchange and transport in regulating the spatial and diurnal variability of CO2 near the surface, and to characterize the well-known ‘signal-to-noise’ problem, i.e. the interference from the biosphere on the interpretation of atmospheric CO2 observations. It is found that differences in the meteorological conditions for different urban areas strongly contribute to the contrast in concentrations. The uncertainty of NEE, as measured by the difference among the three different NEE inputs, has notable impact on regional distribution of CO2 simulated by CMAQ. Larger NEE uncertainty and impact are found over eastern U.S. urban areas than along the western coast. A comparison with tower CO2 measurements at Boulder Atmospheric Observatory (BAO) shows that the CMAQ model using hourly varied and high-resolution CO2 emission from the Vulcan inventory and CarbonTracker optimized NEE reasonably reproduce the observed diurnal profile, whereas switching to different NEE inputs significantly degrades the model performance. Spatial distribution of CO2 is found to correlate with NOx, SO2 and CO, due to their similarity in emission sources and transport processes. These initial results from CMAQ demonstrate the power of a state-of-the art CTM in helping interpret CO2 observations and verify fossil fuel emissions. The ability to simulate CO2 in CMAQ will also facilitate investigations of the utility of traditionally regulated pollutants and other species as tracers to CO2 source attribution.

CMAQ↗

First inverse kinematics measurement of resonances in Be 7 ( α , γ ) C 11 relevant to neutrino-driven wind nucleosynthesis using DRAGON

A possible mechanism to explain the origin of the light $\textit{p}$ nuclei in the Galaxy is the nucleosynthesis in the proton-rich neutrino-driven wind ejecta of core-collapse supernovas via the $\textit{νp}$ process. However, this production scenario is very sensitive to the underlying supernova dynamics and the nuclear physics input. As far as the nuclear uncertainties are concerned, the breakout from the $\textit{pp}$ chains via the 7 Be(α,γ) 11 C reaction has been identified as an important link which can influence the nuclear flow and, therefore, the efficiency of the $\textit{νp}$ process. However, its reaction rate is poorly known over the relevant temperature range, $\textit{T}$ = 1.5–3 GK. We report on the first direct measurement of two resonances of the 7 Be(α,γ) 11 C reaction with previously unknown strengths using an intense radioactive 7 Be beam from the Isotope Separator and Accelerator (ISAC-I) Center facility and the DRAGON recoil separator in inverse kinematics. In this work, we have decreased the 7 Be(α,γ) 11 C reaction rate uncertainty to ≈9.4–10.7% over the relevant temperature region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reweighting Monte Carlo predictions and automated fragmentation variations in Pythia 8

This work reports on a method for uncertainty estimation in simulated collider-event predictions. The method is based on a Monte Carlo-veto algorithm, and extends previous work on uncertainty estimates in parton showers by including uncertainty estimates for the Lund string-fragmentation model. This method is advantageous from the perspective of simulation costs: a single ensemble of generated events can be reinterpreted as though it was obtained using a different set of input parameters, where each event now is accompanied with a corresponding weight. This allows for a robust exploration of the uncertainties arising from the choice of input model parameters, without the need to rerun full simulation pipelines for each input parameter choice. Such explorations are important when determining the sensitivities of precision physics measurements. Accompanying code is available at https://gitlab.com/uchep/mlhad-weights-validation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

MELCOR 2.2 Benchmarks of Peach Bottom NUREG/CR 7155 Uncertainty Analysis

The U.S. Nuclear Regulatory Commission (NRC) performed a first-of-a-kind uncertainty analysis (UA) of the accident progression, radiological releases, and offsite consequences for the State-of- the-Art Reactor Consequence Analyses (SOARCA) of an unmitigated long-term station blackout (LTSBO) severe accident scenario at the Peach Bottom Atomic Power Station. The objective of the UA was to evaluate the robustness of the SOARCA deterministic "best estimate results and conclusions documented in NUREG-1935, and to develop insight into the overall sensitivity of the SOARCA results to uncertainty in key modeling inputs. The study was completed in 2015 and documented in NUREG/CR-7155. Since 2015, two other SOARCA UAs were completed for two pressurized water reactor (PWR) plants. The PWR UAs incrementally updated the approach and methodology, including using the latest release of the MELCOR 2.2 computer code. There were also advances made in the state-of-the-art modeling related to NRC efforts using the Peach Bottom model in NUREG-2206, which provided the technical basis for the containment protection and release reduction rulemaking for boiling water reactors with Mark I and Mark II containments. This report documents the input model changes from the NUREG/CR-7155 study and performs a small number of reference calculations to assess the changes of the new computer code and the model input updates. The objective of the work is to verify whether the updated Peach Bottom MELCOR model and updated version of MELCOR support the conclusions formed in the Peach Bottom SOARCA UA by performing these representative calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Characterization of the Fusion Zone in Laser Beam Welds on SS304L

Designing a welded joint requires consideration of the welding process, component/assembly structure and weld geometry, material composition, properties, and desired joint performance. The designer must also specify a set of acceptance criteria used to evaluate the welds. This study was conducted to support the design team responsible for specifying the spot weld profile by investigating the relationship between weld process parameters and the weld geometry and properties to further inform design optimization. The subjects of this study are two thin stainless steel 304L (SS304L) components joined by Laser Beam Welding (LBW). The components are the Housing, which contains flexible circuitry, and the Cover that protects the circuitry during handling, installation, and under service loads. Figure 1 illustrates coupon versions of the Cover and Housing used to simulate the interface between these components for weld testing. These coupon representations of assembly components were designed to match the material and relevant geometry related to the welded joint while other features not relevant to the welded joint are excluded.The objective of this study was to characterize the sensitivity of the weld geometry and mechanical properties to the input weld parameters to inform further Finite Element Analysis (FEA) efforts in balancing joint strength against assembly deformation caused by the welding process. Characterization was accomplished by varying laser welding parameters around the current design values and measuring weld geometry on the cross-sections. Nanoindentation with a spherical and flat-ended tip was used on polished section cuts of the weld to obtain direct estimates of the stress-strain response of the Fusion Zone (FZ) for comparison to the Base Metal (BM).

42 ENGINEERING↗

Application of Machine Learning Algorithms to Identify Problematic Nuclear Data

In this work we aim to show that Machine learning algorithms are promising tools for the identification of nuclear data that contribute to increased errors in transport simulations. We demonstrate this through an application of a machine learning algorithm (Random Forest) to the Whisper/MCNP6 criticality validation library to identify nuclear data that are associated with an increase of the bias (simulated - experimental $k_{eff}$) in the calculations. Specifically, the $k_{eff}$ sensitivity profiles (w.r.t. nuclear data) of 233 U solution benchmarks are used to predict the bias and Shapley Additive Explanations (SHAP) are used to explain how the sensitivities are related to the predicted bias. The SHAP values can be interpreted as sensitivity coefficients of the machine learning model to the $k_{eff}$ sensitivities which are used to make predictions of bias. Using the SHAP values we can identify specific subsets of nuclear data which have the highest probability of influencing bias. We demonstrate the utility of this method by showing how SHAP values were used to identify an inconsistency in the 19 F inelastic scattering nuclear data. The methodology presented here is not limited to transport problems and can be applied to other simulations if there are experimental measurements to compare against, simulations of those experimental measurements, and the ability to calculate sensitivities of the model output with respect to the data inputs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Soil bacterial and fungal community responses to nitrogen addition across soil depth and microhabitat in an arid shrubland

Arid shrublands are stressful environments, typified by alkaline soils low in organic matter, with biologically-limiting extremes in water availability, temperature, and UV radiation. The widely-spaced plants and interspace biological soil crusts in these regions provide soil nutrients in a localized fashion, creating a mosaic pattern of plant- or crust-associated microhabitats with distinct nutrient composition. With sporadic and limited rainfall, nutrients are primarily retained in the shallow surface soil, patterning biological activity. We examined soil bacterial and fungal community responses to simulated nitrogen (N) deposition in an arid Larrea tridentata-Ambrosia dumosa field experiment in southern Nevada, USA, using high-throughput sequencing of ribosomal RNA genes. To examine potential interactions among the N application, microhabitat and soil depth, we sampled soils associated with shrub canopies and interspace biological crusts at two soil depths (0–0.5 or 0–10 cm) across the N-amendment gradient (0, 7, and 15 kg ha –1 yr –1 ). We hypothesized that localized compositional differences in soil microbiota would constrain the impacts of N addition to a microhabitat distribution that would reflect highly localized geochemical conditions and microbial community composition. The richness and community composition of both bacterial and fungal communities differed significantly by microhabitat and with soil depth in each microhabitat. Only bacterial communities exhibited significant responses to the N addition. Community composition correlated with microhabitat and depth differences in soil geochemical features. Provided the distinct roles of soil bacteria and fungi in major nutrient cycles, the resilience of fungi and sensitivity of bacteria to N amendments suggests that increased N input predicted for many arid ecosystems could shift nutrient cycling toward pathways driven primarily by fungal communities.

54 ENVIRONMENTAL SCIENCES↗

Methodologies for optimal resource allocation to the national space program and new space utilizations. Volume 1: Technical description

The optimal allocation of resources to the national space program over an extended time period requires the solution of a large combinatorial problem in which the program elements are interdependent. The computer model uses an accelerated search technique to solve this problem. The model contains a large number of options selectable by the user to provide flexible input and a broad range of output for use in sensitivity analyses of all entering elements. Examples of these options are budget smoothing under varied appropriation levels, entry of inflation and discount effects, and probabilistic output which provides quantified degrees of certainty that program costs will remain within planned budget. Criteria and related analytic procedures were established for identifying potential new space program directions. Used in combination with the optimal resource allocation model, new space applications can be analyzed in realistic perspective, including the advantage gain from existing space program plant and on-going programs such as the space transportation system.

Source record↗

Gas physics environment of high altitude Jovian entry.

A description of the gas physics environment associated with high-altitude Jupiter entry is reported. Two aspects of environmental analysis are treated. First, limited discussion is devoted to prediction of the heating environment of the probe aeroshell. Prime emphasis is given to the second aspect of the analysis, a description of the gas particle probe interaction phenomena in the rarefied flow regime. Several molecular reemission models have been devised and exercised. Model selection was based upon the desire to bracket the expected reemission characteristics. Consideration was given to examining diffuse reflection as exemplified by the Maxwellian model and lobular, specular reflection as characterized by the Nocilla model. Variation of fundamental input parameters was investigated to determine particle-probe interaction sensitivity. Regardless of neutral species type, results show that the attenuation of the incoming flux of that species by reflected particles is approximately 15 per cent in the stagnation region at 200 km above the cloud tops.

Lefferdo, J. M.↗

The impact of chlorofluoromethane and NO/x/ injections on stratospheric ozone

The chlorofluoromethane problem and the SST problem have been subject to continuing evaluation by the use of one-dimensional stratospheric photochemistry models. The problems to some extent appear to be inversely related to one another through hydroxyl radical reactions; in other words, processes which tend to increase the amount of OH in a model tend to increase the sensitivity of the ozone column to ClX perturbations, while at the same time decreasing the sensitivity to NO(x) perturbations. Evaluations of sensitivity to ClX and NO(x) have varied as new input information has become available and modeling techniques have improved. It is clear that either NO(x) or ClX, or both, will affect the distribution of ozone with altitude.

Stolarski, R. S.↗

Analysis and correlation with theory of rotor lift-limit test data

A wind tunnel test program to define the cruise performance and determine any limitations to lift and propulsive force of a conventional helicopter rotor is described. A 2.96 foot radius model rotor was used. The maximum lift and propulsive force obtainable from an articulated rotor for advance ratios of 0.4 to 0.67, and the blade load growth as the lift approaches the limit are determined. Cruise rotor performance for advance ratios of 0.4 to 0.67 and the sensitivity of the rotor forces and moments to rotor control inputs as the lift limit is approached are established.

Sheffler, M.↗