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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Impact of ferromagnetic electrode length and thickness on Magnetic Tunnel Junction-Based Molecular Spintronic Devices (MTJMSD)

A knowledge gap exists about the impact of variation in length and thickness of ferromagnetic(FM) electrodes on molecular spintronics devices’ ability to manifest molecular coupling impact. Magnetic Tunnel Junction-Based Molecular Spintronic Devices (MTJMSDs) are potential candidates for inventing highly correlated materials and devices and investigating the fundamental science of organic spintronics. This paper reports our experi- mental observations providing the dramatic impact of variation in ferromagnetic electrode length and thickness on paramagnetic molecule-based MTJMSD. Room temperature transport studies were performed to investigate the effect of FM electrode thickness. On the other hand, magnetic force microscopy measurements were con- ducted to understand the effect of FM electrode length extending beyond the molecular junction area, i.e., the site where paramagnetic molecules bridged between two FM. In the strong molecular coupling regime, transport study suggested thickness variation caused ~1000 to million-fold differences in junction conductivity. MFM study revealed near-zero magnetic contrast for pillar-shaped MTJMSD without any extended FM electrode. However, MFM images showed a multitude of microscopic magnetic phases on cross junction shaped MTJMSD where FM electrodes extended beyond the junction area. To understand the intriguing experimental results, we conducted an in-depth theoretical study using Monte Carlo Simulation (MCS) approach. MCS study utilized a Heisenberg atomic model of cross junction shaped MTJMSD to gain insights about room temperature transport and MFM experimental observations of microscopic MTJMSD. To make this study applicable for a wide variety of MTJMSDs, we systematically studied the effect of variation in molecular coupling strength between magnetic molecules and ferromagnetic (FM) electrodes of various dimensions. Here, our theoretical results suggest 25 atoms thick FM electrode show very weak molecular coupling impact and of the same order of experimentally studied FM electrode dimensions showing weak effect. Also, MCS results agree with experimentally observed multiple magnetic phases on lengthy FM electrodes of an MTJMSD.

42 ENGINEERING↗
DOE OSTI2021

Synthesis of magnesiowüstite nanocrystallites embedded in an amorphous silicate matrix via low energy multiple ion implantations

The synthesis process is presented for experimentally simulating modifications in cosmic dust grains using sequential ion implantations or irradiations followed by thermal annealing. Cosmic silicate dust analogues were prepared via implantation of 20–80 keV Fe - , Mg - , and O - ions into commercially available p-type silicon (100)wafers. The as-implanted analogues are amorphous with a Mg/(Fe+Mg) ratio of 0.5 tailored to match theoretical abundances in circumstellar dusts. Before the ion implantations were performed, Monte-Carlo-based ion-solid interaction codes were used to model the dynamic redistribution of the implanted atoms in the silicon substrate. 600 keV helium ion irradiation was performed on one of the samples before thermal annealing. Two samples were thermally annealed at a temperature appropriate for an M-class stellar wind, 1000 K, for 8.3 h in a vacuum chamber with a pressure of 1 x 10 -7 torr. The elemental depth profiles were extracted utilizing Rutherford Backscattering Spectrometry (RBS) in the samples before and after thermal annealing. X-ray diffraction (XRD)analysis was employed for the identification of various phases in crystalline minerals in the annealed analogues. Transmission electron microscopy (TEM) analysis was utilized to identify specific crystal structures. RBS analysis shows redistribution of the implanted Fe, Mg, and O after thermal annealing due to incorporation into the crystal structures for each sample type. XRD patterns along with TEM analysis showed nanocrystalline Mg and Fe oxides with possible incorporation of additional silicate minerals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Several Small or Single Large? Quantifying the Catchment-Wide Performance of On-Site Wastewater Treatment Plants with Inaccurate Sensors

On-site wastewater treatment plants (OSTs) often lack monitoring, resulting in unreliable treatment performance. They thus appear to be a stopgap solution despite their potential contribution to circular water management. Low-maintenance but inaccurate soft sensors are emerging that address this concern. However, how their inaccuracy impacts the catchment-wide treatment performance of a system of many OSTs has not been quantified. Here, we develop a stochastic model to estimate catchment-wide OST performances with a Monte Carlo simulation. In our study, soft sensors with a 70% accuracy improved the treatment performance from 66% of the time functional to 98%. Soft sensors optimized for specificity, indicating the true negative rate, improve the system performance, while sensors optimized for sensitivity, indicating the true positive rate, quantify the treatment performance more accurately. This new insight leads us to suggest programming two soft sensors in practical settings with the same hardware sensor data as input: one soft sensor geared to high specificity for maintenance scheduling and one geared to high sensitivity for performance quantification. Our findings suggest that a maintenance strategy combining inaccurate sensors with appropriate alarm management can vastly improve the mean catchment-wide treatment performance of a system of OSTs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deciphering the Scattering of Mechanically Driven Polymers Using Deep Learning

Here, we present a deep learning approach for analyzing two-dimensional scattering data of semiflexible polymers under external forces. In our framework, scattering functions are compressed into a three-dimensional latent space using a Variational Autoencoder (VAE), and two converter networks establish a bidirectional mapping between the polymer parameters (bending modulus, stretching force, and steady shear) and the scattering functions. The training data are generated using off-lattice Monte Carlo simulations to avoid the orientational bias inherent in lattice models, ensuring robust sampling of polymer conformations. The feasibility of this bidirectional mapping is demonstrated by the organized distribution of polymer parameters in the latent space. By integrating the converter networks with the VAE, we obtain a generator that produces scattering functions from given polymer parameters and an inferrer that directly extracts polymer parameters from scattering data. While the generator can be utilized in a traditional least-squares fitting procedure, the inferrer produces comparable results in a single pass and operates 3 orders of magnitude faster. This approach offers a scalable automated tool for polymer scattering analysis and provides a promising foundation for extending the method to other scattering models, experimental validation, and the study of time-dependent scattering data.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2022

Fabrication, thermal analysis, and heavy ion irradiation resistance of epoxy matrix nanocomposites loaded with silane-functionalized ceria nanoparticles

This paper describes a detailed understanding of how nanofillers function as radiation barriers within the polymer matrix, and how their effectiveness is impacted by factors such as composition, size, loading, surface chemistry, and dispersion. Here, we designed a comprehensive investigation of heavy ion irradiation resistance in epoxy matrix composites loaded with surface-modified ceria nanofillers, utilizing tandem computational and experimental methods to elucidate radiolytic damage processes and relate them to chemical and structural changes observed through thermal analysis, vibrational spectroscopy, and electron microscopy. A detailed mechanistic examination supported by FTIR spectroscopy data identified the bisphenol A moiety as a primary target for degradation reactions. Results of computational modeling by the Stopping Range of Ions in Matter (SRIM) Monte Carlo simulation were in good agreement with damage analysis from surface and cross-sectional SEM imaging. All metrics indicated that ceria nanofillers reduce the damage area in polymer nanocomposites, and that nanofiller loading and homogeneity of dispersion are key to effective damage prevention. The results of this study represent a significant pathway for engineered irradiation tolerance in a diverse array of polymer nanocomposite materials. Numerous areas of materials science can benefit from utilizing this facile and effective method to extend the reliability of polymer materials.

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

Measurement of the A dependence of the ν μ charged-current quasielasticlike cross section as a function of muon and proton kinematics at ⟨E ν ⟩ ∼ 6 GeV

The first simultaneous measurements of the 𝜈 𝜇 quasielasticlike cross section on C, CH, H 2 ⁡O, Fe, and Pb targets as a function of kinematic imbalance variables in the plane transverse to the incoming neutrino direction are presented. These variables combine the muon and proton information to provide a new way to disentangle the effects of the nucleus in quasielasticlike processes. The data were obtained using a wideband 𝜈 𝜇 beam with ⟨E 𝜈 ⟩ ∼ 6 GeV. Cross-section ratios of the different target materials to CH are also shown. These measurements are used to explore the nature of the cross-section 𝐴 scaling, as well as initial and final state interaction effects. Comparisons are made to predictions from a number of commonly used neutrino Monte Carlo event generators. The range of predictions of the different models tends to cover the data but the degree and consistency of the agreement suffers in regions, and on higher 𝐴 targets, where the final state interactions are expected to be more pronounced.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Properties of La 0.7 Ca 0.3 MnO 3 under extreme tensile strain

The complex phase diagram of manganites with simultaneously active spin, charge, orbital, and lattice degrees of freedom continues providing surprises. In a recent groundbreaking experiment, membranes of the perovskite manganite La 0.7 Ca 0.3 MnO 3 (LCMO) deposited on a flexible polymer layer were strained up to 8% [S. S. Hong et al., Science 368, 71 (2020)], much more than achieved by regular strain induced by a rigid substrate. By increasing this strain, a metal-insulator transition was reported. Here we reproduce the results of the experiments using Monte Carlo simulations of the two-orbital double-exchange model including Jahn-Teller distortions at hole density x=1/3. The full phase diagram with varying temperature and Jahn-Teller coupling λ is presented. When the bandwidth W of mobile electrons is reduced, thus when the effective Jahn-Teller coupling λ/W is increased, a metal-insulator transition is found in our simulations, between a ferromagnetic metallic state with uniform charge distribution and an insulator with diagonal charge stripes that retains its ferromagnetic character. In between the hole-rich diagonals, staggered orbital order occurs. We also report resistivity and magnetization measurements alongside with spin correlations and charge structure factors. Our overall conclusions are in agreement with the recent experimental and density functional theory results by Hong et al., and we confirm much earlier ground state predictions of striped ferromagnetic order using energy optimization techniques by T. Hotta et al. [Phys. Rev. Lett. 86, 4922 (2001)]. The experimental observation of one of the states predicted by theory suggests that diagonal stripes could be achieved at other hole densities as well, such as x=1/4, if LCMO membranes with that hole doping were subject to similar strains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Superuniversality from disorder at two-dimensional topological phase transitions

We investigate the effects of quenched randomness on topological quantum phase transitions in strongly interacting two-dimensional systems. We focus first on transitions driven by the condensation of a subset of fractionalized quasiparticles (“anyons”) identified with “electric charge” excitations of a phase with intrinsic topological order. All other anyons have nontrivial mutual statistics with the condensed subset and hence become confined at the anyon condensation transition. Furthermore, using a combination of microscopically exact duality transformations and asymptotically exact real-space renormalization group techniques applied to these two-dimensional disordered gauge theories, we argue that the resulting critical scaling behavior is “superuniversal” across a wide range of such condensation transitions and is controlled by the same infinite-randomness fixed point as that of the 2D random transverse-field Ising model. We validate this claim using large-scale quantum Monte Carlo simulations that allow us to extract zero-temperature critical exponents and correlation functions in (2+1)D disordered interacting systems. We discuss generalizations of these results to a large class of ground-state and excited-state topological transitions in systems with intrinsic topological order as well as those where topological order is either protected or enriched by global symmetries. When the underlying topological order and the symmetry group are Abelian, our results provide prototypes for topological phase transitions between distinct many-body localized phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electron-positron, parton-parton, and photon-photon production of τ -lepton pairs: Anomalous magnetic and electric dipole moments spin effects

Anomalous contributions to the electric and magnetic dipole moments of the τ lepton from new physics scenarios have brought renewed interest in the development of new charge-parity violating signatures in τ -pair production at Belle II energies, and also at higher energies of the Large Hadron Collider and the Future Circular Collider. In this paper, we discuss the effects of spin correlations, including transverse degrees of freedom, in the τ -pair production and decay. These studies include calculating analytical formulas, obtaining numerical results, and building semirealistic observables sensitive to the transverse spin correlations induced by the dipole moments of the τ lepton. The effects of such anomalous contributions to the dipole moments are introduced on top of precision simulations of e − e + → τ − τ + , q q ¯ → τ − τ + and γ γ → τ − τ + processes, involving multibody final states. The τ decays are simulated along with radiative corrections, in particular electroweak box contributions of W W and Z Z exchanges are taken into account. Respective extensions of the Standard Model amplitudes and the reweighting algorithms are implemented into the Monte Carlo, which is used to simulate τ -pair production in e − e + collisions, and the aupinner program, which is used to reweight events with τ pair produced in p p collisions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Fractionalized Excitations Revealed by Entanglement Entropy

Fractionalized excitations develop in many unusual many-body states such as quantum spin liquids, disordered phases that cannot be described using any local order parameter. Also, because these exotic excitations correspond to emergent degrees of freedom, how to probe them and establish their existence is a long-standing challenge. We present a general procedure to reveal the fractionalized excitations using real-space entanglement entropy in critical spin liquids that are particularly relevant to experiments. Furthermore, we show how to use the entanglement entropy to construct the corresponding spinon Fermi surface. Our work defines a new pathway to establish and characterize exotic excitations in novel quantum phases of matter.

2-dimensional systems↗
DOE OSTI2026

mvBayes

SAND2026-16980O mvBayes implements multivariate Bayesian regression using MATLAB and decomposes a multivariate or functional response into components based on a user-specified orthogonal basis. This allows for independent modeling of each component with any chosen univariate Bayesian regression model. This tool includes methods for prediction and visualization, facilitating the evaluation of Bayesian surrogate models through the application of Bayesian theory and Markov Chain Monte Carlo (MCMC) sampling techniques. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Tucker, J. Derek [Sandia National Lab. (SNL-CA), L↗
DOE OSTI2025

Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model: Modeling Archive

This archive contains the data and Python scripts required to reproduce the analyses and figures in the study: Gomez-Velez, J. D., Rathore, S. S., Cohen, M. J., & Painter, S. L. (2025). Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model. Submitted to Water Resources Research. The analysis utilizes the subgrid model Advection Dispersion Equation with Lagrangian Subgrids (ADELS) implemented in the Advanced Terrestrial Simulator (ATS; https://amanzi.github.io/ats/stable/). In this case, the ATS and Amanzi versions are (1) ATS version 1.5.1_f5ba18f8 and (2) Amanzi version 1.6-dev_53444cca4. The repository includes a Jupyter Notebook and the necessary data (Pandas DataFrames stored as pickle files) to generate the figures for the manuscript. Additionally, it contains Python scripts to create ATS input files, run the ATS simulations, and post-process the results. Finally, it provides routines for parameter estimation using the Single-Station Metabolism (SSM) model with the Differential Evolution Adaptive Metropolis (DREAM) Markov Chain Monte Carlo (MCMC) algorithm with ZS enhancements (DREAM-ZS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Initial development of a high temperature life prediction method directly accounting for variability in material properties

This report presents the initial development of a life prediction method that accounts for the variability of the material properties in Grade 91 steel. We account for material variability by both fitting a variable 3-parameter Weibull distribution to the rupture data of the material and by quantifying the effect that the creep deformation properties have on stresses via a Monte Carlo approach. To ensure a reasonable time frame when applying Monte Carlo to the creep stress, we propose a methodology that models the material as an extremely viscous Stokes flow with a non-Newtonian viscosity, therefore solving the problem in a single step. We use the final probabilistic model for structural failure developed in this document to evaluate a Flat Head Vessel designed with the Section III, Division 5 ASME Boiler Pressure and Vessel Code rules with a 100,000h design life and assess its probability of premature failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Knowing our Sun: data fusion for optimizing space weather forecasts

Outline: Space weather: model for predicting solar wind & polarity; Particle filtering: optimization with Monte Carlo; Simulation & observation: twin tests and real data; Back to Earth: independent work in inertial confinement fusion (ICF)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Verification of Boundary Conditions in MCNP

MCNP (Monte Carlo N-Particle) is code that is used to model the transport of different particles. To verify in MCNP6.1 that arrays formed with reflective and periodic boundary conditions produce a neutron multiplication factor (keff) that is equal to or higher than an array of the same size with no boundary conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Rapid Gamma Simulations of TRISO Fuel Elements

As energy demand rises, nuclear energy, particularly from reactors that use tristructural isotropic (TRISO) fuels, has gained attention due to the fuel’s enhanced resistance to radiation damage and high temperatures. This report investigates the modeling capabilities of the Gamma Detector Response and Analysis Software (GADRAS) for TRISO fuels, focusing on the gamma signatures of TRISO particles, which have not been extensively explored. Using the Monte Carlo N-Particle (MCNP) code as a benchmark, we developed both homogeneous and heterogeneous models of TRISO pebbles to analyze gamma spectra. Our findings reveal that the homogeneous and heterogeneous models produced different gamma signatures. Additionally, the GADRAS heterogeneous model significantly reduces computation times compared to MCNP, enabling effective modeling of gamma signatures for safeguards applications. This advancement is essential for the International Atomic Energy Agency (IAEA) in detecting anomalies and potential smuggling attempts in TRISO reactor fuel elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

As-Built Simulation of the High Flux Isotope Reactor

The Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) is an 85 MWt flux trap-type research reactor that supports key research missions, including isotope production, materials irradiation, and neutron scattering. The core consists of an inner and an outer fuel element containing 171 and 369 involute-shaped plates, respectively. The thin fuel plates consist of a U 3 O 8 -Al dispersion fuel (highly enriched), an aluminum-based filler, and aluminum cladding. The fuel meat thickness is varied across the width of the involute plate to reduce thermal flux peaks at the radial edges of the fuel elements. Some deviation from the designed fuel meat shaping is allowed during manufacturing. A homogeneity scan of each fuel plate checks for potential anomalies in the fuel distribution by scanning the surface of the plate and comparing the attenuation of the beam to calibration standards. While typical HFIR simulations use homogenized fuel regions, explicit models of the plates were developed under the Low-Enriched Uranium Conversion Program. These explicit models typically include one inner and one outer fuel plate with nominal fuel distributions, and then the plates are duplicated to fill the space of the corresponding fuel element. Therefore, data extracted from these simulations are limited to azimuthally averaged quantities. To determine the reactivity and physics impacts of an as-built outer fuel element and generate azimuthally dependent data in the element, 369 unique fuel plate models were generated and positioned. This model generates the three-dimensional (i.e., radial–axial–azimuthal) plate power profile, where the azimuthal profile is impacted by features within the adjacent control element region and beryllium reflector. For an as-built model of the outer fuel element, plate-specific homogeneity data, 235 U loading, enrichment, and channel thickness measurements were translated into the model, yielding a much more varied azimuthal power profile encompassed by uncertainty factors in analyses. These models were run with the ORNL-TN and Shift Monte Carlo tools, and they contained upwards of 500,000 cells and 100,000 unique tallies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Status of HEU-Pb in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook

The Department of Energy invests tens of millions of dollars each year to develop the next generation of nuclear engineering modeling & simulation (M&S) tools. These tools are used to analyze advanced reactor designs and the safety of current nuclear operations. As computers become more powerful, we are able to enhance resolution in our calculations. This improved resolution is taking us to a point where the limitations of simulation capability are in the quality of data, including our ability to quantify the uncertainty and sensitivity of the data. In order to model systems of interest with increasing accuracy, the industry must improve key nuclear data measurements. Thus, M&S tools need evaluated and quality-assured experimental data for validation purposes. The International Criticality Safety Benchmark Evaluation Project (ICSBEP) compiles and evaluates experiment data in a handbook that can be used by criticality safety engineers and others to validate computer codes and cross-section libraries at nuclear facilities. Both critical and subcritical experiments are included in the handbook. Figure 1 organizes all the benchmark evaluations that have been performed by the isotope of interest, in this case Pb, and the average neutron energy the system. Compared to other isotopes of interest for nuclear applications, there are few benchmark evaluations for Pb systems. The lack of integral measurements to determine errors in Pb cross-section data has caused the latest nuclear cross-section libraries to over/underestimate k eff compared to experimental results. Therefore, this evaluation fills an important knowledge gap in benchmark evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗