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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Probability Density Function Control for Stochastic Nonlinear Systems using Monte Carlo Simulation

This paper presents a fast implementation framework of output probability density function (PDF) control for a class of stochastic nonlinear systems which are subjected to non-Gaussian noises. The statistical properties of the system outputs can be adjusted by shaping the dynamic output probability density function to track the reference stochastic distribution. However, the dynamic probability density function evolution is very difficult to obtain analytically even if the system model and the stochastic distributions of the noises are known. Motivated by Monte Carlo simulation, the dynamic probability density function can be estimated by sampling data which forms the contribution of this paper. In particular, the sampling points are generated following the stochastic distribution of the noise for each instant. These points go through the system would generate the histogram for system outputs, then the dynamic model can be established based on the dynamic histogram which reflects the randomness and the nonlinear dynamics of the investigated system. Based on the established model, the output probability density function tracking can be achieved and the simulation results and discussions show the effectiveness and benefits of the presented framework.

Zhang, Qichun↗
DOE OSTI2021

Multi-Fidelity Active Subspaces for Wind Farm Uncertainty Quantification

Wind plants operate in stochastic environments characterized by complex turbulent flow dynamics and high-dimensional random variables. A key step in uncertainty quantification studies is sensitivity analysis and dimension reduction that can facilitate the development of surrogate models to be used for forward and inverse propagation or optimization under uncertainty. Prior work has shown active subspaces are an effective tool for identifying important directions in the space of stochastic inputs; however, they have only been applied to single-fidelity wind plant models. In this study, we investigate the efficacy of a multi-fidelity active subspace method for analyzing the uncertainty in wind plant power output. The multi-fidelity active subspace estimator offers the promise of increased accuracy in identifying active subspaces as compared to a single-fidelity estimator for the same computational cost, or a reduction in cost for the same accuracy. This makes the study of uncertainty in larger wind plants and with higher fidelity physics tractable. The multi-fidelity active subspace method is applied to gridded and existing wind plant layouts with single and multiple inflow conditions and its performance for surrogate modeling and uncertainty propagation is compared against a single-fidelity active subspace method. This multi-fidelity approach yields substantial computational speedups of 2x - 3.4x across the test cases along with acceptable accuracy in surrogate modeling and computing statistical moments.

active subspace↗
DOE OSTI2020

Validation of MCNP Critical Benchmark Models of PU-MET-FAST-016

A new centralized repository of high-quality MCNP models of critical benchmark experiments is currently under development at Los Alamos National Laboratory (LANL). The repository is the combined effort of the Nuclear Criticality Safety, Nuclear Data, and Monte Carlo code development/application organizations at LANL. The initial set of benchmark models in the repository are derived from the Whisper Suite provided with MCNP6.2 and are evaluated against the benchmark experiments described in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The goal is to build a single LANL benchmark collection that is up to date with the latest ICSBEP revision, has a formal review and revision process, is contained in an open-source repository, and utilizes new Python tools for improved input and output file review. This paper describes the validation of the models associated with PU-MET-FAST-016, “Flooded 3X3X3 Arrays of 3-kg Plutonium (Pu) Metal Cylinders – Phase I”. In this critical experiment, twenty-seven 3-kg Pu metal cylinders with a height of 46.33 ± 0.15 mm and diameter of 65.25 ± 0.05 mm were arranged in a cubic array. The Pu cylinders were sealed in aluminum cans, and three cans were placed in each of nine perforated aluminum sleeves, as shown in Fig. 1. The nine sleeves were supported by an aluminum frame that allowed for control of the lateral spacing of the cylinders, both in the X and Y directions. The benchmark experiment measured thirteen configurations. Of these configurations, six were modeled using MCNP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Consequences of adsorbate-adsorbate interactions for apparent kinetics of surface catalytic reactions

Lateral adsorbate interactions at catalyst surfaces are known to influence adsorption energies and reaction rates. Lattice-based kinetic Monte Carlo (kMC) simulations are able to capture these influences, but such models are typically parameterized for a specific reaction network and catalyst surface. Here we report kMC simulations to probe the influence of lateral interactions on simulated rates, rate orders, apparent activation energies, and Sabatier plots. We construct a simple, two-step reaction network involving a single adsorbate and rate- limiting diatomic dissociation, employ a generic repulsive lateral interaction model consistent with known adsorbate-adsorbate interactions on metal surfaces, and a rate model consistent with known Brønsted-Evans-Polyani relationships for diatomic dissocations. Furthermore, we juxtapose reaction kinetics over a wide range of reaction conditions and catalyst binding energies, as a function of interaction strength. We find that at a given zero-coverage binding energy and external conditions, adsorbate coverage decreases monotonically with increasing interaction strength, but absolute rates can vary linearly or nonlinearly. Interactions flatten the Sabatier volcano and shift the maximum towards stronger binding. Influences on apparent rate orders and activation energies are modest and are sensitive to interaction-induced adsorbate ordering. Model predictions are sensitive to lattice size effects at high coverages. The results, modelled for a simple reaction system, highlight the generic consequences of lateral interactions and guidance for identifying their signatures in observed kinetics.

kinetic Monte Carlo↗
DOE OSTI2021

Measurement of b -quark fragmentation properties in jets using the decay $B^±$ → $J/ψK^±$ in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

The fragmentation properties of jets containing b-hadrons are studied using charged B mesons in 139 fb -1 of pp collisions at $\sqrt{s}$ = 13 TeV, recorded with the ATLAS detector at the LHC during the period from 2015 to 2018. The B mesons are reconstructed using the decay of B ± into J/ψK ± , with the J/ψ decaying into a pair of muons. Jets are reconstructed using the anti-k t algorithm with radius parameter R = 0.4. The measurement determines the longitudinal and transverse momentum profiles of the reconstructed B hadrons with respect to the axes of the jets to which they are geometrically associated. These distributions are measured in intervals of the jet transverse momentum, ranging from 50 GeV to above 100 GeV. The results are corrected for detector effects and compared with several Monte Carlo predictions using different parton shower and hadronisation models. The results for the longitudinal and transverse profiles provide useful inputs to improve the description of heavy-flavour fragmentation in jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The changing nature of hydroclimatic risks across South Africa

We present results from large ensembles of projected twenty-first century changes in seasonal precipitation and near-surface air temperature for the nation of South Africa. These ensembles are a result of combining Monte Carlo projections from a human-Earth system model of intermediate complexity with pattern-scaled responses from climate models of the Coupled Model Intercomparison Project Phase 5 (CMIP5). These future ensemble scenarios consider a range of global actions to abate emissions through the twenty-first century. We evaluate distributions of surface-air temperature and precipitation change over three sub-national regions: western, central, and eastern South Africa. In all regions, we find that without any emissions or climate targets in place, there is a greater than 50% likelihood that mid-century temperatures will increase threefold over the current climate’s two-standard deviation range of variability. However, scenarios that consider more aggressive climate targets all but eliminate the risk of these salient temperature increases. A preponderance of risk toward decreased precipitation (3 to 4 times higher than increased) exists for western and central South Africa. Strong climate targets abate evolving regional hydroclimatic risks. Under a target to limit global climate warming to 1.5 °C by 2100, the risk of precipitation changes within South Africa toward the end of this century (2065–2074) is commensurate to the risk during the 2030s without any global climate target. Thus, these regional hydroclimate risks over South Africa could be delayed by 30 years and, in doing so, provide invaluable lead-time for national efforts to prepare, fortify, and/or adapt.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Spherical time-encoded radiation imaging simulations

Radiation source localization is important for nuclear nonproliferation and can be obtained using time-encoded imaging systems with unsegmented detectors. A scintillation crystal can be used with a moving coded-aperture mask to vary the detected count rate produced from radiation sources in the far field. The modulation of observed counts over time can be used to reconstruct an image with the known coded-aperture mask pattern. Current time-encoded imaging systems incorporate cylindrical coded-aperture masks and have limits to their fully coded imaging field-of-view. This work focuses on expanding the field-of-view to 4π by using a novel spherical coded-aperture mask. A regular icosahedron is used to approximate a spherical mask. This icosahedron consists of 20 equilateral triangles; the faces of which are each subdivided into four equilateral triangle-shaped voxels which are then projected onto a spherical surface, creating an 80-voxel coded-aperture mask. Furthermore, these polygonal voxels can be made from high-Z materials for gamma-ray modulation and/or low-Z materials for neutron modulation. In this work, we present Monte Carlo N-Particle (MCNP) simulations and simple models programmed in Mathematica to explore image reconstruction capabilities of this 80-voxel coded-aperture mask.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

$f$ 0 (980) production in inelastic pp collisions at s = 5.02 TeV

The measurement of the production of f 0 (980) in inelastic pp collisions at $\sqrt{s}$ = 5.02 TeV is presented. This is the first reported measurement of inclusive f 0 (980) yield at LHC energies. The production is measured at midrapidity, |y| < 0.5, in a wide transverse momentum range, 0 < p T < 16 GeV/c, by reconstructing the resonance in the f 0 (980) → π + π - hadronic decay channel using the ALICE detector. The p T -differential yields are compared to those of pions, protons and $φ$ mesons as well as to predictions from the HERWIG 7.2 QCD-inspired Monte Carlo event generator and calculations from a coalescence model that uses the AMPT model as an input. The ratio of the p T -integrated yield of f 0 (980) relative to pions is compared to measurements in e + e - and pp collisions at lower energies and predictions from statistical hadronisation models and HERWIG 7.2. A mild collision energy dependence of the f 0 (980) to pion production is observed in pp collisions from SPS to LHC energies. All considered models underpredict the p T -integrated 2f 0 (980)/(π + +π - ) ratio. The prediction from the canonical statistical hadronisation model assuming a zero total strangeness content of f 0 (980) is consistent with the data within 1.9σ and is the closest to the data. The results provide an essential reference for future measurements of the particle yield and nuclear modification in p–Pb and Pb–Pb collisions, which have been proposed to be instrumental to probe the elusive nature and quark composition of the f 0 (980) scalar meson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Design principles for energy transfer in the photosystem II supercomplex from kinetic transition networks

Photosystem II (PSII) has the unique ability to perform water-splitting. With light-harvesting complexes, it forms the PSII supercomplex (PSII-SC) which is a functional unit that can perform efficient energy conversion, as well as photoprotection, allowing photosynthetic organisms to adapt to the naturally fluctuating sunlight intensity. Achieving these functions requires a collaborative energy transfer network between all subunits of the PSII-SC. In this work, we perform kinetic analyses and characterise the energy landscape of the PSII-SC with a structure-based energy transfer model. With first passage time analyses and kinetic Monte Carlo simulations, we are able to map out the overall energy transfer network. We also investigate how energy transfer pathways are affected when individual protein complexes are removed from the network, revealing the functional roles of the subunits of the PSII-SC. In addition, we provide a quantitative description of the flat energy landscape of the PSII-SC. We show that it is a unique landscape that produces multiple kinetically relevant pathways, corresponding to a high pathway entropy. These design principles are crucial for balancing efficient energy conversion and photoprotection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Design and simulation of a muon detector to characterize geological overburden

This study presents the design, construction, and simulation of a mobile muon detector tailored for geological overburden characterization. The detector employs plastic scintillator paddles with silicon photomultipliers (SiPMs) and a QuarkNet data acquisition system, offering a portable solution suitable for remote field deployment. The simulator’s modular aluminum frame allows for adjustable geometry and directional sensitivity, while its battery system supports over a week of autonomous operation. Preliminary experimental tests confirmed that its muon flux measurements were consistent with theoretical expectations. A comprehensive simulation framework using Geant4 and CORSIKA was developed to model detector response and overburden effects. Analytical and Monte Carlo methods were used to assess quadrant resolution and infer muon directionality. This work lays the foundation for future overburden mapping and supports the development of reconstruction algorithms for geological applications.

72 - PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

An ultra-fast method for generating synthetic down-scattered neutron data for inertial confinement fusion implosions

In inertial confinement fusion experiments at the National Ignition Facility, asymmetries are probed by a variety of neutron diagnostics, including neutron imaging systems, real-time neutron activation diagnostics (RTNADs), and neutron spectrometers. It is often useful to generate synthetic data based on these diagnostics to validate and tune models. However, current methods of doing so using Monte Carlo particle tracing are time-consuming. In this paper, an ultra-fast method is presented for generating synthetic neutron images, RTNAD data, and spectrometry data using line integrals and 3D convolutions. While it does not contain as much physics as particle tracing codes, it is thousands of times faster and produces nearly identical data. This enables analysis techniques that depend on generating large amounts of synthetic data, which will prove very useful for the study of asymmetries going forward.

Deuterium↗
DOE OSTI2020

Commensuration Effects in Layered Nanoparticle Solids.

We have developed a hierarchical nanoparticle transport simulator (HINTS) and adapted it to study commensuration effects in two classes of nanoparticle (NP) solids: (1) a bilayer NP solid (BNS) with an energy offset and (2) a BNS as part of a field effect transistor (FET). HINTS integrates the ab initio characterization of single NPs with the phonon-assisted tunneling transition model of the NP-NP transitions into a kinetic Monte Carlo based simulation of the charge transport in NP solids. First, we studied a BNS with an interlayer energy offset Delta, possibly caused by a fixed electric field. Our results include the following. First, in the independent energy-offset model, we observed the emergence of commensuration effects when scanning the electron filling factor (FF) across integer values. These commensuration effects were profound as they reduced the mobility by several orders of magnitude. We analyzed these commensuration effects in a five-dimensional parameter space, as a function of the on-site charging energy E-C, energy offset Delta, disorder D, the electron FF, and temperature k(B)T. We demonstrated the complexity of our model by showing that at integer FFs commensuration effects are present in some regions of the parameter space, while they vanish in other regions, thus defining distinct dynamical phases of the model. We determined the phase boundaries between these dynamical phases. Second, using these results as a foundation, we shifted our focus to the experimentally much-studied NP-FETs. NP-FETs are also characterized by an interlayer energy offset Delta, which, in contrast to our first model, is set by the gate voltage V-G and thereby related to the electron FF. We repeated many of our simulations and again demonstrated the emergence of commensuration effects and distinct dynamical phases in these NP-FETs. Notably, the commensuration effects in the NP-FETs showed many similarities to those in the independent energy-offset BNS.

Qu, Luman↗
DOE OSTI2021

Quantifying multipoint ordering in alloys

A central problem in multicomponent lattice systems is to systematically quantify multipoint ordering. Ordering in such systems is often described in terms of pairs, even though this is not sufficient when three-point and higher-order interactions are included in the Hamiltonian. Current models and parameters for multipoint ordering are often only applicable for very specific cases or require approximating a subset of correlated occupational variables on a lattice as being uncorrelated. In this paper, cluster order parameters are introduced to systematically quantify arbitrary multipoint ordering motifs in substitutional systems through direct calculations of normalized cluster probabilities. These parameters can describe multipoint chemical ordering in crystal systems with multiple sublattices, multiple components, and systems with reduced symmetry. These are defined in this paper and applied to quantify four-point chemical ordering motifs in platinum/palladium alloy nanoparticles that are of practical interest to the synthesis of catalytic nanocages. Impacts of chemical ordering on nanocage stability are discussed. It is demonstrated that approximating four-point probabilities from superpositions of lower-order pair probabilities is not sufficient in cases where three- and four-body terms are included in the energy expression. Furthermore, conclusions about the formation mechanisms of nanocages may change significantly when using common pair approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhancement of charge density wave correlations in a Holstein model with an anharmonic phonon potential

The Holstein Hamiltonian describes itinerant electrons whose site density couples to local phonon degrees of freedom. In the single-site limit, at half filling, the electron-phonon coupling results in a double-well structure for the lattice displacement, favoring empty or doubly occupied sites. In two dimensions and on a bipartite lattice in d ≥ 2, intersite hopping causes these doubly occupied and empty sites to alternate in a charge density wave (CDW) pattern when the temperature is lowered. Because a discrete symmetry is broken, this occurs in a conventional second-order transition at finite T cdw . In this paper, we investigate the effect of changing the phonon potential energy to one with an intrinsic double-well structure even in the absence of electron-phonon coupling. While this aids in the initial process of pair formation, the implications for subsequent CDW order are nontrivial. One expects that, when the electron-phonon coupling is too strong, the double wells become deep and the polaron mass gets large, an effect which reduces T cdw . Here, we show here the existence of regions of parameter space where the double-well potential, while aiding local pair formation, does so in a way which also substantially enhances long-range CDW order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

EvoProtGrad (Directed Evolution for Proteins with Gradients) [SWR-23-48]

A Python package for directed evolution on a protein sequence with gradient-based discrete Markov chain monte carlo (MCMC). Users are able to compose custom models that map sequence to function with pretrained models, including protein language models (PLMs), to guide and constrain search. Our package natively integrates with the HuggingFace platform and supports PLMs from transformers. Our MCMC sampler identifies promising amino acids to mutate via model gradients taken with respect to the input (i.e., sensitivity analysis). We allow users to compose their own custom target function for MCMC by leveraging the Product of Experts MCMC paradigm. Each model is an "expert" that contributes its own knowledge about the protein's fitness landscape to the overall target function. The sampler is designed to be more efficient and effective than brute force and random search while maintaining most of the generality and flexibility. Additional information can be found in the related publication: https://iopscience.iop.org/article/10.1088/2632-2153/accacd

Emami, Patrick↗