Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Implementation of a vibrationally linked chemical reaction model for DSMC

A new procedure closely linking dissociation and exchange reactions in air to the vibrational levels of the diatomic molecules has been implemented in both one- and two-dimensional versions of Direct Simulation Monte Carlo (DSMC) programs. The previous modeling of chemical reactions with DSMC was based on the continuum reaction rates for the various possible reactions. The new method is more closely related to the actual physics of dissociation and is more appropriate to the particle nature of DSMC. Two cases are presented: the relaxation to equilibrium of undissociated air initially at 10,000 K, and the axisymmetric calculation of shuttle forebody heating during reentry at 92.35 km and 7500 m/s. Although reaction rates are not used in determining the dissociations or exchange reactions, the new method produces rates which agree astonishingly well with the published rates derived from experiment. The results for gas properties and surface properties also agree well with the results produced by earlier DSMC models, equilibrium air calculations, and experiment.

Carlson, A. B.↗
NASA NTRS1992

Evaluation of thermochemical models for particle and continuum simulations of hypersonic flow

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered that are dominated in turn by vibrational relaxation, weak dissociation, strong dissociation and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Boyd, Iain D.↗
NASA NTRSDate not supplied

Updates to the Predictive Materials Modeling Software Tools

Updates on NASA's efforts to build predictive material models from the micro-scale to the macro-scale are presented. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. The PuMA (microstructure analysis), PATO (macroscale material response), SPARTA (direct simulation Monte Carlo) and ARCHeS (arc heater modeling) codes are described. A range of applications, encompassing calculating effective material properties, study of high shear boundary layer flow over woven materials, new model of silicone-based coatings, new model of mechanical erosion, as well as approach to loose multi-physics coupling, are presented.

material modeling↗
NASA NTRS1980

A radiative transfer model for remote sensing of laser induced fluorescence of phytoplankton in non-homogeneous turbid water

A radiative transfer computer model was developed to characterize the total flux of chlorophyll a fluoresced or backscattered photons when laser radiation is incident on turbid water that contains a non-homogeneous suspension of inorganic sediments and phytoplankton. The radiative transfer model is based on the Monte Carlo technique and assumes that: (1) the aquatic medium can be represented by a stratified concentration profile; and (2) that appropriate optical parameters can be defined for each layer. The model was designed to minimize the required computer resources and run time. Results are presented for an anacystis marinus culture.

Venable, D. D.↗
NASA NTRS1993

Development and application of computational aerothermodynamics flowfield computer codes

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second-approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered. The cases are governed, respectively, by the following: vibrational relaxation; weak dissociation; strong dissociation; and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Venkatapathy, Ethiraj↗
NASA NTRS1995

Magnetic Compton-induced pair cascade model for gamma-ray pulsars

Electrons accelerated to relativistic energies in pulsar magnetospheres will Compton scatter surface thermal emission and nonthermal optical, UV, and soft X-ray emission to gamma-ray energies, thereby initiating a pair cascade through synchrotron radiation and magnetic pair production. This process is proposed as the origin of the high-energy radiation that has been detected from six isolated pulsars. We construct an analytic model of magnetic Compton scattering near the polar cap of isolated pulsar magnetospheres and present approximate analytic derivations for scattered spectra, electron energy-loss rates, and photon luminosities. A Monte Carlo simulation is used to model the pair cascade induced by relativistic electrons scattering photons through the cyclotron resonance. For simplicity, the primary electrons are assumed to be monoenergetic and the nonresonant emission is omitted. Assuming that the angle phi(sub B) between the magnetic and spin axes is approximately equal to the polar-cap angle theta(sub pc), this model can produce both double-peaked and broad single-peaked pulse profiles and account for the trend of harder gamma-ray spectra observed from older pulsars.

Sturner, Steven J.↗
NASA NTRS2013

Nuclear Thermal Rocket Simulation in NPSS

Four nuclear thermal rocket (NTR) models have been created in the Numerical Propulsion System Simulation (NPSS) framework. The models are divided into two categories. One set is based upon the ZrC-graphite composite fuel element and tie tube-style reactor developed during the Nuclear Engine for Rocket Vehicle Application (NERVA) project in the late 1960s and early 1970s. The other reactor set is based upon a W-UO2 ceramic- metallic (CERMET) fuel element. Within each category, a small and a large thrust engine are modeled. The small engine models utilize RL-10 turbomachinery performance maps and have a thrust of approximately 33.4 kN (7,500 lbf ). The large engine models utilize scaled RL-60 turbomachinery performance maps and have a thrust of approximately 111.2 kN (25,000 lbf ). Power deposition profiles for each reactor were obtained from a detailed Monte Carlo N-Particle (MCNP5) model of the reactor cores. Performance factors such as thermodynamic state points, thrust, specific impulse, reactor power level, and maximum fuel temperature are analyzed for each engine design.

Belair, Michael L.↗
NASA NTRS2013

Nuclear Thermal Rocket Simulation in NPSS

Four nuclear thermal rocket (NTR) models have been created in the Numerical Propulsion System Simulation (NPSS) framework. The models are divided into two categories. One set is based upon the ZrC-graphite composite fuel element and tie tube-style reactor developed during the Nuclear Engine for Rocket Vehicle Application (NERVA) project in the late 1960s and early 1970s. The other reactor set is based upon a W-UO2 ceramic-metallic (CERMET) fuel element. Within each category, a small and a large thrust engine are modeled. The small engine models utilize RL-10 turbomachinery performance maps and have a thrust of approximately 33.4 kN (7,500 lbf ). The large engine models utilize scaled RL-60 turbomachinery performance maps and have a thrust of approximately 111.2 kN (25,000 lbf ). Power deposition profiles for each reactor were obtained from a detailed Monte Carlo N-Particle (MCNP5) model of the reactor cores. Performance factors such as thermodynamic state points, thrust, specific impulse, reactor power level, and maximum fuel temperature are analyzed for each engine design.

Belair, Michael L.↗
NASA NTRS2007

Simulating Flights of Future Launch Vehicles and Spacecraft

Marshall Aerospace Vehicle Representation in C (MAVERIC) is a computer program for generic, low-to-high-fidelity simulation of the flight(s) of one or more launch vehicle(s) or spacecraft. MAVERIC is designed to accommodate multi-staged vehicles, powered serially or in parallel, with multiple engines, tanks, and cargo elements. Engines can be of jet or conventional rocket types, using either liquid or solid propellants. MAVERIC includes generic subsystem software models for propulsion systems, mass properties, reaction control systems, aerodynamic properties, guidance systems, and navigation systems. Simulations can be started at points other than liftoff. Also included are guidance-system software models that accommodate the ascent, orbit, coasting, deorbiting, entry, terminal-area-energy-management, approach, and landing phases of flight. Options to use different wind profiles and atmospheres are included. A Monte Carlo capability is provided for modeling dispersions associated with atmospheric effects (including winds), propulsion, navigation, aerodynamics, and mass properties. Failures of engines and other subsystems can be modeled. The program is written in the C programming language, which makes it possible for the program to have high degrees of modularity, reusability, and maintainability, thereby also facilitating modification for modeling new vehicles.

McCarter, James W.↗
NASA NTRS2019

Bayesian Estimation of Earth’s Undiscovered Mineralogical Diversity Using Noninformative Priors

Recently, statistical distributions have been explored to provide estimates of the mineralogical diversity of Earth, and Earth-like planets. In this paper, a Bayesian approach is introduced to estimate Earth’s undiscovered mineralogical diversity. Samples are generated from a posterior distribution of the model parameters using Markov chain Monte Carlo simulations such that estimates and inference are directly obtained. It was previously shown that the mineral species frequency distribution conforms to a generalized inverse Gauss–Poisson (GIGP) large number of rare events model. Even though the model fit was good, the population size estimate obtained by using this model was found to be unreasonably low by mineralogists. In this paper, several zero-truncated, mixed Poisson distributions are fitted and compared, where the Poisson-lognormal distribution is found to provide the best fit. Subsequently, the population size estimates obtained by Bayesian methods are compared to the empirical Bayes estimates. Species accumulation curves are constructed and employed to estimate the population size as a function of sampling size. Finally, the relative abundances, and hence the occurrence probabilities of species in a random sample, are calculated numerically for all mineral species in Earth’s crust using the Poisson-lognormal distribution. These calculations are connected and compared to the calculations obtained in a previous paper using the GIGP model for which mineralogical criteria of an Earth-like planet were given.

Bayesian statistics↗
NASA NTRS2025

Assessment of Model Outcomes Between the Integrated Medical Model (IMM) and the Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT)

The Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT) is a computational model that provides human health and medical risk predictions for crewed spaceflight missions. MEDPRAT utilizes discrete event modeling and dynamic probabilistic simulation to predict critical mission outcomes (total medical events, crew health index, quality time lost, loss of crew life, removal to definitive care), condition occurrences, and resource consumption. Input parameters for MEDPRAT include crew attributes (e.g., sex), types of mission activities (e.g., whether and where crew members perform an extravehicular activity (EVA)), available resources, treatment information, and probability distributions for medical conditions. As an evolution of the Integrated Medical Model (IMM), MEDPRAT provides enhanced capabilities and higher fidelity, and incorporates more appropriate assumptions for long-duration spaceflight. IMM is the currently accepted standard for quantifying spaceflight mission medical risk in NASA operations that uses a probabilistic risk assessment (PRA) approach. MEDPRAT builds on the same logical foundation as IMM but implements the model architecture through highly optimized Monte Carlo sampling methods. An analysis is performed comparing the outputs from IMM with those from MEDPRAT V1.0 and V2.0 for the same reference missions in order to quantify similarities and differences in the model outcomes. The juxtaposition between IMM and MEDPRAT V1.0 and 2.0 shown in this report demonstrates that these two models generate very similar results; where differences in outcomes are shown, these are in accordance with the underlying assumptions and differences in the model architectures. This validation effort further establishes the credibility and reliability of the MEDPRAT software.

Matthew T Prelich↗
NASA NTRSDate not supplied

A Combined Computational, Experimental, and Technology Development Approach to In-Space Laser Manufacturing Maturation at NASA Marshall Space Flight Center

In-space manufacturing (ISM) is emerging as a field vital to continued access and capabilities in the space environment. NASA Marshall Space Flight Center (MSFC) is advancing the frontier of in-space laser manufacturing (ISLM) techniques through work initially focused on maturing laser beam welding (LBW) and laser forming (LF) for use in space. Such techniques proffer the ability to assemble and join structures in space from sheet metal or other stock – extant satellites, in situ resource utilization of Lunar regolith, etc. – by forming to desired shapes and then joining via in-space welding (ISW). ISLM processes are useful for assembly, joining, modification, and repair of structures in free space and on the Lunar surface such as large observatories, antennas, trusses, blast/thermal/radiation shields, pressure vessels, and more. However, these techniques are not yet qualified & certified (Q&C) for regular application in space. It would be prohibitively expensive, laborious, and time-consuming to perform Q&C via traditional experimental approaches as data collection & experimentation in space is resource-intensive. As such, benchmark experiments and focused, properly instrumented technology demonstration efforts in space can collect sufficient data that – when combined with verified computational models in an integrated computational materials engineering (ICME) approach – can validate ICME tools capable of translating more readily obtained ground data to in-space, in situ, computationally informed Q&C of ISLM techniques. Several ISLM projects at MSFC are obtaining the data required to validate ICME tools through both ground and flight experiments. A parabolic flight experiment of LBW under vacuum is manifested for August 2024, including both microgravity and Lunar gravity profiles. This collaboration with the Ohio State University is investigating common aerospace alloys such as 316L stainless steel, 2219 aluminum alloy, and Ti64 titanium alloy. In situ data collection includes videography, thermography, and reference thermocouples to build a thermal model of the welds. This will elucidate the relevant physics when combined with post-flight microstructural examination and mechanical testing. MSFC is also progressing towards a suborbital flight experiment of LBW under vacuum, which could provide reams of data on ISW during sustained, high-quality reduced gravity. The effect of combined thermal (cryogenic and high-temperature) and vacuum exposure on both LBW (NASA-funded) and LF (DARPA-funded) is being investigated through ground experiments. In addition to the copious data collected during these ground experiments, ruggedization of LBW hardware will also be pursued. The datasets from these experiments will be used to validate computational models which will inform future ISLM efforts in an ICME framework. A variety of techniques across lengths scales, from CALPHAD-driven thermodynamics & kinetics to phase field modeling of solidification to kinetic Monte Carlo simulations of grain evolution at the mesoscale, will be employed to accelerate the infusion and eventual Q&C of LBW and LF for use in space. The development of data-driven surrogate models to bridge ground to flight experiments and thereby reduce the need for resource-intensive experiments in space will also be investigated. These ICME techniques, surrogate models, and datasets from ground testing can also be employed to advance manufacturing in terrestrial environments.

in-space welding↗
NASA NTRSDate not supplied

The Containment Assurance Risk Framework of the Mars Sample Return Program

The Mars Sample Return campaign aims at bringing rock and atmospheric samples from Mars to Earth through a series of robotic missions. These missions would collect the samples being cached and deposited on Martian soil by the Perseverance rover, place them in a container, and launch them into Martian orbit for subsequent capture by an orbiter that would bring them back. Given there exists a non-zero probability that the samples contain biological material, precautions are being taken to design systems that would break the chain of contact between Mars and Earth. These include techniques such as sterilization of Martian particles, redundant containment vessels, and a robust reentry capsule capable of accurate landings without a parachute. Requirements exist that the probability of containment not assured of Martian-contaminated material into Earth’s biosphere be less than one in a million. To demonstrate compliance with this strict requirement, a statistical framework was developed to assess the likelihood of containment loss during each sample return phase and make a statement about the total combined mission probability of containment not assured. The work presented here describes this framework, which considers failure modes or fault conditions that can initiate failure sequences ultimately leading to containment not assured. Reliability estimates are generated from databases, design heritage, component specifications, or expert opinion in the form of probability density functions or point estimates and provided as inputs to the mathematical models that simulate the different failure sequences. The probabilistic outputs are then combined following the logic of several fault trees to compute the ultimate probability of containment not assured. Given the multidisciplinary nature of the problem and the different types of mathematical models used, the statistical tools needed for analysis are required to be computationally efficient. While standard Monte Carlo approaches are used for fast models, a multi-fidelity approach to rare event probabilities is proposed for expensive models. In this paradigm, inexpensive low-fidelity models are developed for computational acceleration purposes while the expensive high-fidelity model is kept in the loop to retain accuracy in the results. This work presents an example of end-to-end application of this framework highlighting the computational benefits of a multi-fidelity approach.

Giuseppe Cataldo↗
NASA NTRSDate not supplied

Application of a Bayesian Framework for Plasticity Model Selection

Interpretable Machine Learning (IML) has performed well when tasked with deriving constitutive material models. However, IML has been shown to prefer models that overfit noise in data, which tends to lead to bloat and a decrease in interpretability. Due to these issues, the ability of IML to reliably derive models that fit the data and are both interpretable and generalizable is limited. A method developed recently has shown promise to improve upon traditional IML by using a Bayesian fitness definition for the evolution of free-form models with non-deterministic parameters. This framework was developed for genetic-programming-based symbolic regression(GPSR) and involves model parameter estimation using Sequential Monte Carlo sampling (SMC).The method has demonstrated a reduction in bloat when dealing with noisy data in comparison to conventional GPSR. The results of this framework applied to stress-strain data for copper show models that more effectively predict the experimental data better than was previously shown with GPSR.

plasticity↗
NASA NTRS2010

Development of a Human Motor Model for the Evaluation of an Integrated Alerting and Notification Flight Deck System

A human motor model was developed on the basis of performance data that was collected in a flight simulator. The motor model is under consideration as one component of a virtual pilot model for the evaluation of NextGen crew alerting and notification systems in flight decks. This model may be used in a digital Monte Carlo simulation to compare flight deck layout design alternatives. The virtual pilot model is being developed as part of a NASA project to evaluate multiple crews alerting and notification flight deck configurations. Model parameters were derived from empirical distributions of pilot data collected in a flight simulator experiment. The goal of this model is to simulate pilot motor performance in the approach-to-landing task. The unique challenges associated with modeling the complex dynamics of humans interacting with the cockpit environment are discussed, along with the current state and future direction of the model.

Daiker, Ron↗
NASA NTRS1998

Magellan Aerodynamic Characteristics During the Termination Experiment Including Thruster Plume-Free Stream Interaction

Results are presented on the aerodynamic characteristics of the Magellan spacecraft during the October 1994 Termination Experiment, including the effects of the thruster engine exhaust plumes upon the molecular free stream around the spacecraft and upon the aerodynamics coefficients. As Magellan passed through the Venusian atmosphere, the solar arrays were turned in opposite directions relative to the free stream creating a torque on the spacecraft. The spacecraft control system was programmed to counter the effects of this torque with attitude control engines to maintain an inertially fixed attitude. The orientation and reaction engine telemetry returned from Magellan are used to create a model of the aerodynamic torques. Geometric models of the Magellan spacecraft are analyzed with the aid of both free molecular and Direct Simulation Monte Carlo codes. The simulated aerodynamic torques determined are compared to the measured torques. The Direct Simulation Monte Carlo method is also used to model the attitude engine exhaust plumes, the free stream disturbance caused by these plumes, and the resulting torques acting on the spacecraft compared to no-exhaust plume cases. The effect of the exhaust plumes was found to be sufficiently large that thrust reversal is possible.

Cestero, Francisco J.↗
NASA NTRS2022

Automatic variable selection in ecological niche modeling: A case study using Cassin’s Sparrow (Peucaea cassinii)

MERRA/Max provides a feature selection approach to dimensionality reduction that enables direct use of global climate model outputs in ecological niche modeling. The system accomplishes this reduction through a Monte Carlo optimization in which many independent MaxEnt runs, operating on a species occurrence file and a small set of randomly selected variables in a large collection of variables, converge on an estimate of the top contributing predictors in the larger collection. These top predictors can be viewed as potential candidates in the variable selection step of the ecological niche modeling process. MERRA/Max’s Monte Carlo algorithm operates on files stored in the underlying filesystem, making it scalable to large data sets. Its software components can run as parallel processes in a high-performance cloud computing environment to yield near real-time performance. In tests using Cassin’s Sparrow (Peucaea cassinii) as the target species, MERRA/Max selected a set of predictors from Worldclim’s Bioclim collection of 19 environmental variables that have been shown to be important determinants of the species’ bioclimatic niche. It also selected biologically and ecologically plausible predictors from a more diverse set of 86 environmental variables derived from NASA’s Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) reanalysis, an output product of the Goddard Earth Observing System Version 5 (GEOS-5) modeling system. We believe these results point to a technological approach that could expand the use global climate model outputs in ecological niche modeling, foster exploratory experimentation with otherwise difficult-to-use climate data sets, streamline the modeling process, and, eventually, enable automated bioclimatic modeling as a practical, readily accessible, low-cost, commercial cloud service.

John L. Schnase↗