Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Analysis”

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 199 records · Page 11

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

Developing a Multilingual Auto-coding Interface Control for the MAVERIC-II Dynamics Simulator

Simulation model development in certain high-level languages such as Python, MATLAB, or Simulink are unparalleled by their convenience and rapid turnover time. However, legacy simulation engines often depend on more traditional languages such as FORTRAN or C/C++. The NASA Marshall Aerospace Vehicle Representation in C version II (MAVERIC-II) is a modular, legacy-derived computer program used for high-fidelity, 6 degree-of-freedom (6DOF) simulation for aerospace vehicle flights and analyses of guidance and control performance with built-in mathematical modeling of environmental effects such as wind, atmosphere, and gravity as well as dispersion capability for Monte Carlo analysis. MAVERIC-II is modular in the sense that each component software element of the simulation engine may be supplanted for a higher or lower fidelity version. The design flow of the development of these models is often performed in high-level languages as mentioned previously, which must then be translated into C or C++ code to be integrated into MAVERIC-II. Using principles of model-based design, we propose a unified method of auto-coding and interfacing between several languages and MAVERIC-II, which may be generalized further to any type of 6DOF simulation engine.

Mason Nixon↗

Fitting Strategies of Accretion Column Models and Application to the Broadband Spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT data-set of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies

Philipp Thalhammer↗

Inverse Estimation of Mars 2020 Entry Aeroheating Environments Using MEDLI2 Flight Data

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite collected heating and pressure data during entry into Mars' atmosphere of the Mars 2020 Perseverance rover. MEDLI2 included thermocouples, heat flux sensors, and pressure transducers on both the heatshield and the backshell. This paper covers the inverse estimation of heatshield and backshell surface heating based on the MEDLI2 Instrumented Sensor Plugs (MISPs), a network of thermocouples embedded in thermal protection system plugs across the aeroshell. Monte Carlo analysis was conducted to assess the sensitivity of the surface heat rate and temperature to uncertainties in thermocouple depth and material properties such as density, specific heat capacity, and thermal conductivity. Data from each MISP was also used to estimate the local time of transition from laminar to turbulent flow at each plug location.

Hannah S Alpert↗

A Modern Load Relief Guidance Scheme for Space Launch Vehicles

Launch vehicle load relief algorithms are concerned with realizing a reduction of transient bending moments near maximum dynamic pressure. Traditional approaches to load relief typically use inner-loop acceleration feedback to reduce the wind-induced angle of attack. When implemented in the inner loop, load relief bandwidth is necessarily limited by the achievable stability margins, and when acceleration feedback is employed, by the uncertainty associated with structural modes that couple with the body-mounted accelerometer. The structure of inner loop load relief increases the dimensionality of the flight control gain and filter optimization problem. Most importantly, classical load relief laws do not take advantage of high-rate and high-accuracy GPS-aided inertial velocity data that is readily available from modern strap down IMUs. In this paper, a novel load relief guidance scheme is described that uses direct angle-of-attack feedback in a clever mechanization. The steering commands are determined by examining the wind-perturbed dynamics of a launch vehicle with respect to a gravity turn ascent trajectory. An angle of attack estimate is derived from GPS-aided inertial data and pre-launch range wind measurements, and it is shown that a reduction worst-case rigid-body loads can be realized without requiring air data. The algorithm also includes a high-rate navigation data preprocessing scheme that operates directly on the IMU delta-theta and delta-velocity measurements in order to produce a filtered acceleration estimate at the vehicle center of mass. The outer-loop guidance scheme simplifies the design process for the classical inner-loop autopilot. Algorithm performance is demonstrated using Monte Carlo analysis of a representative liquid booster in a production high fidelity launch vehicle simulation.

NESC↗

Quantifying Uncertainty of Landing and Takeoff Noise for Commercial Supersonic Aircraft

Of the many challenges faced by manufacturers attempting to offer supersonic travel to the public, the uncertainty in predicting the noise of these aircraft in airport operations has an immediate impact. No noise regulation exists in FAA or ICAO for certifying such aircraft, as these organization require solid data, usually from existing aircraft. Manufacturers are taking large risks to design a vehicle not knowing whether it will be allowed to fly. A partial solution to this conundrum is to use physics-based simulations to provide the “data” used to calibrate system-level prediction methods, carefully documenting the uncertainty of the method for application to supersonic aircraft. A close look at the accuracy of empirical prediction methods points to areas where improvements need to be made if noise studies of supersonic aircraft are to be useful. As NASA embarks on a focused research program to improve predictions of noise from the noise-dominant propulsion noise of commercial supersonic aircraft, this paper documents the work done to baseline the uncertainties found in today’s noise prediction methods. A relatively simplistic method was developed, summarizing the error of the empirical methods on a component basis and following their impact on the total aircraft during landing and takeoff operations using Monte Carlo analysis. By this method it is found that current empirical noise prediction methods have an uncertainty of 1.5 EPNdB cumulative for propulsion noise of a representative conventional subsonic passenger aircraft. When applied to likely near-term supersonic commercial aircraft, the uncertainty is 7.6 EPNdB cumulative, a difference that must be reduced if the prediction methods are to guide decision makers.

Airport noise↗

Variance Decomposition of MEDLI2 Reconstructed Heating Using Neural Networks

The Mars Entry, Descent, and Landing Instrumentation (MEDLI2) sensor suite collected data during entry of the Mars 2020 Perseverance rover into Mars’ atmosphere. An inverse estimation of the backshell and heatshield surface aeroheating was performed, using the data from the MEDLI2 Instrumented Sensor Plugs, a network of thermocouples embedded within the thermal protection system across the aeroshell. Monte Carlo analysis was conducted to assess the sensitivity of the surface heat rate, temperature, and heat load to uncertainties in thermocouple depth and material properties. In this paper, a variance decomposition method using Sobol indices was employed to understand the relative contributions of each uncertainty parameter. Performing this analysis using results from the inverse analysis tool FIAT_Opt was found to require incredibly high computation time, and thus machine learning models were trained and evaluated as a surrogate model for FIAT_Opt. This paper demonstrates that machine learning models can be an efficient, accurate alternative to state-of-the-art inverse analysis tools like FIAT_Opt, especially for computationally-expensive processes. Using these models, the sensitivity analysis showed that uncertainties in heat capacity and thermal conductivity were the main drivers for the overall uncertainty in peak reconstructed heating and heat load.

H S Alpert↗

Comparison of Dynamic Stability Testing Techniques with Magnetic Suspension Wind Tunnel

Dynamic stability testing techniques currently utilized at NASA Langley Research Center (LaRC) are conducted in multiple facilities and consists of free flight, forced oscillation, and free-to-oscillate tests. The MIT/NASA/ODU Magnetic Suspension and Balance System (MSBS) has been recommissioned to explore its utility as an additional facility to expand the dynamic stability test capabilities currently available at NASA LaRC. Simulations were created to replicate each current test facility and method as closely as possible. Data collected from the simulated environments was corrupted with replicated noise sources of the different testing environments and then compared to real data collected during tests when such data was available. The corrupted data was then passed through data reduction and System Identification (SID) to estimate the accuracy of the results with the known aerodynamic model that was utilized within the simulation to generate the original data. Magnitudes of noise were varied utilizing Monte Carlo analysis to perform sensitivity analysis of each noise source on the extracted dynamic stability coefficients. Some preliminary results will be presented.

Otoniel Ramirez↗

Comparison of Dynamic Stability Testing Techniques with Magnetic Suspension Wind Tunnel

Dynamic stability testing techniques currently utilized at NASA Langley Research Center (LaRC) are conducted in multiple facilities and consists of free flight, forced oscillation, and free-to-oscillate tests. The MIT/NASA/ODU Magnetic Suspension and Balance System (MSBS) has been recommissioned to explore its utility as an additional facility to expand the dynamic stability test capabilities currently available at NASA LaRC. Simulations were created to replicate each current test facility and method as closely as possible. Data collected from the simulated environments was corrupted with replicated noise sources of the different testing environments and then compared to real data collected during tests when such data was available. The corrupted data was then passed through data reduction and System Identification (SID) to estimate the accuracy of the results with the known aerodynamic model that was utilized within the simulation to generate the original data. Magnitudes of noise were varied utilizing Monte Carlo analysis to perform sensitivity analysis of each noise source on the extracted dynamic stability coefficients. Some preliminary results will be presented.

Otoniel Ramirez↗

Generalized Linear Targeting For Cislunar Flight

An important element of Artemis and NASA’s campaign to explore the Moon is the autonomous onboard two-level targeter (TLT) used during all cislunar flight phases. The function of the TLT is to autonomously recompute the burn targets for the upcoming burn (or multiple burns) in response to navigation and vehicle dispersion providing a solution that meets all of the trajectory constraints. Although the TLT has been utilized previously as a ground-based planning tool, and flown onboard during the Artemis I mission, it’s complexity and iterative nature make is difficult to incorporate into and support rapid analyses such as robust optimal trajectory design applications where speed is essential. In this paper, a set of generalized linear targeting algorithms that mimics many of the properties of the TLT is derived. The generalized algorithms can handle single or multiple impulsive maneuvers, with multiple constraints at multiple fixed or variable times. A linear targeting algorithm for finite burn maneuvers is also derived. The generalized linear targeting algorithms are exceptionally fast and easy to implement in Monte Carlo analysis, linear covariance (LinCov) analysis, and robust optimal trajectory design. Several cislunar flight examples are provided.

Linear Covariance Analysis↗

Machine Learning Algorithms for Alignment Verification of the Roman Space Telescope

The Nancy Grace Roman Telescope is a NASA observatory designed to unravel the secrets of dark energy and dark matter, search for and image exoplanets, and explore many topics in infrared optics. Scheduled to launch no earlier than October 2026, this 2.4 meter aperture telescope has a field of view 100 times greater than the Hubble Space Telescope. The mission is currently in its construction phase, where the telescope and its two instruments will soon be aligned together to ensure proper pupil matching. To help verify this alignment, multiple point sources above the entrance pupil of the telescope will illuminate the optical path through the telescope-instrument system, and shadows of various obstructions in the system will be analyzed using machine learning algorithms to determine the pupil matching error. This presentation discusses the test approach and the machine learning algorithms employed, as well as our uncertainty predictions based on a modeled Monte-Carlo analysis of the test.

Telescope↗

Shape Approximation in Camera Fields of View for Parachute Visibility Applications

The ability to simulate camera line of sight and field of view is a powerful capability for flight mechanics simulations. This paper improves upon an existing camera visibility implementation that leverages geometry to estimate parachute visibility from a singular point by introducing a shape discretization method. The new method approximates a shape as a set of analytical points relative to a target marker within any number of camera fields of view. The method is developed in such a way to be generalized to a myriad of spacecraft applications. This effort focuses on an implementation to simulate parachute canopy visibility within the combined field of view of three cameras but additional applications are discussed. Monte Carlo analysis is used to assess the parachute canopy visibility example and demonstrate its improved performance over the previous method.

Evan Roelke↗

Onboard Autonomous Trajectory Planning for Mars Power Descent

In recent years, there has been an increasing interest in space-qualified processors such as multi-core central processing units and graphics processing units that can withstand the adverse effects of space radiation. These processors can allow parallel programming to perform tasks that typically demand high computational power. One can study guidance schemes that can take advantage of these currently developing processors and provide more robust guidance. Software for Multi-model Autonomous Real-time Trajectories (SMART) guidance can identify robust trajectories by running an onboard Monte Carlo analysis. SMART guidance can take advantage of knowledge updates obtained from the onboard sensors, allowing it to consider the off-nominal cases that it would not typically encounter during the offline trajectory analysis. This work uses the SMART guidance for the powered divert at Mars simulation in Program to Optimize and Simulated Trajectories- II.

Pardha Sai Chadalavada↗

Shape Approximation in Camera Fields of View for Parachute Visibility Applications

The ability to simulate camera line of sight and field of view is a powerful capability for flight mechanics simulations. This paper improves upon an existing camera visibility implementation that leverages geometry to estimate parachute visibility from a singular point by introducing a shape discretization method. The new method approximates a shape as a set of analytical points relative to a target marker within any number of camera fields of view. The method is developed in such a way to be generalized to a myriad of spacecraft applications. This effort focuses on an implementation to simulate parachute canopy visibility within the combined field of view of three cameras but additional applications are discussed. Monte Carlo analysis is used to assess the parachute canopy visibility example and demonstrate its improved performance over the previous method.

Evan Roelke↗

Onboard Autonomous Trajectory Planning for Mars Power Descent

In recent years, there has been an increasing interest in space-qualified processors such as multi-core central processing units and graphics processing units that can withstand the adverse effects of space radiation. These processors can allow parallel programming to perform tasks that typically demand high computational power. One can study guidance schemes that can take advantage of these currently developing processors and provide more robust guidance. Software for Multi-model Autonomous Real-time Trajectories (SMART) guidance can identify robust trajectories by running an onboard Monte Carlo analysis. SMART guidance can take advantage of knowledge updates obtained from the onboard sensors, allowing it to consider the off-nominal cases that it would not typically encounter during the offline trajectory analysis. This work uses the SMART guidance for the powered divert at Mars simulation in Program to Optimize and Simulated Trajectories- II.

Autonomous Planning↗

Validation of Linear Covariance Techniques for Mars Entry, Descent, and Landing Guidance and Navigation Performance Analysis

Current Monte Carlo-based uncertainty analysis methods may require significant computational resources to evaluate the performance of a closed-loop guidance, navigation, and control system. An attractive alternative, particularly during the preliminary and conceptual design phase, is to use linear covariance analysis, which can provide the same statistical information as Monte Carlo methods at a fraction of the computational load. Linear covariance has already been demonstrated in various spaceflight regimes, but only recently has it been applied to atmospheric flight. In this study, a 6-degree-of-freedom formulation of both a linear covariance and Monte Carlo analysis tools are utilized for a Mars entry, descent, and landing scenario which capture both atmospheric and powered flight phases to support precision landing. Comparison of the performance results shows close agreement between the linear covariance and traditional Monte Carlo methods when incorporating an assortment of guidance algorithms and processing a variety of inertial and relative sensor measurements onboard the lander's navigation filter.

James W. Williams↗

Monte Carlo models and analysis of galactic disk gamma-ray burst distributions

Gamma-ray bursts are transient astronomical phenomena which have no quiescent counterparts in any region of the electromagnetic spectrum. Although temporal and spectral properties indicate that these events are likely energetic, their unknown spatial distribution complicates astrophysical interpretation. Monte Carlo samples of gamma-ray burst sources are created which belong to Galactic disk populations. Spatial analysis techniques are used to compare these samples to the observed distribution. From this, both quantitative and qualitative conclusions are drawn concerning allowed luminosity and spatial distributions of the actual sample. Although the Burst and Transient Source Experiment (BATSE) experiment on Gamma Ray Observatory (GRO) will significantly improve knowledge of the gamma-ray burst source spatial characteristics within only a few months of launch, the analysis techniques described herein will not be superceded. Rather, they may be used with BATSE results to obtain detailed information about both the luminosity and spatial distributions of the sources.

Hakkila, Jon↗