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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spin-orbital liquids and insulator-metal transitions on the pyrochlore lattice

The two orbital Hubbard model, with the electrons additionally coupled to a complex magnetic background, arises in the pyrochlore molybdates. The background involves local moments Hund's coupled to the electrons, driving double exchange ferromagnetism, and antiferromagnetic (AF) tendency arising from competing superexchange. The key scales include the Hubbard repulsion and the superexchange, both of which can be tuned in these materials. They control the phase transition from a ferromagnetic metal to a spin glass metal and then to a spin glass (Mott) insulator. We provide a comprehensive description of the ground state of this model using an unrestricted Hartree-Fock scheme implemented via a simulated annealing procedure and establish the metal-insulator transition line for varying Hubbard interaction and superexchange. The electrons see an effective disorder, due to orbital frustration, already in the ferromagnetic phase. The disorder is further enhanced by antiferromagnetic coupling and the resulting magnetic disorder. As a result, increasing AF coupling shifts the metal-insulator transition to lower Hubbard interaction and gives it an additional "Anderson" character. In conclusion, we provide detailed results on the magnetic and orbital correlations, the density of states, and the optical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Determination of the argon spectral function from ( e , e ' p ) data

The E12-14-012 experiment, performed in Jefferson Lab Hall A, has measured the $(e, e'p)$ cross section in parallel kinematics using a natural argon target. Here, we report the full results of the analysis of the data set corresponding to beam energy 2.222 GeV, and spanning the missing momentum and missing energy range $15 \lesssim p_m \lesssim 300$ MeV/c and $12 \lesssim E_m \lesssim 80$ MeV. The reduced cross section, determined as a function of $p_m$ and $E_m$ with $\approx$4\% accuracy, has been fitted using the results of Monte Carlo simulations involving a model spectral function and including the effects of final state interactions. The overall agreement between data and simulations turns out to be quite satisfactory ($\chi^2$/n.d.o.f.=1.9). Furthermore, the resulting spectral function will provide valuable new information, needed for the interpretation of neutrino interactions in liquid argon detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Finite-temperature phase diagram of the Berenstein-Maldacena-Nastase matrix model on the lattice

We investigate the thermal phase structure of the Berenstein-Maldacena-Nastase matrix model using nonperturbative lattice Monte Carlo calculations. Our main analyses span 3 orders of magnitude in the coupling, involving systems with sizes up to 𝑁𝜏 =24 lattice sites and SU⁡(𝑁) gauge groups with 8 ≤𝑁 ≤16. In addition, we carry out extended checks of discretization artifacts for 𝑁𝜏 ≤128 and gauge group SU(4). We find results for the deconfinement temperature that interpolate between the perturbative prediction at weak coupling and the large-𝑁 dual supergravity calculation at strong coupling. While we confirm that the phase transition is first order for strong coupling, it appears to be continuous for weaker couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Colloquium: Eigenvector continuation and projection-based emulators

Eigenvector continuation is a computational method for parametric eigenvalue problems that uses subspace projection with a basis derived from eigenvector snapshots from different parameter sets. It is part of a broader class of subspace-projection techniques called reduced-basis methods. In this Colloquium, the development, theory, and applications of eigenvector continuation and projection-based emulators are presented. In conclusion, the basic concepts are introduced, the underlying theory and convergence properties are discussed, and recent applications for quantum systems and future prospects are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

HEPfit: a code for the combination of indirect and direct constraints on high energy physics models

HEPfit is a flexible open-source tool which, given the Standard Model or any of its extensions, allows to (i) fit the model parameters to a given set of experimental observables; (ii) obtain predictions for observables. HEPfit can be used either in Monte Carlo mode, to perform a Bayesian Markov Chain Monte Carlo analysis of a given model, or as a library, to obtain predictions of observables for a given point in the parameter space of the model, allowing HEPfit to be used in any statistical framework. In the present version, around a thousand observables have been implemented in the Standard Model and in several new physics scenarios. In this paper, we describe the general structure of the code as well as models and observables implemented in the current release.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Measurement of the Lund jet plane in hadronic decays of top quarks and W bosons with the ATLAS detector

The Lund jet plane (LJP) is measured for the first time in $t\bar{t}$ events, using 140 $\textrm{fb}^{-1}$ of $\sqrt{s} = 13$ TeV pp collision data collected with the ATLAS detector at the LHC. The LJP is a two-dimensional observable of the sub-structure of hadronic jets that acts as a proxy for the kinematics of parton showers and hadron formation. The observable is constructed from charged particles and is measured for R = 1.0 anti-$k_t$ jets with transverse momentum above 350 GeV containing the full decay products of either a top quark or a daughter W boson. The other top quark in the event is identified from its decay into a b-quark, an electron or a muon and a neutrino. The measurement is corrected for detector effects and compared with a range of Monte Carlo predictions sensitive to different aspects of the hadronic decays of the heavy particles. In the W-boson-initiated jets, all the predictions are incompatible with the measurement. In the top quark initiated jets, disagreement with all predictions is observed in smaller subregions of the plane, and with a subset of the predictions across the fiducial plane. The measurement could be used to improve the tuning of Monte Carlo generators, for better modelling of hadronic decays of heavy quarks and bosons, or to improve the performance of jet taggers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Validation of the SCEPTRE Boltzmann-CSD Solver

A new Boltzmann-CSD solver has been developed within the SCEPTRE radiation-transport code, based on the 1st-order form of the transport equation, using discontinuous finite elements in space and energy and discrete ordinates in angle. The Boltzmann-CSD solver has been validated against experimental data for electron energy deposition distributions and for electron emission spectra. Comparison of the calculated results with experimental data shows excellent agreement for many of the test configurations and reasonable agreement for other test configurations. The tests have also been modeled with the ITS Monte Carlo code, which also shows excellent to reasonable agreement with the SCEPTRE results and experimental data. The SCEPTRE Boltzmann-CSD solver relies on electron cross sections generated by the legacy CEPXS code, which currently is limited to electron-only Boltzmann-CSD cross sections. Performing full electron-photon radiation transport with the Boltzmann-CSD solver will require further development in the cross section generating code. For the energy- deposition calculations, neglecting photon transport results in at most about 5% overprediction of the energy deposition for high-energy electrons on high-Z targets, and relatively insignificant difference for the other test configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Computational Modeling and Simulation for Nonproliferation: The History (and Present) of Monte Carlo and the MCNP(R) Code at Los Alamos [Slides]

The emergence of the Monte Carlo method as a research tool springs from work done at Los Alamos in the 1940s.Monte Carlo and the MCNP Code have been and continue to be developed at Los Alamos for many decades. From basic science in support of understanding nuclear interaction physics to global security applications, application uses of the code are extensive. Recent R&D projects and code modernization efforts make the MCNP code a great tool for nuclear nonproliferation applications. In collaboration with nuclear safeguards experts, new training has just recently been developed to help new practitioners learn how to use the code for nuclear safeguards applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Optimizing the Heisenberg Vortex Tube for Hydrogen Cooling (Final Technical Report)

Hydrogen utilization at Plug Power sites ranges between 70-93% with the remainder vented and lost to the atmosphere. The goal of this project was to improve hydrogen utilization at Plug Power fulfillment centers via the patented Heisenberg vortex tube (HVT). The HVT combines conventional counterflow vortex tubes with para-orthohydrogen conversion catalyst to allow the cooling of hydrogen flows with no moving parts or external energy inputs. In year 1 of this project the HVT was evaluated to improve hydrogen utilization via the following concepts: 1) increasing liquid hydrogen pump volumetric efficiency by 20% through vapor separation and subcooling of the liquid, 2) reducing liquid hydrogen storage tank boil-off losses by 20% through thermal vapor shielding (TVS), and 3) increase isentropic efficiency of supercritical hydrogen expansion at 40-50 K, from 31% to over 40%. The TVS concept was selected as the most favorable application for further development. While the conceptual analysis was underway, the cryocatalysis hydrogen experiment facility (CHEF) was retrofitted with a new cryocooler, higher pressure condenser tanks, and in-situ fiber-optic Raman probes for ortho-parahydrogen composition analysis at the inlets and outlets of the HVT. To our knowledge this is the first in-situ implementation of cryogenic Raman probes for ortho-parahydrogen analysis. Subsequent testing of catalyzed and non-catalyzed HVT identified a low flow rate, high-conversion efficacy regime well suited for TVS development. Computational fluid dynamics (CFD), Reduced order modeling (ROM), and quantum Monte Carlo atomistic simulations were applied to optimize the design of the HVT implemented for TVS. The simulations and analysis identified a ruthenium-based catalyst as the most optimal for high conversion with little pressure losses and was matched to experimental measurements. These results indicated the design could achieve the 20% boil-off reduction target. An HVT field-trial was designed and constructed by Plug Power and implemented at a customer site. However, budget limitations reduced the amount of catalyst that could be applied to the HVT. With the reduced amount of catalyst, CFD analysis estimated about 3% reduction in boil-off. A similar amount of boil-off reduction was observed during preliminary measurements during tank commissioning. However, the results have a large margin for error and were operated at relatively low tank temperatures when the HVT has the lowest potential for improvement. This project demonstrated the use of para-orthohydrogen catalysis for reducing liquid hydrogen tank boil-off losses. The field trial system is anticipated to provide long-term experimental measurements on tank operational performance that will allow Plug Power to instrument additional tanks in the future. Recommendations for future work include the development of higher activity, lower cost para-orthohydrogen conversion catalyst for use in the HVT. Tank process optimizations could maximize the efficacy of the HVT and further improve hydrogen utilization.

08 HYDROGEN↗
DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Origin of Intense Electron Heating in Relativistic Blast Waves

The modeling of gamma-ray burst afterglow emission bears witness to strong electron heating in the precursor of Weibel-mediated, relativistic collisionless shock waves propagating in unmagnetized electron–ion plasmas. In this Letter, we propose a theoretical model, which describes electron heating via a Joule-like process caused by pitch-angle scattering in the decelerating, self-induced microturbulence and the coherent charge-separation field induced by the difference in inertia between electrons and ions. The emergence of this electric field across the precursor of electron–ion shocks is confirmed by large-scale particle-in-cell (PIC) simulations. Integrating the model using a Monte Carlo-Poisson method, we compare the main observables to the PIC simulations to conclude that the above mechanism can indeed account for the bulk of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Predictive Transient Model of the TREAT-SIRIUS Experiments

To support experiments modeling and transient analysis of the NASA-sponsored SIRIUS experiments in the Transient Reactor Test Facility (TREAT), an innovative predictive transient model has been developed to simulate various transient experiments. The developed model is based on a new approach that utilizes steady-state Monte Carlo calculations along with a surrogate model based on polynomial regression to determine the core reactivity. The reactivity is supplied to an exact point kinetics model that determines the reactor power, and it is coupled to an adiabatic feedback model to determine the fuel average temperature. In this work, the newly developed model is introduced along with validation test results considering SIRIUS-1 experiments at different power levels. The initial test results show a very good agreement with experimental data at different power levels, which indicates that the model can be used to predict the power of other transient tests.

42 ENGINEERING↗
DOE OSTI2022

Sign-free determinant quantum Monte Carlo study of excitonic density orders in a two-orbital Hubbard-Kanamori model

While excitonic instabilities in multiorbital systems recently have come under scrutiny in a variety of transition-metal compounds, understanding emergence of these instabilities from strong electronic interactions has remained a challenge. In this work, we present a sign-problem-free determinant quantum Monte Carlo study of excitonic density orders in a half-filled two-orbital Hubbard-Kanamori model with broken orbital degeneracy, which accounts for the role of Hund's coupling in transition-metal compounds. For strong inverted (negative) Hund's exchange, we find numerical evidence for the emergence of excitonic density order, with competition between anti-ferro-orbital order and Q = (π, π) excitonic density order as a function of orbital splitting and Hund's coupling. While inverted Hund's coupling stabilizes a spin-singlet excitonic density phase for weak orbital splitting, positive Hund's coupling favors a spin-triplet excitonic density phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗