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 631 records · Page 35

Nanoparticle dynamics in the spatial afterglows of nonthermal plasma synthesis reactors

Nonthermal plasma flow tube reactors are industrially scalable systems for the production of nanocrystal (NC) based materials and coatings. One key advantage of nonthermal plasma synthesis is the ability to both synthesize NCs and deposit films in a single reactor, as at the reactor outlet, NCs can be inertially deposited onto a target substrate. The size and morphology of deposited particles can substantially influence the film structure and function. Though NCs are typically near-spherical and monodispersed as-produced in plasma synthesis reactors, NC charge and growth dynamics can be altered substantially when NCs are sampled out of the plasma and through the spatial afterglow region, affecting deposition. Experiments have demonstrated changes of NC size and charge in the spatial afterglow; however, these dynamics remain unexplored and unexplained via theory and simulation. To address this, we developed a constant number Monte Carlo (CNMC) simulation model to examine the mechanisms of NC decharging and growth in the spatial afterglow of plasma flow tube reactors. Collisions between NC and plasma species, diffusive deposition, and electron desorption from NCs are incorporated in the CNMC simulation. The simulation results are specifically compared with previous experiments on Si NCs synthesized from a low pressure Ar-SiH 4 nonthermal plasma reactor. Furthermore, the experiment-model comparison shows that CNMC models can be implemented which accurately model NC size distribution evolution in a spatial afterglow. Simultaneously, results show that improved collision models, energetic species diffusion models, and electron desorption models will be necessary to accurately depict NC dynamics in spatial afterglows.

42 ENGINEERING↗

Analysis of IUE observations of hydrogen in comets

The large body of hydrogen Lyman-alpha observations of cometary comae obtained with the International Ultraviolet Explorer satellite has gone generally unanalyzed because of two main modeling complications. First, the inner comae of many bright (gas productive) comets are often optically thick to solar Lyman-alpha radiation. Second, even in the case of a small comet (low gas production) the large IUE aperture is quite small as compared with the immense size of the hydrogen coma, so an accurate model which properly accounts for the spatial distribution of the coma is required to invert the inferred brightnesses to column densities and finally to H atom production rates. Our Monte Carlo particle trajectory model (MPTM), which for the first time provides the realistic full phase space distribution of H atoms throughout the coma was used as the basis for the analysis of IUE observations of the inner coma. The MCPTM includes the effects of the vectorial ejection of the H atoms upon dissociation of their parent species (H2O and OH) and of their partial collisional thermalization. Both of these effects are crucial to characterize the velocity distribution of the H atoms. A new spherical radiative transfer calculation based on our MCPTM was developed to analyze IUE observations of optically thick H comae. The models were applied to observations of comets P/Giacobini-Zinner and P/Halley.

Combi, Michael R.↗

Analysis of IUE Observations of Hydrogen in Comets

The 15-years worth of hydrogen Lyman-alpha observations of cometary comae obtained with the International Ultraviolet Explorer (IUE) satellite had gone generally unanalyzed because of two main modeling complications. First, the inner comae of many bright (gas productive) comets are often optically thick to solar Lyman-alpha radiation. Second, even in the case of a small comet (low gas production) the large IUE aperture is quite small as compared with the immense size of the hydrogen coma, so an accurate model which properly accounts for the spatial distribution of the coma is required to invert the infrared brightnesses to column densities and finally to H atom production rates. Our Monte Carlo particle trajectory model (MCPTM), which for the first time provides the realistic full phase space distribution of H atoms throughout the coma has been used as the basis for the analysis of IUE observations of the inner coma. The MCPTM includes the effects of the vectorial ejection of the H atoms upon dissociation of their parent species (H2O and OH) and of their partial collisional thermalization. Both of these effects are crucial to characterize the velocity distribution of the H atoms. This combination of the MCPTM and spherical radiative transfer code had already been shown to be successful in understanding the moderately optically thick coma of comet P/Giacobini-Zinner and the coma of comet Halley that varied from being slightly to very optically thick. Both of these comets were observed during solar minimum conditions. Solar activity affects both the photochemistry of water and the solar Lyman-alpha radiation flux. The overall plan of this program here was to concentrate on comets observed by IUE at other time during the solar cycle, most importantly during the two solar maxima of 1980 and 1990. Described herein are the work performed and the results obtained.

Combi, Michael R.↗

Modeling Entanglement-Based Quantum Key Distribution for the NASA Quantum Communications Analysis Suite

One of the most practical, and sought after, applications of quantum mechanics in the field of information science is the use of entanglement distribution to communicate quantum information effectively. Similar to the continued improvements of functional quantum computers over the past decade, advances in demonstrations of entanglement distribution over long distances may enable new applications in aeronautics and space communications. The existing NASA Quantum Communications Analysis Suite (NQCAS) software models such applications, but limited experimental data exists to verify the model’s theoretical results. There is, however, a large body of experimental data in the relevant literature for entanglement-based quantum key distribution (QKD). This paper details a Monte Carlo based QKD model that uses NQCAS input parameters to generate an estimated QKD link budget for verification of NQCAS. The model generates link budget statistics like key rates, error rates, and S values that can then be compared to the experimental values in the literature. Preliminary comparisons show many similarities between the simulated and experimental data, supporting the model’s validity. A verified NQCAS model will inform experimental work conducted in Glenn Research Center’s (GRC) NASA Quantum Metrology Laboratory (NQML), supporting the United States Quantum Initiative and potential NASA missions.

NASA Quantum Communications Analysis Suite↗

A Corrosion Risk Assessment Model for Underground Piping

The Pressure Systems Manager at NASA Ames Research Center (ARC) has embarked on a project to collect data and develop risk assessment models to support risk-informed decision making regarding future inspections of underground pipes at ARC. This paper shows progress in one area of this project - a corrosion risk assessment model for the underground high-pressure air distribution piping system at ARC. It consists of a Corrosion Model of pipe-segments, a Pipe Wrap Protection Model; and a Pipe Stress Model for a pipe segment. A Monte Carlo simulation of the combined models provides a distribution of the failure probabilities. Sensitivity study results show that the model uncertainty, or lack of knowledge, is the dominant contributor to the calculated unreliability of the underground piping system. As a result, the Pressure Systems Manager may consider investing resources specifically focused on reducing these uncertainties. Future work includes completing the data collection effort for the existing ground based pressure systems and applying the risk models to risk-based inspection strategies of the underground pipes at ARC.

Datta, Koushik↗

Apprentice for Event Generator Tuning

Apprentice is a tool developed for event generator tuning. It contains a range of conceptual improvements and extensions over the tuning tool Professor. Its core functionality remains the construction of a multivariate analytic surrogate model to computationally expensive Monte-Carlo event generator predictions. The surrogate model is used for numerical optimization in chi-square minimization and likelihood evaluation. Apprentice also introduces algorithms to automate the selection of observable weights to minimize the effect of mis-modeling in the event generators. We illustrate our improvements for the task of MC-generator tuning and limit setting.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation of the Integrated Medical Model Using Historical Space Flight Data

The Integrated Medical Model (IMM) utilizes Monte Carlo methodologies to predict the occurrence of medical events, utilization of resources, and clinical outcomes during space flight. Real-world data may be used to demonstrate the accuracy of the model. For this analysis, IMM predictions were compared to data from historical shuttle missions, not yet included as model source input. Initial goodness of fit test-ing on International Space Station data suggests that the IMM may overestimate the number of occurrences for three of the 83 medical conditions in the model. The IMM did not underestimate the occurrence of any medical condition. Initial comparisons with shuttle data demonstrate the importance of understanding crew preference (i.e., preferred analgesic) for accurately predicting the utilization of re-sources. The initial analysis demonstrates the validity of the IMM for its intended use and highlights areas for improvement.

Kerstman, Eric L.↗

System identification and model reduction using modulating function techniques

Weighted least squares (WLS) and adaptive weighted least squares (AWLS) algorithms are initiated for continuous-time system identification using Fourier type modulating function techniques. Two stochastic signal models are examined using the mean square properties of the stochastic calculus: an equation error signal model with white noise residuals, and a more realistic white measurement noise signal model. The covariance matrices in each model are shown to be banded and sparse, and a joint likelihood cost function is developed which links the real and imaginary parts of the modulated quantities. The superior performance of above algorithms is demonstrated by comparing them with the LS/MFT and popular predicting error method (PEM) through 200 Monte Carlo simulations. A model reduction problem is formulated with the AWLS/MFT algorithm, and comparisons are made via six examples with a variety of model reduction techniques, including the well-known balanced realization method. Here the AWLS/MFT algorithm manifests higher accuracy in almost all cases, and exhibits its unique flexibility and versatility. Armed with this model reduction, the AWLS/MFT algorithm is extended into MIMO transfer function system identification problems. The impact due to the discrepancy in bandwidths and gains among subsystem is explored through five examples. Finally, as a comprehensive application, the stability derivatives of the longitudinal and lateral dynamics of an F-18 aircraft are identified using physical flight data provided by NASA. A pole-constrained SIMO and MIMO AWLS/MFT algorithm is devised and analyzed. Monte Carlo simulations illustrate its high-noise rejecting properties. Utilizing the flight data, comparisons among different MFT algorithms are tabulated and the AWLS is found to be strongly favored in almost all facets.

Shen, Yan↗

Modified Energy Span Analysis Reveals Heterogeneous Catalytic Kinetics

Mechanistic modeling is a cornerstone of catalyst development generally conducted with microkinetic models or density functional theory-based energy profiles. We extend the energy span model of homogeneous catalysis to heterogeneous systems by introducing the modified energy span analysis (MESA) model by implementing collision theory and gas-phase concentration effects. We determine analytically turnover frequencies, coverages, rate-determining steps, apparent activation energies, and reaction orders in agreement with microkinetic and kinetic Monte Carlo simulations. The model applies to discrete catalyst state systems, including single-atom catalysts, homogeneous systems, and microporous materials. A generalizable rate expression is derived for reactor modeling or mechanistic insights from solely experimentally measured reaction orders. Here, we illustrate MESA on three published mechanisms and reveal unexpected phenomena absent in traditional energy span profiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion model for lightning radiative transfer

A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.

Koshak, William J.↗

Monte Carlo simulation of launchsite winds at Kennedy Space Center

This paper develops and validates an easily implemented model for simulating random horizontal wind profiles over the Kennedy Space Center (KSC) at Cape Canaveral, Florida. The model is intended for use in Monte Carlo launch vehicle simulations of the type employed in mission planning, where the large number of profiles needed for statistical fidelity of such simulation experiments makes the use of actual wind measurements impractical. The model is based on measurements made at KSC and represents vertical correlations by a decaying exponential model which is parameterized via least-squares parameter optimization against the sample data. The validity of the model is evaluated by comparing two Monte Carlo simulations of an asymmetric, heavy-lift launch vehicle. In the first simulation, the measured wind profiles are used, while in the second, the wind profiles are generated using the stochastic model. The simulations indicate that the use of either the measured or simulated wind field results in similar launch vehicle performance.

Queen, Eric M.↗

Light Curves and Spectra from Kilonova Models

Kilonovae are the radioactively powered optical/IR transients following a binary neutron star or neutron star black hole merger. A broad grid (varying ejecta mass, velocity, composition) of 900 2D kilonova models were simulated with Monte Carlo on Grizzly, using MPI+OpenMP. The set of models are being compared to observation AT2017gfo (the electromagnetic counterpart of gravitational wave GW170817), and hint at potential fidelity improvements.

79 ASTRONOMY AND ASTROPHYSICS↗

Vectorized Monte Carlo methods for reactor lattice analysis

Some of the new computational methods and equivalent mathematical representations of physics models used in the MCV code, a vectorized continuous-enery Monte Carlo code for use on the CYBER-205 computer are discussed. While the principal application of MCV is the neutronics analysis of repeating reactor lattices, the new methods used in MCV should be generally useful for vectorizing Monte Carlo for other applications. For background, a brief overview of the vector processing features of the CYBER-205 is included, followed by a discussion of the fundamentals of Monte Carlo vectorization. The physics models used in the MCV vectorized Monte Carlo code are then summarized. The new methods used in scattering analysis are presented along with details of several key, highly specialized computational routines. Finally, speedups relative to CDC-7600 scalar Monte Carlo are discussed.

Brown, F. B.↗

Matrix-Model Simulations Using Quantum Computing, Deep Learning, and Lattice Monte Carlo

Matrix quantum mechanics plays various important roles in theoretical physics, such as a holographic description of quantum black holes, and it underpins the only practical numerical approach to the study of complex high-dimensional supergravity theories. Understanding quantum black holes and the role of entanglement in a holographic setup is of paramount importance for the realization of a quantum theory of gravity. Moreover, a complete numerical understanding of the holographic duality and the emergence of geometric space-time features from microscopic degrees of freedom could pave the way for new discoveries in quantum information science. Euclidean lattice Monte Carlo simulations are the de facto numerical tool for understanding the spectrum of large matrix models and have been used to test the holographic duality. However, they are not tailored to extract dynamical properties or even the quantum wave function of the ground state of matrix models. Quantum computing and deep learning provide potentially useful approaches to study the dynamics of matrix quantum mechanics. If successful in the context of matrix models, these rapidly improving numerical techniques could become the new Swiss army knife of quantum gravity practitioners. In this paper, we perform the first systematic survey for quantum computing and deep-learning approaches to matrix quantum mechanics, comparing them to lattice Monte Carlo simulations. These provide baseline benchmarks before addressing more complicated problems. In particular, we test the performance of each method by calculating the low-energy spectrum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First Measurement of Differential Charged Current Quasielastic-like $\nu_\mu$–Argon Scattering Cross Sections In Kinematic Imbalance Variables With The MicroBooNE Detector

We report the first measurement of flux-integrated multi-differential cross sections for charged-current events with muon neutrinos scattering on argon with solely a muon and a single proton in the final state as a function of kinematic imbalance variables. The measurement was carried out using the Booster Neutrino Beam at Fermi National Accelerator Laboratory within the MicroBooNE Liquid Argon Time Projection Chamber detector with an exposure of 6.79 × 10 20 protons on target. Events were selected to enhance the contribution of charged-current mesonless interactions with one proton detected in the final state. The data discussed here are reported in terms of multidifferential cross sections in kinematic imbalance variables, which are generally sensitive to nuclear effects. The double-differential results in these variables can provide an excellent handle to disentagle specific nuclear aspects not easily isolated via single differential cross sections. Our results pave a path towards identifying regions of the phase-space where future interaction modeling development and Monte Carlo neutrino generator tuning efforts should concentrate. EVENT SELECTION Neutrino oscillation measurements aim to extract neutrino mixing angles, mass differences, the value of the chargeparity violating phase in the lepton sector, and to search for new physics beyond the Standard Model [1, 2]. For that to be achieved, an unprecedented understanding of neutrino-argon interactions is of utmost importance since a growing number of neutrino oscillation experiments employ Liquid Argon Time Projector Chamber (LArTPC) neutrino detectors [3–6]. The accuracy to which these experiments can extract neutrino oscillation parameters requires a good understanding of the neutrino energy. Experimentally, this energy is deduced from the measured kinetic energies of particles that are emitted following the neutrino interaction in the detector. The kinematic properties of such finalstate particles reflect complex dynamics due to nuclear and initial-state effects of the interaction [7]. However, certain categories of nuclear effects can be isolated by variables built specifically to characterize the degeneracy between such effects [8–10]. This note reports cross sections in kinematic variables sensitive to nuclear effects using events with one detected muon with momentum 0.1 < pµ < 1.2 GeV/c, and exactly one proton with 0.3 < pp < 1 GeV/c. This signal definition includes events with any number of protons below 300 MeV/c, neutrons at any momenta, and charged pions with momentum lower than 70 MeV/c. This choice is guided by the fact that their experimental signature of correlated muon-proton pairs is fairly straightforward to reconstruct [11–22]. Such events primarily originate from chargedcurrent (CC) neutrino-nucleon quasielastic (QE) scattering interactions where the neutrino removes a single intact nucleon from the nucleus without producing any additional particles. This definition can also include contributions from interactions that lead to the production of additional particles that are absent from the final state due to nuclear effects, such as pion absorption, or have momenta that are below the experimental detection threshold.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Microcanonical Hamiltonian Monte Carlo

We develop Microcanonical Hamiltonian Monte Carlo (MCHMC), a class of models that follow fixed energy Hamiltonian dynamics, in contrast to Hamiltonian Monte Carlo (HMC), which follows canonical distribution with different energy levels. MCHMC tunes the Hamiltonian function such that the marginal of the uniform distribution on the constant-energy-surface over the momentum variables gives the desired target distribution. We show that MCHMC requires occasional energy-conserving billiard-like momentum bounces for ergodicity, analogous to momentum resampling in HMC. We generalize the concept of bounces to a continuous version with partial direction preserving bounces at every step, which gives energy-conserving underdamped Langevin-like dynamics with non-Gaussian noise (MCLMC). MCHMC and MCLMC exhibit favorable scalings with condition number and dimensionality. We develop an efficient hyperparameter tuning scheme that achieves high performance and consistently outperforms NUTS HMC on several standard benchmark problems, in some cases by orders of magnitude.

97 MATHEMATICS AND COMPUTING↗

Ice Cloud Backscatter Study and Comparison with CALIPSO and MODIS Satellite Data

An invariant imbedding T-matrix (II-TM) method is used to calculate the single-scattering properties of 8-column aggregate ice crystals. The II-TM based backscatter values are compared with those calculated by the improved geometric-optics method (IGOM) to refine the backscattering properties of the ice cloud radiative model used in the MODIS Collection 6 cloud optical property product. The integrated attenuated backscatter-to-cloud optical depth (IAB-ICOD) relation is derived from simulations using a CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite) lidar simulator based on a Monte Carlo radiative transfer model. By comparing the simulation results and co-located CALIPSO and MODIS (Moderate Resolution Imaging Spectroradiometer) observations, the non-uniform zonal distribution of ice clouds over ocean is characterized in terms of a mixture of smooth and rough ice particles. The percentage of the smooth particles is approximately 6 percent and 9 percent for tropical and mid-latitude ice clouds, respectively.

MODIS↗