Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic Modeling and 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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Large-scale Stochastic Optimization for Unit Commitment and Economic Dispatch

Our FY 2012 LDRD research on stochastic optimization with its application in power grid planning and operations has advanced the state-of-the-art of this area in the following aspects: (1) Deploying stochastic unit commitment on large real-world applications, (2) Stochastic optimization coupled with real-time production cost simulation (interleaved), and (3) Real-time stochastic optimization with rolling-horizon look-ahead. In conjunction, a number of computational experiments have been performed to evaluate the capability of the latest stochastic optimization software and explore opportunities for future research. In what follows, we discuss each of these contribution areas, including the results of our computational experiments and their implications on the scalability and solvability of large-scale stochastic optimization models.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1973

A guidance and navigation system for continuous low-thrust vehicles

A midcourse guidance and navigation system for continuous low thrust vehicles was developed. The equinoctial elements are the state variables. Uncertainties are modelled statistically by random vector and stochastic processes. The motion of the vehicle and the measurements are described by nonlinear stochastic differential and difference equations respectively. A minimum time trajectory is defined; equations of motion and measurements are linearized about this trajectory. An exponential cost criterion is constructed and a linear feedback quidance law is derived. An extended Kalman filter is used for state estimation. A short mission using this system is simulated. It is indicated that this system is efficient for short missions, but longer missions require accurate trajectory and ground based measurements.

Jack-Chingtse, C.↗
DOE OSTI2016

SU-G-TeP1-02: Analytical Stopping Power and Range Parameterization for Therapeutic Energy Intervals

Purpose: To develop a simple, analytic parameterization of stopping power and range, which covers a wide energy interval and is applicable to many species of projectile ions and target materials, with less than 15% disagreement in linear stopping power and 1 mm in range. Methods: The new parameterization was required to be analytically integrable from stopping power to range, and continuous across the range interval of 1 µm to 50 cm. The model parameters were determined from stopping power and range data for hydrogen, carbon, iron, and uranium ions incident on water, carbon, aluminum, lead and copper. Stopping power and range data was taken from SRIM. A stochastic minimization algorithm was used to find model parameters, with 10 data points per energy decade. Additionally, fitting was performed with 2 and 26 data points per energy decade to test the model’s robustness to sparse Results: 6 free parameters were sufficient to cover the therapeutic energy range for each projectile ion species (e.g. 1 keV – 300 MeV for protons). The model agrees with stopping power and range data well, with less than 9% relative stopping power difference and 0.5 mm difference in range. As few as, 4 bins per decade were required to achieve comparable fitting results to the full data set. Conclusion: This study successfully demonstrated that a simple analytic function can be used to fit the entire energy interval of therapeutic ion beams of hydrogen and heavier elements. Advantages of this model were the small number (6) of free parameters, and that the model calculates both stopping power and range. Applications of this model include GPU-based dose calculation algorithms and Monte Carlo simulations, where available memory is limited. This work was supported in part by a research agreement between United States Naval Academy and Louisiana State University: Contract No N00189-13-P-0786. In addition, this work was accepted for presentation at the American Nuclear Society Annual Meeting 2016.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRSDate not supplied

Stochastic Verification by Analysis for Autonomous Systems Management Architecture (ASMA)

The Gateway Vehicle Systems Manager (VSM) is the top-level of a distributed, hierarchical software control system. VSM is data-driven and will make decisions related to mission, fault, resource management and vehicle control. These attributes combined with a high degree of autonomy make it susceptible to emergent behavior. In order to achieve the high level of confidence needed in this critical system, the VSM team has developed a multifaceted verification strategy employing traditional verification techniques, simulation, model checking, and runtime verification. Individual algorithms are verified using conventional testing and model checking using assume-guarantee contracts. A discrete event-based simulation approach is being developed to verify timelines. This presentation describes an enhancement to the verification approach using analysis to enhance system robustness by detecting and resolving the potential for emergent behavior. The verification by analysis employs a Software in the Loop (SITL) environment with real flight software executing on emulated processors, simulations of vehicle subsystems, flight dynamics, and human inputs. Since the possible input space and configuration data set are too large for exhaustive testing, a Monte Carlo approach is used to cover feasible scenarios, augmented with corner cases and known higher-risk scenarios. A key problem in using Monte Carlo-based system verification is evaluating test results to ensure that system behavior is correct. The presentation describes the approach the VSM team uses to monitor behavior for compliance with predetermined boundaries and to identify anomalous behavior for further analysis. This presentation describes the multi-level systems approach to verification, and the simulation-based layer that covers the feasible state space: 1. Overview of the Gateway VSM 2. Special challenges due to heterogeneous, hierarchical architecture 3. Modeling and simulation environment using flight software and system simulations 4. Developing input sets to ensure state-space coverage 5. Developing model and data configuration sets to ensure model coverage 6. Interpreting results without predetermined outcomes 7. Lessons learned and future work

Verification and Validation↗
DOE OSTI2026

Cascading economic losses from port disruptions under capacity constrained multimodal freight networks

This study quantifies how throughput disruptions at major seaports cascade through capacity-constrained multimodal freight networks and interregional production systems. We couple an agent-based model (ABM) multimodal freight simulation that resolves rerouting, terminal queueing, and inventory drawdown under binding modal and facility capacities with a multiregional output loss input-output (MRIIM) model that propagates realized delivery shortfalls across regions and sectors. The framework is demonstrated for the Port of Los Angeles using Freight Analysis Framework flows and Bureau of Economic Analysis input-output accounts and is evaluated over a 52-week horizon under deterministic sector targeted shocks and stochastic disruption realizations with uncertain severity and duration. Results indicate nonlinear amplification: realized national losses concentrate in manufacturing and transportation/warehousing even when exogenous port shocks are dispersed, suggesting that congestion spillback and limited short-run substitution can dominate the initial shock allocation. We further evaluate a tabular reinforcement-learning (Q-learning) intervention layer that selects among a small set of implementable system level levers (truck-to-rail and truck-to-barge shift settings) without overriding shipper routing, finding that such interventions reduce total losses for moderate disruptions but yield diminishing returns once substitute modes approach capacity. By linking operational freight behavior to system wide impacts under uncertainty, the proposed ABM-MRIIM pipeline provides a reusable workflow for port disruption stress testing, identification of structurally critical sectors/corridors, and evaluation of resilience interventions under realistic capacity limits.

42 ENGINEERING↗
DOE OSTI2024

Modelling of the electron cyclotron emission burst from a laboratory tokamak plasma with loss-cone maser instability

The maser instability associated with the loss-cone distribution has been widely invoked to explain the radio bursts observed in the astrophysical plasma environment, such as aurora and corona. In the laboratory plasma of a tokamak, events reminiscent of these radio bursts have also been frequently observed as an electron cyclotron emission (ECE) burst in the microwave range (~2f ce near the last closed flux surface) during transient magnetohydrodynamic events. These bursts have a short duration of ~10 μs and display a radiation spectrum corresponding to a radiation temperature T e,rad of over 30 keV while the edge thermal electron temperature T e is only in the range of 1 keV. Suprathermal electrons can be generated through magnetic reconnection, and a loss-cone distribution can be generated through open stochastic field lines in the magnetic mirror of the near-edge region of a tokamak plasma. Radiation modelling shows that a sharp distribution gradient ∂f/∂v ⊥ > 0 at the loss-cone boundary can cause a negative absorption of ECE radiation through the maser instability. The negative absorption then amplifies the radiation so that the microwave intensity is significantly stronger than the thermal value. The significant T e,rad from the simulations suggests the potential role of the loss-cone maser instability in generating the ECE burst in a tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Non-linear dynamic characteristics and optimal control of giant magnetostrictive film subjected to in-plane stochastic excitation

The non-linear dynamic characteristics and optimal control of a giant magnetostrictive film (GMF) subjected to in-plane stochastic excitation were studied. Non-linear differential items were introduced to interpret the hysteretic phenomena of the GMF, and the non-linear dynamic model of the GMF subjected to in-plane stochastic excitation was developed. The stochastic stability was analysed, and the probability density function was obtained. The condition of stochastic Hopf bifurcation and noise-induced chaotic response were determined, and the fractal boundary of the system's safe basin was provided. The reliability function was solved from the backward Kolmogorov equation, and an optimal control strategy was proposed in the stochastic dynamic programming method. Numerical simulation shows that the system stability varies with the parameters, and stochastic Hopf bifurcation and chaos appear in the process; the area of the safe basin decreases when the noise intensifies, and the boundary of the safe basin becomes fractal; the system reliability improved through stochastic optimal control. Finally, the theoretical and numerical results were proved by experiments. The results are helpful in the engineering applications of GMF.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ASPIRES: Airport Shuttle Planning and Improved Routing Event-driven Simulation

Most of the existing traffic simulation packages require significant calibration work to be able to reflect reality. To evaluate special operations including emerging technologies, a microscopic simulation that tracks detailed interactions of all the elements of the traffic systems is usually needed. This type of simulation is usually computationally demanding. This work developed an Airport Shuttle Planning and Improved Routing Event-driven Simulation (ASPIRES) package to simulate and evaluate current, potential, and future airport shuttle operations. Here, the simulation was driven by data and thus did not require much calibration effort. The discrete-event simulation nature of ASPIRES makes the simulation computationally efficient. Simulating 1 day of shuttle operations takes less than 2s. The study site of this work is the Dallas/Fort Worth International Airport in the U.S. The shuttle service that connects the five terminals of the airport and the rental car center was studied. Travel times, dwell times, and passenger arrivals were simulated using empirical distributions derived mainly from real data to capture the stochastic nature of the rental car center shuttle bus operations. Data on bus miles traveled, bus energy consumption, passenger wait times, and passengers left behind at stops were collected to study the trade-off between energy use and passenger experience. Electric bus and on-demand bus operations were also included. The simulation outputs can show passengers statistics at terminals, shuttles statistics, and charging station statistics. ASPIRES cannot be used to model a generic traffic system but is well-suited for fleet systems.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2013

Documentation of a Linked-Population Viability Model to Evaluate Recovery Options for Snake River Fall Chinook Salmon

We developed an age- and life-history based population viability analysis model for linked populations of fall Chinook salmon (LPVA). LPVA was designed to quantify the costs and benefits of re-establishing access to an upstream spawning area using managed reconnection options such as hatchery stocking and translocation. The model incorporates stochasticity in the age at which salmon return to spawn, sex ratio, juvenile life history, and survival rates. Response variables include the proportion of replicate populations that persisted for 80 y (roughly 20 generations), the linked population’s growth rate, λ, time to extinction and summaries of population sizes including average female spawner abundance, final female spawner abundance, and total female abundance. This report attempts to provide a complete model description, including simulation of management options and options for simulating forward with a joint distribution of parameters based on their likelihoods with respect to historical data. Intermediate results illustrating the Bayesian approach to identifying poorly-known parameters are reported here. LPVA simulates reconnection of linked populations by representing management processes such as trapping of adults migrating upstream at two dams and juveniles migrating downstream, and hatchery operations for one or two hatcheries. The purpose of the model is to support efforts to recover the threatened Snake River fall Chinook salmon ESU by evaluating population viability using existing spawning areas with and without a hatchery and to assess the potential costs and benefits of programs to re-establish an upstream spawning area.

Jager, Henriette (Yetta) I.↗
NASA NTRS1977

Set-point control of nonlinear systems

A method of nonlinear stochastic control is applied to the problem of providing set-point control. The procedure is illustrated for a nonlinear second-order plant. The controller design uses a variation of the 'multiple-model adaptive control' approach. Specific problems related to the design are investigated, including learning time and set-point accuracy dependence on the number of approximating models, chord versus tangent approximations, and the a priori Markov switching parameter values. Simulations are provided for different situations under various levels of measurement noise.

Vanlandingham, H. F.↗
NASA NTRS1989

Probalistic structural analysis methods for select space propulsion system structural components

A summary of the status of this five-year project which is now in its third year of research and development is presented. The goal of the project is the development of several methodologies for probabilistic structural modeling. Probabilistic structural modeling consists of stochastic models of material properties, part geometries, boundary conditions, as well as loading conditions. The current presentation focuses on one methodology - coupling of an advanced finite element structural analysis code with probabilistic modeling strategies. The essential algorithm developments for combining the finite element and probabilistic analysis methods are reported. The validity of the resulting probabilistic structural analysis method is confirmed through a series of test problems with exact results based on Monte Carlo simulations. Additionally, the applicability of the method to a Space Propulsion System (a turbine blade) is demonstrated for static stresses.

Cruse, T. A.↗
NASA NTRS1990

Advances and trends in computational structures technology

The major goals of computational structures technology (CST) are outlined, and recent advances in CST are examined. These include computational material modeling, stochastic-based modeling, computational methods for articulated structural dynamics, strategies and numerical algorithms for new computing systems, multidisciplinary analysis and optimization. The role of CST in the future development of structures technology and the multidisciplinary design of future flight vehicles is addressed, and the future directions of CST research in the prediction of failures of structural components, the solution of large-scale structural problems, and quality assessment and control of numerical simulations are discussed.

Noor, A. K.↗
NASA NTRS1993

Fuel injector: Air swirl characterization aerothermal modeling, phase 2, volume 1

A well integrated experimental/analytical investigation was conducted to provide benchmark quality relevant to a prefilming type airblast fuel nozzle and its interaction with the combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM), and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS1993

Fuel Injector: Air swirl characterization aerothermal modeling, phase 2, volume 2

A well integrated experimental/analytical investigation was conducted to provide benchmark quality data relevant to prefilming type airblast fuel nozzle and its interaction with combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) equipment was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM) and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS2010

Absorbed Dose and Dose Equivalent Calculations for Modeling Effective Dose

While in orbit, Astronauts are exposed to a much higher dose of ionizing radiation than when on the ground. It is important to model how shielding designs on spacecraft reduce radiation effective dose pre-flight, and determine whether or not a danger to humans is presented. However, in order to calculate effective dose, dose equivalent calculations are needed. Dose equivalent takes into account an absorbed dose of radiation and the biological effectiveness of ionizing radiation. This is important in preventing long-term, stochastic radiation effects in humans spending time in space. Monte carlo simulations run with the particle transport code FLUKA, give absorbed and equivalent dose data for relevant shielding. The shielding geometry used in the dose calculations is a layered slab design, consisting of aluminum, polyethylene, and water. Water is used to simulate the soft tissues that compose the human body. The results obtained will provide information on how the shielding performs with many thicknesses of each material in the slab. This allows them to be directly applicable to modern spacecraft shielding geometries.

Welton, Andrew↗
NASA NTRSDate not supplied

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗
DOE OSTI2022

Comparison of Lagrangian Superdroplet and Eulerian Double-Moment Spectral Microphysics Schemes in Large-Eddy Simulations of an Isolated Cumulus Congestus Cloud

Abstract Advanced microphysics schemes (such as Eulerian bin and Lagrangian superdroplet) are becoming standard tools for cloud physics research and parameterization development. This study compares a double-moment bin scheme and a Lagrangian superdroplet scheme via large-eddy simulations of nonprecipitating and precipitating cumulus congestus clouds. Cloud water mixing ratio in the bin simulations is reduced compared to the Lagrangian simulations in the upper part of the cloud, likely from numerical diffusion, which is absent in the Lagrangian approach. Greater diffusion in the bin simulations is compensated by more secondary droplet activation (activation above cloud base), leading to similar or somewhat higher droplet number concentrations and smaller mean droplet radius than the Lagrangian simulations for the nonprecipitating case. The bin scheme also produces a significantly larger standard deviation of droplet radius than the superdroplet method, likely due to diffusion associated with the vertical advection of bin variables. However, the spectral width in the bin simulations is insensitive to the grid spacing between 50 and 100 m, suggesting other mechanisms may be compensating for diffusion as the grid spacing is modified. For the precipitating case, larger spectral width in the bin simulations initiates rain earlier and enhances rain development in a positive feedback loop. However, with time, rain formation in the superdroplet simulations catches up to the bin simulations. Offline calculations using the same drop size distributions in both schemes show that the different numerical methods for treating collision–coalescence also contribute to differences in rain formation. The stochastic collision–coalescence in the superdroplet method introduces more variability in drop growth for a given rain mixing ratio.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

Multi-State Physics Models of Aging Passive Components in Probabilistic Risk Assessment

Multi-state Markov modeling has proved to be a promising approach to estimating the reliability of passive components - particularly metallic pipe components - in the context of probabilistic risk assessment (PRA). These models consider the progressive degradation of a component through a series of observable discrete states, such as detectable flaw, leak and rupture. Service data then generally provides the basis for estimating the state transition rates. Research in materials science is producing a growing understanding of the physical phenomena that govern the aging degradation of passive pipe components. As a result, there is an emerging opportunity to incorporate these insights into PRA. This paper describes research conducted under the Risk-Informed Safety Margin Characterization Pathway of the Department of Energy’s Light Water Reactor Sustainability Program. A state transition model is described that addresses aging behavior associated with stress corrosion cracking in ASME Class 1 dissimilar metal welds – a component type relevant to LOCA analysis. The state transition rate estimates are based on physics models of weld degradation rather than service data. The resultant model is found to be non-Markov in that the transition rates are time-inhomogeneous and stochastic. Numerical solutions to the model provide insight into the effect of aging on component reliability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗