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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Enhanced Stationkeeping Maneuver Control Technique for Delta-V Cost Reduction in the Korea Pathfinder Lunar Orbiter

Final document is attached. This paper proposes an enhanced control technique for stationkeeping maneuvers to reduce delta-v costs for the Korea Pathfinder Lunar Orbiter (KPLO). A scheduled circularization control technique exploits patterns in the evolution of the line of apsides and eccentricity to achieve a significant reduction in stationkeeping delta-v costs based on spacecraft requirements. The technique is compared against previous algorithms implemented for maneuver operations of the Lunar Prospector and Lunar Reconnaissance Orbiter (LRO) missions in the USA and KAGUYA in Japan. Through Monte Carlo analysis, the efficacy and robustness of the proposed method are verified, and the technique is shown to meet the operational requirements of KPLO.

Davis, Diane C.↗

NASA NTRS2019

Transfer Trajectory Options for Servicing Sun-Earth-Moon Libration Point Missions

Future missions to the Sun-Earth Libration L1 and L2 regions will require scheduled servicing to maintain hardware and replenish consumables. While there have been statements made by various NASA programs regarding servicing of vehicles at these locations or in Cis-lunar space, a practical transfer study has not been extensively investigated in an operational fashion to determine the impacts of navigation and maneuver errors. This investigation uses dynamical systems and operational models to design transfer trajectories between the Sun-Earth Libration region (QuasiHalo orbit) and the Earth-Moon vicinity (Distant Retrograde Orbit, QuasiHalo Orbit, Halo Orbit, and Near Rectilinear Halo Orbit). We address the total delta V cost of transfers between each pair of locations using a Monte Carlo analysis.

Earth↗

NASA NTRS2019

Modeling and Simulation Techniques for the NASA SLS Service Module Panel Separation Event; from Loosely-Coupled Euler to Fully-Coupled 6-DOF, Time-Accurate, Navier-Stokes Methodologies

An aerodynamic database has been generated for use by the Orion Multi-Purpose Crew Vehicle (MPCV) Program to analyze Service Module (SM) panel jettison from the NASA SLS vehicle. The database is a combination of CFD data for the panel aerodynamic coefficients, and MATLAB code written to query the CFD data. The Cart3D inviscid CFD flow solver was used to generate the panel aerodynamic coefficients for static panel orientations and free stream conditions that can occur during the jettison event. The MATLAB code performs the multivariate interpolation to obtain aerodynamic coefficients. The MATLAB code uses input for SM panel parameters and returns the SM panel aerodynamic force and moment coefficients for use with a Six-Degree-of-Freedom (6-DOF) motion solver to model the jettison event. This paper examines the accuracy of the sequential-static database approach by modeling the panel jettison event with a fully-coupled, time-dependent, viscous, moving-body CFD simulation. The fully-coupled simulation is obtained using the Loci/Chem unstructured Navier-Stokes CFD solver. The results show that the fully-coupled approach agrees well with the loosely-coupled database/6-DOF approach, indicating that unsteady effects are minimal for the panel jettison event. These results suggest that the database/6-DOF approach is sufficient. In addition, this paper presents the development of an uncertainty model for use in Monte Carlo analysis of the panel jettison event. Here viscous CFD simulations are obtained with Loci/Chem and compared to the inviscid CFD forces and moments. An uncertainty model based on model-form error and numerical error is presented.

Hall, Leslie H.↗

NASA NTRS2018

A Component-Level Model of Automatic Dependent Surveillance - Broadcast (ADS-B)

Automatic Dependent Surveillance – Broadcast (ADS-B) is being employed in numerous peer-to-peer initiatives attempting to expand the capacity of the National Airspace System (NAS) or enable mixed operations of manned and unmanned vehicles. Safety assessments of these initiatives rely, in part, on modeling the accuracy of ADS-B in reporting the position and direction of an ownship and surrounding traffic. Frequently, these initiatives utilize a position uncertainty model that applies a reported ADS-B estimation position uncertainty (EPU) value to a Rayleigh distribution and uses a Gauss-Markov random walk to add error to the ADS-B output of a vehicle. This model of ADS-B state error is easy to implement and apply to numerous problems. However, it has a couple of draw-backs. First, the ADS-B state errors are equally probable in all directions. This is a good assumption in situations where aircraft maneuvering is not constrained. However, in situations where the aircraft maneuvering is constrained such as landing, the error distribution is likely to exhibit directionality and the non-directional model may skew results especially when assessing very low probabilities (e.g., 10(exp -9)) of catastrophic encounters. Second, the model does not account for processing latency in the receiving aircraft. NASA Langley Research Center (LaRC) recently examined the feasibility of decreasing the spacing of aircraft on parallel approaches to runways separated by as little as 700 feet. For Monte-Carlo analysis using a high-fidelity simulation of a large transport, LaRC started with a Gauss-Markov model of ADS-B error but then developed a component level model of ADS-B error to increase the fidelity of results.

Madden, Michael M.↗

NASA NTRS2019

Implementation Recommendations and Usage Boundaries for the Two-Dimensional Probability of Collision Calculation

The two-dimensional (2D) probability of collision (𝑃𝑐) estimation method relies on several assumptions that must be satisfied for accurate results. Monte Carlo analysis of ~44,000 conjunctions indicates that 2D-𝑃(sub 𝑐) pro-vides accurate estimates for most typical conjunctions, but occasionally underestimates 𝑃(sub 𝑐) significantly, indicating an assumption violation. A test to detect large-amplitude underestimation inaccuracies can be based on how much “offset-from-TCA” 2D-𝑃(sub 𝑐) values vary during a well-defined time interval bracketing closest approach. The test successfully detects all large-amplitude 2D-𝑃(sub 𝑐) underestimations found to date, but with a high false-alarm rate. The analysis also provides implementation recommendations and usage boundaries for the 2D-𝑃(sub 𝑐) method.

Hall, Doyle T.↗

NASA NTRS2019

Implementation Recommendations and Usage Boundaries for the Two-Dimensional Probability of Collision Calculation

Hall, Doyle T.↗

NASA NTRS2019

Transfer Trajectory Options for Servicing Sun-Earth-Moon Libration Point Missions

Future missions to the Sun-Earth Libration L1 and L2 regions will require scheduled servicing to maintain hardware and replenish consumables. While there have been statements made by various NASA programs regarding servicing of vehicles at these locations or in Cis-lunar space, a practical transfer study has not been extensively investigated in an operational fashion to determine the impacts of navigation and maneuver errors. This investigation uses dynamical systems and operational models to design transfer trajectories between the Sun-Earth Libration region (QuasiHalo orbit) and the Earth-Moon vicinity (Distant Retrograde Orbit, QuasiHalo Orbit, Halo Orbit, and Near Rectilinear Halo Orbit). We address the total ΔV cost of transfers and operational considerations between each pair of locations using a Monte Carlo analysis.

Folta, David C.↗

NASA NTRS2019

Optimization of the Lucy Interplanetary Trajectory via Two-Point Direct Shooting

Lucy is NASA’s next Discovery-class mission and will explore the Trojan asteroids in the Sun-Jupiter L4 and L5 regions. This paper details the design of Lucy’s interplanetary trajectory using a two-point direct shooting transcription, nonlinear programming, and monotonic basin hopping. These techniques are implemented in the Evolutionary Mission Trajectory Generator (EMTG), a trajectory optimization tool developed at NASA Goddard Space Flight Center. We present applications to the baseline trajectory design, Monte Carlo analysis, and operations.

Lucy↗

NASA NTRS2019

Star Tracker Accuracy Improvement and Optimization for Attitude Measurement in Three-Axis

High precision attitude measurement systems obviate the need for the beacon from the receiver making it possible for the spacecraft to beam a laser communications signal to a ground station without the ground station advertising its location. The research presented targets new detection and estimation methods to improve the accuracy in locating stars on a focal plane detector, and an understanding of the effects of changes in the optics design parameters and aberration, including defocus, on the navigation solution itself. This understanding can lead to an optimization of the attitude solution with respect to those optics realm parameter changes. The methodology discussed includes the development of a model of a current star tracker system. Using this model, multiple algorithms are implemented, including a multi-hypothesis method (MHT), to detect and estimate the position of the stars on the focal plane detector. It will be shown that using the MHT for detection and estimation, a greater accuracy can be found for each star estimation from more traditional detection and estimation algorithms. The approach then uses the model to develop statistics of the star tracker and the attitude estimation outputs to understand the accuracy, or variance, of the system's attitude solution. This solution is repeated for a range of defocus aberration, and a lower limit to the variance of the attitude solution is shown. A Cramer Rao lower bound solution is derived for the star tracker system and the results are compared to the Monte Carlo analysis from the model and shown to correlate very well. The approach uses a star image not as a Gaussian spot on the focal plane as done in previous work, and use of an image that includes the effects of aberrations of the optic system, and the effects of under-sampling and noise from the focal plane detector as well. Analysis includes exploring a star tracker's accuracy improvement through the combination of focus error and under-sampling effects alone, possibly contradicting conventional wisdom and approaches.

Lichter, Michael J.↗

NASA NTRS2020

Comparison of Aerocapture Performance Using Bank Control and Direct Force Control with Two Human-Scale Vehicles at Mars

Recent studies of human-scale missions to Mars have included a wide trade space of vehicle configurations and control schemes. Some configurations fly at a low angle of attack with a low L/D,while others fly at a high angle of attack with a mid L/D. Some use bank angle control, while others use direct force control, where the angle of attack and sideslip angle are independently modulated. This paper compares the performance of three vehicle configurations: a low-L/D vehicle using direct force control, a low-L/D vehicle using bank control, and a mid-L/D vehicle using bank control. The reference mission is aerocapture at Mars into a highly elliptical, 1-sol orbit. The trajectories are integrated in three degrees of freedom. All three cases utilize numeric predictor-corrector guidances and emulate control system responses with rate and acceleration limits. The configurations are compared using a Monte Carlo analysis. The robustness of each configuration to increased dispersions is also compared.

Matz, Daniel A.↗

NASA NTRS2019

String Midpoint Monitoring for Li-Ion Batteries

The paper introduces the midpoint monitoring method for detecting cell voltage imbalance in a series string of 8 electrochemical cells. Monte Carlo analysis is performed for normal cell End of Discharge open and closed circuit voltage behavior to assess the extent to which derived guard band voltage limits will detect over discharged cells, provide false positive indication, or fail to detect over discharged cells. Forward work is proposed to refine the battery monitoring method for operational use in 28V lithium ion batteries.

Russell, S.↗

NASA NTRS2018

Detection of Photospheric Features in the Near-Infrared Spectrum of a Class 0 Protostar

We present a near-infrared K-band R ≃ 1500 Keck spectrum of S68N, a Class 0 protostar in the Serpens molecular cloud. The spectrum shows a very red continuum, CO absorption bands, weak or non-existent atomic metal absorptions, and H2 emission lines. The near-IR H2 emission is consistent with excitation in shocks or by X-rays but not by UV radiation. We model the absorption component as a stellar photosphere plus circumstellar continuum emission with wavelength-dependent extinction. A Markov Chain Monte Carlo analysis shows that the most likely model parameters are consistent with a low-temperature, low-gravity photosphere with signifi cant extinction and no more than modest continuum veiling. Its T(eff) ≃ 3260 K effective temperature is similar to that of older, more evolved pre-main-sequence stars, but its surface gravity log g ≃ 2:4 cm/sq.s is approximately 1 dex lower. This implies that the radius of this protostar is a factor of ~ 3 larger than that of 106 yr old T Tauri stars. Its low veiling is consistent with a circumstellar disk having intrinsic near-IR emission that is less than or equal to that of more evolved Class I protostars. Along with the high extinction, this suggests that most of the circumstellar material is in a cold envelope, as expected for a Class 0 protostar. This is the first known detection and analysis of a Class 0 protostar absorption spectrum.

Thomas P Greene↗

NASA NTRS2021

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗

NASA NTRSDate not supplied

Integrated GNC Modifications and Performance Assessments for the Mars Ascent Vehicle, Spin-Stabilized Upper Stage Configuration

The objective of the Mars Ascent Vehicle (MAV), an element of the Mars Sample Return (MSR) campaign, is to successfully transport Martian surface samples from the Martian surface to a stable Low-Martian Orbit. Complete autonomy is required throughout ascent, and orbital insertion is constrained by tight dispersion boundaries. Given mass requirements to merely transport the MAV onto the Martian surface, reducing the total mass of the MAV reduces stress on all other campaign elements. An unguided, spin-stabilized second stage MAV allows for lower overall vehicle mass, at the cost of reduced GNC capability. This paper presents GNC modifications developed and employed to improve the performance of an unguided, spin-stabilized second stage MAV configuration. These updates include a predictive calculation scheme to solve the spin-stabilized attitude necessary for successful orbital insertion, and a spin-stabilized control law for a preliminary low-rate first stage spin-up. Trades were conducted focusing on first and second stage low-rate and high-rate spin rates (and other factors) in MSFC's MANTIS toolkit through 6-degree-of-freedom simulation and Monte Carlo analysis, and are presented here.

Dane Erickson↗

NASA NTRS2020

OSA Spotlight Summary: Designing Miniature X-ray Optics for the SmallSat Lunar Science Mission Concept CubeX

The scientific demand for extended field X-Ray remote sensing for mapping of planetary bodies has produced the CubeSat X-ray Telescope (CubeX) lunar mission concept. Orbiting the Moon at 6,000 km for 1 year sometime during the mid-decade solar maximum, the primary objectives of CubeX are to map lunar surface elemental abundances and to assess feasibility of x-ray pulsar timing navigation. CubeX’s Miniature X-ray Optics (MiXO) design builds on heritage of existing X-ray imaging designs using new electroformed NiCo alloy replication (ENR) to create shells of 1 mm thickness. MiXO will provide 1 arcmin resolution X-Ray imaging for an extended instantaneous field of view, at various energy range sensitivities. Based on the classic Wolter-Type I design (using a paraboloid/hyperboloid grazing mirror pair), mirror parameters were expanded to allow both shell length and focal-plane offsets to vary independently. The design was optimized with a merit function that factored together four sub-merits: PSF sharpness (measured as half-power diameter), consistency across the field, minimized scattering, and maximized effective area. Monte Carlo analysis was used to predict final performance, including surface scatter effects. The demonstrated configurations include a compact 34 shell design with a diameter of 105 mm more efficient at higher energies, and a larger 130 mm/24 shell version with a larger overall effective area at lower energies. The compact design carries more expensive fabrication and assembly cost (with more shells) as well as a more challenging alignment (with denser packing). Both have overall resolution <1’ up to 26’ off-axis and <1.5’ up to 33’ off-axis, whereas a conventional design’s image quality degrades quickly beyond 15’ off-axis. This opens up new possibilities for future SmallSat missions that could use two different configurations for wide energy ranges or midsized instruments using multiple identical modules for larger collecting area.

x-ray↗

NASA NTRSDate not supplied

SLS Integrated Modal Test Uncertainty Quantification using the Hybrid Parametric Variation Method

Uncertainty in structural loading during launch is a significant concern in the development of spacecraft and launch vehicles. Small variations in launch vehicle and payload mode shapes and their interaction can result in significant variation in system loads. In many cases involving large aerospace systems it is difficult, not economical, or impossible to perform a system modal test. However, it is still vital to obtain test results that can be compared with analytical predictions to validate models. Instead, the “Building Block Approach” is used in which system components are tested individually. Component models are correlated and updated to agree as best they can with test results. The Space Launch System consists of a number of components that are assembled into a launch vehicle. Finite element models of the components are developed, reduced to Hurty/Craig-Bampton models and assembled to represent different phases of flight. The only opportunity to obtain modal test data from an assembled Space Launch System will be during the Integrated Modal Test. There is always uncertainty in every model, which flows into uncertainty in predicted system results. Uncertainty Quantification is used to determine statistical bounds on prediction accuracy based on model uncertainty. For the Space Launch System, model uncertainty is at the Hurty/Craig-Bampton component level. Uncertainty in the Hurty/Craig-Bampton components is quantified using the hybrid parametric variation approach that combines parametric and nonparametric uncertainty. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. This type of uncertainty cannot be represented by variations infinite element model input parameters and thus cannot be included in a parametric approach. However, model-form uncertainty can be modeled using a nonparametric approach based on random matrix theory. The hybrid parametric variation method requires the selection of dispersion values for the Hurty/Craig-Bampton fixed-interface eigenvalues, and the Hurty/Craig-Bampton stiffness matrices. Component test/analysis frequency error is used to identify the fixed-interface eigenvalue dispersions, while test/analysis cross-orthogonality is used to identify stiffness dispersion values. The hybrid parametric variation uncertainty quantification approach is applied to the Space Launch System Integrated Modal Test configuration. Monte Carlo analysis is performed, and statistics are determined for modal correlation metrics, frequency response from Integrated Modal Test shakers to selected accelerometers, as well as other metrics for determining how well target modes are excited and identified. If the predicted uncertainty envelopes future Integrated Modal Test results, then there will be increased confidence in the utility of the component-based hybrid parametric variation uncertainty quantification approach.

Uncertainty Quantification↗

NASA NTRSDate not supplied

Developing a Multi-Lingual Autocoding Interface 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. We propose a unified method of autocoding and interfacing between several languages and MAVERIC-II, which may be generalized further to any type of 6dof simulation engine.

Mason Nixon↗

NASA NTRS2020

Detection of Interstellar HC4NC and an Investigation of Isocyanopolyyne Chemistry under TMC-1 Conditions

Abstract We report an astronomical detection of HC4NC for the ﬁrst time in the interstellar medium with the Green Bank Telescope toward the TMC-1 molecular cloud with a minimum signiﬁcance of 10.5σ. The total column density and excitation temperature of HC4NC are determined to be 3.29+8.60 ´ 1011 cm−2 and 6.7+0.3 K, respectively, using -1.20 -0.3 Markov chain Monte Carlo analysis. In addition to HC4NC, HCCNC is distinctly detected whereas no clear detection of HC6NC is made. We propose that the dissociative recombination of the protonated cyanopolyyne, HC5NH+, and the protonated isocyanopolyyne, HC4NCH+, are the main formation mechanisms for HC4NC while its destruction is dominated by reactions with simple ions and atomic carbon. With the proposed chemical networks, the observed abundances of HC4NC and HCCNC are reproduced satisfactorily.

Astrochemistry↗