Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Neural-network based collision operators for the Boltzmann equation

Kinetic gas dynamics in rarefied and moderate-density regimes have complex behavior associated with collisional processes. These processes are generally defined by convolution integrals over a high-dimensional space (as in the Boltzmann operator), or require evaluating complex auxiliary variables (as in Rosenbluth potentials in Fokker-Planck operators) that are challenging to implement and computationally expensive to evaluate. In this work, we develop a data-driven neural network model that augments a simple and inexpensive BGK collision operator with a machine-learned correction term, which improves the fidelity of the simple operator with a small overhead to overall runtime. The composite collision operator has a tunable fidelity and, in this work, is trained using and tested against a direct-simulation Monte-Carlo (DSMC) collision operator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Joint Modeling of Wind Speed and Wind Direction Through a Conditional Approach

Atmospheric near surface wind speed and wind direction play an important role in many applications, ranging from air quality modeling, building design, wind turbine placement to climate change research. It is therefore crucial to accurately estimate the joint probability distribution of wind speed and direction. In this work, we develop a conditional approach to model these two variables, where the joint distribution is decomposed into the product of the marginal distribution of wind direction and the conditional distribution of wind speed given wind direction. To accommodate the circular nature of wind direction, a von Mises mixture model is used; the conditional wind speed distribution is modeled as a directional dependent Weibull distribution via a two-stage estimation procedure, consisting of a directional binned Weibull parameter estimation, followed by a harmonic regression to estimate the dependence of the Weibull parameters on wind direction. A Monte Carlo simulation study indicates that our method outperforms two other approaches in estimation efficiency: one that utilizes periodic spline quantile regression and another that generates data from the commonly used Abe-Ley distribution for cylindrical data. We illustrate our method by using the output from a regional climate model to investigate how the joint distribution of wind speed and direction may change under some future climate scenarios. Our method indicates significant changes in the variation of wind speed with respect to some directions.

17 WIND ENERGY↗
DOE OSTI2021

Adjoint DSMC for nonlinear Boltzmann equation constrained optimization

Applications for kinetic equations such as optimal design and inverse problems often involve finding unknown parameters through gradient-based optimization algorithms. Based on the adjoint-state method, we derive two different frameworks for approximating the gradient of an objective functional constrained by the nonlinear Boltzmann equation. While the forward problem can be solved by the DSMC method, it is difficult to efficiently solve the high-dimensional continuous adjoint equation obtained by the “optimize-then-discretize” approach. This challenge motivates us to propose an adjoint DSMC method following the “discretize-then-optimize” approach for Boltzmann-constrained optimization. We also analyze the properties of the two frameworks and their connections. Here, several numerical examples are presented to demonstrate their accuracy and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enabling the broader adoption of fusion simulation on complex geometry

This project addressed a key barrier to advanced fusion and nuclear simulation: the difficulty of performing high-fidelity Monte Carlo neutronics directly on complex, real-world CAD geometry. Traditional workflows require engineers to rebuild CAD models as simplified constructive solid geometry, a time-consuming and error-prone process that limits design iteration and broader adoption of simulation tools. The goal of this Phase I SBIR was to make CAD-based neutronics practical, accessible, and robust for industrial and research users. During the project, Coreform significantly enhanced the Direct Accelerated Geometry Monte Carlo (DAGMC) workflow and fully integrated it into Coreform Cubit as a first-class capability. Major achievements include optimized material assignment and surface meshing workflows, substantial performance improvements to geometry imprinting and preparation, native export of DAGMC models, and new visualization tools to support OpenMC source definition and lost-particle debugging. Coreform also expanded Cubit’s capabilities as a full OpenMC preprocessor, including the ability to convert OpenMC constructive solid geometry models back into CAD for visualization, multiphysics coupling, and debugging. In collaboration with Argonne National Laboratory, the project delivered comprehensive new DAGMC documentation and training materials, transforming DAGMC from a research-oriented tool into a production-ready workflow. Results were disseminated through tutorials, conference training, and multiple well-attended webinars demonstrating integrated CAD-based neutronics and multiphysics workflows. Overall, this project demonstrated that high-fidelity Monte Carlo simulations can be performed directly on complex CAD geometry, reducing setup time, improving usability, and enabling faster, more informed design decisions for fusion and nuclear energy systems.

42 ENGINEERING↗
DOE OSTI2023

A physics-informed deep learning model of the hot tail runaway electron seed

A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Single molecule tracking of bacterial cell surface cytochromes reveals dynamics that impact long-distance electron transport

Using a series of multiheme cytochromes, the metal-reducing bacterium Shewanella oneidensis MR-1 can perform extracellular electron transfer (EET) to respire redox-active surfaces, including minerals and electrodes outside the cell. While the role of multiheme cytochromes in transporting electrons across the cell wall is well established, these cytochromes were also recently found to facilitate long-distance (micrometer-scale) redox conduction along outer membranes and across multiple cells bridging electrodes. Recent studies proposed that long-distance conduction arises from the interplay of electron hopping and cytochrome diffusion, which allows collisions and electron exchange between cytochromes along membranes. However, the diffusive dynamics of the multiheme cytochromes have never been observed or quantified in vivo, making it difficult to assess their hypothesized contribution to the collision-exchange mechanism. Here, we use quantum dot labeling, total internal reflection fluorescence microscopy, and single-particle tracking to quantify the lateral diffusive dynamics of the outer membrane-associated decaheme cytochromes MtrC and OmcA, two key components of EET in S. oneidensis. We observe confined diffusion behavior for both quantum dot-labeled MtrC and OmcA along cell surfaces (diffusion coefficients D MtrC = 0.0192 ± 0.0018 µm 2 /s, D OmcA = 0.0125 ± 0.0024 µm 2 /s) and the membrane extensions thought to function as bacterial nanowires. We find that these dynamics can trace a path for electron transport via overlap of cytochrome trajectories, consistent with the long-distance conduction mechanism. The measured dynamics inform kinetic Monte Carlo simulations that combine direct electron hopping and redox molecule diffusion, revealing significant electron transport rates along cells and membrane nanowires.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Understanding Reentrance in Frustrated Magnets: The Case of the Er 2 Sn 2 O 7 Pyrochlore

Reentrance, the return of a system from an ordered phase to a previously encountered less-ordered one as a controlled parameter is continuously varied, is a recurring theme found in disparate physical systems, yet its microscopic cause is often not investigated thoroughly. Here, through detailed characterization and theoretical modeling, we uncover the microscopic mechanism behind reentrance in the strongly frustrated pyrochlore antiferromagnet Er 2 Sn 2 O 7 . We use single crystal heat capacity measurements to expose that Er 2 Sn 2 O 7 exhibits multiple instances of reentrance in its magnetic field B vs temperature T phase diagram for magnetic fields along three cubic high symmetry directions. Through classical Monte Carlo simulations, mean field theory, and classical linear spin-wave expansions, we argue that the origins of the multiple occurrences of reentrance observed in Er 2 Sn 2 O 7 are linked to soft modes. These soft modes arise from phase competition and enhance thermal fluctuations that entropically stabilize a specific ordered phase, resulting in an increased transition temperature for certain field values and thus the reentrant behavior. Our work represents a detailed examination into the mechanisms responsible for reentrance in a frustrated magnet and may serve as a template for the interpretation of reentrant phenomena in other physical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deep probabilistic direction prediction in 3D with applications to directional dark matter detectors

Abstract We present the first method to probabilistically predict 3D direction in a deep neural network model. The probabilistic predictions are modeled as a heteroscedastic von Mises-Fisher distribution on the sphere S 2 , giving a simple way to quantify aleatoric uncertainty. This approach generalizes the cosine distance loss which is a special case of our loss function when the uncertainty is assumed to be uniform across samples. We develop approximations required to make the likelihood function and gradient calculations stable. The method is applied to the task of predicting the 3D directions of electrons, the most complex signal in a class of experimental particle physics detectors designed to demonstrate the particle nature of dark matter and study solar neutrinos. Using simulated Monte Carlo data, the initial direction of recoiling electrons is inferred from their tortuous trajectories, as captured by the 3D detectors. For 40 keV electrons in a 70% He 30% CO 2 gas mixture at STP, the new approach achieves a mean cosine distance of 0.104 (26 ∘ ) compared to 0.556 (64 ∘ ) achieved by a non-machine learning algorithm. We show that the model is well-calibrated and accuracy can be increased further by removing samples with high predicted uncertainty. This advancement in probabilistic 3D directional learning could increase the sensitivity of directional dark matter detectors.

Computer Science↗
DOE OSTI2022

Precision measurement of Compton scattering in silicon with a skipper CCD for dark matter detection

Experiments aiming to directly detect dark matter through particle recoils can achieve energy thresholds of $\mathcal{O}(1\,\mathrm{eV})$. In this regime, ionization signals from small-angle Compton scatters of environmental $\gamma$-rays constitute a significant background. Monte Carlo simulations used to build background models have not been experimentally validated at these low energies. We report a precision measurement of Compton scattering on silicon atomic shell electrons down to 23$\,$eV. A skipper charge-coupled device (CCD) with single-electron resolution, developed for the DAMIC-M experiment, was exposed to a $^{241}$Am $\gamma$-ray source over several months. Features associated with the silicon K, L$_{1}$, and L$_{2,3}$-shells are clearly identified, and scattering on valence electrons is detected for the first time below 100$\,$eV. We find that the relativistic impulse approximation for Compton scattering, which is implemented in Monte Carlo simulations commonly used by direct detection experiments, does not reproduce the measured spectrum below 0.5$\,$keV. The data are in better agreement with $ab$$initio$ calculations originally developed for X-ray absorption spectroscopy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization

We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Monte carlo simulations and experiments of all-silica zeolite LTA assembly combining structure directing agents that match cage sizes

In this work, we investigated the influence of organic structure-directing agents (OSDAs) on the formation rates of all-silica zeolite LTA using both simulations and experiments, to shed light on the crystallization process. We compared syntheses using one OSDA with a diameter close to the size of the large cavity in LTA, and two OSDAs of diameters matching the sizes of both the small and large LTA cavities. Reaction-ensemble Monte Carlo (RxMC) simulations predict a speed up of LTA formation using two OSDAs matching the LTA pore sizes; this qualitative result is confirmed by experimental studies of crystallization kinetics, which find a speedup in all-silica LTA crystallization of a factor of 3. Analyses of simulated rings and their Si–O–Si angular energies during RxMC crystallizations show that all ring sizes in the faster crystallization exhibit lower angular energies, on average, than in the slower crystallization, explaining the origin of the speedup through packing effects.

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

Unravelling Magnetic Nanochain Formation in Dispersion for In Vivo Applications

Self-assembly of iron oxide nanoparticles (IONPs) into 1D chains is appealing, because of their biocompatibility and higher mobility compared to 2D/3D assemblies while traversing the circulatory passages and blood vessels for in vivo biomedical applications. In this work, parameters such as size, concentration, composition, and magnetic field, responsible for chain formation of IONPs in a dispersion as opposed to spatially confining substrates, are examined. In particular, the monodisperse 27 nm IONPs synthesized by an extended LaMer mechanism are shown to form chains at 4 mT, which are lengthened with applied field reaching 270 nm at 2.2 T. The chain lengths are completely reversible in field. Using a combination of scattering methods and reverse Monte Carlo simulations the formation of chains is directly visualized. The visualization of real-space IONPs assemblies formed in dispersions presents a novel tool for biomedical researchers. This allows for rapid exploration of the behavior of IONPs in solution in a broad parameter space and unambiguous extraction of the parameters of the equilibrium structures. Additionally, it can be extended to study novel assemblies formed by more complex geometries of IONPs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Statistical emulation of a perturbed basal melt ensemble of an ice sheet model to better quantify Antarctic sea level rise uncertainties

Abstract. Antarctic ice shelves are vulnerable to warming ocean temperatures, and some have already begun thinning in response to increased basal melt rates. Sea level is therefore expected to rise due to Antarctic contributions, but uncertainties in its amount and timing remain largely unquantified. In particular, there is substantial uncertainty in future basal melt rates arising from multi-model differences in thermal forcing and how melt rates depend on that thermal forcing. To facilitate uncertainty quantification in sea level rise projections, we build, validate, and demonstrate projections from a computationally efficient statistical emulator of a high-resolution (4 km) Antarctic ice sheet model, the Community Ice Sheet Model version 2.1. The emulator is trained to a large (500-member) ensemble of 200-year-long 4 km resolution transient ice sheet simulations, whereby regional basal melt rates are perturbed by idealized (yet physically informed) trajectories. The main advantage of our emulation approach is that by sampling a wide range of possible basal melt trajectories, the emulator can be used to (1) produce probabilistic sea level rise projections over much larger Monte Carlo ensembles than are possible by direct numerical simulation alone, thereby providing better statistical characterization of uncertainties, and (2) predict the simulated ice sheet response under differing assumptions about basal melt characteristics as new oceanographic studies are published, without having to run additional numerical ice sheet simulations. As a proof of concept, we propagate uncertainties about future basal melt rate trajectories, derived from regional ocean models, to generate probabilistic sea level rise estimates for 100 and 200 years into the future.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Geometric control of emergent antiferromagnetic order in coupled artificial spin ices

Artificial spin ices (ASIs) composed of coupled nanomagnets offer the possibility to create designer geometrical frustration and manipulate inter-nanomagnet interactions. In particular, by using a dimer motif consisting of two strongly coupled single-domain nano magnets as a building block, we can control and realize intriguing antiferromagnetic physical states in which the magnetic charge is not conserved. Here we create a dimer kagome ASI system within which the antiferromagnetic order is controlled by tuning lattice geometry. Different disordered antiferromagnetic phases and the inter-nanomagnet interactions are visualized directly in demagnetized lattices with different lattice parameters. Monte Carlo simulations establish that the collective ground state consists of disordered antiferromagnetic dimers across the lattice when the intra-dimer interaction is dominant. However, for lattices governed by intra-triad-unit interactions, the ground state exhibits a longrange spin-ordered state in which the vertex magnetic charge is uniform across all triads.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unraveling small-scale defects in irradiated ThO 2 using kinetic Monte Carlo simulations

Point defects and their clusters generated through irradiation can have significant impact on the physical and mechanical properties of materials. However, direct experimental visualization of these small-scale defects using high-resolution scanning transmission electron microscopy remains a challenging task. In this study, using thorium dioxide (ThO 2 ) with the fluorite structure as a model system, we demonstrate the use of ab initio basin-hopping simulations in synergy with object kinetic Monte Carlo simulations as a powerful tool for identifying small defect complexes in irradiated materials. In addition to providing quantitative insights into defect evolution in ThO 2 under irradiation, our study reveals an unexpected role of bound anti-Schottky defect clusters in mediating defect transport. Remarkably, despite their poor thermal stability against dissociation at high temperatures, the transient formation of bound anti-Schottky defects under irradiation and their subsequent migration provide the dominant mechanism for the growth of large interstitial loops that have been experimentally observed in ThO 2 .

36 MATERIALS SCIENCE↗