Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic modelling”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Robust Path Planning and Feedback Design Under Stochastic Uncertainty

Autonomous vehicles require optimal path planning algorithms to achieve mission goals while avoiding obstacles and being robust to uncertainties. The uncertainties arise from exogenous disturbances, modeling errors, and sensor noise, which can be characterized via stochastic models. Previous work defined a notion of robustness in a stochastic setting by using the concept of chance constraints. This requires that mission constraint violation can occur with a probability less than a prescribed value.In this paper we describe a novel method for optimal chance constrained path planning with feedback design. The approach optimizes both the reference trajectory to be followed and the feedback controller used to reject uncertainty. Our method extends recent results in constrained control synthesis based on convex optimization to solve control problems with nonconvex constraints. This extension is essential for path planning problems, which inherently have nonconvex obstacle avoidance constraints. Unlike previous approaches to chance constrained path planning, the new approach optimizes the feedback gain as wellas the reference trajectory.The key idea is to couple a fast, nonconvex solver that does not take into account uncertainty, with existing robust approaches that apply only to convex feasible regions. By alternating between robust and nonrobust solutions, the new algorithm guarantees convergence to a global optimum. We apply the new method to an unmanned aircraft and show simulation results that demonstrate the efficacy of the approach.

autonomuys vehicles↗
DOE OSTI2020

A Radiative Transfer Model for Patchy Landscapes Based on Stochastic Radiative Transfer Theory

The availability of global high-resolution land cover maps provides promising a priori knowledge for characterizing subpixel heterogeneity and improving predictions of directional reflectance of coarse-resolution pixels. Due to mutual shadowing and sheltering effects between the adjacent forest and cropland patches, the spectral nonlinear mixing of patchy ecotones is significant, especially when the sun illuminates the ecotone from the forest side with high solar zenith angle. The spectral linear mixture (SLM) approach leads to overestimation of the bidirectional reflectance factor (BRF) in the red band in the principal plane (PP), with a maximum absolute error (MAE) of 0.0063 and a maximum relative error (MRE) of 52.5%, and to underestimation in the near-infrared band in PP with an MAE of 0.0940 and an MRE of 14.5%. In a scenario with randomly distributed boundary orientations, the overestimation of SLM increases with the degree of fragmentation and the view zenith angle. We propose a Radiative Transfer model for patchy ECotones (RTEC). which improves R² from 0.61 to 0.94 in the red band of Landsat-8 directional reflectance at the validation site. The RTEC model provides an efficient and analytical approach for directional reflectance predictions over heterogeneous patchy landscapes at coarse resolution and will be used for biophysical parameter retrievals [e.g., the leaf area index (LAI)] in future applications.

Zeng, Yelu↗
DOE OSTI2020

Characterization of Partially Observed Epidemics - Application to COVID-19

This report documents a statistical method for the "real-time" characterization of partially observed epidemics. Observations consist of daily counts of symptomatic patients, diagnosed with the disease. Characterization, in this context, refers to estimation of epidemiological parameters that can be used to provide short-term forecasts of the ongoing epidemic, as well as to provide gross information for the time-dependent infection rate. The characterization problem is formulated as a Bayesian inverse problem, and is predicated on a model for the distribution of the incubation period. The model parameters are estimated as distributions using a Markov Chain Monte Carlo (MCMC) method, thus quantifying the uncertainty in the estimates. The method is applied to the COVID-19 pandemic of 2020, using data at the country, provincial (e.g., states) and regional (e.g. county) levels. The epidemiological model includes a stochastic component due to uncertainties in the incubation period. This model-form uncertainty is accommodated by a pseudo-marginal Metropolis-Hastings MCMC sampler, which produces posterior distributions that reflect this uncertainty. We approximate the discrepancy between the data and the epidemiological model using Gaussian and negative binomial error models; the latter was motivated by the over-dispersed count data. For small daily counts we find the performance of the calibrated models to be similar for the two error models. For large daily counts the negative-binomial approximation is numerically unstable unlike the Gaussian error model. Application of the model at the country level (for the United States, Germany, Italy, etc.) generally provided accurate forecasts, as the data consisted of large counts which suppressed the day-to-day variations in the observations. Further, the bulk of the data is sourced over the duration before the relaxation of the curbs on population mixing, and is not confounded by any discernible country-wide second wave of infections. At the state-level, where reporting was poor or which evinced few infections (e.g., New Mexico), the variance in the data posed some, though not insurmountable, difficulties, and forecasts were able to capture the data with large uncertainty bounds. The method was found to be sufficiently sensitive to discern the flattening of the infection and epidemic curve due to shelter-in-place orders after around 90% quantile for the incubation distribution (about 10 days for COVID-19). The proposed model was also used at a regional level to compare the forecasts for the central and north-west regions of New Mexico. Modeling the data for these regions illustrated different disease spread dynamics captured by the model. While in the central region the daily counts peaked in the late April, in the north-west region the ramp-up continued for approximately three more weeks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Analysis of all-frequency variational behavior of the Kirchhoff approximation for a classic surface-scattering model

In testing a stochastic variational principle at high frequencies by using a Kirchhoffean trial function in an idealized model for surface scattering - a randomly embossed plane - we have found not only the predicted high-frequency improvement but also an unexpected low-frequency improvement in the calculated scattering amplitudes. To investigate systematically the all-frequency variational behavior, we consider here the deterministic one-boss case - Rayleigh's classic model whose exact solution is available for comparison - over all wavelengths, polarizations, and configurations of incidence and scattering. We examine analytically in particular the long-wave limit of the variational-Kirchhoff amplitudes; the results demonstrate improvements in both wavelength and angle depedence for horizontal (TM) polarization and some variational improvements for vertical (TE) polarization. This low-frequency behavior in tandem with the foreseen high-frequency improvement leads to good variational-Kirchhoff results through the intermediate resonance-frequency regime for this model.

Bird, J. F.↗
NASA NTRS1996

Indirect Identification of Linear Stochastic Systems with Known Feedback Dynamics

An algorithm is presented for identifying a state-space model of linear stochastic systems operating under known feedback controller. In this algorithm, only the reference input and output of closed-loop data are required. No feedback signal needs to be recorded. The overall closed-loop system dynamics is first identified. Then a recursive formulation is derived to compute the open-loop plant dynamics from the identified closed-loop system dynamics and known feedback controller dynamics. The controller can be a dynamic or constant-gain full-state feedback controller. Numerical simulations and test data of a highly unstable large-gap magnetic suspension system are presented to demonstrate the feasibility of this indirect identification method.

Huang, Jen-Kuang↗
DOE OSTI2022

Surrogate modeling for efficiently, accurately and conservatively estimating measures of risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. Risk measures combine objective probabilities with the subjective values of a decision maker to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical approximations obtained from the training data. These surrogate models limit over-confidence in reliability and safety assessments and produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to a stakeholder’s risk preferences and then present an approach, based on stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. Our surrogate models include biases that permit them to conservatively estimate the target risk measures. We provide theoretical results that show that these biases decay at the same rate as the L 2 error in the surrogate model. Numerical demonstrations confirm that risk-adapted surrogate models do indeed overestimate the target risk measures while converging at the expected rate.

42 ENGINEERING↗
NASA NTRS1990

A Unified Statistical Rain-Attenuation Model for Communication Link Fade Predictions and Optimal Stochastic Fade Control Design Using a Location-Dependent Rain-Statistic Database

A static and dynamic rain-attenuation model is presented which describes the statistics of attenuation on an arbitrarily specified satellite link for any location for which there are long-term rainfall statistics. The model may be used in the design of the optimal stochastic control algorithms to mitigate the effects of attenuation and maintain link reliability. A rain-statistics data base is compiled, which makes it possible to apply the model to any location in the continental U.S. with a resolution of 0-5 degrees in latitude and longitude. The model predictions are compared with experimental observations, showing good agreement.

Manning, Robert M.↗
DOE OSTI2020

Predicting demand for hydrogen station fueling

Full function hydrogen stations are a reality; fuel cell electric vehicle drivers can pull up to commercial fueling stations and receive 3–5 kg in less than 5 min, for an approximately 300-mile range. The demand for hydrogen is increasing, driven by an increase in the fueling of public and private fuel cell vehicles. This study describes the development and value of a model that simulates stochastic future demand at a hydrogen filling station. The predictive hydrogen demand model described in this article is trained from mathematical models constructed from actual hydrogen fill count, amount, and frequency data. Future fill probabilities inform the hour-by-hour demand profile and the station state of either “available, ready to fill” or “available, filling”. For example, a prediction for a station generally dispensing 5,000 kg a week on a Friday afternoon at 4 p.m. is 16 fills, totaling 48.7 kg, with a 0.52 proportion of time spent in “available, filling” state yielding 31 min of filling time. This is a first-of-its kind, published study on predicting future hydrogen demand by the time of day (e.g., hour-by-hour intervals) and day of week. This study can be used for hydrogen station requirements and operation and maintenance strategies and to assess the impact of demand variations and scenarios. Finally, this article presents the current status of hydrogen demand, the model development methods, a set of sample results. Discussion and conclusions concentrate on the value and use of the proposed model.

08 HYDROGEN↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Response to ‘Stochastic and deterministic interpretation of pool models’

Here, we concur with Azizi-Rad et al. (2021) that it is vital to critically evaluate and compare different soil carbon models, and we welcome the opportunity to further describe the unique contribution of the PROMISE model (Waring et al., 2020) to this literature. The PROMISE framework does share many features with established biogeochemical models, as our original manuscript highlighted in Table 1, and our work builds upon model innovations developed by many different groups, including that of Azizi-Rad and colleagues. Yet, the PROMISE framework is distinctive due to where it places mechanistic emphasis, and how these mechanisms are formalized in the mathematical model structure.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Thermal conductivity of heterogeneous mixtures and lunar soils

The theoretical evaluation of the effective thermal conductivity of granular materials is discussed with emphasis upon the heat transport properties of lunar soil. The following types of models are compared: probabilistic, parallel isotherm, stochastic, lunar, and a model based on nonlinear heat flow system synthesis.

Vachon, R. I.↗
NASA NTRS1977

Models for interrupted monitoring of a stochastic process

As computers are added to the cockpit, the pilot's job is changing from of manually flying the aircraft, to one of supervising computers which are doing navigation, guidance and energy management calculations as well as automatically flying the aircraft. In this supervisorial role the pilot must divide his attention between monitoring the aircraft's performance and giving commands to the computer. Normative strategies are developed for tasks where the pilot must interrupt his monitoring of a stochastic process in order to attend to other duties. Results are given as to how characteristics of the stochastic process and the other tasks affect the optimal strategies.

Palmer, E.↗
NASA NTRS1981

A comparison of helium and heavy ion spectra in He/3/ rich solar flares with a model calculation

He isotopes, O and Fe in He(3) rich solar flares are studied during the 1977 to 1979 period with the MPI/UoMd particle experiment on ISEE-1 and ISEE-3. The study revealed that the He(3) spectrum is generally harder than that of He(4), and the O spectrum is harder than that of Fe. The spectra are compared with a stationary model based on stochastic Fermi acceleration in Alfven turbulence including a rigidity dependent particle loss. Model calculations fit the He(3) and He(4) spectra between 0.4 and 4.0 MeV/nucleon, with discrepancies at lower and higher energies, and the O and Fe spectra cannot be fitted simultaneously.

Mobius, E.↗
NASA NTRS1984

On the acceleration of energetic ions in Jupiter's magnetosphere

Several aspects of the problem of high-energy ions in the Jovian magnetosphere are addressed. Voyager observations pertaining to the problem of high-energy ions in the magnetosphere are summarized, and the charge exchange emission of fast neutral sulfur and oxygen atoms and their subsequent recapture by electron impact, charge exchange, and photoionization is considered. Solutions are given to the diffusion equation assuming a source of ions injected with a gyroenergy corresponding to pickup in the middle and outer magnetosphere. It is concluded that no reasonable model parameters exist to produce the required steep spectra of the particle observations with only pickup and adiabatic radial diffusion included. A local acceleration mechanism based on nonadiabatic wave-particle interactions is needed. The assumptions and model predictions of stochastic acceleration by MHD turbulence for the Jovian magnetosphere are described. The model makes a specific correspondence between MHD wave spectrum properties and particle spectrum properties at energies above the Alfven energy.

Barbosa, D. D.↗
DOE OSTI2023

The Construction of Reduced Order Models for the HYBRID Repository

Reduced order models (ROMs) present an opportunity to reduce the computational burden of IES analysis using HYBRID models without significantly sacrificing data accuracy. Other applications of ROMs are digital twins, stochastic optimization studies, model exporting, and protection of proprietary information. In this paper the capability of creating ROMs of the models within the HYBRID repository will be demonstrated. Specifically, a linear regression ROM of a gas turbine and a dynamic mode decomposition with control (DMDc) ROM of a lithium-ion battery are constructed.

97 MATHEMATICS AND COMPUTING↗