Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-principles simulation of light-ion microscopy of graphene

The extreme sensitivity of 2D materials to defects and nanostructure requires precise imaging techniques to verify presence of desirable and absence of undesirable features in the atomic geometry. Helium-ion beams have emerged as a promising materials imaging tool, achieving up to 20 times higher resolution and 10 times larger depth-of-field than conventional or environmental scanning electron microscopes. Here, we offer first-principles theoretical insights to advance ion-beam imaging of atomically thin materials by performing real-time time-dependent density functional theory simulations of single impacts of 10–200 keV light ions in free-standing graphene. Here we predict that detecting electrons emitted from the back of the material (the side from which the ion exits) would result in up to three times higher signal and up to five times higher contrast images, making 2D materials especially compelling targets for ion-beam microscopy. This predicted superiority of exit-side emission likely arises from anisotropic kinetic emission. The charge induced in the graphene equilibrates on a sub-fs time scale, leading to only slight disturbances in the carbon lattice that are unlikely to damage the atomic structure for any of the beam parameters investigated here.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Nutrient Cycling at a Coastal Terrestrial‐Aquatic Interface

Abstract The complex interactions among soil, vegetation, and site hydrologic conditions driven by precipitation and tidal cycles control the biogeochemical transformations and bi‐directional exchange of carbon and nutrients across the terrestrial–aquatic interfaces (TAIs) in coastal regions. This study uses a highly mechanistic model, Advanced Terrestrial Simulator (ATS)‐PFLOTRAN, to explore how these interactions affect exchanges of materials and carbon and nitrogen cycling. We used a transect in the Chesapeake Bay region that spans zones of open water, coastal wetland, transition, and upland forest. We designed several simulation scenarios to parse the effects of the individual controlling factors and the sensitivity of carbon cycling to reaction rate parameters derived from laboratory experiments. Our simulations reveal an active zone for carbon cycling under the transition zones between the wetland and the upland. Evapotranspiration is found to enhance the exchange fluxes between the surface and subsurface domains, resulting in a higher dissolved oxygen concentration in the TAIs. The transport of organic carbon derived from plant leaves and roots provide an additional source of organic carbon needed for the aerobic respiration and denitrification processes in the TAIs. The variability in reaction rate parameters associated with microbial activities is also found to play a dominant role in controlling the heterogeneity and dynamics of the simulated redox conditions. This modeling‐focused exploratory study enabled us to better understand the complex interactions among soil, water and microbes that govern the hydro‐biogeochemical processes at the TAIs, which is an important step toward representing coastal ecosystems in larger‐scale Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Meteorological drivers of the low-cloud radiative feedback pattern effect and its uncertainty

The radiative feedback pattern effect remains a large source of uncertainty for both projections of future trends and interpretations of past trends in global temperature. The pattern effect is defined as the difference in feedbacks between transient and long-term simulations, and past work shows that is primarily attributed to changes in the marine low-cloud radiative feedback. Here we use low cloud meteorological kernels to map out both the primary cloud controlling factors through which changing surface temperature patterns drive changes in low-cloud feedback, as well as the sources of model spread. We find that the pattern effect is almost entirely driven by changes in estimated inversion strength (EIS) in the Southern Hemisphere, particularly in the South East Pacific and Southern Ocean. In both past and future simulations, inter-model spread is primarily caused by model differences in the sensitivity of low clouds to the environmental conditions, rather than differences in the simulated evolution of environmental conditions.

58 GEOSCIENCES↗
DOE OSTI2021

Cyber Physical Grid-Interactive Distributed Energy Resources Control for VPP Dispatch and Regulation: Preprint

This paper presents a cyber-physical algorithm for grid interactive DER control to enable two features of Virtual Power Plants (VPPs) dispatch and grid voltage regulation, considering the communication and security impacts. We first formulate the DER dispatch problem as a real-time, iterative, and grid-interactive DER control problem. Thereafter, we consider a probabilistic traffic model to characterize packet delays and loss in a communication network, and study how the delays enter the process of information exchange among the grid measurement units, local DER controllers and the grid control center that coordinately execute this dispatch algorithm. Finally, we propose a cyber-physical grid-interactive DER control algorithm using the previous message strategy. The tracking and regulation capabilities of this proposed algorithm can be made immune to the asynchrony resulting from the communications network traffic. We carry out simulations to show possible numerical instabilities and sensitivities of the tracking and regulation capabilities on the proposed strategy. Our results exhibit that the uncertainties of the underlying communications infrastructure must be considered in the control algorithm for the VPP tracking and regulation capabilities of any DER in a generic Cyber-Physical System (CPS).

co-simulation↗
DOE OSTI2025

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗
DOE OSTI2023

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗
DOE OSTI2026

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)↗
DOE OSTI2023

Sensitivity Analysis of Genome-Scale Metabolic Flux Prediction

TRIMER, Transcription Regulation Integrated with MEtabolic Regulation, is a genome-scale modeling pipeline targeting at metabolic engineering applications. Using TRIMER, regulated metabolic reactions can be effectively predicted by integrative modeling of metabolic reactions with a Transcription Factor (TF)-gene regulatory network (TRN), which is modeled via a Bayesian network (BN). In this paper, we focus on sensitivity analysis of metabolic flux prediction for uncertainty quantification of BN structures for TRN modeling in TRIMER. We propose a computational strategy to construct the uncertainty class of TRN models based on the inferred regulatory order uncertainty given transcriptomic expression data. With that, we analyze the prediction sensitivity of the TRIMER pipeline for the metabolite yields of interest. The obtained sensitivity analyses can guide Optimal Experimental Design (OED) to help acquire new data that can enhance TRN modeling and achieve specific metabolic engineering objectives, including metabolite yield alterations. Here we have performed small- and large-scale simulated experiments, demonstrating the effectiveness of our developed sensitivity analysis strategy for BN structure learning to quantify the edge importance in terms of metabolic flux prediction uncertainty reduction and its potential to effectively guide OED.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Informing nuclear physics via machine learning methods with differential and integral experiments

Information from differential nuclear-physics experiments and theory is often too uncertain to accurately define nuclear-physics observables such as cross sections or energy spectra. Integral experimental data, representing the applications of these observables, are often more precise but depend simultaneously on too many of them to unambiguously identify issues in the observable with human expert analysis alone. Here, we explore how we can leverage physics knowledge gained from differential experimental data, nuclear theory, integral experiments, and neutron-transport calculations to better understand nuclear-physics observables in the context of the application area represented by integral experiments. We support this task with machine-learning methods to discern trends in a large amount of convoluted data. Differential and integral information was used in an analysis augmented by the random forest and the Shapley additive explanations metric. We chose as an application area one that is represented by criticality measurements and pulsed-sphere neutron-leakage spectra. We show one representative example ( 241 Pu fission observables) where the combination of differential and integral information allowed to resolve issues in data representing these observables. As a starting point, the machine learning (ML) algorithms highlighted several observables as leading potentially to bias in simulating integral experiments. Differential information, paired with sensitivity to integral quantities, allowed us then to pinpoint one specific observable ( 241 Pu fission cross section) as the main driver of bias. The comparison to integral experiments, on the other hand, allowed us to indicate a likely reliable experiment among several discrepant ones for this observables. In other cases (e.g., 239 Pu observables), we were not able to resolve the confounding introduced by integral experiments but instead highlighted the need for targeted new experiments and theory developments to better constrain the nuclear-physics space for the application area represented by integral experiments. We were able to combine information from differential experimental data, nuclear-physics theory, integral experiments, and neutron-transport simulations of the latter experiments with the help of the random forest algorithm and expert judgment. This combination of knowledge allows to improve our description of nuclear-physics observables as applied to a particular application area.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Efficient Distance-based Global Sensitivity Analysis for Terrestrial Ecosystem Modeling

Sensitivity analysis in terrestrial ecosystem modeling is important for understanding controlling processes, guiding model development, and targeting new observations to reduce parameter and prediction uncertainty. Complex and computationally expensive terrestrial ecosystem models (TEM) limit the number of ensemble simulations, requiring sophisticated and efficient methods to analyze sensitivities of multiple model responses to different types of parameter uncertainties. In this study, we propose a distance-based global sensitivity analysis (DGSA) method. DGSA first classifies model response samples into a small set of discrete classes and then calculates the distance between parameter frequency distributions in different classes to measure the parameter sensitivity. The principle is that, if the parameter distribution is the same in each class, then the model response is insensitive to the parameter, while a large difference in the distributions indicates the parameter is influential to the response. Built on this idea, DGSA can be applied to analyze sensitivity of a single and a group of responses to different kinds of parameter uncertainties including continuous, discrete and even stochastic. Besides the main-effect sensitivity from a single parameter, DGSA can also quantify the sensitivity from parameter interactions. Additionally, DGSA is computationally efficient which can use a small number of model evaluations to obtain an accurate and statistically significant result. We applied DGSA to two TEMs, one having eight parameters and three kinds of model responses, and the other having 47 parameters and a long-period response. We demonstrated that DGSA can be used for sensitivity problems with multiple responses and high-dimensional parameters efficiently.

Lu, Dan↗
DOE OSTI2022

Rigorous justification for the space–split sensitivity algorithm to compute linear response in Anosov systems

Ruelle gave a formula for linear response of transitive Anosov diffeomorphisms. Recently, practically computable realizations of Ruelle's formula have emerged that potentially enable sensitivity analysis of certain high-dimensional chaotic numerical simulations encountered in the applied sciences. In this paper, we provide full mathematical justification for the convergence of one such efficient computation, the space–split sensitivity, or S3, algorithm. In S3, Ruelle's formula is computed as a sum of two terms obtained by decomposing the perturbation vector field into a coboundary and a remainder that is parallel to the unstable direction. Such a decomposition results in a splitting of Ruelle's formula that is amenable to efficient computation. Here, we prove the existence of the S3 decomposition and the convergence of the computations of both resulting components of Ruelle's formula.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Coupling Warm Rain With an Eddy Diffusivity/Mass Flux Parameterization: 2. Sensitivities and Comparison to Observations

Abstract This is the second part of the description and validation of the coupling of a warm‐rain microphysics scheme to a stochastic multi‐plume eddy‐diffusivity/mass‐flux parameterization. Part 1 provides the model description and initial validation of the parameterization in a single column model (SCM) against large eddy simulations. This manuscript provides additional comparisons against satellite observations of 500 unique simulation cases, an analysis of model parameter‐based sensitivities in those 500 cases, and an evaluation of the idealized LES case study paradigm for parameterization development. The SCM is forced by a weather reanalysis in the subtropical northeast Pacific region characterized by the stratocumulus‐to‐cumulus cloud transition. 18 scalar outputs are examined to identify the sensitivities to 22 of the model's parameters. The model is highly sensitive to parameters that determine the behavior of convective plumes, which collectively lead to the most sensitivity in 13 of the 18 outputs. Only a few parameters contribute sensitivities that are not highly correlated with those of other parameters. While forcing the model with two canonical field campaign cases results in similar picture of which parameters lead to high sensitivity in model outputs, the sensitivity magnitudes can be overestimated compared to when the model is forced by weather reanalyzes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Verification of Revised and Upcoming Nuclear Data Sensitivity MCNP Features [Poster]

Nuclear data is ubiquitous across nuclear applications. Improved nuclear data means more accurate simulations. Past focus on k eff caused compensating error and areas of unvalidated nuclear data. Sensitivity can be used to optimize experiments to focus on specific areas. Sensitivity capabilities must be expanded and improved to design multivariate experiments focused on nuclear data needs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Calibration of reactive burn and Jones-Wilkins-Lee parameters for simulations of a detonation-driven flow experiment with uncertainty quantification

Here, uncertainties in the explosive-specific parameters of the Jones-Wilkins-Lee (JWL) equation of state (EOS) are carefully considered in hydrodynamic simulations of an explosive experiment to minimize the error in the flow prediction. Experimental data of the leading shock position in the transverse direction over time serves as the prediction metric for quantifying simulation prediction error. The uncertainty quantification technique, global sensitivity analysis, is utilized to determine the JWL parameters to which the transverse shock propagation is most sensitive. A polynomial response surface (PRS) is constructed in the space of the most influential JWL parameters, and the point of minimum error between the experimental data and the PRS yields calibrated JWL parameters for the experimental flow. The simulation results following the parameter calibration show good agreement with the experimental data. It was found that two significant parameters, the heat release per unit mass of reactant Q and JWL model exponent R 1 are strongly related, which makes it difficult to identify accurate values.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Performance assessment of near-fault buildings subjected to physics-based simulated earthquake ground motions with fling step

The effects of the co-seismic static offset (known as fling step) and associated velocity pulses on civil structures have been difficult to study because the static offset is typically removed during the processing of earthquake ground motion records. Simulated ground motions contain fling features and require no processing; therefore, they create new opportunities for representing fling features in seismic hazard analysis and assessing their influence on the seismic demands on near-fault structures. We use physics-based fault rupture simulations to study the characteristics of ground motions with fling step and the sensitivity of the near-fault structural demands to strong fling features. We uncover that simulated ground motions with a large fling step tend to have higher spectral intensity than those without a fling step at the same rupture distance, especially at periods longer than 2 s. As a result, the structural demands on flexible buildings tend to be the most sensitive to the fling features. Statistical analysis suggests that the ground motion spectral shape (represented by spectral accelerations at multiple periods) is—in most cases—a sufficient predictor of the structural demands on near-fault low-rise and mid-rise buildings at locations that are susceptible to strong fling effects. Finally, ground motion record selection experiments reveal that representing the spectral shape features at periods that are most relevant to a given structure may be an effective strategy to reduce the bias in the estimated demands on near-fault long-period structures when the available database of records is considered deficient in fling features.

Fling step↗
DOE OSTI2020

Effects of Cloud Liquid-Phase Microphysical Processes in Mixed-Phase Cumuli Over the Tibetan Plateau

Numerical simulations often overpredict precipitation over the Tibetan Plateau (TP). To examine the factors causing precipitation overprediction, different parameterizations of liquid-phase microphysical processes (accretion, autoconversion, and entrainment-mixing) are implemented into the Morrison microphysics scheme to simulate a TP precipitation event in summer with the Weather Research and Forecasting (WRF) model. The general spatial distribution and temporal trend of precipitation are captured by all simulations, but the precipitation rate is overpredicted. The results from sensitivity experiments suggest that compared to other examined liquid-phase processes, the accretion process is more important in precipitation simulation over the TP region. Further investigation with the Heidke skill scores reveals that accretion parameterization that takes into account the raindrop size produces the most accurate results in terms of the total surface precipitation. This parameterization suppresses spurious accretion and does not produce liquid-phase precipitation until cloud droplets are big enough. Here, it is also confirmed that increasing the model resolution can reduce precipitation overprediction. Results from the case study are confirmed by the use of a one-month simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Droughts in Wind and Solar Power: Assessing Climate Model Simulations for a Net‐Zero Energy Future

Abstract Understanding and predicting “droughts” in wind and solar power availability can help the electric grid operator planning and operation toward deep renewable penetration. We assess climate models' ability to simulate these droughts at different horizontal resolutions, ∼100 and ∼25 km, over Western North America and Texas. We find that these power droughts are associated with the high/low pressure systems. The simulated wind and solar power variabilities and their corresponding droughts during historical periods are more sensitive to the model bias than to the model resolution. Future climate simulations reveal varied future change of these droughts across different regions. Although model resolution does not affect the simulation of historical droughts, it does impact the simulated future changes. This suggests that regional response to future warming can vary considerably in high‐ and low‐resolution models. These insights have important implications for adapting power system planning and operations to the changing climate.

Liu, Xue [Department of Oceanography Texas A&amp,M↗