Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational Monte Carlo”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Aerocapture Performance Analysis for a Neptune-Triton Exploration Mission

A systems analysis has been conducted for a Neptune-Triton Exploration Mission in which aerocapture is used to capture a spacecraft at Neptune. Aerocapture uses aerodynamic drag instead of propulsion to decelerate from the interplanetary approach trajectory to a captured orbit during a single pass through the atmosphere. After capture, propulsion is used to move the spacecraft from the initial captured orbit to the desired science orbit. A preliminary assessment identified that a spacecraft with a lift to drag ratio of 0.8 was required for aerocapture. Performance analyses of the 0.8 L/D vehicle were performed using a high fidelity flight simulation within a Monte Carlo executive to determine mission success statistics. The simulation was the Program to Optimize Simulated Trajectories (POST) modified to include Neptune specific atmospheric and planet models, spacecraft aerodynamic characteristics, and interplanetary trajectory models. To these were added autonomous guidance and pseudo flight controller models. The Monte Carlo analyses incorporated approach trajectory delivery errors, aerodynamic characteristics uncertainties, and atmospheric density variations. Monte Carlo analyses were performed for a reference set of uncertainties and sets of uncertainties modified to produce increased and reduced atmospheric variability. For the reference uncertainties, the 0.8 L/D flatbottom ellipsled vehicle achieves 100% successful capture and has a 99.87 probability of attaining the science orbit with a 360 m/s V budget for apoapsis and periapsis adjustment. Monte Carlo analyses were also performed for a guidance system that modulates both bank angle and angle of attack with the reference set of uncertainties. An alpha and bank modulation guidance system reduces the 99.87 percentile DELTA V 173 m/s (48%) to 187 m/s for the reference set of uncertainties.

Starr, Brett R.↗
NASA NTRS1986

Cometary atmospheres: Modeling the spatial distribution of observed neutral radicals

Progress during the second year of a program of research on the modeling of the spatial distributions of cometary radicals is discussed herein in several major areas. New scale length laws for cometary C2 and CN were determined which explain that the previously-held apparent drop of the C2/CN ratio for large heliocentric distances does not exist and that there is no systematic variation. Monte Carlo particle trajectory model (MCPTM) analysis of sunward and anti-sunward brightness profiles of cometary C2 was completed. This analysis implies a lifetime of 31,000 seconds for the C2 parent and an ejection speed for C2 of approximately 0.5 km/sec upon dissociation from the parent. A systematic reanalysis of published C3 and OH data was begun. Preliminary results find a heliocentric distance dependence for C3 scale lengths with a much larger variation than for C2 and CN. Scale lengths for OH are generally somewhat larger than currently accepted values. The MCPTM was updated to include the coma temperature. Finally, the collaborative effort with the University of Arizona programs has yielded some preliminary CCD images of Comet P/Halley.

Combi, Michael R.↗
NASA NTRS1994

The global resurfacing of Venus

The impact cratering record on Venus is unique among the terrestrial planets. Fully 84% of the craters are in pristine condition, and only 12% are fractured. Remarkably, only 2.5% of the craters and crater-related features are embayed by lava, although intense volcanism and tectonism have affected the entire planet. Furthermore, the spatial and hypsometric distribution of the craters is consistent with a completely random one, including stochastic variations. Monte Carlo simulations of equilibrium resurfacing models result in a minimum of 17 times more embayed craters than observed, or unobserved nonrandom crater distributions for resurfacing areas between 0.03% and 100% of the planet's surface. These models also are not consistent with the number and nonrandom distribution of volcanoes, and the nonrandom distribution of embayed and heavily fractured craters. The constraints imposed by the cratering record strongly indicate that Venus experienced a global resurfacing event about 300 m.y. ago followed by a dramatic reduction of volcanism and tectonism. This global resurfacing event ended abruptly (less than 10 m.y.). The present crater population has accumulated since then and remains largely intact. Thermal history models suggest that similar global resurfacing events probably occured episodically in the past. We show that neither the present level and style of geologic activity nor anything less than global resurfacing could have produced the observed cratering record. The effects of recent geologic activity are much less than those of the earlier global resurfacing event, when the record of all the early heavy bombardment and much of the later light bombardment was erased from the surface by massive volcanism and tectonic activity. Episodic regional resurfacing events that had global effects also occurred on Earth (e.g., the mid-Cretaceous superplume) and probably on Mars. On Mars they may have triggered the catastrophic releases of water that formed the outflow channels.

Strom, Robert G.↗
NASA NTRS2004

Uncertainty Analysis on Heat Transfer Correlations for RP-1 Fuel in Copper Tubing

NASA is studying kerosene (RP-1) for application in Next Generation Launch Technology (NGLT). Accurate heat transfer correlations in narrow passages at high temperatures and pressures are needed. Hydrocarbon fuels, such as RP-1, produce carbon deposition (coke) along the inside of tube walls when heated to high temperatures. A series of tests to measure the heat transfer using RP-1 fuel and examine the coking were performed in NASA Glenn Research Center's Heated Tube Facility. The facility models regenerative cooling by flowing room temperature RP-1 through resistively heated copper tubing. A Regression analysis is performed on the data to determine the heat transfer correlation for Nusselt number as a function of Reynolds and Prandtl numbers. Each measurement and calculation is analyzed to identify sources of uncertainty, including RP-1 property variations. Monte Carlo simulation is used to determine how each uncertainty source propagates through the regression and an overall uncertainty in predicted heat transfer coefficient. The implications of these uncertainties on engine design and ways to minimize existing uncertainties are discussed.

Driscoll, E. A.↗
NASA NTRS2016

Orbit Determination Sensitivity Analysis for the Europa Multiple Flyby Mission Concept

This paper details the results of parametric variations on the notional Europa Mission Concept tour navigation strategy on orbit determination delivery and knowledge errors and the associated statistical ∆V consumption. Approach maneuver placement at encounter minus two and a half days and at encounter minus three and a half days is compared to the reference encounter minus three days location as well as variations to the data cutoff for approach maneuver design, which is baselined one day prior to maneuver execution. An execution error model with fixed and proportional components assigned to the magnitude and pointing directions is applied to each maneuver in a Monte Carlo simulation. Variations on the fixed and proportional contributions to both component errors are simulated and analyzed to show the effect of varying thruster characterization levels. The effects of varying levels of a priori satellite ephemeris errors are also characterized in this study.

Boone, Dylan R.↗
NASA NTRS1989

Monte Carlo simulations of the diurnal variation in seismic detection rate of sporadic meteoroid impacts on the moon

The rate of detection of meteoroid impacts on the moon by the lunar seismic network shows a characteristic diurnal variation. Assuming that these meteoroids have a flux and a preimpact orbital distribution similar to that of fireballs observed by terrestrial camera networks, one can compute the expected diurnal variation for a given set of parameters that describe the seismic wave generation and transmission on the moon. An iterative process to match the theoretical variation with the observed one has led us to the following results: (1) the majority of the detected impact events occur within a closer range of the network than was believed earlier. This results in higher meteoroid flux estimates from lunar seismic data that agree with the terrestrially measured flux. (2) For meteoroid masses smaller than 1000 g, seismic amplitude is approximately proportional to the one-fifth power of the impact speed; for larger masses it is approximately proportional to the eight-fifth power, provided that the terrestrial meteor data used for analysis are not biased. (3) Seismic efficiency of meteoroids smaller than 1000 g is significantly less than that of large meteoroids. (4) Using orbits of fireballs that represent meteorites, we predict that the share of meteorites among the detected impacts is approximately 15 percent assuming that seismic efficiency of the high-density meteorites is the same as that of average meteoroids. A greatly increased seismic efficiency for these high-density objects is not likely.

Oberst, J.↗
NASA NTRS1983

Monte Carlo simulation of the effect of soil moisture variation on the microwave emission from soils

In this paper, results of a Monte Carlo simulation of the effect of noise on the relationship between the microwave emissivity of soil and its moisture content are presented. It is found that whenever the magnitude of the noise for the independent variable, in this case the soil moisture, is increased, both the slope of the regression and the correlation coefficient decrease. In particular, when the noise has magnitude equivalent to a coefficient of variation of 0.25, the slope and correlation coefficient are in good agreement with those obtained from the data of a 21-cm airborne microwave radiometer which was flown over a test site in hand county, South Dakota. The comparison was made using a linear relationship to determine the estimated emissivity from the ground measurements of soil moisture. The linear relationship was derived from a radiative transfer model calculation of the microwave emissivities using realistic soil-moisture profiles. The effect of surface roughness was included in the relationship, and the variability of the surface roughness was also simulated by a Monte Carlo technique.

Mo, T.↗
NASA NTRS1990

Planetary atmosphere representation for Monte-Carlo simulation of aerospace vehicle trajectories

Methods of simulating large scale random atmospheric variations are examined for use in Monte Carlo simulations of hypervelocity vehicles. Random density, temperature, and pressure variations have a significant impact on the performance of such vehicles. Methods for generating random atmosphere profiles are described and compared in the context of a Monte Carlo simulation of the Aeroassist Flight Experiment (AFE). It is seen that these techniques result in statistically significant differences in key simulation point statistics.

Dutton, K.↗
NASA NTRS1978

Three-dimensional radiative exchange factors for nongray, nondiffuse surfaces

A computer program using a variation of the Monte Carlo method was developed to predict radiative exchange factors in three-dimensional configurations with curved surfaces. The model allows for arbitrary emission, absorption, and reflection characteristics and for radiation exchange between surfaces as well as external irradiation through openings. It was shown that exchange factors can be determined efficiently by using the normal Monte Carlo method for closed configurations and/or high surface emissivities, while energy partitioning is vastly superior in the case of open configurations and/or low surface emissivities. Comparison with some experiments performed on the Space Shuttle heat rejection system demonstrates the model's accuracy as well as its superiority over a simplified diffuse-surface analysis, which requires similar amounts of computer time.

Modest, M. F.↗
NASA NTRS1980

Simulation studies of a model of high-density metallic hydrogen

Upper bounds for the ground-state energies of liquid and solid phases of metallic hydrogen and metallic deuterium have been calculated with variational methods and Monte Carlo techniques. At four densities (0.8, 1.2, 1.36, and 1.488) crystalline phases are clearly preferred in the sense that the energy difference, when compared to the liquid, is in excess of the errors inherent in the numerical procedures. At a fifth density (1.6), the energy differences between solid and liquid phases are smaller than these errors.

Mon, K. K.↗
NASA NTRS2004

Impact of Variations on 1-D Flow in Gas Turbine Engines via Monte Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRS2003

A Study of the Impact of Variations on Aerodynamic Flow in Gas Turbine Engines via Monte-Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented with a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, under the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions that motivate this research are: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance. This paper presents an initial attempt to the flow behavior in terms (perturbation) using simulation outputs from a customer-deck model obtained from Pratt&Whitney, (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine. Flow behavior is analyzed for the high pressure compressor and combustion chamber employing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000-sample simulation. Furthermore, the characteristic waves and eave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures; which are clearly the case in the combustion chamber and high-pressure turbine. Finally a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRS1976

Solid-propellant rocket motor internal ballistic performance variation analysis, phase 2

The Monte Carlo method was used to investigate thrust imbalance and its first time derivative throughtout the burning time of pairs of solid rocket motors firing in parallel. Results obtained compare favorably with Titan 3 C flight performance data. Statistical correlations of the thrust imbalance at various times with corresponding nominal trace slopes suggest several alternative methods of predicting thrust imbalance. The effect of circular-perforated grain deformation on internal ballistics is discussed, and a modified design analysis computer program which permits such an evaluation is presented. Comparisons with SRM firings indicate that grain deformation may account for a portion of the so-called scale factor on burning rate between large motors and strand burners or small ballistic test motors. Thermoelastic effects on burning rate are also investigated. Burning surface temperature is calculated by coupling the solid phase energy equation containing a strain rate term with a model of gas phase combustion zone using the Zeldovich-Novozhilov technique. Comparisons of solutions with and without the strain rate term indicate a small but possibly significant effect of the thermoelastic coupling.

Sforzini, R. H.↗
NASA NTRS1992

Shuttle entry aerothermodynamic flight research - The Orbiter Experiments (OEX) Program

Results of the OEX program are summarized with emphasis on the information on entry aerothermodynamic phenomena derived from Space Shuttle operations. The discussion focuses on OEX experiment complement and operational history, freestream environment and vehicle attitude data, aerodynamic force and moment data, aerodynamic surface data, and vehicle configuration data. Attention is also given to orbiter aerodynamic performance, stability and control, high-altitude atmospheric density variability, direct simulation Monte Carlo validation, orbital drag variation, and computational fluid dynamic technique validation.

Throckmorton, David A.↗
NASA NTRS2013

Vehicle Integrated Prognostic Reasoner (VIPR) Metric Report

This document outlines a set of metrics for evaluating the diagnostic and prognostic schemes developed for the Vehicle Integrated Prognostic Reasoner (VIPR), a system-level reasoner that encompasses the multiple levels of large, complex systems such as those for aircraft and spacecraft. VIPR health managers are organized hierarchically and operate together to derive diagnostic and prognostic inferences from symptoms and conditions reported by a set of diagnostic and prognostic monitors. For layered reasoners such as VIPR, the overall performance cannot be evaluated by metrics solely directed toward timely detection and accuracy of estimation of the faults in individual components. Among other factors, overall vehicle reasoner performance is governed by the effectiveness of the communication schemes between monitors and reasoners in the architecture, and the ability to propagate and fuse relevant information to make accurate, consistent, and timely predictions at different levels of the reasoner hierarchy. We outline an extended set of diagnostic and prognostics metrics that can be broadly categorized as evaluation measures for diagnostic coverage, prognostic coverage, accuracy of inferences, latency in making inferences, computational cost, and sensitivity to different fault and degradation conditions. We report metrics from Monte Carlo experiments using two variations of an aircraft reference model that supported both flat and hierarchical reasoning.

Cornhill, Dennis↗
NASA NTRS2007

Determination of Turboprop Reduction Gearbox System Fatigue Life and Reliability

Two computational models to determine the fatigue life and reliability of a commercial turboprop gearbox are compared with each other and with field data. These models are (1) Monte Carlo simulation of randomly selected lives of individual bearings and gears comprising the system and (2) two-parameter Weibull distribution function for bearings and gears comprising the system using strict-series system reliability to combine the calculated individual component lives in the gearbox. The Monte Carlo simulation included the virtual testing of 744,450 gearboxes. Two sets of field data were obtained from 64 gearboxes that were first-run to removal for cause, were refurbished and placed back in service, and then were second-run until removal for cause. A series of equations were empirically developed from the Monte Carlo simulation to determine the statistical variation in predicted life and Weibull slope as a function of the number of gearboxes failed. The resultant L(sub 10) life from the field data was 5,627 hr. From strict-series system reliability, the predicted L(sub 10) life was 774 hr. From the Monte Carlo simulation, the median value for the L(sub 10) gearbox lives equaled 757 hr. Half of the gearbox L(sub 10) lives will be less than this value and the other half more. The resultant L(sub 10) life of the second-run (refurbished) gearboxes was 1,334 hr. The apparent load-life exponent p for the roller bearings is 5.2. Were the bearing lives to be recalculated with a load-life exponent p equal to 5.2, the predicted L(sub 10) life of the gearbox would be equal to the actual life obtained in the field. The component failure distribution of the gearbox from the Monte Carlo simulation was nearly identical to that using the strict-series system reliability analysis, proving the compatibility of these methods.

Zaretsky, Erwin V.↗