Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

What Have We Learned from MERRA about Reanalyses of the Stratosphere?

The Modern-Era Retrospective analysis for Research and Applications (MERRA) reanalysis includes a well-resolved middle atmosphere. The upper boundary of the underlying model is in the mesosphere, near 80km, and the input data streams include the Stratospheric Sounding Unit (SSU) and Advanced Microwave Sounding Unit (AMSU) radiance observations. These two datasets provide observational constraints on the deep-layer thermal structure in approximately the 10-2hPa region, the middle to upper stratosphere, which is above the highest range of most radiosonde ascents. This analysis will focus on the difficulties of producing realistic analyses in the middle to upper stratosphere: these arise largely because of vertical averaging inherent in the AMSU and SSU observations, the sensitivity to model biases in this region, the relative biases among the same channels on different instruments, and the orbital sampling of the satellites (morning or afternoon orbits). These issues will be illustrated with examples from MERRA and enhanced by discussions of potential ways of improving the middle atmosphere in future reanalyses.

Pawson, Steven↗
NASA NTRSDate not supplied

Global Sensitivity Analyses for Test Planning with Black-Box Models for Mars Sample Return

This work describes sensitivity analyses performed on complex black-box models used to support experimental test planning under limited resources in the context of the Mars Sample Return program, which aims at bringing to Earth rock and atmospheric samples from Mars. We develop a systematic workflow that allows the analysts to simultaneously obtain quantitative insights on key drivers of uncertainty, on the direction of impact, and the presence of interactions. We apply novel optimal transport-based global sensitivity measures to tackle the multivariate nature of the output. On the modeling side, we apply multi-fidelity techniques that leverage low-fidelity models to speed up the calculations and make up for the limited amount of high-fidelity samples, while keeping these in the loop for accuracy guarantees. The sensitivity analysis reveals insights useful for the analysts to understand the model's behavior and identify the factors to focus on during testing in order to maximize the value of information extracted from them to ensure mission success when limited resources are available.

Giuseppe Cataldo↗
NASA NTRS2010

Sensitivity Analysis of Launch Vehicle Debris Risk Model

As part of an analysis of the loss of crew risk associated with an ascent abort system for a manned launch vehicle, a model was developed to predict the impact risk of the debris resulting from an explosion of the launch vehicle on the crew module. The model consisted of a debris catalog describing the number, size and imparted velocity of each piece of debris, a method to compute the trajectories of the debris and a method to calculate the impact risk given the abort trajectory of the crew module. The model provided a point estimate of the strike probability as a function of the debris catalog, the time of abort and the delay time between the abort and destruction of the launch vehicle. A study was conducted to determine the sensitivity of the strike probability to the various model input parameters and to develop a response surface model for use in the sensitivity analysis of the overall ascent abort risk model. The results of the sensitivity analysis and the response surface model are presented in this paper.

Gee, Ken↗
NASA NTRS1989

Accurate calculation of control-augmented structural eigenvalue sensitivities using reduced-order models

A method is presented for generating mode shapes for model order reduction in a way that leads to accurate calculation of eigenvalue derivatives and eigenvalues for a class of control augmented structures. The method is based on treating degrees of freedom where control forces act or masses are changed in a manner analogous to that used for boundary degrees of freedom in component mode synthesis. It is especially suited for structures controlled by a small number of actuators and/or tuned by a small number of concentrated masses whose positions are predetermined. A control augmented multispan beam with closely spaced natural frequencies is used for numerical experimentation. A comparison with reduced-order eigenvalue sensitivity calculations based on the normal modes of the structure shows that the method presented produces significant improvements in accuracy.

Livne, Eli↗
DOE OSTI2026

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up↗
DOE OSTI2026

Life cycle analysis of dedicated energy crops for fuel production in the United States

Dedicated energy crops are promising feedstocks to make biofuels including jet fuels. This study applies life cycle analysis (LCA) to estimate direct well-to-wake (WTW) greenhouse gas (GHG) emissions (g CO 2 e/MJ) for jet fuel derived from five energy crops-biomass sorghum, miscanthus, switchgrass, poplar, and willow-via Fischer-Tropsch-to-Jet (FTJ) and Ethanol-to-Jet (ETJ) pathways. The WTW boundary includes direct emissions from biomass production, fuel production, and fuel combustion. The R&D GREET model is expanded to conduct the LCA, using national average biomass yields and farming inputs from the 2023 Billion-Ton Study. In addition, this study estimates emissions from market-mediated effects, including induced land-use change, induced other crop (non-feedstock) production changes, and induced livestock production changes using global economic and emissions factor models. On a per-dry U.S. ton basis, cultivation and harvest emissions are lowest for willow (51,565 g CO 2 e) and highest for biomass sorghum (104,488 g CO 2 e). Per-acre results show similarly high emissions for sorghum and lowest values for poplar and willow. Direct WTW emissions are substantially lower for FTJ (biomass sorghum: 5.5; miscanthus: 10.3; switchgrass: 11.7; poplar: 11.9; and willow: 8.7 g CO 2 e/MJ) than ETJ (33.2; 33.8; 34.8; 36.2; and 31.7 g CO 2 e/MJ, respectively). When market-mediated emissions are included, miscanthus exhibits the lowest total emissions across energy crop pathways. Although results are sensitive to modeling assumptions, they indicate that high-yielding perennial and woody crops, particularly when planted on marginal land, could significantly reduce WTW emissions for bio-jet fuels by combining low direct emissions with soil carbon gains and favorable market-mediated effects.

Billion-Ton Study↗
NASA NTRS1980

The vertical distribution of tropospheric ammonia

A one-dimensional tropospheric photochemical model is used to simulate measured profiles of NH3 obtained with the Infrared Heterodyne Radiometer. The relative roles of homogeneous loss, heterogeneous loss, and vertical eddy transport are discussed in terms of selecting parameters which best fit the measurements. The best fit was obtained for a vertical eddy diffusion coefficient of 200,000/sq cm per sec or greater (corresponding to a characteristic vertical transport time in excess of about 35 days), and a characteristic heterogeneous loss time in excess of 10 days. The characteristic homogeneous chemical loss time was found to be about 40 days at the surface and decreased to about 180 days at 10 km, and not very sensitive to model chemical perturbations. Increased ground-level concentrations of NH3 to about 10 ppb, compared to background surface concentrations of about 1 ppb, were measured several weeks after application of ammonium nitrate fertilizer. This suggests that the volatilization of ammonium nitrate fertilizer is rapid, and an important source of NH3. Because of the characteristic times for the loss mechanisms, synoptic time-scale phenomena may play an important role in determining the tropospheric distribution of NH3 concentrations.

Levine, J. S.↗
NASA NTRS1993

Robust on-off pulse control of flexible space vehicles

The on-off reaction jet control system is often used for attitude and orbital maneuvering of various spacecraft. Future space vehicles such as the orbital transfer vehicles, orbital maneuvering vehicles, and space station will extensively use reaction jets for orbital maneuvering and attitude stabilization. The proposed robust fuel- and time-optimal control algorithm is used for a three-mass spacing model of flexible spacecraft. A fuel-efficient on-off control logic is developed for robust rest-to-rest maneuver of a flexible vehicle with minimum excitation of structural modes. The first part of this report is concerned with the problem of selecting a proper pair of jets for practical trade-offs among the maneuvering time, fuel consumption, structural mode excitation, and performance robustness. A time-optimal control problem subject to parameter robustness constraints is formulated and solved. The second part of this report deals with obtaining parameter insensitive fuel- and time- optimal control inputs by solving a constrained optimization problem subject to robustness constraints. It is shown that sensitivity to modeling errors can be significantly reduced by the proposed, robustified open-loop control approach. The final part of this report deals with sliding mode control design for uncertain flexible structures. The benchmark problem of a flexible structure is used as an example for the feedback sliding mode controller design with bounded control inputs and robustness to parameter variations is investigated.

Wie, Bong↗
DOE OSTI2026

Bayesian Adaptive Polynomial Chaos Expansions

Polynomial chaos expansions (PCEs) are widely used for uncertainty quantification (UQ) tasks, particularly in the applied mathematics community. However, PCE has received comparatively less attention in the statistics literature, and fully Bayesian formulations remain rare—especially with implementations in R. Motivated by the success of adaptive Bayesian machine learning models such as BART, BASS and BPPR, we develop a new fully Bayesian adaptive PCE method with an efficient and accessible R implementation: khaos. Our approach includes a novel proposal distribution that enables data-driven interaction selection and supports a modified g-prior tailored to PCE structure. Through simulation studies and real-world UQ applications, we demonstrate that the Bayesian adaptive PCE provides competitive performance for surrogate modeling, global sensitivity analysis and ordinal regression tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Probabilistic Inference of Low-Surface-Brightness Galaxy Morphological Parameters Using Simulation-Based Inference

Low-surface-brightness galaxies (LSBGs) are diffuse, often dark-matter-dominated systems whose faintness makes their structural parameters difficult to measure reliably in wide-field imaging surveys. Robust parameter inference, including uncertainty quantification, is important for population studies and for comparisons with models of galaxy formation, as future surveys are expected to produce increasingly large samples of diffuse galaxies. In practice, LSBG profile modeling is sensitive to sky- background errors, masking choices, contaminating background sources, and the computational cost of obtaining posterior-level uncertainties for large samples. Motivated by these questions, we develop a simulation-based inference (SBI) framework for estimating posterior distributions of LSBG morphological parameters from simulated galaxy images. Using PyImfit, we generate DES-like single-Sersic profile LSBG images with known position angle, ellipticity, Sersic index, effective surface brightness, and effective radius. We then train a normalizing-flow-based neural posterior estimator using the sbi package to infer these parameters from the simulated images. For isolated simulated galaxies, the SBI posterior recovers the true input parameters, produces posterior predictive residuals consistent with the assumed noise model, and shows good empirical calibration in a DES-motivated test regime. We also compare SBI with PyImfit-based MCMC inference and find broadly comparable posterior constraints, while SBI enables substantially faster posterior sampling after training. Finally, we test robustness to compact background contaminants. A model trained only on isolated galaxies produces undercovered posteriors on contaminated images, whereas training on simulations with variable contaminant positions and fluxes improves calibration across contaminated test sets. These results demonstrate the promise of SBI for scalable, uncertainty-aware LSBG morphology inference, while emphasizing that posterior reliability strongly depends on whether training simulations include relevant observational complications.

Batbayar, Bilguun [U. Chicago (main)]↗
NASA NTRS1995

On treatments of fetch and stability sensitivity in large-area estimates of sensible heat flux over sea ice

Bulk transfer coefficients estimated as a function of atmospheric stability and sea-ice lead width are combined with Arctic meteorological observations and ice thickness data to calculate the sensitivity of turbulent flux estimates to changes in lead width, wind speed, air temperature, and ice thickness for a high-concentration ice pack. These results are considered in terms of bulk transfer parameterizations that use a fixed transfer coefficient r that address atmospheric stability only. On the basis of the fetch-sensitive parameterizations considered here, differences in lead width for widths up to about 200 m can exert a substantial influence on sensible heat transfer coefficients and heat flux from leads under typical Arctic conditions. Fluxes from an open water lead decrease by 34% if fetch increases from 10 m to 100 m. This effect is greatest for open water leads, decreases considerably as leads refreeze, and is negligible for ice thicker than about 0.3 m. If open or newly refrozen leads make up 2% of the ice cover, than an increase in mean fetch from 10 m to 100 m yields a decrease of about 2 W/sq m in areally averaged flux from the ice pack. Calculations using observed and theoretical lead width distributions suggest that parameterizing lead widths in a sea ice model can be done effectively using a single, representative lead width rather than requiring a full distribution of widths. When coupled to the lowest atmospheric boundary layer using a bulk similarity theory model, this sensitivity of heat transfer to fetch results in substantially higher near-surface air temperatures over narrow leads, with equilibrium air temperatures decreasing by about 50% as fetch increases from 10 to 100 m.

Maslanik, J. A.↗
NASA NTRS1982

A preliminary study of the application of HCMM satellite data to define initial and boundary conditions for numerical models: A case study in St. Louis, Missouri

Infrared and visible HCMM data were used to examine the potential application of these data to define initial and boundary conditions for mesoscale numerical models. Various boundary layer models were used to calculate the distribution of the surface heat flux, specific humidity depression (the difference between the specific humidity in the air at approxmately the 10 m level and the specific humidity at the ground), and the eddy vicosity in a 72 km by 72 km area centered about St. Louis, Missouri. Various aspects of the implications of the results on the meteorology of St. Louis are discussed. Overall, the results indicated that a reasonable estimate of the surface heat flux, urban albedo, ground temperature, and specific humidity depression can be obtained using HCMM satellite data. Values of the ground-specific humidity can be obtained if the distribution of the air-specific humidity is available. More research is required in estimating the absolute magnitude of the specific humidity depression because calculations may be sensitive to model parameters.

Vukovich, F. M.↗
DOE OSTI2025

Degrees of rate control and AutoDiff-driven direct sensitivity analysis in heterogeneous catalysis

Despite the wide application and benefits of the degree of rate control (DRC) analysis, several details remain argued, particularly about the conservation of DRCs at transient (TR) and steady-state (SS) conditions, especially for complex reaction networks. This work argues that previous proofs about the conservation properties of DRCs have been incomplete, and we provide new mathematical proofs at TR and SS conditions. In addition, we use both analytical (automatic differentiation) and numerical (finite difference) approaches to compute DRCs for the case study of ethane hydrogenolysis (EH) over Pt(111). This work confirms that at both TR and SS conditions, the sum of all DRCs, i.e., sum of the degrees of kinetic (DKRC) and thermodynamic rate control (DTRC), is conserved at zero. At SS conditions, the sum of DKRC is conserved at 1 while the sum of DTRC is conserved at −1. In corroboration of previous works, we show that the DTRC for any adsorbate at SS is equal to the product of the species coverage and a constant. In contrast, at TR conditions, the individual sums of both DTRC and DKRC are not conserved and can be any real number, with potential implications for the novel field of dynamic catalysis. Finally, we show that the conventional finite difference (FD) approach, only useful at SS, is prone to inaccuracy and very sensitive to the value of the differential change applied. The optimal differential value also varies significantly with system and rate definition. Consequently, we describe and illustrate in this work the application of the automatic differentiation (AD) approach for the more accurate determination of DRCs at both TR and SS conditions.

Automatic differentiation↗
DOE OSTI2026

Response of marine post-frontal clouds to Gulf Stream variability

Understanding how Gulf Stream variation influences cloud morphology is critical for evaluating cloud feedback in the western North Atlantic Ocean and beyond, where mesoscale air-sea interactions dominate. This study investigates the impact of altered mean sea surface temperature (SST) and SST gradients on post-frontal cloud characteristics during cold-air outbreaks, using the Weather Research and Forecasting (WRF) model. Three sensitivity experiments are conducted: a control simulation (default SST), Plus4 (uniform SST increase of 4 K), and Gradplus (SST gradient enhanced by 25 %, centered around mean SST). Results reveal distinctly different responses in boundary layer dynamics and cloud macro-physics. In Plus4, a warmer and moister boundary layer reduces total cloud cover but promotes larger cloud sizes and elongated cloud streets, with diminished liquid water and enhanced ice-phase hydrometeors. Conversely, Gradplus amplifies impacts in the upwind colder SST regions, yielding a drier, colder boundary layer, weaker energy transport, and higher liquid water path but reduced ice water content and cloud lines. Tracer analysis highlights that SST modifications alter airmass sources near cloud tops due to the entrainment of ambient air, with Plus4 amplifying boundary layer contributions to cloud-top regions. These findings underscore the spatially varying effects of SST gradients and mean SST on cloud organization and microphysics, emphasizing the need to resolve ocean-atmosphere coupling in global models to improve the prediction of marine cloud feedback under warming scenarios.

AIR-SEA INTERACTION↗
NASA NTRS1995

The stellar content of 30 doradus derived from spatially integrated ultraviolet spectra: A test of spectral synthesis models

Using the IUE satellite, we have obtained spatially integrated ultraviolet spectra of three areas within the giant H II region 30 Dor in the Large Magellanic Cloud. The spectra correspond to spatial reginswith sizes of 20 sec x 20 sec, 1 min x 1 min, and 3 min x 3 min, all of which are approximately centered on R136. We have performed a spectral synthesis analysis of the spectra of the two larger regions and compared the results with the known stellar content in these regions. The spectral synthesis models are sensitive to the ultraviolet continuum level, the P Cygni profile of the C Iv wavelength 1550 line, the absorption strength of the Si IV wavelength 1400 line, and the emission strength of the He II wavelength 1640 line. The intrinsic continuum levels and the profiles of these stellar wind lines provide constraints on the age and duration of the starburst episode within a region, as well as on the upper curoff mass of the initial mass function. From our analysis we find that the present-day value of the upper cutoff mass in the 1 min x 1 min and 3 min x 3 min regions has a lower limit of approximately 50 solar mass, a result which is in good agreement with several other recent determinations. The age of the starburst episode must be less than approximately 3 Myr, also in agreement with other estimates. Comparison of the observed total numbers of O and W-R stars with those predicted from the various models favors an instantaneous burst of star formation in the regions. However, the differences between the two burst scenarios we investigated (instantaneous and continuous) are small at such a young age, and distinguishing between the two is difficult. We are now confident that these spectral synthesis models can be used to determine the stellar content of more distant star-forming regions.

Vacca, William D.↗
NASA NTRS1991

Sensitivity-based scaling for correlating structural response from different analytical models

This paper presents a sensitivity-based linearly varying scale factor used to reconcile results from simple and refined models for analysis of the same structure. The improved accuracy of the linear scale factor compared to a constant scale factor as well as the commonly used tangent approximation is demonstrated. A wing box structure is used as an example, with displacements, stresses and frequencies correlated. The linear scale factor could permit the use of a simplified model in an optimization procedure during preliminary design to approximate the response given by a refined model over a considerable range of design changes.

Chang, Kwan J.↗
DOE OSTI2025

Regional Oil and gas Aerial Methane Synthesis model (ROAMS) v2.0

The Regional Oil and gas Aerial Methane Synthesis model is a tool to convert the results of wide-area, source-resolved aerial methane remote sensing surveys of oil and natural gas infrastructure in a given region into methane emissions inventories (estimates of the magnitude and breakdown of methane emissions from the surveyed infrastructure). The tool leverages databases of source-resolved methane emissions detected in aerial surveys, aerial survey coverage information (which areas were measured and when), data summarizing surveyed oil and natural gas infrastructure and production (derived from third-party databases), as well as state-of-the-art mechanistic emissions simulation tools to characterize emissions too small for the aerial system to see. The regional methane emissions estimates produced by this tool are much more granular in both space and asset type than common satellite- or flux tower-based regional estimates. Unlike other tools for converting site-level measurements into regional emissions estimates, our unique geostatistical approach integrates aerially measured emissions with limited need for statistical extrapolation, which can be highly sensitive to modeler assumptions. As a result, ROAMS-based estimates of regional methane emissions from oil and gas activity are widely viewed as highly credible, as evidenced by the success of Dr. Sherwin's recent paper in Nature.

Sherwin, Evan [Lawrence Berkeley National Laborato↗
DOE OSTI2025

Regional Oil and gas Aerial Methane Synthesis model (Analytica) (ROAMS Analytica) v1.5.2

The Regional Oil and gas Aerial Methane Synthesis model (Analytica) is a tool to convert the results of wide-area, source-resolved aerial methane remote sensing surveys of oil and natural gas infrastructure in a given region into methane emissions inventories (estimates of the magnitude and breakdown of methane emissions from the surveyed infrastructure). This version is written in the Analytica programming language, and this version accompanies a correction in preparation for submission to Sherwin et al. 2024 (Nature). The tool leverages databases of source-resolved methane emissions detected in aerial surveys, aerial survey coverage information (which areas were measured and when), data summarizing surveyed oil and natural gas infrastructure and production (derived from third-party databases), as well as state-of-the-art mechanistic emissions simulation tools to characterize emissions too small for the aerial system to see. The regional methane emissions estimates produced by this tool are much more granular in both space and asset type than common satellite- or flux tower-based regional estimates. Unlike other tools for converting site-level measurements into regional emissions estimates, our unique geostatistical approach integrates aerially measured emissions with limited need for statistical extrapolation, which can be highly sensitive to modeler assumptions. As a result, ROAMS-based estimates of regional methane emissions from oil and gas activity are widely viewed as highly credible, as evidenced by the success of Dr. Sherwin's recent paper in Nature.

Sherwin, Evan [Lawrence Berkeley National Laborato↗