Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct Simulation 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Dynamic Trajectory Fit to Multisensor Fireball Observations

Meteorites with known orbital origins are key to our understanding of solar system formation and the source of life on Earth. Fireball networks have been developed globally in a unified effort to record and ultimately retrieve these cosmic samples. However, the accuracy of the determined orbit and the likelihood of meteorite recovery depend directly on the accuracy of the chosen meteoroid triangulation method. There are three leading techniques for meteoroid triangulation discussed in the literature: the method of planes, the straight-line least-squares method, and the multiparameter fit method. Here we describe an alternative method to meteoroid triangulation, called the dynamic trajectory fit. This approach uses the meteoroid’s 3D dynamic equations of motion to fit a realistic trajectory directly to multisensor line-of-sight observations. This method has the ability to resolve fragmentation events, fit systematic observatory timing offsets, and determine mass estimates of the meteoroid along its observable trajectory. Through a comprehensive Monte Carlo analysis of over 100,000 trajectory simulations, we find this new method to more accurately estimate meteoroid trajectories of slow entry events (<25 km s{sup −1}) and events observed from low convergence angles (<10°) compared to existing meteoroid triangulation techniques. Additionally, we triangulate an observed fireball event with visible fragmentation using the various triangulation methods to show that the proposed dynamic trajectory fit implementing fragmentation to best match the captured multisensor line-of-sight data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Avoiding ergodicity problems in lattice discretizations of the Hubbard model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. Here, we show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

First-principles calculations↗
DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗
DOE OSTI2020

Critical Hydrogen Coverage Effect on the Hydrogenation of Ethylene Catalyzed by δ-MoC(001): An Ab Initio Thermodynamic and Kinetic Study

Highly selective oxidation of methane to methanol has long been challenging in catalysis. Here, we reveal key steps for the promotion by water when tuning the selectivity of a well-defined CeO 2 /Cu 2 O/Cu(111) catalyst from CO and CO 2 to methanol under a reaction environment with methane, oxygen, and water. Ambient-pressure x-ray photoelectron spectroscopy showed that water added to CH 4 and O 2 led to surface CH 3 O and accelerated methanol production. These results were consistent with density functional theory calculations and kinetic Monte Carlo simulations which showed that water preferentially dissociates over the active Ce ions at the CeO 2 -Cu 2 O/Cu(111) interface. The adsorbed OH blocked O-O bond cleavage that would dehydrogenate CH 3 O to CO and CO 2 , and it directly converted this species to methanol, while O 2 now reoxidized the reduced surface. Water adsorption also displaced the produced methanol into the gas phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enriching 2D transition metal borides via MB XMenes (M = Fe, Co, Ir): Strong correlation and magnetism

Recently, two-dimensional (2D) FeSe-like anti-MXenes (or XMenes), composed of late d-block transition metal M and p-block nonmetal X elements, have been both experimentally and theoretically investigated. Here, we select three 2D borides FeB, CoB and IrB for a deeper investigation by including strong correlation effects, as a fertile ground for understanding and applications. Using a combination of Hubbard corrected first-principles calculations and Monte Carlo simulations, FeB and CoB are found to be ferro- and anti-ferro magnetic, contrasting with the non-magnetic nature of IrB. The metallic FeB XMene monolayer, superior to most of the MXenes or MBenes, exhibits robust ferromagnetism, driven by intertwined direct-exchange and super-exchange interactions between adjacent Fe atoms. The predicted Curie temperature ($T$ C ) of the FeB monolayer via the Heisenberg model reaches an impressive 425 K, with the easy-axis oriented out-of-plane and high magnetic anisotropic energy (MAE). The asymmetry in the spin-resolved transmission spectrum induces a thermal spin current, providing an opportunity for spin filtration. In conclusion, this novel 2D FeB material is expected to hold great promise as an information storage medium and find applications in emerging spintronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Water-promoted interfacial pathways in methane oxidation to methanol on a CeO 2 -Cu 2 O catalyst

Highly selective oxidation of methane to methanol has long been challenging in catalysis. Here, we reveal key steps for the pro­motion of this reaction by water when tuning the selectivity of a well-defined CeO 2 /Cu 2 O/Cu(111) catalyst from carbon monoxide and carbon dioxide to methanol under a reaction environment with methane, oxygen, and water. Ambient-pressure x-ray photoelectron spectroscopy showed that water added to methane and oxygen led to surface methoxy groups and accelerated methanol production. These results were consistent with density functional theory calculations and kinetic Monte Carlo simulations, which showed that water preferentially dissociates over the active cerium ions at the CeO 2 –Cu 2 O/Cu(111) interface. The adsorbed hydroxyl species blocked O-O bond cleavage that would dehydrogenate methoxy groups to carbon monoxide and carbon dioxide, and it directly converted this species to methanol, while oxygen reoxidized the reduced surface. Water adsorption also displaced the produced methanol into the gas phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Measurements of Nuclear Effects and the $\bar \nu_\mu +H \to \mu^+ + n$ Cross Section in MINERνA with Neutron Tagging

MINERνA, or Main INjector ExpeRiment for ν-A, at Fermilab, is an experimentdedicated to the study of neutrino-nucleus interactions in the GeV regime. Itsgoal is to illustrate the interplay between hadronic and nuclear physics and measure intranuclear dynamics crucial for the present and future neutrino oscillationmeasurements. We first measure a set of variables sensitive to how Monte Carlo(MC) simulations of neutrino-nucleus interactions implement binding energy andthen move on to measure the antineutrino CCQE cross section on the hydrogentargets in MINERνA’s CH detector. We have developed a method to preferentially select events on the hydrogen by comparing outgoing neutrons’ directions totheoretical neutron directions assuming two-body interactions. We measured thecross section, extracted the axial form factor, and performed a Z-expansion fit. Weobserve larger values in the axial form factor at high Q2than current best fits. Finally, we show a preliminary selection of events with both protons and neutrons toinvestigate nuclear processes responsible for producing these final states.

Cai, Tejin↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2019

Data for secondary-electron production from ion precipitation at Jupiter III: Target and projectile processes in H{sup +}, H, and H{sup −} + H{sub 2} collisions

To extend the range of data required for modeling the secondary-electron production from ion precipitation into the upper atmosphere of Jupiter, inelastic processes for collisions of 1 keV to 25 MeV H{sup +}, H, and H{sup −} with H{sub 2} are considered. As in other work treating the dominant heavy-ion species of magnetospheric origin, O and S ions (Schultz et al., 2017, 2019; Gharibnejad et al., 2019) the classical trajectory Monte Carlo method is employed to describe the secondary-electron-producing channels (single and double ionization, transfer ionization, and single and double stripping) as well as the other inelastic channels (single and double charge transfer and projectile and target excitation) required to model the energy loss and charge state evolution of the precipitating ions in their passage through the atmospheric gas. Data is described and tabulated both as directly obtained from these calculations and normalized to widely accepted recommended values (Hunter et al., 1990) for channels for which recommendations exist. As in the previous work, the overall accuracy and completeness of the data presented is verified by use of a Monte Carlo ion-transport simulation to obtain the stopping power and ion-fraction populations as a function of impact energy in comparison with accepted values. The addition of the present data to models of secondary-electron production in Jupiter’s atmosphere improves such model’s ability to interpret in situ observations of the precipitating ions’ effect by the spacecraft Juno as well as enhancing the physical reality of models of the coupling of the Jovian magnetosphere, ionosphere, and atmosphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Methods, systems, and apparatuses for calculating global fluence for neutron and photon monte carlo transport using expected value estimators

Global fluence estimators may be calculated on accelerators and processors for neutron and photon Monte Carlo transport. Monte Carlo random walk simulation may be performed on the processors and the calculation of a Volumetric-Ray-Casting (VRC) estimator may be offloaded to the accelerators. The VRC estimator may modify an expected-value estimator to extend a pseudo-particle ray along the direction of the emitted particle from source and collision event through not only the event volume, but also through all volumes that describe the problem geometry. Additionally, many pseudo-particle rays may be sampled per event, rather than just a single pseudo-particle ray per event, in order to provide more complete angular coverage.

Sweezy, Jeremy Ed↗
DOE OSTI2024

Laboratory realization of relativistic pair-plasma beams

Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black-hole and neutron-star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with electron-positron pairs. Their role in the dynamics of such environments is in many cases believed to be fundamental, but their behavior differs significantly from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies, which are rather limited. We present the first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN’s Super Proton Synchrotron (SPS) accelerator. Monte Carlo simulations agree well with the experimental data and show that the characteristic scales necessary for collective plasma behavior, such as the Debye length and the collisionless skin depth, are exceeded by the measured size of the produced pair beams. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Diagnosing hot-spot symmetry in surrogate ignition experiments via secondary DT-neutron spectroscopy at the NIF

The directional energy spectrum of neutrons generated from the in-flight fusion reaction of 1-MeV tritons contains information about the hot-spot symmetry. The National Ignition Facility (NIF) fields Symmetry Capsule (Symcap) implosions, which have historically measured the symmetry of the radiation, drive by measuring the hot-spot shape via x-ray self-emission. Symcaps are used to tune the hot-spot symmetry for ignition experiments at the NIF. This work shows the relationship between directional secondary DT-n spectra and x-ray imaging data for a large database of Symcap implosions. A correlation is observed between the relative widths of the DT-n spectra measured with nTOFs and the shape measured with x-ray imaging. A Monte Carlo model, which computes the directional secondary DT-n spectrum, is used to interpret the results. A comparison of the x-ray and secondary DT-n data with the Monte Carlo model indicates that 56% of the variance between the two datasets is explained by a P2 asymmetry. More advanced simulations using HYDRA suggest that the unaccounted variance is due to P1 and P4 asymmetries present in the hot spot. The comparison of secondary DT-n data and x-ray imaging data to the modeling shows the DT-n data contain important information that supplements current P2 measurements and contain new information about the P1 asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Neural MUSE Analysis

Researchers at Oak Ridge National Laboratory (ORNL) created data as part of the MUSE (Multi-Agency Urban Search Experiment Detector and Algorithm Test Bed) project simulating illicit nuclear materials located in various buildings along a road. In the simulation, a truck containing a radiation detector drives down the road gathering listmode data (counting the and energy of incident gamma radiation). Building materials, source shielding, driving speed, truck direction, truck location on the road, source type, and source placement are all varied between runs of the data set. This data was created using deterministic neutron transport and Monte Carlo methods through a combination of SCALE, MAVRIC, MCNP, and GADRAS. As part of a follow-on NA-22 project, two Kaggle competitions were created to determine the best algorithms for finding and identifying gamma sources in this simulated urban environment. The winning algorithm was neural network-based and had a test accuracy of 76.4% accuracy for source identification. This work seeks to build upon this work and improve the results through the application of novel machine learning techniques. As a first step, the data was classified by a simple Convolutional Neural Network (CNN) To accomplish this, the data was first preprocessed into “waterfall plots.” These plots are composed of energy vs count plots that are stacked vertically to show progression in time. The horizontal axis indicating the particle energy incorporated user defined bin spacing with options for in linear-, logarithmic-, square root-, and user-spaced bins. The z or color dimension showed the number of counts corresponding the energy-time combination. This data was then used to generate more data, by generating a local estimate of the mean of the distribution for a bin and then randomly re-sampling that bin from a Poisson distribution. Once all of this data was generated, it was fed into a well-known CNN architecture, ResNet50. The output layer of this model was removed and replaced with layers corresponding to the shape desired isotope outputs. The provided training data was used to train the classifier and the remaining testing data was used to evaluate the model. Results are soon to be forthcoming.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Improving Detector Systematic Uncertainties Through Data-Driven Machine Learning

Detector simulation in liquid argon time projection chambers (LArTPCs) is a constant challenge. In particular the modeling of electrons response on wires is highly nontrivial. However, new machine learning techniques exist which can be leveraged to ameliorate these concerns. We present a novel methodology to attempt to learn from cosmic muon data in the ICARUS detector how reconstructed wire signals are influenced by features of the hits such as location, particle direction, angle relative to the wire plane, etc. A model can then be generated which can apply the learned mapping to Monte Carlo events to create a more data-like simulation sample. By creating such a sample we expect to reduce the systematic uncertainties at ICARUS due to our detector modeling.

Hausner, Harry [Fermilab] (ORCID:0000000188932280)↗
DOE OSTI2026

Improving Detector Systematic Uncertainties Through Data-Driven Machine Learning

Detector simulation in liquid argon time projection chambers (LArTPCs) is a constant challenge. In particular the modeling of electrons response on wires is highly nontrivial. However, new machine learning techniques exist which can be leveraged to ameliorate these concerns. We present a novel methodology to attempt to learn from cosmic muon data in the ICARUS detector how reconstructed wire signals are influenced by features of the hits such as location, particle direction, angle relative to the wire plane, etc. A model can then be generated which can apply the learned mapping to Monte Carlo events to create a more data-like simulation sample. By creating such a sample we expect to reduce the systematic uncertainties at ICARUS due to our detector modeling.

Hausner, Harry [Fermilab] (ORCID:0000000188932280)↗
DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Clear-sky detection for PV degradation analysis using multiple regression

A method is presented to detect clear-sky periods for plane-of-array irradiance time-averaged data that is based on the algorithm originally described by Reno and Hansen. Here we show this new method improves the state-of-the-art by providing accurate detection at longer data averaging intervals. Moreover, our new method detects clear periods in plane-of-array data, which is novel. The new method is developed by applying a Design of Experiment approach to optimize the parameters used in the Reno method, and Monte Carlo simulations are used to understand the robustness of the found parameters. Clear-sky detection accuracy is compared among four methods: the Reno method, the default clear-sky filter in RdTools, the Ellis method, and the method outlined in this work, using a hand-labeled two-year data set of 1-min plane-of-array irradiance for a fixed tilt system. The RdTools clear-sky filter is marred by excessive false positives. The other methods all perform well at 1-min data intervals; the method developed here provides more accurate detection at longer data averaging intervals. We show that the parameters are directly linked to the data frequency in the hope that these input variables may not have to be optimized for every data frequency and location. However, only a single fixed system in one location was carefully examined. Finally, we illustrate how accurate determination of clear-sky conditions helps to eliminate data noise and bias in the assessment of long-term performance of PV plants.

14 SOLAR ENERGY↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗