Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A study on the scale dependence of mixing indices for Eulerian multiphase models

Abstract Mixing can vary based on the scale at which the system is observed, and a mixing index that can capture the features at different length scales is desirable. In this article, we analyze the scale dependence of the mixing indices developed for Eulerian multiphase models. Relevant length scales are distinguished by filtering solid fraction fields. The scale‐dependence study is first done on manufactured fields of solid fraction to assess the performance of the mixing indices. The study is extended to a two‐dimensional CFD simulation of the segregation of a bidisperse gas–solid mixture. The local mixing index performs well in capturing the spatial variation of mixing at different scales. The scale dependence of two global mixing indices is considered in the study, where the state of mixing is defined based on statistical measures. We demonstrate that the choice of measures influences the sensitivity of mixing indices to mixing at different scales.

Nagawkar, Barlev R.↗
NASA NTRS2014

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GCR measurements is used by deep space mission planners for the certification of microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n - 20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determined the GCR flux at a given time applying an emperical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized more than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried outward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so emperical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.↗
NASA NTRS2014

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GR measurements is used by deep space mission planners for the certification of micro-electronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n-20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determines the GCR flux at a given time applying an empirical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized for than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried onward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so empirical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on the gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.↗
NASA NTRS1987

Effectiveness of back-to-back testing

Three models of back-to-back testing processes are described. Two models treat the case where there is no intercomponent failure dependence. The third model describes the more realistic case where there is correlation among the failure probabilities of the functionally equivalent components. The theory indicates that back-to-back testing can, under the right conditions, provide a considerable gain in software reliability. The models are used to analyze the data obtained in a fault-tolerant software experiment. It is shown that the expected gain is indeed achieved, and exceeded, provided the intercomponent failure dependence is sufficiently small. However, even with the relatively high correlation the use of several functionally equivalent components coupled with back-to-back testing may provide a considerable reliability gain. Implications of this finding are that the multiversion software development is a feasible and cost effective approach to providing highly reliable software components intended for fault-tolerant software systems, on condition that special attention is directed at early detection and elimination of correlated faults.

Vouk, Mladen A.↗
DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling delayed thermal runaway in nitric acid-soaked cat litter mixed with radioactive waste

Thermal ignition of radioactive waste within a 55-gallon drum was simulated by using a pressure-dependent waste decomposition model (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) calibrated with data from full-scale drum experiments (Parker et al. in The thermolytic response of a surrogate RNS waste mixture at the drum scale. Los Alamos National Laboratory Report LA-UR-16-21760, 2016) and validated with experiments from multiple laboratories (Hobbs et al. in Thermal analysis of aged nitric acid-soaked kitty litter in TRU waste drums-23370.WM2023 Conference, Phoenix, AZ, 2023). The acceleration of nitric acid chemistry reacting with an organic cat litter leading to thermal ignition was likely triggered by a restricted vent in the drum. Here, we address whether the form of the rate equation in (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) is sufficient to extrapolate thermal ignition within aged drums of similar content that have been stored in Texas for over nine years by investigating four different reaction rate forms for waste decomposition. A critical reaction rate reduction analysis is performed on each of these models to determine if delayed thermal runaway within vented aged waste is possible after nine years. We found that a pressure-dependent first-order rate expression not only predicted the accidental ignition of the waste drum, but the form also matches multiple experiments from different laboratories. Even though the waste composition decreases over time, the model predicts that acceleration leading to thermal runaway is possible if the waste is confined, even after 9 years. In conclusion, waste containing oxidizers such as nitric acid should not be mixed with organic adsorbents, especially if the waste is confined.

cookoff↗
NASA NTRS1971

The Effects of Parachute System Mass and Suspension-Line Elastic Properties on the LADT #3 Viking Parachute Inflation Load

Analytical calculations have considered the effects of 1) varying parachute system mass, 2) suspension-line damping, and 3) alternate suspension-line force-elongation data on the canopy force history. Results indicate the canopy force on the LADT #3 parachute did not substantially exceed the recorded vehicle force reading and that the above factors can have significant effects on the canopy force history. Analytical calculations have considered the effects of i) varying parachute system mass, 2) suspension line damping, and 3) different suspension-line force-elongation data on the canopy force history. Based on the results of this study the following conclusions are drawn: Specifically, 1. At the LADT #3 failure time of 1.70 seconds, the canopy force ranged anywhere from 15.7% below to 2.4% above the vehicle force depending upon the model and data used. Therefore, the canopy force did not substantially exceed the recorded vehicle force reading. 2. At a predicted full inflation time of 1.80 seconds the canopy force would be greater than the vehicle force by from 1.1% to 10.6%, again depending upon the model and data used. Generally, 3. At low altitudes, enclosed and apparent air mass can significantly effect the canopy force calculated and should, therefore, not be neglected. 4. The canopy force calculations are sensitive to decelerator physical properties. In this case changes in the damping and/or force-elongation characteristics produced significant changes in the canopy force histories. Accurate prediction of canopy force histories requires accurate inputs in these areas.

Talay, Theodore A.↗
NASA NTRS1981

The possible effect of 11 year solar cycle on the thermal structure of the stratosphere and ozone distribution investigated in a 2-D model

A two-dimensional time dependent radiative-photochemical-dynamical model of the stratosphere which takes into account the feedback effect of temperature on ozone concentration is used to investigate the possible changes in the seasonal and latitudinal variations of temperature, ozone, and other minor constituents associated with the UV flux variations from the solar minimum to the solar maximum. The transport model determines the zonal wind, non-zonal eddy transports, and eddy momentum fluxes, while the photochemical system includes the important reactions affecting the concentrations of ozone and other relevant trace constituents in the stratosphere. Experiments are carried out by changing the position of the sun every day and evaluating the daily average photodissociation integrals for the relevant minor constituents and ozone heating. Values are found for the ambient stratosphere and a stratosphere perturbed by solar cycle variations, and a rough comparison is made with observed ozone and temperature trends.

Vupputuri, R. K. R.↗
DOE OSTI2026

A Multi-Region SEIR Model Incorporating Inter-County Mobility and Time-Dependent Transmission Dynamics: Application to COVID-19 Disease Outbreak Data in North Carolina.

Classical infectious disease compartmental models typically do not incorporate spatial heterogeneity or mobility. We develop a multi-region susceptible-exposed-infected-recovered (SEIR) model in which disease dynamics are coupled via inter-region mobility and the transmission rate is both region and time dependent. We calibrate the model using rolling averages of daily COVID-19 data in all 100 North Carolina counties. Mobility parameters are prescribed using daily inter-county commuter data. The number of transmission rate parameters is substantially reduced by hypothesizing that the dynamics correlate with county-level population density. Parameter estimation is carried out using several objective functions with error terms at different scales. An additive combination of least squares error at the county-level and the state-level, along with a quadratic transmission rate polynomial, yields the lowest overall error at both spatial scales. The calibrated model is used to simulate regional effects of perturbing disease transmission rates in adjacent counties and to illustrate effects of the state’s mobility infrastructure on disease dynamics and spread for a new disease outbreak.

COVID-19 modeling↗
NASA NTRS1988

Data dependent systems methodology for lumped mass modeling of structures

Limitations of the frequency domain methods in analyzing structura1 vibrations has created an awareness of the comparative merits of the time domain methods. Although time domain methods would be ideal for modeling large precisions space systems, the popular methods based on fitting theoretical response to actual data by least squares are too sensitive to noise and require too much data to be suitable for orbiting space crafts. This paper briefly reviews the theory and illustrative applications of a time domain methodology called Data Dependent Systems (DDS) that eliminates these limitations. Simulation results are presented to demonstrate a better than 4-place accuracy in the identifications of all system parameters, both modal (frequencies, damping ratios, and mode shapes) and physical (mass, stiffness, and damping matrices).

Pandit, Sudhakar M.↗
NASA NTRS2002

Hydrodynamic Models of Line-Driven Accretion Disk Winds III: Local Ionization Equilibrium

We present time-dependent numerical hydrodynamic models of line-driven accretion disk winds in cataclysmic variable systems and calculate wind mass-loss rates and terminal velocities. The models are 2.5-dimensional, include an energy balance condition with radiative heating and cooling processes, and includes local ionization equilibrium introducing time dependence and spatial dependence on the line radiation force parameters. The radiation field is assumed to originate in an optically thick accretion disk. Wind ion populations are calculated under the assumption that local ionization equilibrium is determined by photoionization and radiative recombination, similar to a photoionized nebula. We find a steady wind flowing from the accretion disk. Radiative heating tends to maintain the temperature in the higher density wind regions near the disk surface, rather than cooling adiabatically. For a disk luminosity L (sub disk) = solar luminosity, white dwarf mass M(sub wd) = 0.6 solar mass, and white dwarf radii R(sub wd) = 0.01 solar radius, we obtain a wind mass-loss rate of M(sub wind) = 4 x 10(exp -12) solar mass yr(exp -1) and a terminal velocity of approximately 3000 km per second. These results confirm the general velocity and density structures found in our earlier constant ionization equilibrium adiabatic CV wind models. Further we establish here 2.5D numerical models that can be extended to QSO/AGN winds where the local ionization equilibrium will play a crucial role in the overall dynamics.

Pereyra, Nicolas Antonio↗
DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗
DOE OSTI2025

First measurement of symmetric cumulants of hexagonal flow harmonics in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Correlations between event-by-event fluctuations of anisotropic flow harmonics are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV, as recorded by the ALICE detector at the LHC. This study presents correlations up to the hexagonal flow harmonic 𝑣 6 , which was measured for the first time. The magnitudes of these higher-order correlations are found to vary as a function of collision centrality and harmonic order. These measurements are compared to viscous hydrodynamic model calculations with EKRT initial conditions and to the iEBE-VISHNU model with T R ⁢ENTo initial conditions. The observed discrepancies between the data and the model calculations vary depending on the harmonic combinations. Due to the sensitivity of model parameters estimated with Bayesian analyses to these higher-order observables, the results presented in this work provide new and independent constraints on the initial conditions and transport properties in theoretical models used to describe the system created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Towards Realistic Models of Heterogeneous Catalysis: Simulations of Oxidation Catalysis from First Principles

Among the most significant developments in heterogeneous catalysis in the last 20 years is the emergence of microkinetic models, often parameterized from density functional theory (DFT) calculations, used to quantify observed catalyst performance and to guide the discovery of new catalytic materials. While DFT directly reports binding energies and elementary step activation energies, the free energies that enter into microkinetic models must be computed from additional approximations. Frequently these free energy approximations assume ideal behavior—that adsorbates do not interact with one another, or that adsorbates are immobile or vibrate harmonically about a binding site. These assumptions can and do have an impact on predicted catalyst performance and potentially even on predicted trends. Our work relates to three categories of non-idealities, explored in the context of nitrogen oxidation and reduction catalysis on metal surfaces. One component relates to adsorbate translational free energy, the contribution least well described by conventional models, in which we develop modeling approaches that improve accuracy with limited increase in computational expense. A second component relates adsorbate-adsorbate interactions, in which we develop benchmark on-lattice interaction models and compare kinetic predictions with conventional mean-field models, in an effort to develop more robust coverage-dependent mean field modeling approaches. The last relates the most fundamental assumption of all—that of energy equipartition—in an effort to rationalize and guide plasma-enhanced catalysis.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1977

Evaluation of turbulence models for three primary types of shock-separated boundary layers

Zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence eddy viscosity models were used in computing three basic types of shock-separated boundary-layer flows. The three basic types of shock boundary-layer interaction discussed are: (1) a normal shock wave at transonic speeds, (2) a compression corner shock at supersonic speeds, and (3) an incident oblique shock at hypersonic speeds. The models tested are simple, unmodified models used extensively for incompressible, unseparated flows. A comparison of computed and measured results for the compressible, separated flows described herein indicates that model performance is dependent on flow configuration with no distinct superiority of one model over the other for all three flow configurations.

Coakley, T. J.↗
NASA NTRS1980

Evaluation and comparison of tropical analyses during DST-5 and DST-6

Data systems tests were conducted to assess the adequacy of the global data base for numerical analysis and forecasting and in the process to determine the impact of meteorological satellite data. The results of these tests indicate that the satellite data impact is model and season dependent but definitely dependent on the method of data assimilation and the numerical model used to produce the analyses or initial conditions.

Sadler, J. C.↗
NASA NTRS1989

An entrance region friction factor model applied to annular seal analysis - Theory versus experiment for smooth and honeycomb seals

A friction factor model is developed for the entrance-region of a duct. The model is used in an annular gas seal analysis similar to Nelson's (1984). Predictions of the analysis are compared to experimental results for a smooth-stator/smooth-rotor seal and three honeycomb-stator/smooth-rotor seals. The model predicts a leakage and direct damping well. The model overpredicts the dependence of cross-coupled stiffness on fluid prerotation. The model predicts direct stiffness poorly.

Elrod, D.↗
NASA NTRS1995

Differenced Range Versus Integrated Doppler (DRVID) ionospheric analysis of metric tracking in the Tracking and Data Relay Satellite System (TDRSS)

The Differenced Range (DR) Versus Integrated Doppler (ID) (DRVID) method exploits the opposition of high-frequency signal versus phase retardation by plasma media to obtain information about the plasma's corruption of simultaneous range and Doppler spacecraft tracking measurements. Thus, DR Plus ID (DRPID) is an observable independent of plasma refraction, while actual DRVID (DR minus ID) measures the time variation of the path electron content independently of spacecraft motion. The DRVID principle has been known since 1961. It has been used to observe interplanetary plasmas, is implemented in Deep Space Network tracking hardware, and has recently been applied to single-frequency Global Positioning System user navigation This paper discusses exploration at the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) of DRVID synthesized from simultaneous two-way range and Doppler tracking for low Earth-orbiting missions supported by the Tracking and Data Relay Satellite System (TDRSS) The paper presents comparisons of actual DR and ID residuals and relates those comparisons to predictions of the Bent model. The complications due to the pilot tone influence on relayed Doppler measurements are considered. Further use of DRVID to evaluate ionospheric models is discussed, as is use of DRPID in reducing dependence on ionospheric modeling in orbit determination.

Radomski, M. S.↗