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 469 records · Page 26

Probabilistic Design of a Wind Tunnel Model to Match the Response of a Full-Scale Aircraft

approach is presented for carrying out the reliability-based design of a plate-like wing that is part of a wind tunnel model. The goal is to design the wind tunnel model to match the stiffness characteristics of the wing box of a flight vehicle while satisfying strength-based risk/reliability requirements that prevents damage to the wind tunnel model and fixtures. The flight vehicle is a modified F/A-18 aircraft. The design problem is solved using reliability-based optimization techniques. The objective function to be minimized is the difference between the displacements of the wind tunnel model and the corresponding displacements of the flight vehicle. The design variables control the thickness distribution of the wind tunnel model. Displacements of the wind tunnel model change with the thickness distribution, while displacements of the flight vehicle are a set of fixed data. The only constraint imposed is that the probability of failure is less than a specified value. Failure is assumed to occur if the stress caused by aerodynamic pressure loading is greater than the specified strength allowable. Two uncertain quantities are considered: the allowable stress and the thickness distribution of the wind tunnel model. Reliability is calculated using Monte Carlo simulation with response surfaces that provide approximate values of stresses. The response surface equations are, in turn, computed from finite element analyses of the wind tunnel model at specified design points. Because the response surface approximations were fit over a small region centered about the current design, the response surfaces were refit periodically as the design variables changed. Coarse-grained parallelism was used to simultaneously perform multiple finite element analyses. Studies carried out in this paper demonstrate that this scheme of using moving response surfaces and coarse-grained computational parallelism reduce the execution time of the Monte Carlo simulation enough to make the design problem tractable. The results of the reliability-based designs performed in this paper show that large decreases in the probability of stress-based failure can be realized with only small sacrifices in the ability of the wind tunnel model to represent the displacements of the full-scale vehicle.

Mason, Brian H.↗

Pterodactyl: Integrated Control Design for Precision Targeting of Deployable Entry Vehicles

Deployable Entry Vehicles (DEVs) enable in-situ scientific exploration at destinations with atmospheres across the solar system. Because they stow in a compact form and deploy only when ready to enter the atmosphere, DEVs relax the volume constraint imposed by rigid aeroshells. This work seeks to do for a DEV what the Wright Brothers did to propel modern day aviation: develop the guidance and control (G&C) methods that will make maneuvering and precision landing of DEV a reality. The Pterodactyl project objective is to deliver an integrated G&C methodology for a DEV, based on a detailed analysis that utilizes a Multi-disciplinary, Design, Analysis and Optimization (MDAO) framework. The current state-of-the-art for blunt body entry, G&C is rooted in the precision landing of vehicles such as Mars Science Laboratory (MSL) and Apollo, which used a propulsive reaction control system (RCS) to steer. Recent research has taken a particular interest in non-propulsive control for DEVs, including direct force control (angle of attack modulation via control surfaces or mass movement) and drag modulation (discrete change in ballistic coefficient). Using the MDAO framework that includes a guidance and control model to explore multiple control concepts for a DEV will shed light on the best design approach for these vehicles. In Pterodactyl, we will complete this study for a novel DEV concept, and then we will fabricate a functional prototype to help validate the design. The project is expected to down-select to a final control architecture by the end of 2018, and complete fabrication of the prototype by the end of 2019.The DEV chosen for detailed study in this project is the Adaptable Deployable Entry and Placement Technology (ADEPT). ADEPT uses a revolutionary 3D-woven carbon fabric that is foldable, can serve as primary structure, and can survive the extreme heating environment of atmospheric entry. The specific configuration of ADEPT under investigation is called Lifting Nano-ADEPT (LNA). LNA is designed for secondary payloads missions that require precision landing either for scientific objectives at a target destination or for payload recovery at Earth.The MDAO framework being created through this research, called COBRA-Pt (Composite Beam Roll-Up Solar Array-Prototype), will combine three critical elements of the system design: a guidance algorithm with Monte Carlo, a parametric control model, and vehicle geometry details. Novel control models being studied are deployable aerodynamic surfaces as well as shape morphing. These concepts will be compared at the system level with a more traditional propulsive RCS by comparing several key performance parameters. Upon completion of the design study, a functional prototype of LNA will be fabricated that will include the integration of guidance software and relevant control actuators. We expect this study will provide critical data that could feed into the development of an Earth-based flight test of LNA. The COBRA-Pt framework will provide a modular system by which to study any DEV concept in any atmosphere.

Smith, Brandon↗

Performance Risk Model (PRisM) Proof-of-Concept: An Operational Decision Support Tool to Predict Crew Performance in Space from Available Performance Tests

The Crew Health and Performance-Probabilistic Risk Assessment (CHP-PRA) team at NASA Glenn Research Center has developed a range of tools to evaluate astronaut health during spaceflight and to optimize the medical set required for missions. Among these, the Performance Risk Model (PRisM) represents a novel advancement, extending CHP-PRA’s focus beyond medical systems into the domain of human performance. Such tool could be pivotal in optimizing astronaut capabilities during space travel, thereby enhancing overall mission success. PRisM leverages data from well-established performance assessments conducted during previous crewed space missions and Earth analogs to strategically predict outcomes for planned tasks, even when direct performance testing has not been conducted for those specific tasks. To evaluate performance, PRisM references the various metrics outlined in NASA-STD-3001 and employs a methodology to integrate different performance scales. This framework analyzes the contribution of various human system task categories (HSTCs) to task execution and compares these contributions to the HSTCs’ involvement in other known tests. The model further employs a Monte Carlo simulation to sample performance scores from their distribution in operationally relevant tests such as those in Mulavara et al. (2018) and, by leveraging similarities in HSTC involvement, transfers this knowledge to actual mission tasks, such as those outlined in the "Human Exploration of Mars: Preliminary List of Crew Tasks”. The current PRisM proof-of-concept includes analyses of the impact of exercise and specific medications on astronauts’ performance, with provisions to incorporate additional countermeasures as data becomes available. Furthermore, the tool is customizable to include any system necessary to fully encompass the domain of human systems and can be adapted to evaluate performance for any spaceflight activities as requested by operational stakeholders. PRisM has the potential to assist the Human Research Program in exploring the capabilities trade space for optimized crew performance.

performance modeling↗

Testing Models of Resonant Compton Scattering in X-Ray Pulsars

Over the performance period covered by the grant, the principal investigator modified a Monte Carlo Compton scattering code to model the propagation of x-rays through the magnetosphere of accreting neutron stars. These modifications were made to enable the author to compare the observations of x-ray pulsars to theoretical models of the system. The original code was designed to study relativistic plasmas with one of two geometries: a plane parallel plasma with a differential relativistic bulk velocity, and a static spherically symmetric plasma.- This code did not treat gravitational bending or bulk motion in the magnetosphere of a neutron star. Under the grant, the author incorporated code to trace light paths in a Schwarzschild metric. The code was modified to keep track of the photon polarization during propagati on. The investigator also modified the code so that bulk motion in an axisymmetric system is treated properly. An approximate treatment for resonant Compton scattering was added to the code. Finally, code was added that creates model observables that can be compared to observations, such as projected x-ray emission maps and energy-dependent light curves. Comparison to observations is now commencing.

Brainerd, Jerome J.↗

Simulation of Vacuum Chamber Pressure Distribution with Surrogate Modeling and Uncertainty Quantification

A major challenge in understanding differences in electric propulsion performance in ground tests and in space operations concerns the pressure distribution within the test vacuum chamber. The chamber backpressure is much higher than experienced in space, modifying thruster performance and plume dynamics. Numerical simulation is a key element to determining the background conditions in non-ideal vacuum chamber environments. An important parameter for the accurate simulation of chamber backpressure is the probability that an atom will stick to a cryogenic panel used to pump away the plume gases. This quantity can be used to model vacuum pumps in particle-based kinetic numerical methods. In this work, a three-dimensional direct simulation Monte Carlo code is used to model neutral xenon atoms flowing from the anode of the H9 Hall Effect Thruster within the University of Michigan’s Large Vacuum Test Facility. Simulated pressures are compared with ion gauge pressure measurements to infer the effective sticking coefficient of the chamber’s vacuum pumps. A pressure predicting surrogate model is developed for inference of pump sticking coefficients and for uncertainty quantification. This information enables accurate and useful kinetic simulations of electric propulsion thruster plasma plumes in vacuum chambers.

DSMC↗

A study of two statistical methods as applied to shuttle solid rocket booster expenditures

The state probability technique and the Monte Carlo technique are applied to finding shuttle solid rocket booster expenditure statistics. For a given attrition rate per launch, the probable number of boosters needed for a given mission of 440 launches is calculated. Several cases are considered, including the elimination of the booster after a maximum of 20 consecutive launches. Also considered is the case where the booster is composed of replaceable components with independent attrition rates. A simple cost analysis is carried out to indicate the number of boosters to build initially, depending on booster costs. Two statistical methods were applied in the analysis: (1) state probability method which consists of defining an appropriate state space for the outcome of the random trials, and (2) model simulation method or the Monte Carlo technique. It was found that the model simulation method was easier to formulate while the state probability method required less computing time and was more accurate.

Perlmutter, M.↗

Chain pooling to minimize prediction error in subset regression

Monte Carlo studies using population models intended to represent response surface applications are reported. Simulated experiments were generated by adding pseudo random normally distributed errors to population values to generate observations. Model equations were fitted to the observations and the decision procedure was used to delete terms. Comparison of values predicted by the reduced models with the true population values enabled the identification of deletion strategies that are approximately optimal for minimizing prediction errors.

Holms, A. G.↗

A general cratering-history model and its implications for the lunar highlands

Through analysis of a large number of Monte Carlo and Markov Chain simulations, a model for determining crater accumulation and crater obliteration histories has been derived. The model generally applies to populations of large craters. It predicts that the following relationships hold for subequilibrium-density crater populations: (1) the more negative the production function's exponent, alpha (N near D super alpha) the lower the crater density at which the population size-frequency distribution will significantly depart from its production function; (2) the more negative the production function's exponent, the less obliteration a crater population will sustain after a set number of impacts. Application of the model to the lunar highlands implies (1) the production function for the large craters is highly structured, resembling the observed size-frequency distribution and not the function N near D to the -2; (2) even the densely cratered highlands have not attained crater saturation or equilibrium. Direct simulations of the highlands' crater population supports the model's implications.

Woronow, A.↗

Verification of reflectance models in turbid waters

Inherent optical parameters of very turbid waters were used to evaluate existing water reflectance models. Measured upwelling radiance spectra and Monte Carlo simulations of the radiative transfer equations were compared with results from models based upon two flow, quasi-single scattering, augmented isotropic scattering, and power series approximation. Each model was evaluated for three separate components of upwelling radiance: (1) direct sunlight; (2) diffuse skylight; and (3) internally reflected light. Limitations of existing water reflectance models as applied to turbid waters and possible applications to the extraction of water constituent information are discussed.

Tanis, F. J.↗

Luminosity correlations in quasars

Simulations are conducted with and without flux thresholds in an investigation of quasar luminosity correlations by means of a Monte Carlo analysis, for various model distributions of quasars in X-rays and optical luminosity. For the case where the X-ray photons are primary, an anticorrelation between X-ray-to-optical luminosity ratio and optical luminosity arises as a natural consequence which resembles observations. The low optical luminosities of X-ray selected quasars can be understood as a consequence of the same effect, and similar conclusions may hold if the X-ray and optical luminosities are determined independently by a third parameter, although they do not hold if the optical photons are primary. The importance of such considerations is demonstrated through a reanalysis of the published X-ray-to-optical flux ratios for the 3CR sample.

Chanan, G. A.↗

Sensitivity of airborne fluorosensor measurements to linear vertical gradients in chlorophyll concentration

A semianalytic Monte Carlo radiative transfer simulation model for airborne laser fluorosensors has been extended to investigate the effects of inhomogeneities in the vertical distribution of phytoplankton concentrations in clear seawater. Simulation results for linearly varying step concentrations of chlorophyll are presented. The results indicate that statistically significant differences can be seen under certain conditions in the water Raman-normalized fluorescence signals between nonhomogeneous and homogeneous cases. A statistical test has been used to establish ranges of surface concentrations and/or verticl gradients in which calibration by surface samples would by inappropriate, and the results are discussed.

Venable, D. D.↗

Universal amplitude ratios and the interfacial tension near consolute points of binary liquid mixtures

The interfacial tension of the binary liquid mixtures triethylamine/heavy water, triethylamine/water, and cyclohexane/methanol near their consolute temperatures is determined experimentally by measuring the densities of the phases and the capillary lengths. The results are presented in tables and graphs and compared with published data obtained by other methods and with theoretical predictions. The findings are shown to be in good agreement with previous experiments, confirming the discrepancy between temperature-independent universal amplitude ratios calculated from such data and those predicted by renormalization-group calculations or Monte Carlo simulations of Ising models.

Moldover, M. R.↗

Neutral cometary atmospheres. V - C2 and CN in comets

The radial Haser scale lengths for the production and destruction of C2 and CN are calculated on the basis of brightness profiles for comets Bennett 1970 II and Kohoutek 1973 XII, applying the Monte Carlo particle-trajectory model of Combi and Delsemme (1980). Findings reported include C2/CN ratio independent of heliocentric distance, C2 ejection velocity from its parent molecule about 0.5 km/s, and 1-AU photochemical lifetime 31,000 s for the C2 parent molecule and 120,000 s for C2.

Combi, M. R.↗

Evaluation of concepts and methods for extraction of flow parameters from laser transit anemometry data

A Monte-Carlo type simulation system was developed to test concepts and associated algorithms proposed by Mayo and Smart (1984) for the extraction of two-dimensional flow parameters from laser transit anemometry data sets. The Monte-Carlo LTA simulation program models a flow field as well as the sample volume geometry of the LTA system, thereby permitting simulation of two-dimensional data acquisition in a three-dimensional flow field. The results of these tests on the proposed processing concepts suggest that the instrument is capable of measuring mean velocities to less than +0.15 percent error and flow angles to less then +0.3 percent for turbulence intensities of up to 15.0 percent.

Humphreys, W. M., Jr.↗

An Exploratory Exercise in Taguchi Analysis of Design Parameters: Application to a Shuttle-to-space Station Automated Approach Control System

The chief goals of the summer project have been twofold - first, for my host group and myself to learn as much of the working details of Taguchi analysis as possible in the time allotted, and, secondly, to apply the methodology to a design problem with the intention of establishing a preliminary set of near-optimal (in the sense of producing a desired response) design parameter values from among a large number of candidate factor combinations. The selected problem is concerned with determining design factor settings for an automated approach program which is to have the capability of guiding the Shuttle into the docking port of the Space Station under controlled conditions so as to meet and/or optimize certain target criteria. The candidate design parameters under study were glide path (i.e., approach) angle, path intercept and approach gains, and minimum impulse bit mode (a parameter which defines how Shuttle jets shall be fired). Several performance criteria were of concern: terminal relative velocity at the instant the two spacecraft are mated; docking offset; number of Shuttle jet firings in certain specified directions (of interest due to possible plume impingement on the Station's solar arrays), and total RCS (a measure of the energy expended in performing the approach/docking maneuver). In the material discussed here, we have focused on single performance criteria - total RCS. An analysis of the possibility of employing a multiobjective function composed of a weighted sum of the various individual criteria has been undertaken, but is, at this writing, incomplete. Results from the Taguchi statistical analysis indicate that only three of the original four posited factors are significant in affecting RCS response. A comparison of model simulation output (via Monte Carlo) with predictions based on estimated factor effects inferred through the Taguchi experiment array data suggested acceptable or close agreement between the two except at the predicted optimum point, where a difference outside a rule-of-thumb bound was observed. We have concluded that there is most likely an interaction effect not provided for in the original orthogonal array selected as the basis for our experimental design. However, we feel that the data indicates that this interaction is a mild one and that inclusion of its effect will not alter the location of the optimum.

Deal, Don E.↗

Dynamic Stability of Uncertain Laminated Beams Under Subtangential Loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. In this research, we focus on the deterministic and probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. Thus a shear-deformable laminated beam element, including warping effects, is derived to study the deterministic and probabilistic response of laminated beams. This twenty-one degrees of freedom element can be used for solving both static and dynamic problems. In the first-order shear deformable model used here we have employed a more accurate method to obtain the transverse shear correction factor. The dynamic version of the principle of virtual work for laminated composites is expressed in its nondimensional form and the element tangent stiffness and mass matrices are obtained using analytical integration The stability is studied by giving the structure a small disturbance about an equilibrium configuration, and observing if the resulting response remains small. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: Exact Monte Carlo Simulation, Sensitivity Based Monte Carlo Simulation, and Probabilistic FEA. These methods were integrated into the developed finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems, as well as to study the stability analysis, using the dynamic criterion.

Goyal, Vijay K.↗

Predictions of secondary neutrons and their importance to radiation effects inside the International Space Station

As part of a study funded by NASA MSFC to assess thecontribution of secondary particles in producing radiation damage to optoelectronics devices located on the International Space Station (IS), Monte Carlo calculations have been made to predict secondary spectra vs. shielding inside ISS modules and in electronics boxes attached on the truss (Armstrong and Colborn, 1998). The calculations take into account secondary neutron, proton, and charged pion production from the ambient galactic cosmic-ray (GCR) proton, trapped proton, and neutron albedo environments. Comparisons of the predicted neutron spectra with measurments made on the Mir space station and other spacecraft have also been made (Armstrong and Colborn, 1998). In this paper, some initial results from folding the predicted neutron spectrum inside ISS modules from Armstrong and Colborn (1998) with several types of radiation effects response functions related to electronics damage and astronaut-dose are given. These results provide an estimate of the practical importance of neutrons compared to protons in assessing radiation effects for the ISS. Also, the important neutron energy ranges for producing these effects have been estimated, which provides guidance for onboard neutron measurement requirements.

Neutrons↗