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 379 records · Page 21

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low earth orbit for scientific, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with 8 Teledyne dosimeters over varying shielding effectiveness. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 –AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat 3 unit (U) models. Using NOVICE model radiation analysis coding, the shielding effectiveness's, based on a generic CubeSat 3U structure with 4 electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cmto4.44 cm (0.550 g/cm2–12.0 g/cm2) thick aluminum. For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID)increases by nearly a magnitude for thin-walled aluminum 0.550g/cm2-0.686 g/cm2(0.204 cm –0.254 cm) typical CubeSat structures. The shielding effectiveness by NOVICE Sigma estimates, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from Earth's magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast, to a thin walled generic 3U CubeSat, Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm2from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the spacecraft’s mission value versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned paper is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low Earth orbit for science, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low Earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with eight Teledyne µDosimeters over varying shielding effectivenesses. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 – AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat three unit (3U) models. Using NOVICE model radiation analysis coding, the shielding effectivenesses, based on a generic CubeSat 3U structure with four electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cm to 4.44 cm (0.550 g/cm 2 to 12.0 g/cm 2 ). For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID) increases by nearly a magnitude for thin-walled aluminum 0.550-g/cm 2 - 0.686-g/cm 2 (0.204-cm – 0.254-cm) typical CubeSat structures. The shielding effectiveness estimated by NOVICE Sigma, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one-year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from the Earth’s magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast to a thin-walled generic 3U CubeSat, the Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm 2 from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the mission value of the spacecraft versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗

Uncertainty Analysis of the NASA Glenn 8x6 Supersonic Wind Tunnel

This paper presents methods and results of a detailed measurement uncertainty analysis that was performed for the 8- by 6-foot Supersonic Wind Tunnel located at the NASA Glenn Research Center. The statistical methods and engineering judgments used to estimate elemental uncertainties are described. The Monte Carlo method of propagating uncertainty was selected to determine the uncertainty of calculated variables of interest. A detailed description of the Monte Carlo method as applied for this analysis is provided. Detailed uncertainty results for the uncertainty in average free stream Mach number as well as other variables of interest are provided. All results are presented as random (variation in observed values about a true value), systematic (potential offset between observed and true value), and total (random and systematic combined) uncertainty. The largest sources contributing to uncertainty are determined and potential improvement opportunities for the facility are investigated.

Error↗

Analysis of large solid propellant rocket engine exhaust plumes using the direct simulation Monte Carlo method

A new solution procedure has been developed to analyze the flowfield properties in the vicinity of the Inertial Upper Stage/Spacecraft during the 1st stage (SRMI) burn. Continuum methods are used to compute the nozzle flow and the exhaust plume flowfield as far as the boundary where the breakdown of translational equilibrium leaves these methods invalid. The Direct Simulation Monte Carlo (DSMC) method is applied everywhere beyond this breakdown boundary. The flowfield distributions of density, velocity, temperature, relative abundance, surface flux density, and pressure are discussed for each species for 2 sets of boundary conditions: vacuum and freestream. The interaction of the exhaust plume and the freestream with the spacecraft and the 2-stream direct interaction are discussed. The results show that the low density, high velocity, counter flowing free-stream substantially modifies the flowfield properties and the flux density incident on the spacecraft. A freestream bow shock is observed in the data, located forward of the high density region of the exhaust plume into which the freestream gas does not penetrate. The total flux density incident on the spacecraft, integrated over the SRM1 burn interval is estimated to be of the order of 10 to the 22nd per sq m (about 1000 atomic layers).

Hueser, J. E.↗

Using Importance Sampling Monte Carlo to Analyze Aircraft Takeoff and Landing

Aircraft have achieved high levels of reliability, so the probability of undesirable dynamics is very low. The low probability of a bad takeoff or landing challenges nondeterministic methods to quantify the uncertainty of these improbable events. In the present work, a Monte Carlo method is defined to more effectively quantify these unlikely events that is suitable for uncertainty analysis of a mature design with many uncertainty parameters. The proposed Monte Carlo method has been applied to the takeoff and landing of the X-59 Quiet SuperSonic Technology (QueSST). Takeoff and landing simulation profiles were defined using a combination of industry standards and piloted simulation experience. The proposed method is shown to be able to quantify events with significantly fewer samples than would be required with conventional Monte Carlo. These results improve understanding of the sensitivity of the takeoff and landing dynamics to inform further studies and test planning.

Jeffrey Ouellette↗

Using Importance Sampling Monte Carlo to Analyze Aircraft Takeoff and Landing

Aircraft have achieved high levels of reliability, so the probability of undesirable dynamics is very low. The low probability of a bad takeoff or landing challenges nondeterministic methods to quantify the uncertainty of these improbable events. In the present work, a Monte Carlo method is defined to more effectively quantify these unlikely events that is suitable for uncertainty analysis of a mature design with many uncertainty parameters. The proposed Monte Carlo method has been applied to the takeoff and landing of the X-59 Quiet SuperSonic Technology (QueSST). Takeoff and landing simulation profiles were defined using a combination of industry standards and piloted simulation experience. The proposed method is shown to be able to quantify events with significantly fewer samples than would be required with conventional Monte Carlo. These results improve understanding of the sensitivity of the takeoff and landing dynamics to inform further studies and test planning.

Jeffrey Ouellette↗

Three dimensional particle simulation of high altitude rocket plumes

The interaction of two nozzles exhausting into vacuum generates a complex three-dimensional shock structure. The shock structure and resulting plume flowfield is characterized by the nozzle separation distance. For the appropriate range of penetration Knudsen numbers, the analysis of this shock structure can be suitably accomplished through a Monte Carlo simulation. This paper describes the application of a general three-dimensional Monte Carlo simulation on the Connection Machine CM-2 to the analysis of the plume self-interaction shock in the near field. Results are presented for two cases, corresponding to a small and a large nozzle separation distance. The results correctly reproduce the expected flow features and demonstrate the ability of this method to properly simulate the start of the plume self-interaction shock. This has significance not only for allowing analysis of the self-interacting plume in the near field, but also for allowing the subsequent simulation of the far field flow through the use of a continuation downstream exit boundary.

Dagum, Leonardo↗

Application of Monte Carlo Technique for Determining Maneuvering Loads from Statistical Information on Airplane Motions

Results of a statistical analysis of horizontal-tail loads on a fighter airplane are presented. The data were obtained from a number of operational training missions with flight at altitudes up to about 50,000 feet and at Mach numbers up to 1.22. The analysis was performed to determine the feasibility of calculating horizontal-tail load from data on the flight conditions and airplane motions. In the analysis the calculated loads are compared with the measured loads for the different types of missions performed. The loads were calculated by two methods: a direct approach and a Monte Carlo technique. The procedures used and some of the problems associated with the data analysis are discussed. frequencies of occurrence of tail loads of given magnitudes are derived from statistical information on the flight quantities. In the direct method, a time history of tail load is calculated from time-history measurements of the flight quantities. The Monte Carlo method could be useful for extending loads information for design of prospective airplanes . For the Monte Carlo method, the The results indicate that the accuracy of loads, regardless of the method used for calculation, is largely dependent on the knowledge of the pertinent airplane aerodynamic characteristics and center-of-gravity location. In addition, reliable Monte Carlo results require an adequate sample of statistical data and a knowledge of the more important statistical dependencies between the various flight conditions and airplane motions.

Hamer, Harold A.↗

Uncertainty Analysis of the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel 2019 Characterization Test 14-Foot, 5.8 Percent Test Section Porosity

This paper presents methods and results of a measurement uncertainty analysis that was performed for the 8- by 6-foot Supersonic Wind Tunnel located at the NASA Glenn Research Center. The uncertainty analysis is based on data collected during a characterization test following major structural modifications to the facility and upgrades to the facility’s data and control systems. The statistical methods and engineering judgments used to estimate elemental uncertainties are described in this report. The Monte Carlo method of propagating uncertainty was selected to determine the uncertainty of calculated variables of interest. A description of the Monte Carlo method as applied for this analysis is provided. The primary variable of interest (VOI) for this facility is free stream Mach number. In addition to determining the uncertainty in Mach number, the uncertainty in free stream values of static pressure, total pressure, dynamic pressure, total temperature, static temperature, Reynolds number, and air speed were also calculated. Uncertainty results are presented as random (unpredictable variation in repeated measurements), systematic (potential offset between observed and true value), and total (random and systematic combined) uncertainty for Mach number. Systematic uncertainty results are presented for the other VOIs. Individual uncertainty sources are presented both dimensionally and as percent contributions to uncertainty in all VOIs, to aid in the identification of primary uncertainty sources.

Uncertainty Analysis↗

Uncertainty Analysis of NASA Glenn's 8- by 6-Foot Supersonic Wind Tunnel

An analysis was performed to determine the measurement uncertainty of the Mach Number of the 8- by 6-foot Supersonic Wind Tunnel at the NASA Glenn Research Center. This paper details the analysis process used, including methods for handling limited data and complicated data correlations. Due to the complexity of the equations used, a Monte Carlo Method was utilized for this uncertainty analysis. A summary of the findings are presented as pertains to understanding what the uncertainties are, how they impact various research tests in the facility, and methods of reducing the uncertainties in the future.

uncertainty↗

Assessing GEO and LEO Repeating Conjunctions Using High Fidelity Brute Force Monte Carlo Simulations

Probability of collision (P(sub c)) estimates for Earth-orbiting satellites typically assume a temporally-isolated conjunction event. However, under certain conditions two objects may experience multiple high-risk close approach events over the course of hours or days. In these repeating conjunction cases, the P(sub c) accumulates as each successive encounter occurs. The NASA Conjunction Assessment Risk Analysis team has updated its “brute force Monte Carlo” (BFMC) software to estimate such accumulating P(sub c) values for repeating conjunctions. This study describes the updated BFMC algorithm and discusses the implications for conjunction risk assessment.

Baars, Luis↗

Assessing Geo and Leo Repeating Conjunctions Using High Fidelity Brute Force Monte Carlo Simulations

Probability of collision (P(sub c)) estimates for Earth-orbiting satellites typically assume a temporally-isolated conjunction event. However, under certain conditions two objects may experience multiple high-risk close approach events over the course of hours or days. In these repeating conjunction cases, the P(sub c) accumulates as each successive encounter occurs. The NASA Conjunction Assessment Risk Analysis team has updated its “brute force Monte Carlo” (BFMC) software to estimate such accumulating P(sub c) values for repeating conjunctions. This study describes the updated BFMC algorithm and discusses the implications for conjunction risk assessment.

Baars, Luis↗

Stochastic Verification by Analysis for Autonomous Systems Management Architecture (ASMA)

The Gateway Vehicle Systems Manager (VSM) is the top-level of a distributed, hierarchical software control system. VSM is data-driven and will make decisions related to mission, fault, resource management and vehicle control. These attributes combined with a high degree of autonomy make it susceptible to emergent behavior. In order to achieve the high level of confidence needed in this critical system, the VSM team has developed a multifaceted verification strategy employing traditional verification techniques, simulation, model checking, and runtime verification. Individual algorithms are verified using conventional testing and model checking using assume-guarantee contracts. A discrete event-based simulation approach is being developed to verify timelines. This presentation describes an enhancement to the verification approach using analysis to enhance system robustness by detecting and resolving the potential for emergent behavior. The verification by analysis employs a Software in the Loop (SITL) environment with real flight software executing on emulated processors, simulations of vehicle subsystems, flight dynamics, and human inputs. Since the possible input space and configuration data set are too large for exhaustive testing, a Monte Carlo approach is used to cover feasible scenarios, augmented with corner cases and known higher-risk scenarios. A key problem in using Monte Carlo-based system verification is evaluating test results to ensure that system behavior is correct. The presentation describes the approach the VSM team uses to monitor behavior for compliance with predetermined boundaries and to identify anomalous behavior for further analysis. This presentation describes the multi-level systems approach to verification, and the simulation-based layer that covers the feasible state space: 1. Overview of the Gateway VSM 2. Special challenges due to heterogeneous, hierarchical architecture 3. Modeling and simulation environment using flight software and system simulations 4. Developing input sets to ensure state-space coverage 5. Developing model and data configuration sets to ensure model coverage 6. Interpreting results without predetermined outcomes 7. Lessons learned and future work

Verification and Validation↗

Estimating the Rate of Occurrence of Renal Stones in Astronauts

Changes in urine chemistry, during and post flight, potentially increases the risk of renal stones in astronauts. Although much is known about the effects of space flight on urine chemistry, no inflight incidence of renal stones in US astronauts exists and the question "How much does this risk change with space flight?" remains difficult to accurately quantify. In this discussion, we tackle this question utilizing a combination of deterministic and probabilistic modeling that implements the physics behind free stone growth and agglomeration, speciation of urine chemistry and published observations of population renal stone incidences to estimate changes in the rate of renal stone presentation. The modeling process utilizes a Population Balance Equation based model developed in the companion IWS abstract by Kassemi et al. (2016) to evaluate the maximum growth and agglomeration potential from a specified set of urine chemistry values. Changes in renal stone occurrence rates are obtained from this model in a probabilistic simulation that interrogates the range of possible urine chemistries using Monte Carlo techniques. Subsequently, each randomly sampled urine chemistry undergoes speciation analysis using the well-established Joint Expert Speciation System (JESS) code to calculate critical values, such as ionic strength and relative supersaturation. The Kassemi model utilizes this information to predict the mean and maximum stone size. We close the assessment loop by using a transfer function that estimates the rate of stone formation from combining the relative supersaturation and both the mean and maximum free stone growth sizes. The transfer function is established by a simulation analysis which combines population stone formation rates and Poisson regression. Training this transfer function requires using the output of the aforementioned assessment steps with inputs from known non-stone-former and known stone-former urine chemistries. Established in a Monte Carlo system, the entire renal stone analysis model produces a probability distribution of the stone formation rate and an expected uncertainty in the estimate. The utility of this analysis will be demonstrated by showing the change in renal stone occurrence predicted by this method using urine chemistry distributions published in Whitson et al. 2009. A comparison to the model predictions to previous assessments of renal stone risk will be used to illustrate initial validation of the model.

gravitational physiology↗

Investigation of the external flow analysis for density measurements at high altitude

A Monte Carlo program was developed for modeling the flow field around the space shuttle in the vicinity of the shuttle upper atmosphere mass spectrometer experiment. The operation of the EXTERNAL code is summarized. Issues associated with geometric modeling of the shuttle nose region and the modeling of intermolecular collisions including rotational energy exchange are discussed as well as a preliminary analysis of vibrational excitation and dissociation effects. The selection of trial runs is described and the parameters used for them is justified. The original version and the modified INTERNAL code for the entrance problem are reviewed. The code listing is included.

Bienkowski, G. K.↗