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 1,171 records · Page 65

Aerocapture and aeromaneuvering at Mars

This paper presents the results of a guidance analysis study for aerocapture and aeromaneuvering at Mars. A 4-state, near-optimal, hybrid guidance algorithm is derived for applications throughout the paper. It is based on the critical assumption that Loh's function remains constant between successive guidance updates. Besides the guidance algorithm, the paper presents a performance evaluation of guidance derived from estimated or 'navigated' spacecraft state, provided by an onboard navigation subsystem for real-time applications on Mars. For the performance evaluation, the spacecraft state is augmented by such parameters as the reference atmospheric density and scale-height, the banking angle and angle of attack for spacecraft control and the accelerometer measurement model errors. Finally, the paper examines also, the influence of random uncertainties in atmospheric density on aeroassisted maneuvers, by a Monte Carlo type of simulation.

Vijayaraghavan, A.↗

Functional nanocomposites prepared by self-assembly and polymerization of diacetylene surfactants and silicic acid

Conjugated polymer/silica nanocomposites with hexagonal, cubic, or lamellar mesoscopic order were synthesized by self-assembly using polymerizable amphiphilic diacetylene molecules as both structure-directing agents and monomers. The self-assembly procedure is rapid and incorporates the organic monomers uniformly within a highly ordered, inorganic environment. By tailoring the size of the oligo(ethylene glycol) headgroup of the diacetylene-containing surfactant, we varied the resulting self-assembled mesophases of the composite material. The nanostructured inorganic host altered the diacetylene polymerization behavior, and the resulting nanocomposites show unique thermo-, mechano-, and solvatochromic properties. Polymerization of the incorporated surfactants resulted in polydiacetylene (PDA)/silica nanocomposites that were optically transparent and mechanically robust. Molecular modeling and quantum calculations and (13)C spin-lattice relaxation times (T(1)) of the PDA/silica nanocomposites indicated that the surfactant monomers can be uniformly organized into precise spatial arrangements prior to polymerization. Nanoindentation and gas transport experiments showed that these nanocomposite films have increased hardness and reduced permeability as compared to pure PDA. Our work demonstrates polymerizable surfactant/silica self-assembly to be an efficient, general approach to the formation of nanostructured conjugated polymers. The nanostructured inorganic framework serves to protect, stabilize, and orient the polymer, mediate its performance, and provide sufficient mechanical and chemical stability to enable integration of conjugated polymers into devices and microsystems.

NASA Discipline Life Support Systems↗

An analysis of energy deposition in a tissue equivalent proportional counter onboard the space shuttle

An improved prediction for space radiations in the lower earth orbits measured by the shuttle TEPC is obtained when energy loss straggling and chord length distribution of the detector are considered. A generalized analytic model is used to describe the energy deposition of direct ion interaction events in a micron-size detector. The transport calculation accounting for the shuttle configuration is accomplished by using a new version of HZETRN that has been extensively verified with laboratory and flight data. The agreement of predicted and measured lineal energy spectra is within 70% for the region above 2 keV/micrometer but within a factor of 2.3 underpredicted for the region below this value. The inclusion of indirect delta ray events in the model is needed before possible causes for the underprediction below 2 keV/micrometer can be assessed.

Flight Experiment↗

Launch Vehicle Propulsion Design with Multiple Selection Criteria

The approach and techniques described herein define an optimization and evaluation approach for a liquid hydrogen/liquid oxygen single-stage-to-orbit system. The method uses Monte Carlo simulations, genetic algorithm solvers, a propulsion thermo-chemical code, power series regression curves for historical data, and statistical models in order to optimize a vehicle system. The system, including parameters for engine chamber pressure, area ratio, and oxidizer/fuel ratio, was modeled and optimized to determine the best design for seven separate design weight and cost cases by varying design and technology parameters. Significant model results show that a 53% increase in Design, Development, Test and Evaluation cost results in a 67% reduction in Gross Liftoff Weight. Other key findings show the sensitivity of propulsion parameters, technology factors, and cost factors and how these parameters differ when cost and weight are optimized separately. Each of the three key propulsion parameters; chamber pressure, area ratio, and oxidizer/fuel ratio, are optimized in the seven design cases and results are plotted to show impacts to engine mass and overall vehicle mass.

Shelton, Joey D.↗

MARS Science Laboratory Post-Landing Location Estimation Using Post2 Trajectory Simulation

The Mars Science Laboratory (MSL) Curiosity rover landed safely on Mars August 5th, 2012 at 10:32 PDT, Earth Received Time. Immediately following touchdown confirmation, best estimates of position were calculated to assist in determining official MSL locations during entry, descent and landing (EDL). Additionally, estimated balance mass impact locations were provided and used to assess how predicted locations compared to actual locations. For MSL, the Program to Optimize Simulated Trajectories II (POST2) was the primary trajectory simulation tool used to predict and assess EDL performance from cruise stage separation through rover touchdown and descent stage impact. This POST2 simulation was used during MSL operations for EDL trajectory analyses in support of maneuver decisions and imaging MSL during EDL. This paper presents the simulation methodology used and results of pre/post-landing MSL location estimates and associated imagery from Mars Reconnaissance Orbiter s (MRO) High Resolution Imaging Science Experiment (HiRISE) camera. To generate these estimates, the MSL POST2 simulation nominal and Monte Carlo data, flight telemetry from onboard navigation, relay orbiter positions from MRO and Mars Odyssey and HiRISE generated digital elevation models (DEM) were utilized. A comparison of predicted rover and balance mass location estimations against actual locations are also presented.

Davis, J. L.↗

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

processing↗

Uncertainty Aware Structural Topology Optimization Via a Stochastic Reduced Order Model Approach

This work presents a stochastic reduced order modeling strategy for the quantification and propagation of uncertainties in topology optimization. Uncertainty aware optimization problems can be computationally complex due to the substantial number of model evaluations that are necessary to accurately quantify and propagate uncertainties. This computational complexity is greatly magnified if a high-fidelity, physics-based numerical model is used for the topology optimization calculations. Stochastic reduced order model (SROM) methods are applied here to effectively 1) alleviate the prohibitive computational cost associated with an uncertainty aware topology optimization problem; and 2) quantify and propagate the inherent uncertainties due to design imperfections. A generic SROM framework that transforms the uncertainty aware, stochastic topology optimization problem into a deterministic optimization problem that relies only on independent calls to a deterministic numerical model is presented. This approach facilitates the use of existing optimization and modeling tools to accurately solve the uncertainty aware topology optimization problems in a fraction of the computational demand required by Monte Carlo methods. Finally, an example in structural topology optimization is presented to demonstrate the effectiveness of the proposed uncertainty aware structural topology optimization approach.

Aguilo, Miguel A.↗

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.↗

Near Rectilinear Halo Orbit Determination with Simulated DSN Observations

This paper presents the results of a high-fidelity simulation of spacecraft orbit determination in a near rectilinear halo orbit (NRHO). Others in the literature have examined this problem with linear covariance analysis, but the highly-nonlinear dynamics of this orbit challenge the assumptions underlying such analyses. The present work builds on similar analysis performed by other authors to contribute a fuller understanding of the operational requirements for NRHO navigation. The present work serves as a check to the assumptions of previous studies and an independent verification of those results. The results from the literature are extended by quantifying the space of orbital states from which a spacecraft with given control authority can safely return to the nominal path. Spacecraft state uncertainty estimates are evaluated as a function of time. Simulated range and range-rate measurements with the Deep Space Network (DSN) ground stations are used to model orbit determination accuracy. Orbit maintenance maneuvers are performed using both short-horizon and long-horizon stationkeeping targeting. Monte Carlo analysis of orbit determination and stationkeeping is performed. This paper quantifies the achievable state uncertainty with deep space network (DSN)-only range and range-rate observations. This paper also addresses requirements on the frequency of DSN observation periods and correlates ground contact frequency with navigation accuracy. The results of several related studies are presented and discussed: the effect of missing ground station passes, the effect of missing stationkeeping maneuvers, the sensitivity of the spacecraft state estimate to realistic error sources, and stationkeeping propellant budget.

Nathan L Parrish↗

Comparing Landsat-7 ETM+ and ASTER Imageries to Estimate Daily Evapotranspiration Within a Mediterranean Vineyard Watershed

We compared the capabilities of Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imageries for mapping daily evapotranspiration (ET) within a Mediterranean vineyard watershed. We used Landsat and ASTER data simultaneously collected on four dates in 2007 and 2008, along with the simplified surface energy balance index (S-SEBI) model. We used previously ground-validated good quality ASTER estimates as reference, and we analyzed the differences with Landsat retrievals in light of the instrumental factors and methodology. Although Landsat and ASTER retrievals of S-SEBI inputs were different, estimates of daily ET from the two imageries were similar. This is ascribed to the S-SEBI spatial differencing in temperature, and opens the path for using historical Landsat time series over vineyards.

Temperature measurement↗

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'propagate-burn-propagate' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons between it and three methods using modified linearized covariance propagation are made. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. In the second method, the maneuver is not modeled but the uncertainty in its magnitude is accounted for by the inclusion of a process noise matrix. In the third method, which is essentially a hybrid of the first two, the nominal portion of the maneuver is included via the state transition matrix while a process noise matrix is used to account for the magnitude uncertainty. Since this method also correctly accounts for the presence of the maneuver in the nominal orbit, it is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, Pc. However, applications for the two other methods exist and are briefly discussed. Despite the fact that the process model ('propagate-burn-propagate') that was studied was very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗

Application of a Monte Carlo accuracy assessment tool to TDRS and GPS

In support of a NASA study on the application of radio interferometry to satellite orbit determination, MITRE developed a simulation tool for assessing interferometric tracking accuracy. Initially, the tool was applied to the problem of determining optimal interferometric station siting for orbit determination of the Tracking and Data Relay Satellite (TDRS). Subsequently, the Orbit Determination Accuracy Estimator (ODAE) was expanded to model the general batch maximum likelihood orbit determination algorithms of the Goddard Trajectory Determination System (GTDS) with measurement types including not only group and phase delay from radio interferometry, but also range, range rate, angular measurements, and satellite-to-satellite measurements. The user of ODAE specifies the statistical properties of error sources, including inherent observable imprecision, atmospheric delays, station location uncertainty, and measurement biases. Upon Monte Carlo simulation of the orbit determination process, ODAE calculates the statistical properties of the error in the satellite state vector and any other parameters for which a solution was obtained in the orbit determination. This paper presents results from ODAE application to two different problems: (1)determination of optimal geometry for interferometirc tracking of TDRS, and (2) expected orbit determination accuracy for Global Positioning System (GPS) tracking of low-earth orbit (LEO) satellites. Conclusions about optimal ground station locations for TDRS orbit determination by radio interferometry are presented, and the feasibility of GPS-based tracking for IRIDIUM, a LEO mobile satellite communications (MOBILSATCOM) system, is demonstrated.

Pavloff, Michael S.↗

Development of Monte Carlo Capability for Orion Parachute Simulations

Parachute test programs employ Monte Carlo simulation techniques to plan testing and make critical decisions related to parachute loads, rate-of-descent, or other parameters. This paper describes the development and use of a MATLAB-based Monte Carlo tool for three parachute drop test simulations currently used by NASA. The Decelerator System Simulation (DSS) is a legacy 6 Degree-of-Freedom (DOF) simulation used to predict parachute loads and descent trajectories. The Decelerator System Simulation Application (DSSA) is a 6-DOF simulation that is well suited for modeling aircraft extraction and descent of pallet-like test vehicles. The Drop Test Vehicle Simulation (DTVSim) is a 2-DOF trajectory simulation that is convenient for quick turn-around analysis tasks. These three tools have significantly different software architectures and do not share common input files or output data structures. Separate Monte Carlo tools were initially developed for each simulation. A recently-developed simulation output structure enables the use of the more sophisticated DSSA Monte Carlo tool with any of the core-simulations. The task of configuring the inputs for the nominal simulation is left to the existing tools. Once the nominal simulation is configured, the Monte Carlo tool perturbs the input set according to dispersion rules created by the analyst. These rules define the statistical distribution and parameters to be applied to each simulation input. Individual dispersed parameters are combined to create a dispersed set of simulation inputs. The Monte Carlo tool repeatedly executes the core-simulation with the dispersed inputs and stores the results for analysis. The analyst may define conditions on one or more output parameters at which to collect data slices. The tool provides a versatile interface for reviewing output of large Monte Carlo data sets while preserving the capability for detailed examination of individual dispersed trajectories. The Monte Carlo tool described in this paper has proven useful in planning several Crew Exploration Vehicle parachute tests.

Moore, James W.↗

Confidence bounds on structural reliability

Different approaches for quantifying physical, statistical, and model uncertainties associated with the distribution parameters which are aimed at determining structural reliability are described. Confidence intervals on the distribution parameters of the input random variables are estimated using four algorithms to evaluate uncertainty of the response. Design intervals are evaluated using either Monte Carlo simulation or an iterative approach. A first order approach can be used to compute a first approximation of the design interval, but its accuracy is not satisfactory. The regression approach which combines the iterative approach with Monte Carlo simulation is capable of providing good results if the performance function can be accurately represented using regression analysis. It is concluded that the design interval-based approach seems to be quite general and takes into account distribution and model uncertainties.

Mehta, S. R.↗

Spacecraft Radiation Analysis

The radiation interface in spacecrafts using radioisotope thermoelectric generators is studied. A Monte Carlo analysis of the radiation field that includes scattered radiation effects, produced neutron and gamma photon isoflux contours as functions of distance from the RTG center line. It is shown that the photon flux is significantly depressed in the RTG axial direction because of selfshielding. Total flux values are determined by converting the uncollided flux values into an equivalent RTG surface source and then performing a Monte Carlo analysis for each specific dose point. Energy distributions of the particle spectra completely define the radiation interface for a spacecraft model.

Harris, D. W.↗

Dissociation rate of bromine diatomics in an argon heat bath

The evolution of a collection of 300 K bromine diatomics embedded in a heat bath of argon atoms at 1800 K was studied by computer, and a dissociation-rate constant for the reaction Br2 + BR + Ar yields Br + Ar was determined. Previously published probability distributions for energy and angular momentum transfers in classical three-dimensional Br2-Ar collisions were used in conjunction with a newly developed Monte Carlo scheme for this purpose. Results are compared with experimental shock-tube data and the predictions of several other theoretical models. A departure from equilibrium is obtained which is significantly greater than that predicted by any of these other theories.

Razner, R.↗

Ablation effects in oxygen-lead fragmentation at 2.1 GeV/nucleon

The mechanism of particle evaporation was used to examine ablation effects in the fragmentation of 2.1 GeV/nucleon oxygen nuclei by lead targets. Following the initial abrasion process, the excited projectile prefragment is assumed to statistically decay in a manner analogous to that of a compound nucleus. The decay probabilities for the various particle emission channels are calculated by using the EVAP-4 Monte Carlo computer program. The input excitation energy spectrum for the prefragment is estimated from the geometric ""clean cut'' abrasion-ablation model. Isotope production cross sections are calculated and compared with experimental data and with the predictions from the standard geometric abrasion-ablation fragmentation model.

Townsend, L. W.↗