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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Monte Carlo analysis of the Titan III/Transfer Orbit Stage guidance system for the Mars Observer mission

An important part of space launch vehicle mission planning for a planetary mission is the integrated analysis of guidance and performance dispersions for both booster and upper stage vehicles. For the Mars Observer mission, an integrated trajectory analysis was used to maximize the scientific payload and to minimize injection errors by optimizing the energy management of both vehicles. This was accomplished by designing the Titan III booster vehicle to inject into a hyperbolic departure plane, and the Transfer Orbit Stage (TOS) to correct any booster dispersions. An integrated Monte Carlo analysis of the performance and guidance dispersions of both vehicles provided sensitivities, an evaluation of their guidance schemes and an injection error covariance matrix. The polynomial guidance schemes used for the Titan III variable flight azimuth computations and the TOS solid rocket motor ignition time and burn direction derivations accounted for a wide variation of launch times, performance dispersions, and target conditions. The Mars Observer spacecraft was launched on 25 September 1992 on the Titan III/TOS vehicle. The post flight analysis indicated that a near perfect park orbit injection was achieved, followed by a trans-Mars injection with less than 2sigma errors.

Bell, Stephen C.↗

DOE OSTI2023

Fitting a deep generative hadronization model

Hadronization is a critical step in the simulation of high-energy particle and nuclear physics experiments. As there is no first principles understanding of this process, physically-inspired hadronization models have a large number of parameters that are fit to data. Deep generative models are a natural replacement for classical techniques, since they are more flexible and may be able to improve the overall precision. Proof of principle studies have shown how to use neural networks to emulate specific hadronization when trained using the inputs and outputs of classical methods. However, these approaches will not work with data, where we do not have a matching between observed hadrons and partons. In this paper, we develop a protocol for fitting a deep generative hadronization model in a realistic setting, where we only have access to a set of hadrons in data. Our approach uses a variation of a Generative Adversarial Network with a permutation invariant discriminator. We find that this setup is able to match the hadronization model in Herwig with multiple sets of parameters. This work represents a significant step forward in a longer term program to develop, train, and integrate machine learning-based hadronization models into parton shower Monte Carlo programs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Variations of the physical parameters of the blazar Mrk 421 based on analysis of the spectral energy distributions

Abstract We report on the variations of the physical parameters of the jet observed in the blazar Mrk 421, and discuss the origin of X-ray flares in the jet, based on analysis of several spectral energy distributions (SEDs). The SEDs are modeled using the one-zone synchrotron self-Compton model, its parameters determined using a Markov chain Monte Carlo method. The lack of data at TeV energies means many of the parameters cannot be uniquely determined and are correlated. These are studied in detail. We find that the optimal solution can be uniquely determined only when we apply a constraint to one of four parameters: the magnetic field (B), the Doppler factor, the size of the emitting region, and the normalization factor of the electron energy distribution. We used 31 sets of SEDs from 2009 to 2014 with optical–UV data observed with UVOT/Swift and the Kanata telescope, X-ray data with XRT/Swift, and γ-ray data with the Fermi Large Area Telescope. The result of our SED analysis suggests that, in the X-ray faint state, the emission occurs in a relatively small area (∼1016 cm) with a relatively strong magnetic field (B ∼ 10−1 G). The X-ray bright state shows a tendency opposite to that of the faint state, that is, a large emitting area (∼1018 cm), probably downstream of the jet, and a weak magnetic field (B ∼ 10−3 G). The high X-ray flux was due to an increase in the maximum energy of electrons. On the other hand, the presence of two kinds of emitting areas implies that the one-zone model is unsuitable for reproducing at least part of the observed SEDs.

Yamada, Yurika↗

NASA NTRS1997

Probabilistic Micromechanics and Macromechanics for Ceramic Matrix Composites

The properties of ceramic matrix composites (CMC's) are known to display a considerable amount of scatter due to variations in fiber/matrix properties, interphase properties, interphase bonding, amount of matrix voids, and many geometry- or fabrication-related parameters, such as ply thickness and ply orientation. This paper summarizes preliminary studies in which formal probabilistic descriptions of the material-behavior- and fabrication-related parameters were incorporated into micromechanics and macromechanics for CMC'S. In this process two existing methodologies, namely CMC micromechanics and macromechanics analysis and a fast probability integration (FPI) technique are synergistically coupled to obtain the probabilistic composite behavior or response. Preliminary results in the form of cumulative probability distributions and information on the probability sensitivities of the response to primitive variables for a unidirectional silicon carbide/reaction-bonded silicon nitride (SiC/RBSN) CMC are presented. The cumulative distribution functions are computed for composite moduli, thermal expansion coefficients, thermal conductivities, and longitudinal tensile strength at room temperature. The variations in the constituent properties that directly affect these composite properties are accounted for via assumed probabilistic distributions. Collectively, the results show that the present technique provides valuable information about the composite properties and sensitivity factors, which is useful to design or test engineers. Furthermore, the present methodology is computationally more efficient than a standard Monte-Carlo simulation technique; and the agreement between the two solutions is excellent, as shown via select examples.

Murthy, Pappu L. N.↗

DOE OSTI2021

Dynamic life-cycle carbon analysis for fast pyrolysis biofuel produced from pine residues: implications of carbon temporal effects

Abstract Background Woody biomass has been considered as a promising feedstock for biofuel production via thermochemical conversion technologies such as fast pyrolysis. Extensive Life Cycle Assessment studies have been completed to evaluate the carbon intensity of woody biomass-derived biofuels via fast pyrolysis. However, most studies assumed that woody biomass such as forest residues is a carbon–neutral feedstock like annual crops, despite a distinctive timeframe it takes to grow woody biomass. Besides, few studies have investigated the impacts of forest dynamics and the temporal effects of carbon on the overall carbon intensity of woody-derived biofuels. This study addressed such gaps by developing a life-cycle carbon analysis framework integrating dynamic modeling for forest and biorefinery systems with a time-based discounted Global Warming Potential (GWP) method developed in this work. The framework analyzed dynamic carbon and energy flows of a supply chain for biofuel production from pine residues via fast pyrolysis. Results The mean carbon intensity of biofuel given by Monte Carlo simulation across three pine growth cases ranges from 40.8–41.2 g CO 2 e MJ −1 (static method) to 51.0–65.2 g CO 2 e MJ −1 (using the time-based discounted GWP method) when combusting biochar for energy recovery. If biochar is utilized as soil amendment, the carbon intensity reduces to 19.0–19.7 g CO 2 e MJ −1 (static method) and 29.6–43.4 g CO 2 e MJ −1 in the time-based method. Forest growth and yields (controlled by forest management strategies) show more significant impacts on biofuel carbon intensity when the temporal effect of carbon is taken into consideration. Variation in forest operations and management (e.g., energy consumption of thinning and harvesting), on the other hand, has little impact on the biofuel carbon intensity. Conclusions The carbon temporal effect, particularly the time lag of carbon sequestration during pine growth, has direct impacts on the carbon intensity of biofuels produced from pine residues from a stand-level pine growth and management point of view. The carbon implications are also significantly impacted by the assumptions of biochar end-of-life cases and forest management strategies.

09 BIOMASS FUELS↗

NASA NTRS2005

Primal and Dual Integrated Force Methods Used for Stochastic Analysis

At the NASA Glenn Research Center, the primal and dual integrated force methods are being extended for the stochastic analysis of structures. The stochastic simulation can be used to quantify the consequence of scatter in stress and displacement response because of a specified variation in input parameters such as load (mechanical, thermal, and support settling loads), material properties (strength, modulus, density, etc.), and sizing design variables (depth, thickness, etc.). All the parameters are modeled as random variables with given probability distributions, means, and covariances. The stochastic response is formulated through a quadratic perturbation theory, and it is verified through a Monte Carlo simulation.

Patnaik, Surya N.↗

NASA NTRSDate not supplied

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

NASA NTRSDate not supplied

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

microstructure↗

NASA NTRS1997

Probabilistic Mesomechanical Fatigue Model

A probabilistic mesomechanical fatigue life model is proposed to link the microstructural material heterogeneities to the statistical scatter in the macrostructural response. The macrostructure is modeled as an ensemble of microelements. Cracks nucleation within the microelements and grow from the microelements to final fracture. Variations of the microelement properties are defined using statistical parameters. A micromechanical slip band decohesion model is used to determine the crack nucleation life and size. A crack tip opening displacement model is used to determine the small crack growth life and size. Paris law is used to determine the long crack growth life. The models are combined in a Monte Carlo simulation to determine the statistical distribution of total fatigue life for the macrostructure. The modeled response is compared to trends in experimental observations from the literature.

Tryon, Robert G.↗

NASA NTRS2003

The Development of Design Tools for Fault Tolerant Quantum Dot Cellular Automata Based Logic

We are developing software to explore the fault tolerance of quantum dot cellular automata gate architectures in the presence of manufacturing variations and device defects. The Topology Optimization Methodology using Applied Statistics (TOMAS) framework extends the capabilities of the A Quantum Interconnected Network Array Simulator (AQUINAS) by adding front-end and back-end software and creating an environment that integrates all of these components. The front-end tools establish all simulation parameters, configure the simulation system, automate the Monte Carlo generation of simulation files, and execute the simulation of these files. The back-end tools perform automated data parsing, statistical analysis and report generation.

Armstrong, Curtis D.↗

DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗

NASA NTRS2008

DSN Array Simulator

The DSN Array Simulator (wherein 'DSN' signifies NASA's Deep Space Network) is an updated version of software previously denoted the DSN Receive Array Technology Assessment Simulation. This software (see figure) is used for computational modeling of a proposed DSN facility comprising user-defined arrays of antennas and transmitting and receiving equipment for microwave communication with spacecraft on interplanetary missions. The simulation includes variations in spacecraft tracked and communication demand changes for up to several decades of future operation. Such modeling is performed to estimate facility performance, evaluate requirements that govern facility design, and evaluate proposed improvements in hardware and/or software. The updated version of this software affords enhanced capability for characterizing facility performance against user-defined mission sets. The software includes a Monte Carlo simulation component that enables rapid generation of key mission-set metrics (e.g., numbers of links, data rates, and date volumes), and statistical distributions thereof as functions of time. The updated version also offers expanded capability for mixed-asset network modeling--for example, for running scenarios that involve user-definable mixtures of antennas having different diameters (in contradistinction to a fixed number of antennas having the same fixed diameter). The improved version also affords greater simulation fidelity, sufficient for validation by comparison with actual DSN operations and analytically predictable performance metrics.

Tikidjian, Raffi↗

NASA NTRS1973

Reflection of a shock wave from a thermally accommodating wall - Molecular simulation.

Reflection of a plane shock wave from a wall has been simulated on a microscopic scale using a direct simulation Monte Carlo technique of the type developed by Bird. A monatomic gas model representing argon was used to describe the fluid medium and a simple one-parameter accommodation coefficient model was used to describe the gas-surface interaction. The influence of surface accommodation was studied parametrically by varying the accommodation coefficient from zero to one. Results are presented showing the temporal variations of flow field density, and mass, momentum, and energy fluxes to the wall during the shock wave reflection process. The energy flux was used to determine the wall temperature history. Comparisons with experiment are found to be satisfactory where data are available.

Deiwert, G. S.↗

NASA NTRS2007

The UDF05 Follow-up of the HUDF: I. The Faint-End Slope of the Lyman-Break Galaxy Population at zeta approx. 5

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP, and NICMOS (Fll0W, Fl60W) imaging of three fields, two of which coincide with the NICP1-4 NICMOS parallel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity Function (LF) of z approximately 5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based on a V - i, i - z selection criterion, we identify a sample of 101 and 133 candidate z approximately 5 galaxies down to z(sub 850) = 28.5 and 29.25 magnitudes in the NICP12 field and in the HUDF, respectively. Using an extensive set of Monte Carlo simulations we derive corrections for observational biases and selection effects, and construct the rest-frame 1400 Angstroms LBG LF over the range M(sub 1400) = [-22.2, -17.1], i.e. down to approximately 0.04 L(sub *) at z = 5. We show that: (i) Different assumptions for the SED distribution of the LBG population, dust properties and intergalactic absorption result in a 25% variation in the number density of LBGs at z = 5 (ii) Under consistent assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense in z = 5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; (iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical parameters, and has a value of alpha approximately -1.6, similar to the faint-end slope of the LF that has been measured for LBGs at z = 3 and z = 6. Our study therefore supports no variation in the faint-end of the LBG LF over the whole redshift range z = 3 to z = 6. The comparison with theoretical predictions suggests that (a,) the majority of the stars in the z = 5 LBG population are produced with a Top-Heavy IMF in merger-driven starbursts, and that (b) possibly, either the fraction of stellar mass produced in starburst, or the fraction of high mass stars in the bursts is increased towards the bright end of the LF.

Oesch, P. A.↗

NASA NTRS1995

Bidirectional reflectance of oceanic waters: A comparison of modeled and measured upward radiance fields

The bidirectional reflectance of oceanic waters is conveniently described in a normalized way by forming the ratio of the upwelling irradiance E(sub u) to any upwelling radiance L(sub u)(theta prime, phi). This ratio, Q (theta prime, theta(sub 0), (phi(sub 0) - phi), where theta prime, phi are the nadir and azimuth angles for the upward radiance and theta(sub 0), phi(sub 0) are the zenith and azimuth angles of the Sun, has been determined from measurements at sea and computed via Monte Carlo simulations using the inherent optical properties measured in the field and appropriate boundary conditions (clear sky, no wind, varying Sun angle). Experimental ad computed Q values are in excellent agreement. This successful comparison confirms the importance of the bidirectional character of ocean reflectance, already pointed out from a purely numerical approach without field validation, and corroborates the extended range of the Q variations. The later point is of importance when interpreting the marine signals detected by an ocean color satellite-borne sensor. The validation is extended by considering the historical data for the radiance distributions in Lake Pend Oreille determined at various depths. The closure issue in ocean optics is examined by solving the direct problem of radiative transfer and through a model-data comparison in terms of radiance field.

Morel, Andre↗

NASA NTRS2003

Estimation and Simulation of Slow Crack Growth Parameters from Constant Stress Rate Data

Closed form, approximate functions for estimating the variances and degrees-of-freedom associated with the slow crack growth parameters n, D, B, and A(sup *) as measured using constant stress rate ('dynamic fatigue') testing were derived by using propagation of errors. Estimates made with the resulting functions and slow crack growth data for a sapphire window were compared to the results of Monte Carlo simulations. The functions for estimation of the variances of the parameters were derived both with and without logarithmic transformation of the initial slow crack growth equations. The transformation was performed to make the functions both more linear and more normal. Comparison of the Monte Carlo results and the closed form expressions derived with propagation of errors indicated that linearization is not required for good estimates of the variances of parameters n and D by the propagation of errors method. However, good estimates variances of the parameters B and A(sup *) could only be made when the starting slow crack growth equation was transformed and the coefficients of variation of the input parameters were not too large. This was partially a result of the skewered distributions of B and A(sup *). Parametric variation of the input parameters was used to determine an acceptable range for using closed form approximate equations derived from propagation of errors.

Salem, Jonathan A.↗

NASA NTRS1986

Bipropellant rocket exhaust plume analysis on the Galileo spacecraft

This paper describes efforts to quantify the contaminant flow field produced by 10 N thrust bipropellant rocket engines used on the Galileo spacecraft. The prediction of the composition of the rocket exhaust by conventional techniques is found to be inadequate to explain experimental observations of contaminant deposition on moderately cold (200 K) surfaces. It is hypothesized that low volatility contaminants are formed by chemical reactions which occur on the surfaces. The flow field calculations performed using the direct simulation Monte Carlo method give the expected result that the use of line-of-sight plume shields may have very little effect on the flux of vapor phase contaminant species to a surface, especially if the plume shields are located so close to the engine that the interaction of the plume with the shield is in the transition flow regime. It is shown that significant variations in the exhaust plume composition caused by nonequilibrium effects in the flow field lead to very low concentrations of species which have high molecular weights in the more rarefied regions of the flow field. Recommendations for the design of spacecraft plume shields and further work are made.

Guernsey, C. S.↗

NASA NTRS2002

A Simplified, Closed-Form Method for Screening Spacecraft Orbital Heating Variations

A closed-form analytical technique has been developed to screen orbital average heating variations as a function of beta angle, altitude, surface area, and surface optical properties. Using planetary view factor equations for surfaces parallel-to and normal-to the local vertical, a cylindrical umbral shadow approximation, and a simplified albedo flux model, heating rate equations are formulated and then integrated to obtain orbital average heating. The results are compared to detailed analytical predictions using Monte Carlo integration and an assessment of error is presented.

Rickman, S. L.↗