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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗
NASA NTRS1996

Computer Modeling of Non-Isothermal Crystallization

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homogeneous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, size-dependent growth rates, and surface crystallization are also included. Model predictions are compared with experimental measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the peak in the steady-state rate is demonstrated.

Kelton, K. F.↗
DOE OSTI2021

Computational Approach for Rapidly Predicting Temperature–Dependent Polymer Solubilities Using Molecular–Scale Models

One promising approach to recycle multicomponent plastic waste (e. g., multilayer plastic films) is selective dissolution. Selective dissolution is a solvent-mediated process in which differences in polymer solubility in a carefully chosen solvent system are exploited to recover a target polymer. Here, a computational approach was developed that rapidly predicts temperature-dependent polymer solubilities to guide the design of solvent systems for solvent-mediated polymer recycling. Polymer conformations were obtained from molecular dynamics simulations by modeling the polymer as a short oligomer and then used as input to the conductor-like screening model for real solvents (COSMO-RS) for solubility predictions. Using polyethylene (PE) and ethylene vinyl alcohol (EVOH) as representative polymers, the effect of simulation parameters was systematically studied, and predicted solubilities were found to be in good agreement with experimental measurements. Furthermore, the applicability of the approach was demonstrated by identifying selective solvents for PE and EVOH dissolution from a library of 524 solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗
NASA NTRS1991

Biogeochemical modeling at mass extinction boundaries

The causes of major mass extinctions is a subject of considerable interest to those concerned with the history and evolution of life on earth. The primary objectives of the proposed plan of research are: (1) to develop quantitative time-dependent biogeochemical cycle models, coupled with an ocean atmosphere in order to improve the understanding of global scale physical, chemical, and biological processes that control the distribution of elements important for life at times of mass extinctions; and (2) to develop a comprehensive data base of the best available geochemical, isotopic, and other relevant geologic data from sections across mass extinction boundaries. These data will be used to constrain and test the biogeochemical model. These modeling experiments should prove useful in: (1) determining the possible cause(s) of the environmental changes seen at bio-event boundaries; (2) identifying and quantifying little-known feedbacks among the oceans, atmosphere, and biosphere; and (3) providing additional insights into the possible responses of the earth system to perturbations of various timescales. One of the best known mass extinction events marks the Cretaceous/Tertiary (K/T) boundary (66 Myr ago). Data from the K/T boundary are used here to constrain a newly developed time-dependent biogeochemical cycle model that is designed to study transient behavior of the earth system. Model results predict significant fluctuations in ocean alkalinity, atmospheric CO2, and global temperatures caused by extinction of calcareous plankton and reduction in the sedimentation rates of pelagic carbonates and organic carbon. Oxygen-isotome and other paleoclimatic data from K/T time provide some evidence that such climatic fluctuations may have occurred, but stabilizing feedbacks may have acted to reduce the ocean alkalinity and carbon dioxide fluctuations.

Rampino, M. R.↗
DOE OSTI2020

Helicity-dependent extension of the McLerran–Venugopalan model

We construct a generalization of the McLerran-Venugopalan (MV) model including helicity effects for a longitudinally polarized target (a proton or a large nucleus). The extended MV model can serve as the initial condition for the helicity generalization of the JIMWLK evolution equation constructed in our previous paper, as well as for calculation of helicity-dependent observables in the quasi-classical approximation to QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2025

Time-Dependent Satellite Explosion Probabilities for Long-Term Orbital Debris Environment Modeling

On-orbit accidental explosions are a significant contributor to the growth of the orbital debris population. Certain types of satellites have been seen to exhibit different behaviors and timelines of explosions. Recent assessments of on-orbit explosions for three categories of satellites – Russian Sistema Obespecheniya Zapuska (SOZ, Proton 4th stage attitude and ullage motors) units, spacecraft, and rocket bodies – suggest that on-orbit explosions can be modeled using a continuous, time-dependent probability of explosion. This paper discusses the methodology for developing a time-dependent probability of explosion as a function of the satellite’s orbital lifetime based on historical explosions. The SOZ units are found to be best modeled by a modified Gaussian probability distribution, with the peak probability of explosion occurring around 10.4 years on-orbit. Spacecraft and rocket bodies both exhibit exponentially decaying explosion probabilities, with rocket bodies showing fast-, medium-, and slow-time modes of decay. Total cumulative probabilities of explosion are approximately 57% for SOZ units, 4% for spacecraft, and 2% for rocket bodies. Examples of implementation in a simulation of the orbital debris environment using NASA's long-term evolutionary model are given, showing the overall explosion behavior over two centuries is significantly different from that using a previously implemented constant-probability type-dependent explosion rate model and provides a better representation of the historical explosion record.

Alyssa P Manis↗
DOE OSTI2025

Prevention of resistive wall tearing mode major disruptions with feedback

Resistive wall tearing modes (RWTMs) can cause major disruptions. A signature of RWTMs is that the rational surface is sufficiently close to the wall to interact with it. For (m,n) = (2,1) modes, a RWTM requires normalized minor radius of the rational surface ρ q2 ≥ 0.75, which can also be expressed as q 75 ≤ 2. Major disruptions can occur when the criterion is satisfied. This is confirmed in simulations and theory and in a DIII-D locked mode disruption database. The q 75 < 2 criterion is valid at high β as well as at low β. A very important feature of RWTMs is that they can be feedback stabilized. If the ρ q2 criterion is not satisfied, or if the wall is ideally conducting, then the mode does not produce a major disruption, although it can produce a minor disruption. Feedback, or rotation of the mode at the wall by complex feedback, can emulate an ideal wall, preventing major disruptions. The ρ q2 criterion depends weakly on the wall radius. A simple geometric model of its dependence on wall radius is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

The longitude dependence of the dayside F region trough - A detailed model-observation comparison

The nighttime main F-region trough extends into the sunlit afternoon sector. This trough feature exhibits both a strong magnetic-activity dependence and a longitude (UT) dependence. Whalen (1987), using IGY ionosonde data, showed that both of these effects are readly extracted from f0F2 observations. This study shows that the longitude effect is the same as that contained in the Utah State University time-dependent ionospheric model. It arises from the offset of the geomagnetic axis from the geographic axis. The magnetic-activity dependency is associated with the westward convection in the afternoon sector. It is also contained in the ionospheric model via the empirical magnetospheric convection model.

Sojka, J. J.↗
DOE OSTI2020

DFT and microkinetic comparison of Pt, Pd and Rh-catalyzed ammonia oxidation

Ammonia oxidation is the heart of the Ostwald process and is important in emissions control. Catalytic behaviors are a function of conditions and are observed to vary across the platinum group metals (PGMs) Pt, Pd, Rh. Here, we combine density functional theory computations and microkinetic modeling to rationalize these dependencies. We compute reactions over model (2 1 1) and (1 1 1) surfaces of PGMs. Binding energies are similar on Pd and Pt and generally greater on Rh, while activation energies vary across all metals. Rates on (2 1 1) surfaces are greater than (1 1 1) surfaces. Additionally, the stepped Pt is most active and stepped Rh most selective to N 2 at ammonia slip conditions, while at Ostwald process conditions, stepped Pd is most active and stepped Pt most selective to NO. Degree of rate and selectivity control analysis provides insights into the reactions limiting performance of PGMs. Both activation barriers and surface coverages influence rates and selectivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Chemistry↗
NASA NTRS1980

Clarification process: Resolution of decision-problem conditions

A model of a general process which occurs in both decisionmaking and problem-solving tasks is presented. It is called the clarification model and is highly dependent on information flow. The model addresses the possible constraints of individual indifferences and experience in achieving success in resolving decision-problem conditions. As indicated, the application of the clarification process model is only necessary for certain classes of the basic decision-problem condition. With less complex decision problem conditions, certain phases of the model may be omitted. The model may be applied across a wide range of decision problem conditions. The model consists of two major components: (1) the five-phase prescriptive sequence (based on previous approaches to both concepts) and (2) the information manipulation function (which draws upon current ideas in the areas of information processing, computer programming, memory, and thinking). The two components are linked together to provide a structure that assists in understanding the process of resolving problems and making decisions.

Dieterly, D. L.↗
NASA NTRS1993

Decay of equatorial ring current ions and associated aeronomical consequences

The decay of the major ion species which constitute the ring current is studied by solving the time evolution of their distribution functions during the recovery phase of a moderate geomagnetic storm. In this work, only equatorially mirroring particles are considered. Particles are assumed to move subject to E x B and gradient drifts. They also experience loses along their drift paths. Two loss mechanisms are considered: charge exchange with neutral hydrogen atoms and Coulomb collisions with thermal plasma in the plasmasphere. Thermal plasma densities are calculated with a plasmaspheric model employing a time-dependent convection electric field model. The drift-loss model successfully reproduces a number of important and observable features in the distribution function. Charge exchange is found to be the major loss mechanism for the ring current ions; however the important effects of Coulomb collisions on both the ring current and thermal populations are also presented. The model predicts the formation of a low-energy (less than 500 eV) ion population as a result of energy degradation caused by Coulomb collision of the ring current ions with the plasmaspheric electrons; this population may be one source of the low-energy ions observed during active and quiet periods in the inner magnetosphere. The energy transferred to plasmaspheric electrons through Coulomb collisions with ring current ions is believed to be the energy source for the electron temperature enhancement and the associated 6300 A (stable auroral red (SAR) arc) emission in the subauroral region. The calculated energy deposition rate is sufficient to produce a subauroral electron temperature enhancement and SAR arc emissions that are consistent with observations of these quantities during moderate magnetic activity levels.

Fok, M.-C.↗
NASA NTRS1996

Studies of Disks Around the Sun and Other Stars

We are conducting research designed to enhance our understanding of the evolution and detectability of comet clouds and disks. This area holds promise for also improving our understanding of outer solar system formation, the bombardment history of the planets, the transport of volatiles and organics from the outer solar system to the inner planets, and to the ultimate fate of comet clouds around the Sun and other stars. According to 'standard' theory, both the Kuiper Disk and the Oort Cloud are (at least in part) natural products of the planetary accumulation stage of solar system formation. One expects such assemblages to be a common attribute of other solar systems. Therefore, searches for comet disks and clouds orbiting other stars offer a new method for inferring the presence of planetary systems. This two-element program consists modeling collisions in the Kuiper Disk and the dust disks around other stars. The modeling effort focuses on moving from our simple, first-generation, Kuiper disk collision rate model, to a time-dependent, second-generation model that incorporates physical collisions, velocity evolution, dynamical erosion, and various dust transport mechanisms. This second generation model will be used to study the evolution of surface mass density and the object-size spectrum in the disk. The observational effort focuses on obtaining submm/mm-wave flux density measurements of 25-30 IR excess stars in order to better constrain the masses, spatial extents and structure of their dust ensembles.

Stern, S. Alan↗
NASA NTRS1997

Studies of Disks Around the Sun and Other Stars

This is a NASA Origins of Solar Systems research program, and this NASA Headquarters grant has now been transferred to a new grant at NASA GSFC (NAG5-4082). Thus the need for this 'Final Report' on a project that is not, in fact, complete. We are conducting research designed to enhance our understanding of the evolution and detectability of comet clouds and disks. This area holds promise for also improving our understanding of outer solar system formation, the bombardment history of the planets, the transport of volatiles and organics from the outer solar system to the inner planets, and to the ultimate fate of comet clouds around the Sun and other stars. According to "standard" theory, both the Kuiper Belt and the Oort Cloud are (at least in part) natural products of the planetary accumulation stage of solar system formation. One expects such assemblages to be a common attribute of other solar systems. Our program consists of modeling collisions in the Kuiper Belt and the dust disks around other stars. The modeling effort focuses on moving from our simple, first-generation, Kuiper Belt collision rate model, to a time-dependent, second-generation model that incorporates physical collisions, velocity evolution, dynamical erosion, and various dust transport mechanisms. This second generation model is to be used to study the evolution of surface mass density and the object-size spectrum in the disk.

Stern, S. Alan↗
NASA NTRS1995

Plasmasphere Modeling with Ring Current Heating

Coulomb collisions between ring current ions and the thermal plasma in the plasmasphere will heat the plasmaspheric electrons and ions. During a storm such heating would lead to significant changes in the temperature and density of the thermal plasma. This was modeled using a time- dependent, one-stream hydrodynamic model for plasmaspheric flows, in which the model flux tube is connected to the ionosphere. The model simultaneously solves the coupled continuity, momentum, and energy equations of a two-ion (H(+) and O(+) quasineutral, currentless plasma. Heating rates due to collisions with ring current ions were calculated along the field line using a kinetic ring current model. First, diurnally reproducible results were found assuming only photoelectron heating of the thermal electrons. Then results were found with heating of the H(+) ions by the ring current during the recovery phase of a magnetic storm.

Guiter, S. M.↗
NASA NTRS2006

A dynamic fault tree model of a propulsion system

We present a dynamic fault tree model of the benchmark propulsion system, and solve it using Galileo. Dynamic fault trees (DFT) extend traditional static fault trees with special gates to model spares and other sequence dependencies. Galileo solves DFT models using a judicious combination of automatically generated Markov and Binary Decision Diagram models. Galileo easily handles the complexities exhibited by the benchmark problem. In particular, Galileo is designed to model phased mission systems.

propulsion system↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗