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 325 records · Page 18

Investigation of Magmatic Activities on Early Mars Using Igneous Mineral Chemistry in Gale Crater, Mars

One objective of rover missions is exploring the geological context of the surroundings. Over the years, igneous petrology and sedimentology have been disconnected, the first investigating magmatic processes and volcanic activities, and the second seeking environmental conditions in the past and assessing the habitability of the planet. Although different, one is related to the other: igneous rocks are altered and broken down, leading to the formation of sedimentary rocks, which can in turn be used to back out the nature of their magmatic source. The Curiosity rover that landed in the 3.7 Gyr old impact crater Gale is traveling through sedimentary rocks. About fifty float rocks have been observed, and several of them with ambiguous texture and composition have been classified as igneous or sedimentary depending on studies such as Jake_M. The composition of several unambiguous igneous rocks has been analyzed [4- 6] but their heterogeneity at a larger instrumental (measurement size < 2 cm) scale prevents the measurement of a bulk composition as performed on Earth. An original approach avoiding these two last issues is to consider igneous mineral chemistry analyzed within igneous and sedimentary rocks to assess magmatic processes that could have formed them. Most Curiosity data are used to explore ancient environmental conditions, and a significant number of compositional analyses are under-explored for constraining magmatic activities. We will present how we can make use of sedimentary data for investigating igneous processes in the vicinity of Gale crater. Geological Context: We focus on the first 750 martian days, corresponding to measurements in a coherent lacustrine sedimentary unit called Bradbury, because all sedimentary rocks were sourced from the same watershed and appear to have a consistent source with minimal alteration [2-3]. Igneous detrital minerals including feldspar and pyroxene, are observed in sedimentary rocks. Monte Carlo models showed that minimal cation loss is observed based on the composition of all Bradbury rocks, implying negligible weathering [3]. Although clay minerals are detected in few rocks [7], chemical compositions of rocks can be explained by a mixture of primary igneous minerals [3]. Variation of composition within Bradbury rocks can be explained by mineral sorting and one distinct source component. While a common magmatic source is suggested, Bradbury sediments likely come from several volcanic eruptions from a single magmatic chamber [9- 10]. The occurrence of alkali minerals like sanidine and K-rich rocks throughout Bradbury supports the presence of a potassic component, likely trachytic, while plagioclase and a mafic composition suggest a basaltic component [8-9]. Instruments: Mineral chemistry can be estimated by three instruments onboard Curiosity. The CheMin instrument enables detection of mineral assemblages using X-ray diffraction (XRD). Using Rietveld refinement, each mineral is identified according to their 1D XRD pattern [11]. Note that distinction between pyroxene minerals is challenging with the CheMin instrument due to overlapping peaks on XRD patterns and low angular resolution of the instrument [12]. Then, using least square regression and optimization algorithms based on unit-cell parameters, mineral chemistry has been estimated by [11]. Plagioclase composition has been estimated using the NaAlSi3O8- CaAl2Si2O8 system and alkali feldspar is based on the NaAlSi3O8-KAlSi3O8 system (stars in Fig.1). Two mudstone samples (John Klein and Cumberland) and one sandstone sample (Windjana) were analyzed by CheMin at Yellowknife Bay and Kimberley, respectively. Figure 1. Ternary diagrams of feldspars (top) and pyroxene (bottom) quadrilateral. Stars correspond to CheMin composition and the gray patches to ChemCam composition. The colored dots are the composition of feldspar and pyroxene that crystallized during fractional crystallization at FMQ+1 of a melt extracted at distinct melting degree during the adiabatic ascent of a primitive mantle composition, without any water (left panels) and with 0.5 wt.% of water (right panels) at distinct pressure. The ChemCam instrument enables the analysis of the chemical compositions of rocks at hundreds of micrometer scale (350-550 μm) using laser induced breakdown spectroscopy (LIBS), which may provide the composition of minerals when they are larger than the beam spot (>550 μm) [13]. Within >5000 LIBS points, we performed a typical stoichiometric filtering allowing us to distinguish 56 feldspar and 10 pyroxene mineral compositions (grey patches in Fig. 1). Finally, the Alpha Particle X-ray Spectrometer (APXS) analyzes the composition of rocks with a 1.6 cm diameter spot size. Monte Carlo mass balance modeling allowed [3] to decipher a feldspar range varying between An30 and An40 (Fig. 1). Discussion: Although there could be a more complex history and other ways to form the whole compositional range of igneous minerals analyzed within the Bradbury formation, we are presenting here simple magmatic pathways commonly occurring on Earth using the thermodynamical softwares pMELTS and rhyoliteMELTS [14]. The objective is to find reasonable igneous processes that produce minerals that parallel the compositions of feldspar and pyroxene analyzed by the Curiosity rover. As commonly observed for mid-ocean ridge basalts, the adiabatic ascent of a primitive mantle composition [15] partially melting at 2 GPa has been modeled, followed by the extraction of a liquid at distinct degrees of partial melting, which undergoes fractional crystallization at an oxygen fugacity +1 log unit above the fayalite-magnetite-quartz (FMQ) buffer within the crust (0.02-0.4 GPa) with H2O = 0- 0.5 wt. %. These latter conditions correspond to those recorded within igneous clasts from the Noachian martian breccia NWA 7034 and paired and within Gale igneous rocks (colored dots in Fig. 1) [16-17]. To check the reliability of these 2-step models, we also tested fractional crystallization at similar conditions (FMQ+1; P=0.02-0.4 GPa; H2O = 0-0.5 wt. %) of starting compositions corresponding to that of magmas with distinct melting degrees obtained from isobaric experiments at 2 GPa [18]. Mineral compositions obtained from both models are similar. As shown on Fig. 1, the whole range of observed feldspar compositions cannot be reproduced by fractionation of one magma only. Indeed, while alkali feldspar and Na-plagioclase likely crystallized from fractional crystallization of a low-degree melt (here <15%), plagioclase and pyroxene can only be formed by fractional crystallization of a higher degree melt (here >19%). The corresponding liquid descent lines are broadly in agreement with compositions estimated by ChemCam corresponding to float igneous rocks (Fig. 2) [4-6]. Figure 2. Silica versus alkali content. Lines show the liquid lines of descent from magmas with distinct degrees of melting. Gray patches represents the composition of Gale igneous rocks [4-6]. Trachytic to rhyolitic magmas crystallize alkali feldspar, and andesite to dacite magmas likely form plagioclase. Therefore, at least two starting magmas at distinct melting degrees, which could easily come from a single mantle source, are necessary to explain the whole compositional range of feldspar and pyroxene analyzed within Bradbury rocks. Conclusion: Because rocks from the Bradbury formation are likely originating from the same magmatic source with minimal weathering as supported by several studies using different approaches, igneous mineral chemistry analyzed by CheMin and ChemCam allows us to back out reasonable magmatic pathways that could have crystallized them. Fractional crystallization of at least two starting magmas originating from distinct melting degrees of a single mantle source can explain the whole range of feldspar and pyroxene composition. Both alkaline and sub-alkaline liquids can be produced, with compositions corresponding to those of the igneous rocks analyzed by ChemCam within the Bradbury formation, highlighting the complexity of Mars magmatism.

Payre, V.↗

Probabilistic Assessment and Uncertainty Analysis of CO2 Storage Capacity of the Morrow B Sandstone—Farnsworth Field Unit

This paper presents probabilistic methods to estimate the quantity of carbon dioxide (CO2) that can be stored in a mature oil reservoir and analyzes the uncertainties associated with the estimation. This work uses data from the Farnsworth Field Unit (FWU), Ochiltree County, Texas, which is currently undergoing a tertiary recovery process. The input parameters are determined from seismic, core, and fluid analyses. The results of the estimation of the CO2 storage capacity of the reservoir are presented with both expectation curve and log probability plot. The expectation curve provides a range of possible outcomes such as the P90, P50, and P10. The deterministic value is calculated as the statistical mean of the storage capacity. The coefficient of variation and the uncertainty index, P10/P90, is used to analyze the overall uncertainty of the estimations. A relative impact plot is developed to analyze the sensitivity of the input parameters towards the total uncertainty and compared with Monte Carlo. In comparison to the Monte Carlo method, the results are practically the same. The probabilistic technique presented in this paper can be applied in different geological settings as well as other engineering applications.

03 NATURAL GAS↗

Degree and direction of polarization of multiple scattered light. I - Homogeneous cloud layers.

It is shown that the Monte Carlo method can provide useful information on the polarization and its direction for homogeneous layers corresponding to the haze C, nimbostratus, and ice crystals models. Results for various optical thicknesses and two solar zenith angles show the variation with these parameters. In most cases the direction of polarization with respect to the meridian plane (which contains the outgoing photon) does not change appreciably with the optical thickness of the scattering layer, so that the conclusion that the direction of polarization is determined mainly by the direction of the incident and scattered photons, obtained for single scattering, should apply approximately when multiple scattering is taken into account.

Kattawar, G. W.↗

Development and computer implementation of design/analysis techniques for multilayered composite structures. Probabilistic fiber composite micromechanics

Probabilistic composite micromechanics methods are developed that simulate expected uncertainties in unidirectional fiber composite properties. These methods are in the form of computational procedures using Monte Carlo simulation. The variables in which uncertainties are accounted for include constituent and void volume ratios, constituent elastic properties and strengths, and fiber misalignment. A graphite/epoxy unidirectional composite (ply) is studied to demonstrate fiber composite material property variations induced by random changes expected at the material micro level. Regression results are presented to show the relative correlation between predictor and response variables in the study. These computational procedures make possible a formal description of anticipated random processes at the intraply level, and the related effects of these on composite properties.

Stock, Thomas A.↗

Probabilistic Fiber Composite Micromechanics

Probabilistic composite micromechanics methods are developed that simulate expected uncertainties in unidirectional fiber composite properties. These methods are in the form of computational procedures using Monte Carlo simulation. The variables in which uncertainties are accounted for include constituent and void volume ratios, constituent elastic properties and strengths, and fiber misalignment. A graphite/epoxy unidirectional composite (ply) is studied to demonstrate fiber composite material property variations induced by random changes expected at the material micro level. Regression results are presented to show the relative correlation between predictor and response variables in the study. These computational procedures make possible a formal description of anticipated random processes at the intra-ply level, and the related effects of these on composite properties.

Stock, Thomas A.↗

Quantum Zeno Monte Carlo for computing observables

The recent development of logical quantum processors marks a pivotal transition from the noisy intermediate-scale quantum (NISQ) era to the fault-tolerant quantum computing (FTQC) era. These devices have the potential to address classically challenging problems with polynomial computational time using quantum properties. However, they remain susceptible to noise, necessitating noise resilient algorithms. We introduce Quantum Zeno Monte Carlo (QZMC), a classical-quantum hybrid algorithm that demonstrates resilience to device noise and Trotter errors while showing polynomial computational cost for a gapped system. QZMC computes static and dynamic properties without requiring initial state overlap or variational parameters, offering reduced quantum circuit depth.

Han, Mancheon [Korea Institute for Advanced Study ↗

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a Spectral Model Based on Charge Transport for the Swift/BAT 32K CdZnTe Detector Array

The properties of 32K CdZnTe (4 x 4 sq mm large, 2 mm thick) detectors have been studied in the pre-flight calibration of the Burst Alert Telescope (BAT) on-board the Swift Gamma-ray Burst Explorer (scheduled for launch in November 2004). In order to understand the energy response of the BAT CdZnTe array, we first quantify the mobility-lifetime (mu tau) products of carriers in individual CdZnTe detectors, which produce a position dependency in the charge induction efficiency and results in a low energy tail in the energy spectrum. Based on a new method utilizing (57)Co spectra obtained at different bias voltages, the mu tau for electrons ranges from 5.0 x 10(exp -4) to 1.0 x 10(exp -2) sq cm/V while the mu tau for holes ranges from 1.3 x 10(exp -5 to 1.8 x 10(exp -4) sq cm/V. We find that this wide distribution of mu tau products explains the large diversity in spectral shapes between CdZnTe detectors well. We also find that the variation of mu tau products can be attributed to the difference of crystal ingots or manufacturing harness. We utilize the 32K sets of extracted mu tau products to develop a spectral model of the detector. In combination with Monte Carlo simulations, we can construct a spectral model for any photon energy or any incident angle.

Sato, Goro↗

Some Interesting Applications of Probabilistic Techiques in Structural Dynamic Analysis of Rocket Engines

Numerical and Analytical methods developed to determine damage accumulation in specific engine components when speed variation included. Dither Life Ratio shown to be well over factor of 2 for specific example. Steady‐State assumption shown to be accurate for most turbopump cases, allowing rapid calculation of DLR. If hot‐fire speed data unknown, Monte Carlo method developed that uses speed statistics for similar engines. Application of techniques allow analyst to reduce both uncertainty and excess conservatism. High values of DLR could allow previously unacceptable part to pass HCF criteria without redesign. Given benefit and ease of implementation, recommend that any finite life turbomachine component analysis adopt these techniques. Probability Values calculated, compared, and evaluated for several industry‐proposed methods for combining random and harmonic loads. Two new excel macros written to calculate combined load for any specific probability level. Closed form Curve fits generated for widely used 3(sigma) and 2(sigma) probability levels. For design of lightweight aerospace components, obtaining accurate, reproducible, statistically meaningful answer critical.

Brown, Andrew M.↗

Stochastic robustness of linear control systems

A simple numerical procedure for estimating the stochastic robustness of a linear, time-invariant system is described. Monte Carlo evaluation of the system's eigenvalues allows the probability of instability and the related stochastic root locus to be estimated. This definition of robustness is an alternative to existing deterministic definitions that address both structured and unstructured parameter variations directly. This analysis approach treats not only Gaussian parameter uncertainties but non-Gaussian cases, including uncertain-but-bounded variations. Trivial extensions of the procedure admit alternate discriminants to be considered. Thus, the probabilities that stipulated degrees of instability will be exceeded or that closed-loop roots will leave desirable regions also can be estimated. Results are particularly amenable to graphical presentation.

Stengel, Robert F.↗

Stochastic robustness of linear time-invariant control systems

A simple numerical procedure for estimating the stochastic robustness of a linear, time-invariant system is described. Monte Carlo evaluation of the system's eigenvalues allows the probability of instability and the related stochastic root locus to be estimated. This definition of robustness is an alternative to existing deterministic definitions that address both structured and unstructured parameter variations directly. This analysis approach treats not only Gaussian parameter uncertainties but also non-Gaussian cases, including uncertain-but-bounded variations. Trivial extensions of the procedure permit alternate discriminants to be considered. Thus the probabilities that stipulated degrees of instability will be exceeded or that closed-loop roots will leave desirable regions can also be estimated. Results are particularly amenable to graphical presentation.

Stengel, Robert F.↗

Fatigue Reliability of Gas Turbine Engine Structures

The results of an investigation are described for fatigue reliability in engine structures. The description consists of two parts. Part 1 is for method development. Part 2 is a specific case study. In Part 1, the essential concepts and practical approaches to damage tolerance design in the gas turbine industry are summarized. These have evolved over the years in response to flight safety certification requirements. The effect of Non-Destructive Evaluation (NDE) methods on these methods is also reviewed. Assessment methods based on probabilistic fracture mechanics, with regard to both crack initiation and crack growth, are outlined. Limit state modeling techniques from structural reliability theory are shown to be appropriate for application to this problem, for both individual failure mode and system-level assessment. In Part 2, the results of a case study for the high pressure turbine of a turboprop engine are described. The response surface approach is used to construct a fatigue performance function. This performance function is used with the First Order Reliability Method (FORM) to determine the probability of failure and the sensitivity of the fatigue life to the engine parameters for the first stage disk rim of the two stage turbine. A hybrid combination of regression and Monte Carlo simulation is to use incorporate time dependent random variables. System reliability is used to determine the system probability of failure, and the sensitivity of the system fatigue life to the engine parameters of the high pressure turbine. 'ne variation in the primary hot gas and secondary cooling air, the uncertainty of the complex mission loading, and the scatter in the material data are considered.

Cruse, Thomas A.↗

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗

Cometary Activity Begins at Kuiper Belt Distances: Evidence from C/2017 K2

We study the development of activity in the incoming long-period comet C/2017 K2 over the heliocentric distance range 9 ≲ r {sub H} ≲ 16 au. The comet continues to be characterized by a coma of submillimeter-sized and larger particles ejected at low velocity. In a fixed co-moving volume around the nucleus we find that the scattering cross section of the coma, C, is related to the heliocentric distance by a power law, C∝r{sub H}{sup -s}, with heliocentric index s = 1.14 ± 0.05. This dependence is significantly weaker than the r {sub H} {sup -2} variation of the insolation as a result of two effects. These are, first, the heliocentric dependence of the dust velocity and, second, a lag effect due to very slow-moving particles ejected long before the observations were taken. A Monte Carlo model of the photometry shows that dust production beginning at r {sub H} ~ 35 au is needed to match the measured heliocentric index, with only a slight dependence on the particle size distribution. Mass-loss rates in dust at 10 au are of order 10{sup 3} kg s{sup -1}, while loss rates in gas may be much smaller, depending on the unknown dust to gas ratio. Consequently, the ratio of the nongravitational acceleration to the local solar gravity, α', may, depending on the nucleus size, attain values of ~10{sup -7} ≲ α' ≲ 10{sup -5}, comparable to values found in short-period comets at much smaller distances. Nongravitational acceleration in C/2017 K2 and similarly distant comets, while presently unmeasured, may limit the accuracy with which we can infer the properties of the Oort cloud from the orbits of long-period comets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Physical characteristics of cometary dust from dynamical studies - A review

Progress made in the determination of the physical characteristics of cometary dust particles from studies of dust tail dynamics is reviewed. Applications of the combined dynamical photometric approach of Finson and Probstein (1968) to studies of cometary tails exhibiting continuous light intensity variations are discussed, with attention given to determinations of the particle-size-related distribution function of the solar radiation pressure exerted on the particles, the contribution of comets to the interplanetary dust, calculations of dust ejection rates and a Monte Carlo approach to the analysis of dust tails. Investigations of dust streamers and striae, which are believed to be related to comet outbursts entailing brief but sharp enhancements of dust production, are then reviewed, with particular attention given to observations of Comet West 1976 VI. Finally, the question of cometary particle type is addressed, and it is pointed out that the presence of submicron absorbing particles in the striae of Comet West is not incompatible with the presence of micron-size dielectric particles in the inner coma.

Sekanina, Z.↗

Localized Interactions in Neutrino Simulations

The Deep Underground Neutrino Experiment (DUNE) requires precise modeling of neutrino--nucleus interactions to achieve its neutrino-oscillation measurement goals. GENIE, the Monte Carlo event generator used by DUNE, exhibits a known discrepancy with MicroBooNE measurements of transverse kinematic imbalance (TKI): the data display a larger high-TKI tail while maintaining a peak similar to that predicted by the baseline GENIE model. Previous variations of final-state interaction (FSI) strength affected both the peak and the tail and therefore did not resolve the discrepancy. This work investigates whether inconsistencies between local and global treatments of intranuclear physics contribute to the observed mismodeling. The GENIE FSI routines were modified to use the particle position when generating nucleons, thereby introducing a local-density treatment in the hA and hN intranuclear models for their 2018 and 2025 implementations. Meson-exchange-current (MEC) localization was also tested for the hN 2018 model. Comparisons of the intranuclear scattering-center momentum and its radial dependence confirm that the localization was implemented as intended. Localization produces modest changes in the proton momentum spectrum, primarily at low momentum, but only a small change in the TKI distribution between approximately \SI{0.25}{\giga\electronvolt} and \SI{0.40}{\giga\electronvolt}. These changes are insufficient to account for the discrepancy with MicroBooNE data. Although a localized treatment improves the internal consistency of the GENIE model, the origin of the TKI discrepancy remains unresolved.

Bulla, Braden [Unlisted, US, IL]↗

Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion

Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). Furthermore, the Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental study of shock-wave reflection from a thermally accommodating wall.

Shock-tube experiments have been conducted to study the nonequilibrium gas-surface interaction which occurs when a thick shock wave in argon reflects from a coplanar, heat-conducting wall. Fast-response instrumentation was used to monitor variations in temperature and normal stress on the surface of the shock-tube end wall during and immediately following reflection of the incident shock wave. The laboratory observations are compared with computer predictions obtained by Deiwert using the direct-simulation Monte Carlo method, and excellent agreement is obtained when a suitable average thermal accommodation coefficient is chosen for the wall surface.

Hanson, R. K.↗