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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

X-ray microscope for solidification studies

This report covers the second 6 month period for the year March 1, 1994 to February 28, 1995. The material outlined in this semi-annual report continues from the previous semi-annual report. The Fein Focus Inc. x-ray source was delivered in September and coincides with the beginning of the second 6 month effort. As a result, and as outlined in the statement of work, this period was dedicated to the evaluation, testing and calibration of the x-ray source. In addition, in this period the modeling effort was continued and extended by the Tiger series of Monte-Carlo simulation programs for photon and electron interactions with materials obtained from the Oak Ridge RISC Library. Some further calculations were also made with the absorption model.

Kaukler, William↗
NASA NTRS1995

Comparisons of the Maxwell and CLL gas/surface interaction models using DSMC

The behavior of two different models of gas-surface interactions is studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic forces and heat transfer between the Maxwell and the Cercignani-Lampis-Lord (CLL) models for flat plate configurations at freestream conditions corresponding to a 140 km orbit around Venus. The size of the flat plate represents one of the solar panels on the Magellan spacecraft, and the freestream conditions correspond to those experienced during aerobraking maneuvers. Results are presented for both a single flat plate and a two-plate configuration as a function of angle of attack and gas-surface accommodation coefficients. The two-plate system is not representative of the Magellan geometry but is studied to explore possible experiments that might be used to differentiate between the two gas-surface interaction models. The Maxwell and CLL models produce qualitatively similar results for the aerodynamic forces and heat transfer on a single flat plate. However, the flow fields produced with the two models are qualitatively different for both the single-plate and two-plate calculations. These differences in the flowfield lead to predictions of the angle of attack for maximum heat transfer in a two plate configuration that are distinctly different for the two gas-surface interactions models.

Hedahl, Marc O.↗
NASA NTRS1992

The statistical analysis of circadian phase and amplitude in constant-routine core-temperature data

Accurate estimation of the phases and amplitude of the endogenous circadian pacemaker from constant-routine core-temperature series is crucial for making inferences about the properties of the human biological clock from data collected under this protocol. This paper presents a set of statistical methods based on a harmonic-regression-plus-correlated-noise model for estimating the phases and the amplitude of the endogenous circadian pacemaker from constant-routine core-temperature data. The methods include a Bayesian Monte Carlo procedure for computing the uncertainty in these circadian functions. We illustrate the techniques with a detailed study of a single subject's core-temperature series and describe their relationship to other statistical methods for circadian data analysis. In our laboratory, these methods have been successfully used to analyze more than 300 constant routines and provide a highly reliable means of extracting phase and amplitude information from core-temperature data.

NASA Discipline Regulatory Physiology↗
NASA NTRS2018

6DOF Testing of the SLS Inertial Navigation Unit

The Navigation System on the NASA Space Launch System (SLS) Block 1 vehicle performs initial alignment of the Inertial Navigation System (INS) navigation frame through gyrocompass alignment (GCA). Because the navigation architecture for the SLS Block 1 vehicle is a purely inertial system, the accuracy of the achieved orbit relative to mission requirements is very sensitive to initial alignment accuracy. The assessment of this sensitivity and many others via simulation is a part of the SLS Model-Based Design and Model-Based Requirements approach. As a part of the aforementioned, 6DOF Monte Carlo simulation is used in large part to develop and demonstrate verification of program requirements. To facilitate this and the GN&C flight software design process, an SLS-Program-controlled Design Math Model (DMM) of the SLS INS was developed by the SLS Navigation Team. The SLS INS model implements all of the key functions of the hardware-namely, GCA, inertial navigation, and FDIR (Fault Detection, Isolation, and Recovery)-in support of SLS GN&C design requirements verification. Despite the strong sensitivity to initial alignment, GCA accuracy requirements were not verified by test due to program cost and schedule constraints. Instead, the system relies upon assessments performed using the SLS INS model. In order to verify SLS program requirements by analysis, the SLS INS model is verified and validated against flight hardware. In lieu of direct testing of GCA accuracy in support of requirement verification, the SLS Navigation Team proposed and conducted an engineering test to, among other things, validate the GCA performance and overall behavior of the SLS INS model through comparison with test data. This paper will detail dynamic hardware testing of the SLS INS, conducted by the SLS Navigation Team at Marshall Space Flight Center's 6DOF Table Facility, in support of GCA performance characterization and INS model validation. A 6-DOF motion platform was used to produce 6DOF pad twist and sway dynamics while a simulated SLS flight computer communicated with the INS. Tests conducted include an evaluation of GCA algorithm robustness to increasingly dynamic pad environments, an examination of GCA algorithm stability and accuracy over long durations, and a long-duration static test to gather enough data for Allan Variance analysis. Test setup, execution, and data analysis will be discussed, including analysis performed in support of SLS INS model validation.

Geohagan, Kevin↗
NASA NTRS1993

PDF approach for compressible turbulent reacting flows

The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.

Hsu, A. T.↗
NASA NTRS2004

Radiolytic Model for Chemical Composition of Europa's Atmosphere and Surface

The overall objective of the present effort is to produce models for major and selected minor components of Europa s neutral atmosphere in 1-D versus altitude and in 2-D versus altitude and longitude or latitude. A 3-D model versus all three coordinates (alt, long, lat) will be studied but development on this is at present limited by computing facilities available to the investigation team. In this first year we have focused on 1-D modeling with Co-I Valery Shematovich s Direct Simulation Monte Carlo (DSMC) code for water group species (H2O, O2, O, OH) and on 2-D with Co-I Mau Wong's version of a similar code for O2, O, CO, CO2, and Na. Surface source rates of H2O and O2 from sputtering and radiolysis are used in the 1-D model, while observations for CO2 at the Europa surface and Na detected in a neutral cloud ejected from Europa are used, along with the O2 sputtering rate, to constrain source rates in the 2-D version. With these separate approaches we are investigating a range of processes important to eventual implementation of a comprehensive 3-D atmospheric model which could be used to understand present observations and develop science requirements for future observations, e.g. from Earth and in Europa orbit. Within the second year we expect to merge the full water group calculations into the 2-D version of the DSMC code which can then be extended to 3-D, pending availability of computing resources. Another important goal in the second year would be the inclusion of sulk and its more volatile oxides (SO, SO2).

Cooper, John F.↗
NASA NTRSDate not supplied

SWIPE: Spectral Water Inversion Processor and Emulator

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will be discussing the progress made developing SWIPE: Spectral Water Inversion Processor and Emulator. SWIPE is a platform for advanced modeling of coastal and inland aquatic habitats. The goal is create a comprehensive and cohesive system to leverage recent advancements in computation and machine learning to develop a synthetic training ground for sensitivity studies and algorithm development. The four principal facets of SWIPE include: 1. Advanced two-layer coated sphere bio-optical modeling and GPU radiative transfer modeling, 2. Big Data involving massive synthetic spectral libraries of optical properties of various global aquatic particles, surface reflectance, and top-of-atmosphere reflectance, all at hyperspectral resolution leveraging high-end computing systems at NASA Ames Research Center, 3. Deep Learning for algorithm development for water quality inversion of concentrations of common biogeophysical variables as well as optics, full uncertainty characterization by water type, and forward emulation, and lastly, 4. Image Processing for application of developed retrieval algorithms for both hyperspectral and multispectral sensors with experimental corrections for global adjacency, noise, sunglint, and benthic reflectance. This presentation will demonstrate the Equivalent Algal Populations (EAP) two-layer coated sphere scattering model which has been used develop spectral libraries of hyperspectral inherent optical properties of roughly 80 species of phytoplankton, covering 15 different classes and nine taxonomic functional types. The EAP model was also used to derive spectral properties of 10 different non-algal particle functional types. Examples of how the SMART-G (Speed-up Monte-carlo Advanced Radiative Transfer using GPU) radiative transfer code is used to model optically complex aquatic signals will be presented and discussed in the context of creating a massive synthetic database which can leverage the full power of next generation machine learning techniques and high end computing for water quality inversion. We will discuss our active investigation in things like appropriate model architectures, dimensionality reduction techniques such as PCA and autoencoders, uncertainty quantification and abstaining, and which variables actually benefit most from hyperspectral information versus multispectral resolution. We are also curious about questions relating to cost/benefit analysis in terms of computation resources, neural network complexity, and data volumes. Answers to these questions will hopefully elaborate on cost efficiency for potential future sensor design considerations.

SWIPE↗
NASA NTRS1991

Fluid hydrogen at high density - Pressure dissociation

A model for the Helmholtz free energy of fluid hydrogen at high density and high temperature is developed. This model aims at describing both pressure and temperature dissociation and ionization and bears directly on equations of state of partially ionized plasmas, as encountered in astrophysical situations and high-pressure experiments. This paper focuses on a mixture of hydrogen atoms and molecules and is devoted to the study of the phenomenon of pressure dissociation at finite temperatures. In the present model, the strong interactions are described with realistic potentials and are computed with a modified Weeks-Chandler-Andersen fluid perturbation theory that reproduces Monte Carlo simulations to better than 3 percent. Theoretical Hugoniot curves derived from the model are in excellent agreement with experimental data.

Saumon, Didier↗
NASA NTRS1991

Bremsstrahlung and energetic electrons in supernovae

Explicit Monte Carlo calculations are presented for the rate of Compton scattering and for the spectrum recoil electrons per gram at discrete depths within models of supernovae, with emphasis on the bremsstrahlung emissivity within supernova interiors. The spectrum of the bremsstrahlung emissivity is calculated at E to the -1.3 for a wide range of supernova models. The radial and temporal dependence of that emissivity is calculated for two models of SN 1987A and for one model of a Type Ia and of a He core explosion model of a Type Ib. The spectrum emerging from the surface, which is then evaluated by Monte Carlo techniques, was found to be sensitive to the photoelectric opacity and inverted by it, i.e., the spectral luminosity of the bremsstrahlung component increases slowly between 1 and 30 keV even though the emissivity declines with energy within that energy band. The bremsstrahlung luminosity was found to dominate that of the primary scattered gammas for E of less than 20 keV.

Clayton, Donald D.↗
NASA NTRS2016

Characterizing Epistemic Uncertainty for Launch Vehicle Designs

NASA Probabilistic Risk Assessment (PRA) has the task of estimating the aleatory (randomness) and epistemic (lack of knowledge) uncertainty of launch vehicle loss of mission and crew risk, and communicating the results. Launch vehicles are complex engineered systems designed with sophisticated subsystems that are built to work together to accomplish mission success. Some of these systems or subsystems are in the form of heritage equipment, while some have never been previously launched. For these cases, characterizing the epistemic uncertainty is of foremost importance, and it is anticipated that the epistemic uncertainty of a modified launch vehicle design versus a design of well understood heritage equipment would be greater. For reasons that will be discussed, standard uncertainty propagation methods using Monte Carlo simulation produce counter intuitive results, and significantly underestimate epistemic uncertainty for launch vehicle models. Furthermore, standard PRA methods, such as Uncertainty-Importance analyses used to identify components that are significant contributors to uncertainty, are rendered obsolete, since sensitivity to uncertainty changes are not reflected in propagation of uncertainty using Monte Carlo methods. This paper provides a basis of the uncertainty underestimation for complex systems and especially, due to nuances of launch vehicle logic, for launch vehicles. It then suggests several alternative methods for estimating uncertainty and provides examples of estimation results. Lastly, the paper describes how to implement an Uncertainty-Importance analysis using one alternative approach, describes the results, and suggests ways to reduce epistemic uncertainty by focusing on additional data or testing of selected components.

Novack, Steven D.↗
NASA NTRS2016

Characterizing Epistemic Uncertainty for Launch Vehicle Designs

NASA Probabilistic Risk Assessment (PRA) has the task of estimating the aleatory (randomness) and epistemic (lack of knowledge) uncertainty of launch vehicle loss of mission and crew risk and communicating the results. Launch vehicles are complex engineered systems designed with sophisticated subsystems that are built to work together to accomplish mission success. Some of these systems or subsystems are in the form of heritage equipment, while some have never been previously launched. For these cases, characterizing the epistemic uncertainty is of foremost importance, and it is anticipated that the epistemic uncertainty of a modified launch vehicle design versus a design of well understood heritage equipment would be greater. For reasons that will be discussed, standard uncertainty propagation methods using Monte Carlo simulation produce counter intuitive results and significantly underestimate epistemic uncertainty for launch vehicle models. Furthermore, standard PRA methods such as Uncertainty-Importance analyses used to identify components that are significant contributors to uncertainty are rendered obsolete since sensitivity to uncertainty changes are not reflected in propagation of uncertainty using Monte Carlo methods.This paper provides a basis of the uncertainty underestimation for complex systems and especially, due to nuances of launch vehicle logic, for launch vehicles. It then suggests several alternative methods for estimating uncertainty and provides examples of estimation results. Lastly, the paper shows how to implement an Uncertainty-Importance analysis using one alternative approach, describes the results, and suggests ways to reduce epistemic uncertainty by focusing on additional data or testing of selected components.

Novack, Steven D.↗
NASA NTRS2018

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗
NASA NTRS2016

Influence of Population Variation of Physiological Parameters in Computational Models of Space Physiology

The earliest manifestations of Visual Impairment and Intracranial Pressure (VIIP) syndrome become evident after months of spaceflight and include a variety of ophthalmic changes, including posterior globe flattening and distension of the optic nerve sheath. Prevailing evidence links the occurrence of VIIP to the cephalic fluid shift induced by microgravity and the subsequent pressure changes around the optic nerve and eye. Deducing the etiology of VIIP is challenging due to the wide range of physiological parameters that may be influenced by spaceflight and are required to address a realistic spectrum of physiological responses. Here, we report on the application of an efficient approach to interrogating physiological parameter space through computational modeling. Specifically, we assess the influence of uncertainty in input parameters for two models of VIIP syndrome: a lumped-parameter model (LPM) of the cardiovascular and central nervous systems, and a finite-element model (FEM) of the posterior eye, optic nerve head (ONH) and optic nerve sheath. Methods: To investigate the parameter space in each model, we employed Latin hypercube sampling partial rank correlation coefficient (LHSPRCC) strategies. LHS techniques outperform Monte Carlo approaches by enforcing efficient sampling across the entire range of all parameters. The PRCC method estimates the sensitivity of model outputs to these parameters while adjusting for the linear effects of all other inputs. The LPM analysis addressed uncertainties in 42 physiological parameters, such as initial compartmental volume and nominal compartment percentage of total cardiac output in the supine state, while the FEM evaluated the effects on biomechanical strain from uncertainties in 23 material and pressure parameters for the ocular anatomy. Results and Conclusion: The LPM analysis identified several key factors including high sensitivity to the initial fluid distribution. The FEM study found that intraocular pressure and intracranial pressure had dominant impact on the peak strains in the ONH and retro-laminar optic nerve, respectively; optic nerve and lamina cribrosa stiffness were also important. This investigation illustrates the ability of LHSPRCC to identify the most influential physiological parameters, which must therefore be well-characterized to produce the most accurate numerical results.

biomechanics↗
NASA NTRS2019

Investigation of Direct Force Control for Planetary Aerocapture at Neptune

In this work, a direct force control numerical predictor-corrector guidance architecture is developed to enable Neptune aerocapture using flight-heritage blunt body aeroshells. A linear aerodynamics model is formulated for a Mars Science Laboratory-derived aeroshell. The application of calculus of variations shows that the optimal angle of attack and side-slip angle control laws are bang-bang. A closed-loop numerical predictor-corrector direct force control guidance algorithm is developed and numerically simulated using the Program to Optimize Simulated Trajectories II. The Monte Carlo simulated trajectories are demonstrated to be robust to the modeled dispersions in aerodynamics, atmospheric density, and entry state. An aerocapture technology trade study demonstrates that blunt body direct force control aerocapture enables similar performance as slender body bank angle control but halves the peak g-loading.

Deshmukh, Rohan G.↗
NASA NTRS1976

Statistical mechanics of light elements at high pressure. IV - A model free energy for the metallic phase

A large quantity of data on the thermodynamic properties of hydrogen-helium metallic liquids have been obtained in extended computer calculations in which a Monte Carlo code essentially identical to that described by Hubbard (1972) was used. A model free energy for metallic hydrogen with a relatively small mass fraction of helium is discussed, taking into account the definition of variables, a procedure for choosing the free energy, values for the fitting parameters, and the evaluation of the entropy constants. Possibilities concerning a use of the obtained data in studies of the interiors of the outer planets are briefly considered.

Dewitt, H. E.↗
NASA NTRS1989

'Hypersonic slip flows' and issues on extending continuum model beyond the Navier-Stokes level

Results of a direct simulation Monte Carlo method for a hypersonic flow about a flat plate at a 40 deg angle of attack have been compared with corresponding results from a theory for fully viscous shock layers (FVSLs). Using the 13-moment equations for a Maxwell gas, it is demonstrated that nonequilibrium and equilibrium FVSL flows can be correlated. With the exception of the pressure density, most of the flow properties along a streamline (including heat flux and shear and normal stresses) are correctly predicted to leading order by the Navier-Stokes model.

Cheng, H. K.↗
NASA NTRS1986

The Oort cloud in transition

The evolution of theoretical and empirical models of the Oort cloud (OC) since it was first proposed by Oort in 1950 is traced, and the main features of current models are discussed, in a general review. Consideration is given to work on the classical OC (Monte Carlo simulations of OC evolution, population and mass estimates, and OC perturbation by passing stars and giant molecular clouds), models of a massive inner OC (simulations of planetesimal-swarm evolution in the Uranus-Neptune zone and IRAS observations of circumstellar dust shells), evidence for random and/or periodic comet showers, and the possible role of the Galactic missing mass. The current OC model comprises an almost spherical outer (10,000-100,000-AU) cloud of mass 7-8 earth mass and population (1.4-2.3) x 10 to the 12th, and a flat disklike inner (40-10,000 AU) cloud of mass 100-200 earth mass and population (1-10) x 10 to the 13th.

Weissman, P. R.↗
NASA NTRS1974

The economic value of remote sensing of earth resources from space: An ERTS overview and the value of continuity of service. Volume 3: Intensive use of living resources, agriculture. Part 2: Distribution effects

The results of an investigation of the value of improving information for forecasting future crop harvests are described. A theoretical model is developed to calculate the value of increased speed of availablitiy of that information. The analysis of U.S. domestic wheat consumption was implemented. New estimates of a demand function for wheat and of a cost of storage function were involved, along with a Monte Carlo simulation for the wheat spot and future markets and a model of market determinations of wheat inventories. Results are shown to depend critically on the accuracy of current and proposed measurement techniques.

Bradford, D. F.↗