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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

42 ENGINEERING↗
NASA NTRSDate not supplied

Turbulence Effects on Shaped Booms: Propagation Simulations Using KZKFourier

Upcoming X-59 aircraft flight tests as part of NASA’s Quesst Mission are expected to occur in a range of atmospheric conditions. Sensitivity of ground waveform acoustic metrics to turbulent perturbations during propagation complicates the determination of noise levels. A series of simulations through turbulence was executed using the KZKFourier model of Stout et al. to develop a database of results for expanding the functionality of NASA tools for estimating turbulence effects on shaped-boom ground waveforms. For 45 cases covering a seven-factor design space, multiple realizations of turbulence were generated to characterize statistical turbulence effects on levels of six acoustic metrics. An individual simulation in the study produced over one thousand waveforms across a virtual microphone array. Different measures were used to evaluate refinement of results with increasing number of simulations. Data suggest that mean effects and variability were most strongly influenced by propagation distance and velocity fluctuation intensity. Although local increases in acoustic metrics were common, the overall average result in all cases was a reduction in metric levels. Sensitivity of acoustic metrics varied, with mean reductions of up to 2 dB in Perceived Level occurring across the 45 cases.

turbulence↗
DOE OSTI2025

Clumping index estimation with 30°-tilted cameras in row crops: Evaluation of methods and segment size effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Furthermore, applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

30°-tilted camera↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗
DOE OSTI2026

Catalytic resonance theory for the kinetic signatures of multiple wavelength photocatalysis

Multi-step surface chemistry promoted with one or two wavelengths of light exhibited distinct kinetic signatures indicative of the light-adsorbate interactions and the number of photon-sensitive elementary steps. In this work, kinetic Monte Carlo simulations identified the kinetic response of general surface mechanisms to variation in the per-site photon flux of one or two wavelengths of incident light. Photocatalytic rates were described via a non-dimensional photon flux, identifying multiple kinetic regimes with unique degrees of rate control. Under photon-controlled kinetic conditions, maximum quantum yield and turnover frequency were obtained under constant illumination, while pulsing of one or more light sources was shown to exhibit slower rates and less efficient photocatalytic promotion due to the inherent dynamic nature of light, which comprises a stream of photons. Simulations provided distinct kinetic signatures in catalytic rates and Arrhenius kinetics for specific light-surface-adsorbate interactions corresponding to photocatalytic promotion of specific steps in surface chemistry.

catalytic resonance↗
DOE OSTI2025

Rising atmospheric moisture escalates the future impact of atmospheric rivers in the Antarctic climate system

Abstract Despite their relative rarity, atmospheric rivers are key contributors to the surface mass balance of Antarctica. However, the future role of atmospheric rivers in modulating Antarctic sea-level contributions is a major area of uncertainty. Here, we leverage high-resolution climate simulations to show that Antarctic atmospheric rivers are highly sensitive to future increases in atmospheric moisture, leading to a doubling of atmospheric river frequencies and 2.5 × increase in precipitation from 2066–2100 under present-day thresholds for atmospheric river detection. However, future precipitation impacts are critically dependent on the detection threshold: accounting for moisture increases in the threshold produces smaller, regional changes in atmospheric river frequency, primarily resulting from an eastward shift in the polar jet maximum wind speeds. Our results underscore the importance of using large ensembles to quantify Antarctic atmospheric river responses to variability in projected moisture, which may not be captured when using only a few ensemble members.

Environmental Sciences & Ecology↗
DOE OSTI2023

Effects of the heat source on the steady-state transport in gradient-driven global gyrokinetic simulations

Heat sources are routinely used in gradient-driven global gyrokinetic simulations to prevent profile relaxation and ensure steady-state turbulence and transport. To have any predictive capability, gradient-driven simulations need to be insensitive to the particular choice of a source. Several forms of the heat source are investigated using gyrokinetic δf particle-in-cell simulations. The ion heat flux is found to be sensitive to the source strength when the source rate is below a fraction (∼1/10) of the linear growth rate.

Physics↗
DOE OSTI2025

Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators

We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high k modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

New Approach Methodology for Assessing Inhalation Risks of a Contact Respiratory Cytotoxicant: Computational Fluid Dynamics-Based Aerosol Dosimetry Modeling for Cross-Species and In Vitro Comparisons

Regulatory agencies are considering alternative approaches to assessing inhalation toxicity that utilizes in vitro studies with human cells and in silico modeling in lieu of additional animal studies. In support of this goal, computational fluid-particle dynamics models were developed to estimate site-specific deposition of inhaled aerosols containing the fungicide, chlorothalonil, in the rat and human for comparisons to prior rat inhalation studies and new human in vitro studies. Under bioassay conditions, the deposition was predicted to be greatest at the front of the rat nose followed by the anterior transitional epithelium and larynx corresponding to regions most sensitive to local contact irritation and cytotoxicity. For humans, simulations of aerosol deposition covering potential occupational or residential exposures (1–50 µm diameter) were conducted using nasal and oral breathing. Aerosols in the 1–5 µm range readily penetrated the deep region of the human lung following both oral and nasal breathing. Under actual use conditions (aerosol formulations >10 µm), the majority of deposited doses were in the upper conducting airways. Beyond the nose or mouth, the greatest deposition in the pharynx, larynx, trachea, and bronchi was predicted for aerosols in the 10–20 µm size range. Only small amounts of aerosols >20 µm penetrated past the pharyngeal region. Using the ICRP clearance model, local retained tissue dose metrics including maximal concentrations and areas under the curve were calculated for each airway region following repeated occupational exposures. These results are directly comparable with benchmark doses from in vitro toxicity studies in human cells leading to estimated human equivalent concentrations that reduce the reliance on animals for risk assessments.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Surface acoustic waves and lattice vibrations in two-dimensional Ti 3 ⁢C 2 ⁢𝑇 𝑥 (𝑇 =O, F) MXene films

We investigated surface acoustic wave (SAW) propagation and lattice vibrations in two-dimensional (2D) titanium carbide (Ti 3 ⁢C 2 ⁢𝑇 x ) MXene films as a function of surface termination and layer stacking, using atomistic simulations. We found that SAW propagation velocity is highly sensitive to both single-layer properties and interlayer bonding. Surface terminations significantly modulate wave behavior, with oxygen and fluorine terminations producing distinct effects on wave propagation, with oxygen-terminated monolayers exhibiting 20% higher wave speeds than fluorine counterparts due to strengthened intralayer bonds. Key observations include the transition from one to two layers causing wave speed variations, and the development of interlayer modes that generate more dispersed lattice vibrations. As the film layer thickness increases, SAW propagation becomes predominantly confined to the upper surface, with coherence of vibrational modes diminishing in multilayer structures. These findings suggest MXene terminations and layer stacking are crucial parameters for controlling SAW behavior, offering promising avenues for novel acoustic wave device applications.

2-dimensional systems↗
DOE OSTI2020

High-temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. Finally, this parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The physiological basis for estimating photosynthesis from Chl>em>a fluorescence

We report solar-induced Chl fluorescence (SIF) offers the potential to curb large uncertainties in the estimation of photosynthesis across biomes and climates, and at different spatiotemporal scales. However, it remains unclear how SIF should be used to mechanistically estimate photosynthesis. In this study, we built a quantitative framework for the estimation of photosynthesis, based on a mechanistic light reaction model with the Chla fluorescence of Photosystem II (SIFPSII) as an input (MLR-SIF). Utilizing 29 C 3 and C 4 plant species that are representative of major plant biomes across the globe, we confirmed the validity of this framework at the leaf level. The MLR-SIF model is capable of accurately reproducing photosynthesis for all C 3 and C 4 species under diverse light, temperature, and CO 2 conditions. We further tested the robustness of the MLR-SIF model using Monte Carlo simulations, and found that photosynthesis estimates were much less sensitive to parameter uncertainties relative to the conventional Farquhar, von Caemmerer, Berry (FvCB) model because of the additional independent information contained in SIF PSII . Once inferred from direct observables of SIF, SIF PSII provides ‘parameter savings’ to the MLR-SIF model, compared to the mechanistically equivalent FvCB model, and thus avoids the uncertainties arising as a result of imperfect model parameterization. Our findings set the stage for future efforts to employ SIF mechanistically to improve photosynthesis estimates across a variety of scales, functional groups, and environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Shock Tube and Flame Speed Measurements of 2,4,4-Trimethyl-1-Pentene: A Co-Optima Biofuel

Abstract The combustion of 2,4,4-trimethyl-1-pentene (diisobutylene, C8H16), which is a biofuel and a component of surrogate fuels, is examined in this work. Carbon monoxide time–histories and ignition delay times are collected behind reflected shock waves utilizing a shock tube and mid-infrared laser absorption spectroscopy. Measurements were obtained near 10 atm pressure during stoichiometric oxidation of 0.15%C8H16/O2/Ar. Simulated results from chemical kinetic models are provided, and sensitivity analyses are used to discuss differences between models for both ignition delay times and carbon monoxide formation. In addition, laminar burning speeds are obtained at 1 atm, 428 K, and equivalence ratios, phi, between 0.91 and 1.52 inside a spherical chamber facility. Measured burning speeds are found to be less than that of ethanol over the equivalence ratio span. Burning speed measurements are compared to predictions of chemical kinetic mechanisms and are in agreement for the richest conditions; however, at lean conditions, the model predicts a far slower-burning speed. The maximum burning speed occurs at an equivalence ratio of 1.08 with a magnitude of 0.70 m/s. The current work provides the crucial experimental data needed for assessing the feasibility of this biofuel and for the development of future combustion chemical kinetics models.

Energy & Fuels↗
DOE OSTI2022

Analytical Sensitivity Analysis of a Spent Nuclear Fuel Cask

Here, we report nuclear science and engineering is a field increasingly dominated by computational studies resulting from increasingly powerful computational tools. As a result, analytical studies, which previously pioneered nuclear engineering, are increasingly viewed as secondary or unnecessary. However, analytical solutions to reduced-fidelity models can provide important information concerning the underlying physics of a problem and aid in guiding computational studies. Similarly, there is increased interest in sensitivity analysis studies. These studies commonly use computational tools. However, providing a complementary sensitivity study of relevant analytical models can lead to a deeper analysis of a problem. This work provides the analytical sensitivity analysis of the one-dimensional (1D) cylindrical mono-energetic neutron diffusion equation using the forward sensitivity analysis procedure (FSAP) developed by Cacuci. Further, these results are applied to a reduced-fidelity model of a spent nuclear fuel cask, demonstrating how computational analysis might be improved with a complementary analytic sensitivity analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Sea Ice and Cloud Processes Mediating Compensation between Atmospheric and Oceanic Meridional Heat Transports across the CMIP6 Preindustrial Control Experiment

Abstract Bjerknes compensation (BJC) refers to the anticorrelation observed between atmospheric and oceanic heat transport (AHT/OHT) variability, particularly on decadal to longer time scales that may be important to the predictability of the climate system. This study investigates the spread in BJC across fully coupled simulations of phase 6 of the Coupled Model Intercomparison Project (CMIP6) and critical processes (particularly related to sea ice and clouds) that may contribute to that spread. BJC on decadal to longer time scales is confirmed across all the simulations evaluated, and it is strongest in the Northern Hemisphere (NH) between 60° and 70°N. At these latitudes, BJC appears to be primarily driven by the exchange of turbulent fluxes (sensible and latent) in the Greenland, Iceland, and Barents Seas. Metrics to break down how sea ice and clouds uniquely modify the radiative balance of the polar atmosphere during anomalous OHT events are presented. These metrics quantify the impacts of sea ice and clouds on surface and top of atmosphere (latent, sensible, longwave, and shortwave radiative) energy fluxes. Cloud responses tend to counter the clear sky impacts over the Marginal Ice Zone (MIZ). It is further shown that the degree of BJC present in a simulation at high latitudes is heavily influenced by the sensitivity of the sea ice to OHT, which is most influential over the MIZ. These results are qualitatively robust across models and explain the intermodel spread in NH BJC in the preindustrial control experiment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

On Hydrodynamics Generated from Electrostatic Discharges: A Heat Source Perspective

A parameter space map of a heat source depositing energy into a cylindrical spark channel modeling the hydrodynamic effects of an electrostatic discharge event is explored. The four parameters of spatial weight, λ R , rise time, τ r , fall time, τ f , and energy per unit length deposited in to the system, E 0 /h, were varied to observe resulting hydrodynamic phenomena. The energy budget con-sequences of the associated hydrodynamic flow are discussed in terms of the remaining electrical energy available for the joule heating of sensitive energetic materials (or “victims”) are discussed,where hydrodynamic simulations are validated by experimental data. Future results and conclusions will be added to this report as they are produced.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Streamlining the Coupling of BISON and Dakota Through the NEAMS Workbench

Metallic nuclear fuels for use in advanced reactors are an active area of research and development. Robust, accurate metallic fuel performance models are necessary for the design, analysis, and licensing of such reactors. However, metallic fuel performance models require additional development; they are not as mature as uranium dioxide fuel performance models. To support further metallic fuel development, Oak Ridge National Laboratory and the University of Florida have streamlined the coupling of the BISON fuel performance code with Design Analysis Kit for Optimization and Terascale Applications (Dakota) statistical analysis tool through the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench. This work included performing three different sensitivity analyses on metallic nuclear fuel models in BISON. The analyses examined were a general model of the IFR-1 experiment, the X430 experiment T654 pin, and the X430 experiment T651 pin. The results suggest that BISON and Dakota can be integrated through NEAMS Workbench to perform sensitivity and uncertainty analyses and visualize the results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗