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 721 records · Page 40

Collisional Evolution of the Enceladus Neutral Cloud

Water vapor ejected from Saturn's small moon Enceladus easily escapes its meager gravity to form a Saturn-encircling cloud with a low collision rate. Observations show that the cloud is quite broad in the radial direction, and we show here that collisions, though quite rare, may be largely responsible for this radial spreading. We modeled this cloud using the Direct Simulation Monte Carlo method, as fluid methods would be inappropriate for such a tenuous gas.

magnetic fields↗

Deep Mapping of Small Solar System Bodies with Galactic Cosmic Ray Secondary Particle Showers

Galactic cosmic rays rain steadily from all directions onto asteroids and comets. The interaction of these high-energy ions produces a cascade of secondary particles, including muons, which can penetrate the solid interiors of small solar system bodies. Muons, which are produced in abundance in Earth's atmosphere, have been used to image large structures on Earth, including the Great Pyramids and the interior of volcanoes. In this study, we demonstrate that the transmitted flux of muons is sensitive to the interior density structure of asteroids and comets, less than a few hundred meters in diameter. Muonography has the potential to fill a critical gap in our knowledge of the deep interiors of small bodies, providing information needed for planetary defense, in situ resource utilization, and planetary science. We use Monte Carlo codes (MCNPX and FLUKA), which accurately model galactic cosmic ray showers, to explore systematic variations in the production of muons in solid surfaces. Results of these calculations confirm the scaling of muon production in Earth's atmosphere to solid regolith materials, as predicted by a simple, semi-empirical model. Muons are primarily produced in the top meter of the regoliths of asteroids and comets. Their rate of production is over three orders of magnitude lower than in Earth's atmosphere and depends strongly on regolith density. In practice, the use of muonography to characterize the interiors of small solar system bodies must overcome their low rate of production and their dependence on regolith density, which can vary over the surface of asteroids and comets. We show that interior contrast can be resolved using a muon telescope (hodoscope) with about 1 sq m aperture with integration times ranging from hours to weeks. Design concepts for a practical hodoscope that could be deployed in situ or on an orbiting spacecraft, are described. Regolith density within the top meter of an asteroid can be determined from radar observations. A concept for a pilot mission that combines remote radar measurements with in situ muonography of a near-Earth asteroid is presented. Perceived challenges and next steps for the development of the concept are described.

Prettyman, T. H.↗

Desensitized Optimal Filtering and Sensor Fusion Toolkit

Analytical Mechanics Associates, Inc., has developed a software toolkit that filters and processes navigational data from multiple sensor sources. A key component of the toolkit is a trajectory optimization technique that reduces the sensitivity of Kalman filters with respect to model parameter uncertainties. The sensor fusion toolkit also integrates recent advances in adaptive Kalman and sigma-point filters for non-Gaussian problems with error statistics. This Phase II effort provides new filtering and sensor fusion techniques in a convenient package that can be used as a stand-alone application for ground support and/or onboard use. Its modular architecture enables ready integration with existing tools. A suite of sensor models and noise distribution as well as Monte Carlo analysis capability are included to enable statistical performance evaluations.

Karlgaard, Christopher D.↗

An Analytic Formulation of Ejecta Distributions over Airless Bodies

With recent plans to revisit the Moon by robotic and crewed spacecraft, there has been a renewed interest in understanding the ejecta environment on the Moon and other airless bodies. Meteoroid and asteroid impacts can excavate large amounts of material from the surface and, above an airless body, can send this material long distances on (essentially) ballistic orbits. This ejecta, while typically traveling slower than the impactor, can nevertheless achieve high enough speeds to endanger surface operations. Accurate knowledge of this phenomenon is necessary in order to design appropriate shielding for human activities, both for activities on the surface of the Moon and for orbiters in near-lunar space. This phenomenon is also important in the transport of particles above other Solar System objects, such as Jovian satellites, where this ejecta creates a kind of ever-present “halo” of particles around the gravitating body. While Monte Carlo techniques have been successfully used to model this environment, there are useful analytic expressions, developed for use in modeling the meteoroid environment, that can be used to describe this environment as well. Such analytic tools can shed light on the altitude, velocity, and directionality of these ejecta environments.

Mark Matney↗

Assistive Detect and Avoid for Pilots in the Cockpit

Aircraft not receiving radar services rely on see and avoid and radio coordination via Common Traffic Advisory Frequencies to remain well clear of each other and avoid mid-air collisions. Radio coordination is usually performed in the vicinity of non-towered airports whereas non-radar services en-route operations rely solely on see and avoid. This paper presents the results of a simulation study of the effectiveness of assistive detect and avoid technologies when used to enhance pilots’ ability to see and avoid nearby traffic. Three different experimental conditions are modeled, representing “unaided see and avoid”, “see and avoid with traffic advisories”, and “see and avoid with assistive detect and avoid technology”. The effectiveness of see and avoid is evaluated using a set of head-on, crossing, and overtaking encounter scenarios and a model of visual acquisition embedded in a Monte Carlo simulation. The effectiveness of assistive detect and avoid is estimated for the same encounter scenarios. A prototype system for detect and avoid and a summary of results are presented. Preliminary results strongly suggest that assistive detect and avoid could greatly enhance the capabilities of flight crews to avoid traffic and remain well clear.

collision, detect and avoid, resolution, well clea↗

Assistive Detect and Avoid for Pilots in the Cockpit

Aircraft not receiving radar services rely on see and avoid and radio coordination via Common Traffic Advisory Frequencies to remain well clear of each other and avoid mid-air collisions. Radio coordination is usually performed in the vicinity of non-towered airports whereas non-radar services en-route operations rely solely on see and avoid. This paper presents the results of a simulation study of the effectiveness of assistive detect and avoid technologies when used to enhance pilots’ ability to see and avoid nearby traffic. Three different experimental conditions are modeled, representing “unaided see and avoid”, “see and avoid with traffic advisories”, and “see and avoid with assistive detect and avoid technology”. The effectiveness of see and avoid is evaluated using a set of head-on, crossing, and overtaking encounter scenarios and a model of visual acquisition embedded in a Monte Carlo simulation. The effectiveness of assistive detect and avoid is estimated for the same encounter scenarios. A prototype system for detect and avoid and a summary of results are presented. Preliminary results strongly suggest that assistive detect and avoid could greatly enhance the capabilities of flight crews to avoid traffic and remain well clear.

collision↗

Can Remotely Sensed Snow Disappearance Explain Seasonal Water Supply?

Understanding the relationship between remotely sensed snow disappearance and seasonal water supply may become vital in coming years to supplement limited ground based, in situ measurements of snow in a changing climate. For the period 2001–2019, we investigated the relationship between satellite derived Day of Snow Disappearance (DSD)—the date at which snow has completely disappeared—and the seasonal water supply, i.e., the April—July total streamflow volume, for 15 snow dominated basins across the western U.S. A Monte Carlo framework was applied, using linear regression models to evaluate the predictive skill—defined here as a model’s ability to accurately predict seasonal flow volumes—of varied predictors, including DSD and in situ snow water equivalent (SWE), across a range of spring forecast dates. In all basins there is a statistically significant relationship between mean DSD and seasonal water supply (p ≤ 0.05), with mean DSD explaining roughly half of the variance. Satellite-based model skill improves later in the forecast season, surpassing the skill of in-situ-based (SWE) models in skill in 10 of the 15 basins by the latest forecast date. We found little to no correlation between model error and basin characteristics such as elevation and the ratio of snow water equivalent to total precipitation. Despite a relatively short data record, this exploratory analysis shows promise for improving seasonal water supply prediction, in particular for snow dominated basins lacking in situ observations.

snow remote sensing↗

Minimum length Pb/SCIN detector for efficient cosmic ray identification

A study was made of the performance of a minimal length cosmic ray shower detector that would be light enough for space flight and would provide efficient identification of positions and protons. Cosmic ray positions are mainly produced in the decay chain of: Pion yields Muon yields Positron and they provide a measure of the matter density traversed by primary protons. Present positron flux measurements are consistent with the Leaky Box and Halo models for sources of cosmic rays. Abundant protons in the space environment are a significant source of background that would wash out the positron signal. Protons and positrons produced very distictive showers of particles when they enter matter; many studies have been published on their behavior on large calorimeter detectors. The challenge is to determine the minimal material necessary (minimal calorimeter depth) for positive particles identification. The primary instrument for the investigation is the Monte Carlo code GEANT, a library of programs from CERN that can be used to model experimental geometry, detector responses and particle interaction processes. The use of the Monte Carlo approach is crucial since statistical fluctuations in shower shape are significant. Studies conducted during the 1988 summer program showed that straightforward approaches to the problem achieved 85 to 90 percent correct identification, but left a residue of 10 to 15 percent misidentified particles. This percentage improved to a few percent when multiple shower-cut criteria were applied to the data. This summer, the same study was extended to employ several physical and statistical methods of identifying response of the calorimeter and the efficiency of the optimal shower cuts to off-normal incidence particle was determined.

Snyder, H. David↗

A Model-Based, Bayesian Solution for Characterization of Complex Damage Scenarios in Aerospace Composite Structures

Ultrasonic damage detection and characterization is commonly used in nondestructive evaluation (NDE) of aerospace composite components. In recent years there has been an increased development of guided wave based methods. In real materials and structures, these dispersive waves result in complicated behavior in the presence of complex damage scenarios. Model-based characterization methods utilize accurate three dimensional finite element models (FEMs) of guided wave interaction with realistic damage scenarios to aid in defect identification and classification. This work describes an inverse solution for realistic composite damage characterization by comparing the wavenumber-frequency spectra of experimental and simulated ultrasonic inspections. The composite laminate material properties are first verified through a Bayesian solution (Markov chain Monte Carlo), enabling uncertainty quantification surrounding the characterization. A study is undertaken to assess the efficacy of the proposed damage model and comparative metrics between the experimental and simulated output. The FEM is then parameterized with a damage model capable of describing the typical complex damage created by impact events in composites. The damage is characterized through a transdimensional Markov chain Monte Carlo solution, enabling a flexible damage model capable of adapting to the complex damage geometry investigated here. The posterior probability distributions of the individual delamination petals as well as the overall envelope of the damage site are determined.

H. Reed↗

Study of multi-dimensional radiative energy transfer in molecular gases

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical arrow band model with an exponential-tailed inverse intensity distribution. Consideration of spectral correlation results in some distinguishing features of the Monte Carlo formulations. Validation of the Monte Carlo formulations has been conducted by comparing results of this method with other solutions. Extension of a one-dimensional problem to a multi-dimensional problem requires some special treatments in the Monte Carlo analysis. Use of different assumptions results in different sets of Monte Carlo formulations. The nongray narrow band formulations provide the most accurate results.

Liu, Jiwen↗

The injection and acceleration of particles in oblique shocks - A unified Monte Carlo description

Important problems are discussed concerning the effects of shock geometry on the injection of thermal particles and the process of standard Fermi particle acceleration in general. In pursuit of a self-consistent model that simultaneously describes shock structure and particle acceleration in shocks of arbitrary obliquity, a Monte Carlo simulation developed for parallel shocks is generalized to address oblique geometries. Attention is given to initial results concerning the ways in which injection and acceleration efficiency varies with obliquity and Mach number.

Baring, Matthew G.↗

The OpenSE Cookbook: A Practical, Recipe Based Collection of Patterns, Procedures, and Best Practices for Executable Systems Engineering for the Thirty Meter Telescope

The OpenSE Cookbook is an open-sourced collection of patterns, procedures, and best practices targeted for systems engineers who seek guidance on applying model-based and executable systems engineering (MBSE) using SysML. Its content has emerged from the system level modeling effort on the European Framework Program 6 (FP6) and the Thirty Meter Telescope (TMT). The TMT MBSE approach applied the Executable Systems Engineering Method (ESEM) and the open-source Engineering Environment (OpenMBEE) to specify, analyze, and verify requirements of TMT’s Alignment and Phasing System (APS) and the Narrow Field Infrared Adaptive Optics System (NFIRAOS). In these applications, implicit dependencies are made explicit in a formal model through the use of ESEM, OpenMBEE, and SysML modeling constructs. The value proposition for applying this MBSE approach was to establish precise requirements and fine-grained traceability to system designs, and to verify key requirements beginning early in development. The integration of ESEM and the OpenMBEE tooling infrastructure (providing linked-data and web-operability) is a significant added value for the MBSE approach. The APS is responsible for the overall pre-adaptive optics wavefront quality, using starlight to measure wavefront errors and align the TMT optics. In the formally integrated and executable SysML model, simulations are performed to analyze the impact of changed requirements and verify specified constraints for various operational scenarios. The APS team used several modeling patterns to capture information such as the requirements, the operational scenarios, involved subsystems and their interaction points, the estimated or required time durations, and the mass and power consumption. Adaptive optics systems are designed to sense real-time atmospheric turbulence and correct the telescope’s optical beam to remove its effect. The system model for the adaptive optics operational modes was developed to capture sequence behaviors and operational scenarios to run Monte-Carlo simulations for verifying acquisition time, observing efficiency, and operational behavior requirements. The model is particularly useful for investigating the effect of parallelization, identifying interface issues, and re-ordering sequence acquisition tasks. A former version of the Cookbook (which is now updated to MBSE challenges, goals, and lessons learned) included modeling guidelines and conventions for all system aspects, hierarchy levels, and views, which were developed during for the Active Phasing Experiment (APE), an opto-mechatronical system technology demonstrator for the Extremely Large Telescope (ELT). The Cookbook utilizes the above mentioned system models as real-world case-studies to demonstrate and document the applications of the recipes, providing also instructional examples and addressing the available tooling support. The Cookbook is accompanied by a number of SysML models and aodel libraries which facilitate model authoring and maintenance. The Cookbook covers the different aspects of Systems Engineering such as management of Requirements, Design (behavior and structure), Interfaces, Interdisciplinary Integration, Analysis, Trade Studies, and Technical Resources. This paper presents the background, motivation, architecture, and highlights some key content of the Cookbook. For example, interface management, error budget management, requirements verification, Monte Carlo driven analysis, and timing analysis of operational scenarios. The paper discusses how the capabilities of OpenMBEE contributed significantly to the adoption of executable systems engineering.

Brower, Eric↗

Stress Rupture Life Reliability Measures for Composite Overwrapped Pressure Vessels

Composite Overwrapped Pressure Vessels (COPVs) are often used for storing pressurant gases onboard spacecraft. Kevlar (DuPont), glass, carbon and other more recent fibers have all been used as overwraps. Due to the fact that overwraps are subjected to sustained loads for an extended period during a mission, stress rupture failure is a major concern. It is therefore important to ascertain the reliability of these vessels by analysis, since the testing of each flight design cannot be completed on a practical time scale. The present paper examines specifically a Weibull statistics based stress rupture model and considers the various uncertainties associated with the model parameters. The paper also examines several reliability estimate measures that would be of use for the purpose of recertification and for qualifying flight worthiness of these vessels. Specifically, deterministic values for a point estimate, mean estimate and 90/95 percent confidence estimates of the reliability are all examined for a typical flight quality vessel under constant stress. The mean and the 90/95 percent confidence estimates are computed using Monte-Carlo simulation techniques by assuming distribution statistics of model parameters based also on simulation and on the available data, especially the sample sizes represented in the data. The data for the stress rupture model are obtained from the Lawrence Livermore National Laboratories (LLNL) stress rupture testing program, carried out for the past 35 years. Deterministic as well as probabilistic sensitivities are examined.

Murthy, Pappu L. N.↗

Application of an approach and landing system model to the space shuttle orbiter vehicle

An analytic model and procedure are described which can be used to estimate probability distributions of touchdown conditions, without utilizing a Monte Carlo simulation. Example applications are presented, including one-step and two-step flare strategies for the space shuttle orbiter. The computed probability distributions for the two-step flare are compared with the corresponding histograms from an independent Monte Carlo simulation of the same situation. The computed distributions are very good fits to the Monte-Carlo-generated histograms for pertinent longitudinal and lateral variables at three points (including touchdown) along the vehicle's final approach trajectory. The model, thus verified, has more inherent power for assessing the effects of system element changes than conventional Monte Carlo techniques.

Johnson, W. A.↗

A Monte Carlo study of parameters affecting computer simulations of crater saturation density

Computer models of cratered surfaces often use inputs of uncertain nature and importance. This work evaluates the sensitivity of the resulting crater-saturation estimates to the input parameters, principally applicable to the study of craters upward from 8 km diameter. In order of decreasing importance, crater saturation simulations are found to be sensitive to: (1) the dynamic range of crater diameters used, (2) the effectiveness of ejecta-blanket obliteration assumed, and (3) the number of points taken to describe the crater rim. The size of the largest crater in proportion to the size of the simulated surface has no effect on the results when the edges of the simulated surface are correctly treated and craters are not counted simply by integers. Craters should be counted by their fractions lying within the simulated area. A similar procedure is recommended when gathering crater size-density data from images.

Woronow, A.↗

Program Models A Laser Beam Focused In An Aerosol Spray

Monte Carlo analysis performed on packets of light. Program for Analysis of Laser Beam Focused Within Aerosol Spray (FLSPRY) developed for theoretical analysis of propagation of laser pulse optically focused within aerosol spray. Applied for example, to analyze laser ignition arrangement in which focused laser pulse used to ignite liquid aerosol fuel spray. Scattering and absorption of laser light by individual aerosol droplets evaluated by use of electromagnetic Lorenz-Mie theory. Written in FORTRAN 77 for both UNIX-based computers and DEC VAX-series computers. VAX version of program (LEW-16051). UNIX version (LEW-16065).

Barton, J. P.↗

Updating Mars-GRAM to Increase the Accuracy of Sensitivity Studies at Large Optical Depths

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an 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). During the Mars Science Laboratory (MSL) site selection process, it was discovered that Mars-GRAM, when used for sensitivity studies for MapYear=0 and large optical depth values such as tau=3, is less than realistic. From the surface to 80 km altitude, Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). MGCM results that were used for Mars-GRAM with MapYear set to 0 were from a MGCM run with a fixed value of tau=3 for the entire year at all locations. This has resulted in an imprecise atmospheric density at all altitudes. As a preliminary fix to this pressure-density problem, density factor values were determined for tau=0.3, 1 and 3 that will adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear 0 with Thermal Emission Spectrometer (TES) observations for MapYears 1 and 2 at comparable dust loading. Currently, these density factors are fixed values for all latitudes and Ls. Results will be presented from work being done to derive better multipliers by including variation with latitude and/or Ls by comparison of MapYear 0 output directly against TES limb data. The addition of these more precise density factors to Mars-GRAM 2005 Release 1.4 will improve the results of the sensitivity studies done for large optical depths.

Justh, Hiliary L.↗

Multiple Rings of Millimeter Dust Emission in the HD 15115 Debris Disk

We present observations of the HD 15115 debris disk from the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm that capture this intriguing system with the highest resolution (0."6 or 29 au) at millimeter wavelengths to date. This new ALMA image shows evidence for two rings in the disk separated by a cleared gap. By fitting models directly to the observed visibilities within a Markov Chain Monte Carlo framework, we are able to characterize the millimeter continuum emission and place robust constraints on the disk structure and geometry. In the best-fit model of a power-law disk with a Gaussian gap, the disk inner and outer edges are at 43.9 ± 5.8 au (0."89 ± 0."12) and 92.2 ± 2.4 au (1."88 ± 0."49), respectively, with a gap located at 58.9 ± 4.5 au (1."2 ± 0."10) with a fractional depth of 0.88 ± 0.10 and a width of 13.8 ± 5.6 au (0."28 ± 0."11). Because we do not see any evidence at millimeter wavelengths for the dramatic east–west asymmetry seen in scattered light, we conclude that this feature most likely results from a mechanism that only affects small grains. Using dynamical modeling and our constraints on the gap properties, we are able to estimate a mass for the possible planet sculpting the gap to be 0.16 ± 0.06 M(sub Jup).

MacGregor, Meredith A.↗