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 703 records · Page 39

Utilizing Mars Global Reference Atmospheric Model (Mars-GRAM 2005) to Evaluate Entry Probe Mission Sites

Engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL)1. Traditional Mars-GRAM options for representing the mean atmosphere along entry corridors include: a) TES Mapping Years 1 and 2, with Mars-GRAM data coming from MGCM model results driven by observed TES dust optical depth; and b) TES Mapping Year 0, with user-controlled dust optical depth and Mars-GRAM data interpolated from MGCM model results driven by selected values of globally-uniform dust optical depth. From the surface to 80 km altitude, Mars-GRAM is based on NASA Ames Mars General Circulation Model (MGCM). Mars-GRAM and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA), with altitudes referenced to the MOLA areoid, or constant potential surface. Mars-GRAM 2005 has been validated2 against Radio Science data, and both nadir and limb data from the Thermal Emission Spectrometer (TES)

Justh, Hilary L.↗

Discrete Event Simulation-Based Timeline Validation Using R2U2

The Gateway Vehicle Systems Manager (VSM), the top-level software control system in a distributed, hierarchical Autonomous System Management Architecture is, like most modern spacecraft software control systems, heavily data-driven. For example, schedules (timelines) will be developed on the ground and, due to the high degree of autonomy, contain complex procedures involving conditional branching, variable timing, and resource contention resolution. In order to verify that an uploaded timeline will function correctly, it is necessary to explore the feasible set of possible executions. While it is possible to test a timeline using a mission simulation, the complexity of the system and duration of a timeline limits the number of trials and therefore the test coverage. To address this problem, the VSM team is using a discrete event system model that can rapidly generate from a timeline sets of event sequences using Monte Carlo techniques. To achieve rapid and trustworthy checking of the event sequences, we use an offline version of the runtime model checking tool R2U2. This presentation describes the approach the VSM team is using to implement the discrete event simulation and evaluate event sequences using R2U2. The presentation will discuss: 1. Description of the timelines by VSM in the context of VSM operations 2. Expansion of a timeline into a sequence of atomic events 3. Adjustment, in the Monte Carlo environment, of an event sequence to account for uncertainty, external events, and failures 4. Definition of R2U2 input and mission-time linear temporal logic files 5. Generation and use of R2U2 verdict sequences 6. Lessons learned and future work

Verification↗

Investigation of Next-Generation Earth Radiation Budget Radiometry

The current effort addresses two issues important to the research conducted by the Thermal Radiation Group at Virginia Tech. The first research topic involves the development of a method which can properly model the diffraction of radiation as it enters an instrument aperture. The second topic involves the study of a potential next-generation space-borne radiometric instrument concept. Presented are multiple modeling efforts to describe the diffraction of monochromatic radiant energy passing through an aperture for use in the Monte-Carlo ray-trace environment. Described in detail is a deterministic model based upon Heisenberg's uncertainty principle and the particle theory of light. This method is applicable to either Fraunhofer or Fresnel diffraction situations, but is incapable of predicting the secondary fringes in a diffraction pattern. Also presented is a second diffraction model, based on the Huygens-Fresnel principle with a correcting obliquity factor. This model is useful for predicting Fraunhofer diffraction, and can predict the secondary fringes because it keeps track of phase. NASA is planning for the next-generation of instruments to follow CERES (Clouds and the Earth's Radiant Energy System), an instrument which measures components of the Earth's radiant energy budget in three spectral bands. A potential next-generation concept involves modification of the current CERES instrument to measure in a larger number of wavelength bands. This increased spectral partitioning would be achieved by the addition of filters and detectors to the current CERES geometry. The capacity of the CERES telescope to serve for this purpose is addressed in this thesis.

Coffey, Katherine L.↗

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.↗

Component-level Performance and Mass Sensitivity Analysis of NEP MW-class Power System

Nuclear electric propulsion (NEP) is a promising option towards enabling missions to Mars and is an area of interest for NASA’s Space Nuclear Propulsion project. This project is currently investigating technology development opportunities for an NEP vehicle. Physics-based modeling can be used in the early stages of technology development to gain understanding of the effects of technology and performance assumptions on the system performance and mass. This information can then inform technology maturation planning for near term development. By using a Brayton power conversion model and vehicle mass model for megawatt class NEP applications, a sensitivity analysis is performed to assess the impact of individual components’ performance on the power conversion system performance and system mass. A Monte Carlo simulation is also used to determine the variability in system mass based on uncertainty within the modeling parameters. The sensitivity analysis shows a high sensitivity to power conversion inlet temperature, compressor inlet temperature, and recuperator performance. A Monte Carlo analysis suggests a range of -10% to +15% for a 90% confidence interval on system mass based on the uncertainties in the model inputs.

Nuclear electric propulsion↗

Nuclear electromagnetic cascades from alpha particles incident on an iron absorber

Mean values and fluctuations of the nuclear-electromagnetic cascade development resulting from alpha particles incident on an iron absorber was determined for the primary energy range 10 to 300 GeV/nucleon. This has been accomplished by using a three-dimensional Monte Carlo simulation of the cascade process. The model used was first adjusted until it gave predictions for a small ionization spectrometer which agreed with measurements obtained by exposing the same spectrometer to 10, 20.5, and 28 GeV/c protons at an accelerator. The calcuations were then modified to apply to high energy alpha particles. The results show that measurements of the energies of alpha particles incident on an iron spectrometer with six interaction lengths total depth can be made with accuracies ranging from approximately 5% at 300 GeV/nucleon to approximately 16% at 10 GeV/nucleon.

Jones, W. V.↗

Production of high-energy gamma rays by cosmic ray interactions in the atmosphere and lunar surface

The production of gamma rays above 10 MeV by the interaction of cosmic rays in the atmosphere and the lunar surface is simulated using Monte Carlo methods. The calculation incorporates a new model of high-energy nucleon-nucleus interactions based on empirical fits to inclusive cross sections for the production of pions and nucleons. The atmospheric gamma ray flux is calculated as a function of direction, energy, and atmospheric depth. These calculations are compared with observations from balloons and from the SAS 2 satellite. Estimates of the flux of earth albedo electrons produced by cosmic ray interactions are presented. The lunar gamma ray albedo is calculated and compared with an upper limit based on SAS 2 measurements.

Morris, D. J.↗

Impact of flying qualities on mission effectiveness for helicopter air combat, volume 1

A computer simulation to investigate the impact of flying qualities on mission effectiveness is described. The objective of the study was to relate the effects of flying qualities, such as precision of flight path control and pilot workload, to the ability of a single Scout helicopter, or helicopter team, to accomplish a specified anti-armor mission successfully. The model of the actual engagement is a Monte Carlo simulation that has the capability to assess the effects of helicopter characteristics, numbers, tactics and weaponization on the force's ability to accomplish a specific mission against a specified threat as a function of realistic tactical factors. A key feature of this program is a simulation of micro-terrain features and their effects on detection, exposure, and masking for nap-of-the-earth (NOE) flight.

Harris, T. M.↗

Direct simulation of rotational and vibrational nonequilibrium

The ways in which energy transfer is calculated in the Direct Simulation Monte Carlo method are presented. An energy exchange model that deals with translational and rotational modes is described. A model for simulating the transfer of energy between the translational and vibrational modes is presented as well.

Boyd, Iain D.↗

A correction algorithm for particle size distribution measurements made with the forward-scattering spectrometer probe

A correction algorithm for evaluating the particle size distribution measurements of atmospheric aerosols obtained with a forward-scattering spectrometer probe (FSSP) is examined. A model based on Poisson statistics is employed to calculate the average diameter and rms width of the particle size distribution. The dead time and coincidence errors in the measured number density are estimated. The model generated data are compared with a Monte Carlo simulation of the FSSP operation. It is observed that the correlation between the actual and measured size distribution is nonlinear. It is noted that the algorithm permits more accurate calculation of the average diameter and rms width of the distribution compared to uncorrected measured quantities.

Lock, James A.↗

The solar ionisation rate deduced from Ulysses measurements and its implications to interplanetary Lyman alpha-intensity

We have developed a computer code which can be used to study 3-dimensional and time-dependent effects of the solar cycle on the interplanetary (IP) hydrogen distribution. The code is based on the inverted Monte Carlo simulation. In this work we have modelled the temporal behaviour of the solar ionisation rate. We have assumed that during the most of the time of the solar cycle there is an anisotopic latitudinal structure but right at the solar maximum the anisotropy disappears. The effects of this behaviour will be discussed both in regard to the IP hydrogen distribution and IP Lyman a a-intensity.

Summanen, T.↗

Prototype Conflict Alerting Logic for Free Flight

This paper discusses the development of a prototype alerting system for a conceptual Free Flight environment. The concept assumes that datalink between aircraft is available and that conflicts are primarily resolved on the flight deck. Four alert stages are generated depending on the likelihood of a conflict. If the conflict is not resolved by the flight crews, Air Traffic Control is notified to take over separation authority. The alerting logic is based on probabilistic analysis through modeling of aircraft sensor and trajectory uncertainties. Monte Carlo simulations were used over a range of encounter situations to determine conflict probability. The four alert stages were then defined based on probability of conflict and on the number of avoidance maneuvers available to the flight crew. Preliminary results from numerical evaluations and from a piloted simulator study at NASA Ames Research Center are summarized.

Yang, Lee C.↗

Evaluation of the Effect of Water Vapor On the Performance of NASA's NMRO Catalysts for Carbon Monoxide Oxidation

The Noble Metal Reducible Oxide (NMRO) catalysts for the low temperature oxidation of carbon monoxide were developed by NASA for the reoxidation of carbon monoxide which forms by the dissociation of carbon dioxide during the operation of sealed carbon dioxide lasers. The NMRO catalyst, which consists of a noble metal in conjunction with a reducible metal oxide, was evaluated under conditions that will be encountered in a CO2 laser operation, namely temperatures in the range 298 to 373 K and no significant reaction gas components other than CO, CO2 and O2. The NMRO catalysts may have significant potential for spin-off applications such as the prevention of carbon monoxide build-up in closed spaces as in space vehicle cabins or submarines, and the elimination of the cold start-up problem of automobile exhaust catalysts. The most significant difference in the conditions of these possible future applications is the high moisture content of the gases to be processed. Lack of understanding of the effects of water vapor and high temperature on catalyst activity and operation for extended periods are currently the main stumbling blocks for the transfer of this NASA technology to be used for commercial purposes. In the original proposal the following objectives were stated: To obtain experimental data on the adsorption, desorption and reaction characteristics of CO and CO2 the catalysts under high moisture conditions; to collect evidence on the presence of carbonate and hydroxyl surface species and their involvement in the CO oxidation mechanism; and to model the reaction system using a Monte-Carlo simulation to gain insight on the various steps involved. After the work has commenced the NASA technical monitor Mr. David Scheyer informed us that there was increased interest in the possible use of the NMRO catalysts as automobile exhaust catalysts and therefore NASA wanted to know whether the catalysts can operate at high temperatures as well as with high moisture gases. At that meeting it was decided that investigation of the high temperature performance of the NMRO catalysts should be given priority and replace the Monte-Carlo simulation objective. As a result, the modified objectives of the investigation were taken as the investigation of the high-temperature activity of the NMRO catalysts, and the effect of water vapor on the performance of these catalysts.

Akyurtlu, Ates↗

Estimates of M2 Tidal Energy Dissipation from TOPEX/Poseidon Altimeter Data

Most of the tidal energy dissipation in the ocean occurs in shallow seas, as has long been recognized. However, recent work has suggested that a significant fraction of the dissipation, perhaps 1 TW or more, occurs in the deep ocean. This paper builds further evidence for that conclusion. More than 6 years of data from the TOPEX/Poseidon satellite altimeter are used to map the tidal dissipation rate throughout the world ocean. The dissipation rate is estimated as a balance between the rate of working by tidal forces and the energy flux divergence, computed using currents derived by least squares fitting of the altimeter data and the shallow water equations. Such calculations require dynamical assumptions, in particular about the nature of dissipation. To assess sensitivity of dissipation estimates to input assumptions, a large suite of tidal inversions based on a wide range of drag parameterizations and employing both real and synthetic altimeter data are compared. These experiments and Monte Carlo error fields from a generalized inverse model are used to establish error uncertainties for the dissipation estimates. Owing to the tight constraints on tidal elevation fields provided by the altimeter, area integrals of the energy balance are remarkably insensitive to required dynamical assumptions. Tidal energy dissipation is estimated for all major shallow seas (excluding individual polar seas) and compared with previous model and data-based estimates. Dissipation in the open ocean is significantly tnhanced around major bathymetric features, in a manner consistent with simple theories the generation of baroclinic tides.

Egbert, Gary D.↗

State Estimation of International Space Station Centrifuge Rotor With Incomplete Knowledge of Disturbance Inputs

This thesis develops a state estimation algorithm for the Centrifuge Rotor (CR) system where only relative measurements are available with limited knowledge of both rotor imbalance disturbances and International Space Station (ISS) thruster disturbances. A Kalman filter is applied to a plant model augmented with sinusoidal disturbance states used to model both the effect of the rotor imbalance and the 155 thrusters on the CR relative motion measurement. The sinusoidal disturbance states compensate for the lack of the availability of plant inputs for use in the Kalman filter. Testing confirms that complete disturbance modeling is necessary to ensure reliable estimation. Further testing goes on to show that increased estimator operational bandwidth can be achieved through the expansion of the disturbance model within the filter dynamics. In addition, Monte Carlo analysis shows the varying levels of robustness against defined plant/filter uncertainty variations.

SINUSOIDAL DISTURBANCE STATES↗

Methodolgy For Evaluation Of Technology Impacts In Space Electric Power Systems

The Analysis and Management branch of the Power and Propulsion Office at NASA Glenn Research Center is responsible for performing complex analyses of the space power and In-Space propulsion products developed by GRC. This work quantifies the benefits of the advanced technologies to support on-going advocacy efforts. The Power and Propulsion Office is committed to understanding how the advancement in space technologies could benefit future NASA missions. They support many diverse projects and missions throughout NASA as well as industry and academia. The area of work that we are concentrating on is space technology investment strategies. Our goal is to develop a Monte-Carlo based tool to investigate technology impacts in space electric power systems. The framework is being developed at this stage, which will be used to set up a computer simulation of a space electric power system (EPS). The outcome is expected to be a probabilistic assessment of critical technologies and potential development issues. We are developing methods for integrating existing spreadsheet-based tools into the simulation tool. Also, work is being done on defining interface protocols to enable rapid integration of future tools. Monte Carlo-based simulation programs for statistical modeling of the EPS Model. I decided to learn and evaluate Palisade's @Risk and Risk Optimizer software, and utilize it's capabilities for the Electric Power System (EPS) model. I also looked at similar software packages (JMP, SPSS, Crystal Ball, VenSim, Analytica) available from other suppliers and evaluated them. The second task was to develop the framework for the tool, in which we had to define technology characteristics using weighing factors and probability distributions. Also we had to define the simulation space and add hard and soft constraints to the model. The third task is to incorporate (preliminary) cost factors into the model. A final task is developing a cross-platform solution of this framework.

Holda, Julie↗

Analysis of a Hypergolic Propellant Explosion During Processing of Launch Vehicles in the VAB

NASA is developing launch vehicles to support missions to Low Earth Orbit (LEO), the moon and deep space. Whether manned or unmanned, the vehicle components will likely be integrated in the Vehicle Assembly Building (VAB) at Kennedy Space Center (KSC) and typically include a fueled spacecraft (SC) that sits on top of one or more stages. The processing of a fueled SC involves hazardous operations when it is brought into the VAB Transfer Aisle and lifted a significant height for mating with lower stages. Accidents resulting from these hazardous operations could impact unrelated personnel working in buildings adjacent to the VAB. Safe separation distances based on the DOD Explosives Standards Quantity-Distance (Q-D) approach result in large IBD arcs. This paper presents site-specific air blast and fragmentation hazard analyses for comparison with the Q-D arcs as well as consequence and risk analyses to provide added information for the decision maker. A new physics-based fragmentation model is presented that includes: a) the development of a primary fragment list (which defines the fragment characteristics) associated with a hypergolic propellant explosion, b) a description of a 3D fragment bounce model, c) the results of probabilistic Monte-Carlo simulations (that include uncertainties in the fragment characteristics) to determine: i) the hazardous fragment density distance, ii) the expected number of wall/roof impacts and penetrations to over 40 buildings adjacent to the VAB, and iii) the risk to building occupants.

Chrostowski, Jon D.↗