Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discretization”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electron-Proton Scattering Event Generation using Structured Tokenization

Recent work such as Omnijet-$\alpha$ has demonstrated that effective tokenization combined with transformer-based architectures can produce effective foundation models for jet physics. While tokenization may help models capture generalizable event characteristics, it also introduces discretization errors that may compromise the precision required for downstream physics analyses. As the number and complexity of the particle features grow, these errors are likely to grow proportionally. In this study, we investigate new tokenization strategies to improve the application of generative transformer models to \textsc{Pythia8} simulations of electron-proton scattering at the Electron-Ion Collider. Specifically, we propose a feature-based structured tokenization approach that utilizes multiple tokens per particle, improving expressivity, while reducing the total number of unique tokens needed. We evaluate this method against grid-based binning, K-means clustering, and vector-quantized variational auto-encoders on the event simulations. Our results show that feature-based structured tokenization reduces discretization error, leading to more accurate generative modeling of particle-level events.

Goldenberg, Steven [Thomas Jefferson National Acce↗
DOE OSTI2025

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High-fidelity Pebble Bed Reactor Depletion Based on Pebble Tracking Transport in Griffin

The pebble tracking transport (PTT) method is a high-fidelity, heterogeneous deterministic transport technique for pebble bed reactor analysis. It discretizes the broad-group neutron transport equation in space and angle with the discontinuous finite element and the discrete ordinates method, and utilizes various solving techniques, including mesh sweeping and diffusion acceleration, to provide pebble- wise reaction rates. This work presents the extension of the PTT method to enable fuel depletion capability in the Griffin code. We discuss the implementation details of the PTT-based high-fidelity depletion where isotope inventory of all individual pebbles is tracked through pre-determined pebble flow paths in the core. The implementation is verified with a generic pebble bed reactor model. Some preliminary equilibrium core results are included. Future works are also discussed.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

In situ irradiation of spent nuclear fuels

To improve the economics of commercial nuclear reactors, nuclear vendors and utilities are seeking approval for increased burnup usage of the existing nuclear fleet. This is especially critical for meeting the clean energy initiative by the United States Government, Department of Energy (DOE). However, one of the key challenges the nuclear industry faces in this regard is that the regions exposed to high burnup and low temperatures exhibit a fine-grained microstructure with large bubbles known as high-burnup structure (HBS) [1]. The formation of HBS has been correlated to the diminished performance of the reactor, as well as fuel fragmentation and pulverization during transient and accidental conditions [2]. Therefore, it is paramount to understand the mechanisms for HBS formation along with its impact on the properties and performance of nuclear fuels. While existing programs, such as Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Fuel Campaign (AFC) are focusing on evaluating the performance impact of HBS, the physical mechanisms contributing to HBS formation are still not fully understood. In addition, having predictive capabilities and sound understanding of the microstructural evolution of nuclear fuel is essential for accelerated development, qualification, and deployment of new nuclear materials and novel reactor designs for advanced nuclear reactors. There is a lack of consensus among researchers regarding the mechanisms leading to such restructuring observed in HBS. Grain subdivision due to polygonization versus recrystallization, continuous versus discrete recrystallization occurring in tandem or conjunction, etc., have been proposed and debated. In general, it is hypothesized that defect accumulation and dislocation interaction within the grains cause the realignment of dislocations into grain boundaries, leading to the new subgrain formation, which over time transforms into new grains. However, due to the lack of transient data, the importance of fission rate, irradiation, thermal, and stress history of the fuel on the restructuring could not be assessed. In situ microstructural evolution under various irradiation conditions is desired to bridge this gap. Alternatively, phase-field-based models have been developed to capture HBS formation via discrete recrystallization utilizing the classical nucleation approach [3–5]. However, in these models, the grain nucleation criteria are often defined based on empirical relations for burnup and fission gas density leading to dislocation density change. A mechanistic approach to capture the dislocation interaction with the microstructural features leading to grain subdivision is lacking.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Autonomous Operations for Advanced Reactors Utilizing Supervisory Control

Automation is a critical tenet of reactor plant operations as reliance on nuclear energy increases. Nuclear power plants require a large workforce which does not scale with output; that is, the cost per megawatt increases as reactor output becomes smaller. The economic viability of advanced reactors, particularly small modular reactors (SMRs) and microreactors, requires a significantly reduced onsite workforce. The logical solution is establishing a systematic process of elimination of reliance on human operators, and to the extent possible, replacing these actions with automated functions. In this paper, we propose a method for such transformation to establish a robust technical basis to enable transition to autonomy. Our method is based on finite state automata (FSA)—also known as finite state machines (FSMs). Relying on this method allows us to exploit the rich set of mathematical proofs available in the field of regular languages. FSA are one of the mathematical tools to model discrete event systems (DES). These properties are applied to produce an automated startup controller for the Massachusetts Institute of Technology Research Reactor (MITR). The startup procedure is captured in terms of discrete changes from one state to another while an independent supervisory control system directs the sequence of states and alerts a human in the event of an abnormal operation. First, the design and behavior of the MITR rod control system were modeled in Simulink. Then, the startup procedure was applied to the rod control system and the DES performed a startup by procedurally withdrawing rods to the subcritical position. The simulation also stops rod motion in response to an uncontrollable event and restarts rod motion once the event has been cleared.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2026

Improvements to In-Canal Assay of High Specific Activity 60Co at the ATR

Improvements to the quality of measurements in the assaying of High Specific Activity (HSA) cobalt 60 produced at the Advanced Test Reactor at the Idaho National Laboratory are presented. Such assays are conducted regularly between cycles at the ATR, and due to the radiological safety concerns, they are performed within the ATR canal. A standard ion chamber detector was placed in a collimated assay tube within the canal. Another tube, parallel to the detector, was used to position various cobalt 60 standards and targets. The collimated port between these tubes allows for the targets to be pulled across via a winch system. Previous measurements were performed at discrete 1.5-inch steps using a hand-cranked winch system. Charge was collected and integrated for 10 seconds at each discrete step during a scan. Two substantial improvements to the assay include an electronically controlled, motor-driven winch system and a software controlled continuous charge collection. Combined, these changes improved the quality of measurements as well as streamline the assay process. Further examination of these charge profiles indicated a reasonable amount of photon scattering, or shine, before a line-of-sight is established. This created inaccuracies in the reported activity of the cobalt targets. A series of MCNP-6.3 calculations were performed in an effort to quantify and correct for this shine effect. The results of the calculations and how they were employed to a create shine correction factor for the in-canal cobalt assay at the ATR are presented. The overall bias showed an overestimation of the HSA cobalt 60 activity by around 7-8%.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2026

Residual symmetries and scalar multiplet vacuum alignment in non-Abelian flavour models

We demonstrate that, upon minimizing a renormalizable, single-scalar potential invariant under a non-Abelian symmetry, special orientations in the associated vacuum alignment of the scalar multiplet correspond to the preservation of a discrete residual flavour symmetry in the broken phase of the theory. Conversely, we show that these special scalar alignments are perturbed when additional Lagrangian operators (e.g. renormalizable, multi-flavon operators and/or effective, higher-dimensional operators) are present that break said residual symmetry, leading to a vacuum reorientation and phenomenological consequences. We therefore construct a one-to-one correspondence principle between broken residual symmetries and vacuum alignment corrections, providing a mechanism to identify (and correct) a subtle but persistent form of phenomenologically relevant fine-tuning embedded in — but often ignored by — most successful non-Abelian flavour models. We first establish this correspondence in a set of toy models based on the S4 permutation symmetry, and then apply the lessons learned to the more realistic A4 Altarelli-Feruglio and ∆(27) Universal Texture Zero models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗
DOE OSTI2025

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗
DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2025

Toucan: A performance portable, scalable implementation of the DECA algorithm

In the field of additive manufacturing (AM), cellular automata (CA) is extensively used to simulate microstructural evolution during solidification. However, while traditional CA approaches are relatively fast, they still require a substantial number of time steps, are limited to moderate volumes, and are relatively difficult to improve through parallelism due to the highly localized nature of the solidification front. Here, to address these issues of time to solution and load balancing, we introduce Toucan, a parallel, performance-portable, and scalable code written in C++ with the Kokkos library that leverages the discrete event inspired cellular automata (DECA) algorithm to perform parallel-in-time (PinT) grain growth simulations. Toucan effectively mitigates load balancing issues by distributing the computational workload more evenly across processors, enhancing scalability and efficiency. We conduct both strong and weak scaling studies on up to 64 GPUs on the Frontier supercomputer, demonstrating that Toucan significantly outperforms the current state-of-the-art, time-stepped CA code, ExaCA, on both single and multi-GPU simulations. Even in AM-specific weak scaling scenarios, Toucan maintains near-ideal scaling, in contrast to the linear increase observed with ExaCA due to the moving laser raster pattern. This study highlights Toucan’s potential to transform microstructural simulations in AM by radically improving both efficiency and scalability over existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unified cohesive zone model (UCZM) for fracturing and fragmenting solids

Here, a Unified Cohesive Zone Model (UCZM), which inherits most of the advantages while overcoming the shortcomings of existing Cohesive Zone Models (CZMs), is proposed. Similar to the traditional extrinsic CZM approach, UCZM dynamically inserts the cohesive elements into the system based on local material states (e.g., stress, strain). However, the transition from continua to discontinua is smoothly achieved, thereby eliminating the “time-discontinuous” issue seen in the extrinsic CZM. Moreover, within the novel UCZM framework, the point of transition from continua to discontinua is controllable through the introduction of crack initialization criteria. As a result, the UCZM allows any material models (e.g., elastic, plastic, damage models) for continuum solids and for discrete fracture behavior to work together. In essence, both an enhanced extrinsic cohesive zone model and an intrinsic cohesive zone model can be represented by the proposed unified model. The proposed UCZM has been verified through different numerical examples. The work demonstrates that the UCZM is a highly effective approach for modeling fracture and fragmentation processes in solids.

42 ENGINEERING↗
DOE OSTI2024

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗
DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Effect of 3D-printed surface textures on wear mechanism in 3-body abrasion of soil

This study systematically investigates the enhancement of wear resistance in 3D printed surface textures through both experimental and theoretical approaches. Three distinct surface morphologies (Smooth Surface, Surface with uniformly distributed Pits, and Surface with uniformly distributed Bumps) were fabricated using High-Impact Polystyrene, where the meso-scale textures were precisely controlled through the 3D printing process. Wear behavior was evaluated using a 3-body wear tester in an abrasive particle environment, analyzing the influence of surface textures under various operating conditions. Systematic wear tests revealed that optimally designed surface textures achieved a remarkable 77 % reduction in wear compared to the worst-performing sample. The wear mechanisms were comprehensively characterized through weight loss measurements, Scanning Electron Microscopy (SEM), and Energy Dispersive Spectroscopy (EDS) analyses, elucidating the surface morphology changes and their interaction with wear particles. Notably, the study identified how the geometric characteristics of surface textures influence the movement of wear particles and the distribution of contact stresses. Discrete element method simulations corroborated the experimental findings, providing theoretical validation for the enhanced wear resistance of the optimal structure. The high correlation between simulated wear patterns and experimental results validates the reliability of the proposed design methodology. In conclusion, these results demonstrate that 3D printed surface texturing offers a cost-effective and scalable approach to significantly improve wear resistance in engineering applications, presenting a practical alternative to conventional, high-cost surface engineering methods.

3-body abrasion↗
DOE OSTI2026

Counting topological interface modes using simplicial characteristic classes

A computational approach for predicting the number of topological interface modes (TIMs) in hermitian systems using the spectral flow—monopole correspondence is presented. The number of TIMs is determined by calculating the Chern number of a complex line bundle of local polarisation vectors over a phase space sphere surrounding a Weyl point. The Chern number is computed by constructing the simplicial first Chern class of a discrete vector bundle on a simplicial mesh. This approach is gauge invariant, derivative free, structure preserving, and robust to noise. The algorithm is shown to reproduce the expected number of TIMs for the case of equatorial fluid waves and the topological Langmuir cyclotron wave. The possibility of using this algorithm to analyse experimental measurements of bulk wave polarisations and predict the associated number of TIMs is explored in a synthetic example.

discrete vector bundles↗
DOE OSTI2024

Evolution of Efimov states

The Efimov phenomenon manifests itself as an emergent discrete scaling symmetry in the quantum three-body problem. In the unitarity limit, it leads to an infinite tower of three-body bound states with energies forming a geometric sequence. Here in this work, we study the evolution of these so-called Efimov states using relativistic scattering theory. We identify them as poles of the three-particle S matrix and trace their trajectories in the complex energy plane as they evolve from virtual states through bound states to resonances. We dial the scattering parameters toward the unitarity limit and observe the emergence of the universal scaling of energies and couplings - a behavior known from the non-relativistic case. Interestingly, we find that Efimov resonances follow unusual, cyclic trajectories accumulating at the three-body threshold and then disappear at some values of the two-body scattering length. We propose a partial resolution to this ?missing states? problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗