Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete element modeling”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

MOOSE framework enhancements for meshing reactor geometries

MOOSE is an open-source, parallel finite element framework designed to permit rapid development of robust multi-physics modeling capabilities. Under the DOE-NEAMS program, numerous solvers have been developed utilizing the open-source MOOSE framework for multiphysics reactor analysis. These solvers require input finite element meshes representing the discretized geometry. Typically, reactor analysts turn to licensed external tools for creation of reactor geometry meshes. Recently, enhancements have been added to the MOOSE framework to mesh common reactor geometries and improve MOOSE-based application user workflows. Support for hexagonal pins, assemblies, and cores has been added, and Cartesian support has been extended. Options for modeling static and rotating control drums within a hexagonal assembly are now available. Pin, assembly, and plane regions can be identified through automatically applied tags on the mesh called 'reporting IDs' for easier post- processing of physics results. An external open-source triangle routine has been leveraged within MOOSE to mesh core periphery zones. A set of reactor geometry builder routines further streamline the construction of hexagonal and Cartesian cores and include the ability to assign materials to regions during mesh generation. The new meshing routines are available through the MOOSE framework in the open-source 'Reactor' module, and the resulting directly within MOOSE-based applications or exported as Exodus II files for use in other finite element solvers. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2024

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗
DOE OSTI2024

Investigation of thermal damage in explosive bridgewire detonators via discrete element method simulations

Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. In conclusion, we present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

cohesive powders↗
DOE OSTI2024

Performance Portable Graphics Processing Unit Acceleration of a High-Order Finite Element Multiphysics Application

The Lawrence Livermore National Laboratory (LLNL) will soon have in place the El Capitan exascale supercomputer, based on advanced micro devices (AMD) graphics processing units (GPUs). As part of a multiyear effort under the National Nuclear Security Administration (NNSA) Advanced Simulation and Computing (ASC) program, we have been developing marbl, a next generation, performance portable multiphysics application based on high-order finite elements. In previous years, we successfully ported the Arbitrary Lagrangian–Eulerian (ALE), multimaterial, compressible flow capabilities of marbl to nvidia GPUs as described in Vargas et al. Here, in this paper, we describe our ongoing effort in extending marbl's GPU capabilities with additional physics, including multigroup radiation diffusion and thermonuclear burn for high energy density physics (HEDP) and fusion modeling. We also describe how our portability abstraction approach based on the raja Portability Suite and the mfem finite element discretization library has enabled us to achieve high performance on AMD based GPUs with minimal effort in hardware-specific porting. Throughout this work, we highlight numerical and algorithmic developments that were required to achieve GPU performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

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

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sierra/SolidMechanics 5.2 Capabilities in Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.2 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

42 ENGINEERING↗
DOE OSTI2022

Sierra/SolidMechanics 5.4 Capabilities in Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.4 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Sierra/SolidMechanics 5.6 Capabilities in Development

This user's guide documents capabilities in Sierra/SolidMechanics which remain "in-development" and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.6 User's Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sierra/SolidMechanics 5.8 In-Development Manual

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.8 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

42 ENGINEERING↗
DOE OSTI2022

Sierra/SolidMechanics 5.10 In-Development Manual

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.10 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

42 ENGINEERING↗
DOE OSTI2023

Sierra/SolidMechanics 5.16 Capabilities in Development Manual

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.16 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Accelerating computational modeling and design of high-entropy alloys

High-entropy alloys, with N elements and compositions {$c_{ν = 1,N}$} in competing crystal structures, have large design spaces for unique chemical and mechanical properties. In this work, to enable computational design, we use a metaheuristic hybrid Cuckoo search (CS) to construct alloy configurational models on the fly that have targeted atomic site and pair probabilities on arbitrary crystal lattices, given by supercell random approximates (SCRAPs) with S sites. Our Hybrid CS permits efficient global solutions for large, discrete combinatorial optimization that scale linearly in a number of parallel processors, and linearly in sites S for SCRAPs. For example, a four-element, 128-site SCRAP is found in seconds—a more than 13,000-fold reduction over current strategies. Our method thus enables computational alloy design that is currently impractical. We qualify the models and showcase application to real alloys with targeted atomic short-range order. Being problem-agnostic, our Hybrid CS offers potential applications in diverse fields.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Benchmarking numerical methods for impact and cratering applications

Large scale computational models are important for studying impact cratering events that are prevalent both on Earth and, more broadly, in this solar system. To address these problems, models must reliably account for both large length scales (e.g., kilometers) and relatively long time scales (hundreds of seconds). This work benchmarks two such approaches, a more traditional hydrodynamics approach and a finite-discrete element method (FDEM), for impact cratering applications. Both 2D and 3D results are discussed for two different impact velocities, 5 km/s and 20 km/s, striking normal to the target and, for 3D simulations, 45° from vertical. In addition, comparisons to previously published data are presented. Finally, differences in how these methods model damage are discussed. Ultimately, both approaches show successful modeling of several different impact scenarios.

58 GEOSCIENCES↗
DOE OSTI2024

Modeling the effect of wind speed and direction shear on utility‐scale wind turbine power production

Abstract Wind speed and direction variations across the rotor affect power production. As utility‐scale turbines extend higher into the atmospheric boundary layer (ABL) with larger rotor diameters and hub heights, they increasingly encounter more complex wind speed and direction variations. We assess three models for power production that account for wind speed and direction shear. Two are based on actuator disc representations, and the third is a blade element representation. We also evaluate the predictions from a standard power curve model that has no knowledge of wind shear. The predictions from each model, driven by wind profile measurements from a profiling LiDAR, are compared to concurrent power measurements from an adjacent utility‐scale wind turbine. In the field measurements of the utility‐scale turbine, discrete combinations of speed and direction shear induce changes in power production of −19% to +34% relative to the turbine power curve for a given hub height wind speed. Positive speed shear generally corresponds to over‐performance and increasing magnitudes of direction shear to greater under‐performance, relative to the power curve. Overall, the blade element model produces both higher correlation and lower error relative to the other models, but its quantitative accuracy depends on induction and controller sub‐models. To further assess the influence of complex, non‐monotonic wind profiles, we also drive the models with best‐fit power law wind speed profiles and linear wind direction profiles. These idealized inputs produce qualitative and quantitative differences in power predictions from each model, demonstrating that time‐varying, non‐monotonic wind shear affects wind power production.

Energy & Fuels↗