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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Monte Carlo simulation methodology for characterizing the energy response function and resolving power of a digital Flat Panel Detector (FPD) assembly

The use of high-resolution, digital flat panel detectors (FPDs) makes quantitative assessments of radiographic and 3D tomographic representations feasible. However, the energy response function and spatial resolving power of the detector must be known. Here, we present a methodology for characterizing an FPD assembly using Monte Carlo simulations. The detector assembly modeled here was fielded in the 3-Ring Flash X-ray Imaging System designed and built at Lawrence Livermore National Laboratory (LLNL), which employs fifteen Carestream ® HPX-DR3543 FPDs enclosed in protective housings. Detailed Monte Carlo (MCNP6) simulations of an as-fielded detector assembly were performed to estimate its energy response function in terms of analog-to-digital unit (ADU) “counts” per pixel per incident x-ray as a function of energy up to ∼ 500 keV and its theoretical (idealized) resolving power. The practical resolving power of the FPD assembly in a realistic scenario was then estimated by simulating radiographic images of a variant of the United Kingdom (U.K.) “Kaleidoscope” resolution target.

Physics↗
DOE OSTI2023

X-ray Inspection Model Validation with Physical Dosimetry

A non-functional printed circuit board assembly was developed using ULTIboard and a simple computational transport model was created for use with the Monte Carlo software MCNP. EBT 3 film radiography was used to compare with Monte Carlo simulation to validate the board representation. The computational transport model was then revised to include connection pins, solder balls, highly attenuating internal structures, and copper trace distribution. Results from the revised model were compared to the EBT 3 films to observe improvements to dose profiles. It was found that this method was useful in verifying the placement of components, as the dose profiles were observed to follow the same trends. The experiment was then repeated using XRQA 2 films to achieve the same level of contrast with 1% the dose of EBT 3 films. It was found that high contrast may be achieved using these films to identify major issues with the model geometry, at a cost of dose profile accuracy. A second validation method was applied to the model using 37 CaF 2 thermoluminescent dosimeters (TLDs). TLD measurements were compared with the simplified and complex transport models to identify the features that have the greatest impact on simulation accuracy. The TLD calibration to CaF 2 was found to be accurate within 5.6%, while calibration to dose in Si was found to be accurate within 4.7%. Finally, it was observed that the accurate representation of solder balls and proper modeling of highly attenuating internal structures had the greatest impact on simulation accuracy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Modeling the spectral modification of lower hybrid wave in the presence of drift-wave type density fluctuation in the scrape-off-layer of the EAST tokamak

In this work, the spectrum change of the lower hybrid (LH) waves caused by a low-frequency density fluctuation in the scrape-off-layer is studied by applying the wave scattering model developed by Bonoli and Ott [Phys. Fluids 25, 359 (1982)] via a Monte Carlo method. A ray-tracing model specific to this purpose is developed to evaluate the probability distribution of both the poloidal refractive index (N θ ) and the parallel refractive index (N ∥ ) of the LH wave at the last closed flux surface (LCFS).The probability distributions of N ∥ and N θ at the LCFS are studied using the EAST parameters as a function of wave frequency, the initial N ∥ , and the polar injection position, which may influence the lower hybrid current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Ghent Hybrid model in NuWro: a new neutrino single-pion production model in the GeV regime

Neutrino-induced single-pion production constitutes an essential interaction channel in modern neutrino oscillation experiments, with its products building up a significant fraction of the observable hadronic final states. Frameworks of oscillation analyses strongly rely on Monte Carlo neutrino event generators, which provide theoretical predictions of neutrino interactions on nuclear targets. Thus, it is crucial to integrate state-of-the-art single-pion production models with Monte Carlo simulations to prepare for the upcoming systematics-dominated landscape of neutrino measurements. In this work, we present the implementation of the Ghent Hybrid model for neutrino-induced single-pion production in the NuWro Monte Carlo event generator. The interaction dynamics includes coherently-added contributions from nucleon resonances and a non-resonant background, merged into the pythia branching predictions in the deep-inelastic regime, as instrumented by NuWro. This neutrino-nucleon interaction model is fully incorporated into the nuclear framework of the generator, allowing it to account for the influence of both initial- and final-state nuclear medium effects. We compare the predictions of this integrated implementation with recent pion production data from accelerator-based neutrino experiments. The results of the novel model show improved agreement of the generator predictions with the data and point to the significance of the refined treatment of the description of pion-production processes beyond the ∆ region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Collinear drop

We introduce collinear drop jet substructure observables, which are unaffected by contributions from collinear radiation, and systematically probe soft radiation within jets. These observables can be designed to be either sensitive or insensitive to process-dependent soft radiation originating from outside the jet. Such collinear drop observables can be exploited as variables to distinguish quark, gluon, and color neutral initiated jets, for testing predictions for perturbative soft radiation in Monte Carlo simulations, for assessing models and universality for hadronization corrections, for examining the efficiency of pileup subtraction methods, and for any other application that leaves an imprint on soft radiation. We discuss examples of collinear drop observables that are based both on clustering and on jet shapes. Using the soft-collinear effective theory we derive factorization expressions for collinear drop observables from QCD jets, and carry out a resummation of logarithmically enhanced contributions at next-to-leading-logarithmic order. We also identify an infinite class of collinear drop observables for which the leading double logarithms are absent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Fitting a deep generative hadronization model

Hadronization is a critical step in the simulation of high-energy particle and nuclear physics experiments. As there is no first principles understanding of this process, physically-inspired hadronization models have a large number of parameters that are fit to data. Deep generative models are a natural replacement for classical techniques, since they are more flexible and may be able to improve the overall precision. Proof of principle studies have shown how to use neural networks to emulate specific hadronization when trained using the inputs and outputs of classical methods. However, these approaches will not work with data, where we do not have a matching between observed hadrons and partons. In this paper, we develop a protocol for fitting a deep generative hadronization model in a realistic setting, where we only have access to a set of hadrons in data. Our approach uses a variation of a Generative Adversarial Network with a permutation invariant discriminator. We find that this setup is able to match the hadronization model in Herwig with multiple sets of parameters. This work represents a significant step forward in a longer term program to develop, train, and integrate machine learning-based hadronization models into parton shower Monte Carlo programs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Anomaly Detection in Materials Digital Twins with Multiscale ICME for Additive Manufacturing

Detecting anomaly in fatigue and fracture experimental materials science is an interesting yet challenging topic. The reasons are threefold. First, the anomalous microstructure feature that gives rise to structural failure is small, sometimes in the order of 10 -7 of the interrogated volume. This, in turn, results in a highly imbalanced classification problem in machine learning (ML). Second, the consequence is high, in the sense that the test specimen is destructed in such case. Third, the convolution between microstructure stochasticity and the small probability of void nucleation, growth, and coalescence makes failure and fracture a hard-to-predict and challenging problem in materials science due to its irreproducibility, even experimentally. In this paper, we developed a materials digital twin and applied anomaly detection methods to detect voids and anomaly in additive manufacturing (AM). The materials digital twin is driven by two integrated computational materials engineering (ICME) models, which are kinetic Monte Carlo (kMC) and crystal plasticity finite element method (CPFEM). In conclusion, we demonstrated that by using anomaly detection, it is possible to detect voids and other defects in materials digital twin, which paves way for future research in integrating materials digital twin with its physical counterpart.

ICME↗
DOE OSTI2023

The importance of capturing power system operational details in resource adequacy assessments

Traditional methods for assessing the resource adequacy (RA) of a power system are becoming obsolete due to emerging trends such as the increasing deployment of variable renewable energy and storage. Consequently, analysts are recommending that RA be assessed using a Monte Carlo simulation approach that models chronological power system operations over many instances of possible operating conditions. However, this approach is necessarily more complex and computationally demanding, which is an obstacle to real-world implementation. Here, in this study, we investigate which operational details of power systems are important to capture in order to accurately evaluate a system's RA, versus details that add complexity but do not meaningfully affect RA results. To do so, we develop a probabilistic RA assessment framework by adapting an existing production cost model and apply it to a case study based on the IEEE Reliability Test System. Our results indicate that multi-year data, storage dispatch, and transmission limits are key details to incorporate. Accurate RA results can be obtained using non-economic dispatch strategies as long as they are coordinated with detailed operational strategies. We also demonstrate how popular expectation-based RA metrics can mask important differences in the characteristics of loss of load events.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

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

Cluster-Spin-Glass Magnetic Behavior and Morphology in the Coordination Polymer Alloys Fe y Co 1– y BTT

We report the synthesis of a series of pseudo-1D coordination polymer (CP) materials with the formula Fe y Co 1–y BTT (BTT = 1,3,5-benzenetrithiolate). These materials were structurally characterized by PXRD Rietveld, EXAFS, and PDF analyses, revealing that the CP superstructure enables a continuous and isomorphous alloy between the two homometallic compounds. Lower Fe loadings exhibit emergent spin glass magnetic behavior, such as memory effects and composition-dependent spin glass response time constants ranging from 6.9 × 10 –9 s to 1.8 × 10 –6 s. These data are consistent with the formation of spin clusters within the lattice. The magnetic behavior in these materials was modeled via replica exchange Monte Carlo simulation, which provides a good match for the experimentally measured spin glassing and magnetic phase transitions. Furthermore, these findings underscore how the rigid superstructure of CP and MOF scaffolds can enable the systematic tuning of physical properties, such as the spin glass behavior described here.

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

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

Abstract The origin of the pseudogap behavior, found in many high- T c superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a “gap-filling" rather than “gap-closing" behavior, similar to the one observed in cuprate superconductors. Our results present direct evidence of the pseudogap state, driven by superconducting fluctuations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

DC Hall coefficient of the strongly correlated Hubbard model

The Hall coefficient is related to the effective carrier density and Fermi surface topology in non-interacting and weakly interacting systems. In strongly correlated systems, the relation between the Hall coefficient and single-particle properties is less clear. Clarifying this relation would give insight into the nature of transport in strongly correlated materials that lack well-formed quasiparticles. In this work, we investigate the DC Hall coefficient of the Hubbard model using determinant quantum Monte Carlo in conjunction with a recently developed expansion of magneto-transport coefficients in terms of thermodynamic susceptibilities. At leading order in the expansion, we observe a change of sign in the Hall coefficient as a function of temperature and interaction strength, which we relate to a change in the topology of the apparent Fermi surface. We also combine our Hall coefficient results with optical conductivity values to evaluate the Hall angle, as well as effective mobility and effective mass based on Drude theory of metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A cyclic electrochemical strategy to produce acetylene from CO 2 , CH 4 , or alternative carbon sources

This paper presents the first combined measurement of the double-differential muon neutrino and antineutrino charged-current cross sections with no pions in the final state on hydrocarbon at the off-axis near detector of the T2K experiment. The data analyzed in this work comprise 5.8 × 10$^{20}$ and 6.3 × 10$^{20}$ protons on target in neutrino and antineutrino mode respectively, at a beam energy peak of 0.6 GeV. Using the two measured cross sections, the sum, difference, and asymmetry were calculated with the aim of better understanding the nuclear effects involved in such interactions. The extracted measurements have been compared with the prediction from different Monte Carlo generators and theoretical models showing that the difference between the two cross sections have interesting sensitivity to nuclear effects.

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

Role of hot electrons in shock ignition constrained by experiment at the National Ignition Facility

Shock ignition is a scheme for direct drive inertial confinement fusion that offers the potential for high gain with the current generation of laser facility; however, the benefits are thought to be dependent on the use of low adiabat implosions without laser–plasma instabilities reducing drive and generating hot electrons. A National Ignition Facility direct drive solid target experiment was used to calibrate a 3D Monte Carlo hot-electron model for 2D radiation-hydrodynamic simulations of a shock ignition implosion. The α = 2.5 adiabat implosion was calculated to suffer a 35% peak areal density decrease when the hot electron population with temperature T ℎ = 55 keV and energy E ℎ = 13 kJ was added to the simulation. Optimizing the pulse shape can recover ~1/3 of the peak areal density lost due to a change in shock timing. Here, despite the harmful impact of laser–plasma instabilities, the simulations indicate shock ignition as a viable method to improve performance and broaden the design space of near ignition high adiabat implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

CADIS and FW-CADIS Variance Reduction in Gamma Transport for Predicting Prompt Forensics Signatures

The goal of prompt nuclear forensics is to determine the characteristics of a nuclear detonation based on the signatures available almost immediately after the explosion. An important characteristic is the reaction time history (RTH), a measure of the device’s rate of neutron multiplication. The RTH can be estimated by observation of the gamma radiation emitted from the detonation, which can be detected directly or observed indirectly as Teller light. Gamma transport simulations used to predict these radiation fields are often modeled stochastically using the Monte Carlo N-Particle (MCNP) code, which can be a computationally demanding task due to the number of particle histories needed to achieve statistical convergence. In an attempt to improve the efficiency of these calculations, we evaluate two variance reduction techniques: Consistent Adjoint-Driven Importance Sampling (CADIS) and Forward-Weighted Consistent Adjoint-Driven Importance Sampling (FW-CADIS). These methods use a deterministically calculated adjoint flux to create weight windows and source biasing that guide MCNP sampling. We study the utility of CADIS and FW-CADIS for their use in MCNP gamma transport for nuclear forensics prediction simulations. Furthermore, the results demonstrate that both CADIS and FW-CADIS improve the accuracy for forensics-focused simulations, with CADIS being most beneficial in direct detection and FW-CADIS being ideal for computing a global Teller light source.

CADIS↗
DOE OSTI2021

ASAS J071404+7004.3 – a close, bright nova-like cataclysmic variable with gusty winds

ABSTRACT Despite being bright (V ≃ 11.8) and nearby (d = 212 pc) ASAS J071404+7004.3 has only recently been identified as a nova-like cataclysmic variable. We present time-resolved optical spectroscopy obtained at the Isaac Newton and the Hiltner and McGraw-Hill Telescopes, together with Swift X-ray and ultraviolet observations. We combined these with TESS photometry and find a period of 3.28 h and a mass transfer rate of $4\!-\!9\times 10^{-9}\, {\mathrm{M_{\odot }\, yr}^{-1}}$. Historical photometry shows at least one low state establishing the system as a VY Scl star. Our high-cadence spectroscopy also revealed rapidly changing winds emanating from the accretion disc. We have modelled these using the Monte Carlo python code and shown that all the emission lines could emanate from the wind – which would explain the lack of double-peaked lines in such systems. In passing, we discuss the effect of variability on the position of cataclysmic variables in the Gaia Hertzsprung–Russell diagram.

Inight, K. (ORCID:0000000222002416)↗