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 1,621 records · Page 90

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Putting multipoint kinetics to the test: preliminary analysis of the MUSE-4 experiment with APOLLO3

Modeling neutron kinetics is a challenging task and over the years, several approaches have been developed offering different tradeoffs between precision and cost. Probably the most popular one is point kinetics (PK), which, despite being very cheap, provides good results for many applications. One of the limitations of PK is that it is not able to model effects that happen before the prompt jump time scale. In practice, applications, where those time scales are interesting, are rare and those for which experimental data are available are even rarer. One of those cases is the MUSE-4 experiment for which the focus is on the very fast behavior and where we have no thermal feedback. As we have shown in a recent paper, a simple extension of PK, known as multipoint kinetics (MPK), can overcome this limitation. In this paper, which is intended as a preliminary study, we show how a simple MPK model can qualitatively reproduce the fast behavior that we cannot model with PK. Moreover, we will show how to exploit its low cost and compact form to model also slow effects as the precursors built up, behaviours that are not affordable to model with more expensive techniques, such as kinetic Monte Carlo simulations. For the periodic pulse transients with no feedbacks we will also provide an analytical solution in closed form. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CO2-SCREEN

CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN) is a tool which was developed by the Department of Energy’s National Energy Technology Laboratory (DOE-NETL) to screen saline formations for prospective CO2 storage resources. CO2-SCREEN applies U.S. DOE methods and equations for estimating prospective CO2 storage resources for saline formations. This provides a dependable method for calculating prospective CO2 storage resources which allows for consistent comparison of results between different research efforts. CO2-SCREEN latest version (v5.0) is a python application that utilizes a Java front end user interface. The 3.0 version consists of an Excel spreadsheet containing geologic inputs and outputs, linked to a GoldSim Player model that calculates prospective CO2 storage resources via Monte Carlo simulation. https://www.youtube.com/watch?v=lhakk-HYfOI

CO2 Storage Resource,Prospective Storage,Screening↗

Enhanced Visualization Capabilities in OpenMC: Report of Activities under a DOE GAIN Voucher

This work was performed under a Cooperative Research and Development Agreement (CRADA) through the U.S. Department of Energy’s Gateway for Accelerated Innovations in Nuclear (GAIN) program in partnership with Oklo Inc., the industry participant. The broad goal of this project was to improve the visualization capabilities of OpenMC related to model development and simulation results. OpenMC is a Monte Carlo particle transport code that is developed and maintained by a broad community of contributors, of which Argonne National Laboratory plays a leading role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of Charged-Current $ν_μ$ and $ν_e$ Interactions Towards a Search for Low-Energy $ν_Ε$ Excess Using Wire-Cell In Microboone

This technote summarizes the existing work in searching for $ν_e$ low-energy excess (eLEE) in MicroBooNE Booster Neutrino Beam (BNB) data stream based on the Wire-Cell event reconstruction paradigm. The charged-current $ν_μ$ and $ν_e$ events are selected from the 5.3e19 POT open data from the BNB beam, 6.37e20 POT far sideband from the BNB beam data, and 2.10e20 POT data from Neutrinos at the Main Injector (NuMI) beam. The charged-current $ν_e$ selection results from the BNB data that are sensitive to the eLEE search are not included. Various comparisons between data and Monte Carlo predictions are performed to validate the overall model and demonstrate the power of the analysis techniques. Physics sensitivities in terms of the exclusion and the discovery potential are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of charged-current $\nu_{\mu}$ and $\nu_e$ interactions with wire-cell in MicroBooNE towards a search for low-energy $\nu_e$ excess

This technote summarizes the existing work in searching for $\nu_e$ low-energy excess (eLEE) in MicroBooNE Booster Neutrino Beam (BNB) data stream based on the Wire-Cell event reconstruction paradigm. The charged-current $\nu_{\mu}$ and $\nu_e$ events are selected from the 5.3e19 POT data from the BNB beam and 2.06e20 POT data from Neutrinos at the Main Injector (NuMI) beam. The charged-current $\nu_e$ selection results from the BNB data that are sensitive to the eLEE search are not included. Various comparisons between data and Monte Carlo predictions are performed to validate the overall model and demonstrate the power of the analysis techniques.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

Symmetry Breaking in an Extended-O(2) Model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we consider an extended-O(2) model that differs from the ordinary O(2) model by an explicit symmetry breaking term. Its coupling allows to smoothly interpolate between the O(2) model (zero coupling) and a $q$-state clock model (infinite coupling). In the latter case, a $q$-state clock model can also be defined for non-integer values of $q$. Thus, such a limit can also be considered as an analytic continuation of an ordinary $q$-state clock model to non-integer $q$. The phase diagram of the extended-O(2) model in the infinite coupling limit was established in our previous work, where it was shown that for non-integer $q$, there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we investigate the model at finite values of the coupling using Monte Carlo and tensor methods. The results may be relevant for configurable Rydberg-atom arrays.

Hostetler, Leon↗

Computational Optimization of A 3D Printed Collimator

This contribution describes the computational methodology behind an optimization procedure for a scattered beam collimator. The workflow includes producing a file that can be manufactured via additive methods. A conical collimator, optimized for neutron diffraction experiments in a high pressure clamp cell, is presented as an example. In such a case the scattering from the sample is much smaller than that of the pressure cell. Monte Carlo Ray tracing in MCViNE was used to model scattering from a Si powder sample and the cell. A collimator was inserted into the simulation and the number and size of channels were optimized to maximize the rejection of the parasitic signal coming from the complex sample environment. Constraints, provided by the additive manufacturing process as well as a specific neutron diffractometer, were also included in the optimization. The source code and the tutorials are available in c3dp (Islam (2019)).

74 ATOMIC AND MOLECULAR PHYSICS↗