Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Granular flow”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Pseudo-viscous modeling of transport in dense granular flows for thermal energy storage applications

Dense, granular flows were examined to effectively capture and model bulk viscous properties in thin packed beds. A modified Couette cell with particle image velocimetry was used to experimentally determine pseudo-viscosity properties of four particulate media with varying morphologies: (1) iron oxide-coated SiO 2 particles, (2) CARBOBEAD CP30-60 particles, (3) CARBOBEAD CP40-100 particles, and (4) Al 2 O 3 beads. The pseudo-viscosity functions were fitted using a power law to correlate the measured shear stress as a function of measured shear rate. The pseudo-viscous functions were used as inputs to computation fluid dynamics models for a single-phase viscous fluid to predict granular flow profiles. Steady-state free surface velocity profiles at angular velocities <7 rad/s predicted by the model were in good agreement with the experimental particle image velocimetry measurements, resulting in Pearson correlation coefficients of 0.97 for iron-oxide coated SiO 2 particles and 0.95 for CP30-60 particles. As a result, this alternative approach to measuring pseudo-viscous properties under shearing and modeling bulk transport behavior of granular flow using computation fluid dynamics model offered significant reduction in computational load compared to discrete element methods.

14 SOLAR ENERGY↗
DOE OSTI2025

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗
DOE OSTI2024

Impact of Janssen effect on thermal transport in granular flow

Using a modulated photothermal radiometry (MPR) technique capable of measuring thermal transport in particle beds, we find asymptotically increasing effective thermal conductivity and decreasing near-wall thermal resistance along gravity-driven downward granular flows through a 1-meter-long narrow vertical channel owing to the Janssen effect. The Janssen effect is confirmed by similar asymptotic trends from particle bed apparent mass measurements as well as separate MPR heat transfer measurements on stationary particle beds under compression at different heights of the channel with and without pressure screening. Furthermore, local heat transfer coefficient along the 1-meter-long channel was modeled based on measured spatially resolved thermal conductivity and near-wall thermal resistance due to the Janssen effect. This work reveals for the first time the impact of the over century-old Janssen effect on heat transfer in granular media. Furthermore, the results from this work can lead to a better understanding of heat transfer in granular flow.

14 SOLAR ENERGY↗
DOE OSTI2025

Dense granular flows with MFIX-Exa

This report extends the linear spring dashpot collision model the discrete element method available in MFIX-Exa to include static a static tangential friction force. Additionally, two rolling friction models frequently used in the literature are also implemented. The governing equations are provided with an emphasis on the new terms. The new model is validated by comparison to existing experimental data of single particle oblique collisions. The model is then tested on three dense granular flow problems: the formation of static piles, the discharge from a flat-bottom hopper and the self-induced granular Rayleigh-Taylor instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Continuum Correlations from CFD-DEM Modeling of Conduction Heat Transfer in Granular Flows

Heat transfer between a surface and flowing particles is analyzed to improve the accuracy of continuum models for wall-to-bed heat transfer in a fluidized bed. Discrete element modeling (DEM) is used to model a fluidized bed heat exchanger where heat enters the system through a heated wall. The DEM heat transfer predictions are validated against published experimental work (Brewster et al., 2024) with less than 15% error. In previous work by Morris et al. (2015), a continuum model was developed using data from high-fidelity DEM simulations of chute flows. In the current study, the continuum model is extended and validated for fluidized beds. The sensitivity of the continuum heat transfer model parameters, which was not quantified in previous studies, is also investigated. It is observed that for a given particle with specific properties, e.g. the particle size, roughness, and conduction lens radius, the continuum correlation developed for heat transfer from a heated boundary to the particle bed depends mainly on the solid fraction or porosity of the particle bed for a given fluid. The new continuum heat transfer model is then validated over a wide range of superficial velocities via comparisons to both discrete element and experimental data. It is shown that this correlation is valid for a large range of particle flow conditions from chute flows to fluidized beds with less than 10% error as compared to DEM predictions.

14 SOLAR ENERGY↗
DOE OSTI2024

Heat transfer coefficients of moving particle beds from flow-dependent thermal conductivity and near-wall resistance

Accurate determination of heat transfer coefficients for flowing packed particle beds is essential to the design of particle heat exchangers and other thermal and thermochemical equipment. While such dense granular flows mostly fall into the well-known plug-flow regime, the discrete nature of granular materials alters the thermal transport processes in both the near-wall and bulk regions of flowing particle beds from their stationary counterparts. As a result, heat transfer correlations based on the stationary particle bed thermal conductivity could be inadequate for flowing particles in a heat exchanger. Most earlier works have achieved a reasonable agreement with experiments by treating granular heat transfer media as a plug-flow continuum with a near-wall thermal resistance in series. However, the thermal conductivity values of the continuum were often obtained from measurements on stationary beds owing to the difficulty of flowing bed measurements. In this work, it was found that the properties of a stationary bed are highly sensitive to the method of particle packing and there is a decrease in the particle bed thermal conductivity and increase in the near-wall thermal resistance, measured as an effective air gap thickness, on the onset of particle flow. These variations in thermal conductivity of stationary and flowing particle beds can lead to errors in heat transfer coefficient calculations. Therefore, the heat transfer coefficients for granular flows were calculated using experimentally determined flowing particle bed thermal conductivity and near-wall air gap for ceramic particles – CARBO CP 40/100 (mean diameter = 275 µm), HSP 40/70 (404 µm) and HSP 16/30 (956 µm); at velocities of 5–15 mm·s –1 ; and temperatures of 300–650 °C. The thermal conductivity and air gap values for CP 40/100 and HSP 40/70 were further used to calculate heat transfer coefficients across different particle bed temperatures and velocities for different parallel-plate heat exchanger dimensions. Furthermore, these calculations, which show good agreement with measured HTC values reported in literature, can be used as a guide for heat exchanger designs. Graphical abstract

14 SOLAR ENERGY↗
DOE OSTI2026

Reverse segregation and self-organization in inclined chute flows of bidisperse granular mixtures

In the usual segregation scenario for stable inclined chute flows of bidisperse mixtures of fine and coarse spherical particles, coarse particles rise toward the free surface, forming a coarse-rich region atop the flowing pile. Beyond a threshold coarse-to-fine diameter ratio of approximately 4, conversely, the weight of the coarse particles exceeds the segregation driving forces, causing individual coarse particles to sink within the pile and producing a reversed segregation state. However, an understanding of the collective evolution of the pile structure is still lacking when the particle diameter ratio exceeds 4 and the coarse-particle mass fraction is appreciable. To explore this broadly bidisperse limit, we perform discrete element method simulations considering mean particle diameter ratios of up to 8 and coarse-particle mass fractions spanning 0.1 to 0.9. The steady-state flow profiles reveal several intriguing behaviors that depend on the diameter ratio and mass fraction. These include a previously identified transition from usual to reverse segregation and a newfound tendency to self-organize into alternating coarse- and fine-rich particle layers stacked along the shear gradient direction, with layer thickness dictated by the coarse-particle diameter. A fuller understanding of segregation at this scale could pave the way for enhanced mixing or demixing techniques at the commercial scale.

granular flow↗
DOE OSTI2024

Triboplasma assisted chemical conversion in granular systems: a semi-analytic model

We present a semi-analytic model for a novel plasma-assisted chemical conversion pathway using triboplasmas generated in granular flows. Triboelectric charge relaxation is a well known phenomena where the potential generated from contact charging of particles exceeds the breakdown voltage of the background gas. In this work, we extend the triboelectric charge relaxation theory to include non equilibrium plasma energy and particle balance equations to predict the formation of dissociated and excited species that act as precursors to chemical conversion, for example in plasma-assisted ammonia synthesis. Our example case study with nitrogen background gas and teflon/aluminum tribomaterial system yielded high excited nitrogen species densities per collision that are comparable to current plasma-assisted conversion pathways. We also present a regime diagram for various gases where Paschen breakdown parameters are used to determine whether triboplasmas can be formed for a given effective work-function difference between two materials. Our sensitivity studies indicate particle velocity, particle radius, solids fraction and space charge effects play a critical role in overall plasma densities and excited species production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Machine Learning Aided Modeling of Granular Materials: A Review

Artificial intelligence (AI) has become a buzzy word since Google’s AlphaGo beat a world champion in 2017. In the past five years, machine learning as a subset of the broader category of AI has obtained considerable attention in the research community of granular materials. This work offers a detailed review of the recent advances in machine learning-aided studies of granular materials from the particle-particle interaction at the grain level to the macroscopic simulations of granular flow. This work will start with the application of machine learning in the microscopic particle-particle interaction and associated contact models. Then, different neural networks for learning the constitutive behaviour of granular materials will be reviewed and compared. Finally, the macroscopic simulations of practical engineering or boundary value problems based on the combination of neural networks and numerical methods are discussed. We hope readers will have a clear idea of the development of machine learning-aided modelling of granular materials via this comprehensive review work.

42 ENGINEERING↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2025

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗
DOE OSTI2024

Implementation of the Glued Sphere Discrete Element Model for Non-Spherical Particles in MFiX Software

To enhance solver capabilities, simulation flexibility and model validation within the MFiX software, the U.S. Department of Energy (DOE) is funding efforts to develop and integrate the glued-sphere discrete element method into the latest version of MFiX as a dedicated computational module. The glued-sphere discrete element method is a numerical technique to depict the behavior of non-spherical particles in granular flows or particulate systems by representing them as a collection of component spheres. These spheres are bonded together to approximate the shape and mechanical/chemical properties of a more complex particle. The method effectively reuses the existing sphere-sphere collision algorithm, interphase momentum and heat transfer calculations utilized in the traditional discrete element method, extending these capabilities to non-spherical particles. Additionally, this method explicitly resolves intra-particle temperature and species distributions. The MFiX glued-sphere computational module includes tools for generating glued sphere configurations, a dedicated solver, and visualization capabilities in post-processing. More specifically within the computational module, collision detection and calculations were first performed on component spheres and then mapped onto non-spherical particles. The linear spring-dashpot model was utilized to simulate the sphere-sphere interactions.

Ke, Renjie↗
DOE OSTI2025

Using Solid Particles as Heat Transfer Fluid-Use in CSP Plants (CRADA Final Report)

NREL assisted the BRIDGE project led by the University of Colorado (CU), Boulder, on developing granular flow and heat transfer simulation tool to serve for an enclosed particle-receiver design. NREL will provide the particle receiver design configuration and operating condition to assist the model setup and may provide certain testing results for model verifications. NREL will attend project reviewing meeting and provide application back-ground including concentrating solar energy technology to collaborate with CU researchers.

14 SOLAR ENERGY↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2024

Self-assembled reconfigurable pump architectures via magnetic colloidal swarms

Self-assembled swarms of interactive active units, which are adaptive and dynamically reconfigurable to accommodate different functionalities, represent a promising platform for the development of next-generation robotics. Here, we utilize the emergent collective behavior of active magnetic colloids confined in quasi-two-dimensional arrays of overlapping wells to demonstrate the self-organization of a colloidal swarm into a dynamic pump architecture capable of controlled transport of passive cargo particles. This dynamic architecture provides a global unidirectional looping flow pattern along the entire length of the system. We show that the flow direction of the dynamic swarm-based pump can be externally controlled by a phase shift of a driving magnetic field energizing the swarm. The experimental observations are supported by computational modeling based on phenomenological coarse-grained particle dynamics coupled to shallow-water Navier-Stokes hydrodynamics. In conclusion, our findings demonstrate how the emergent collective behavior of a swarm can be orchestrated into a desired functionality by exploiting the interplay between activity and confinement potentials.

36 MATERIALS SCIENCE↗