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 37 records · Page 2

Granular bed consolidation, creep, and armoring under subcritical fluid flow

In this work, we show that a freshly sedimented granular bed settles and creeps forward over extended periods of time under an applied hydrodynamic shear stress which is below the critical value for bedload transport. The rearrangements are found to last over a timescale which is millions of times the sedimentation timescale of a grain in the fluid. Compaction occurs uniformly throughout the bed, but creep is observed to decay exponentially with depth and decreases over time. The granular volume fraction in the bed is found to increase logarithmically, saturating at the random close packing value Φ rcp ≈ 0.64, while the surface roughness is observed to remain essentially unchanged. We demonstrate that an increasingly higher shear stress is required to erode the bed, after a subcritical shear is applied, which results in an increase in its volume fraction. Thus, we find that bed armoring occurs due to a deep shear-induced relaxation of the bed toward the volume fraction associated with the glass transition.

58 GEOSCIENCES↗

Continuum modeling of rate-dependent granular flows in SPH

In this paper, we discuss a constitutive law for modeling rate-dependent granular flows that has been implemented in smoothed particle hydrodynamics (SPH). We model granular materials using a viscoplastic constitutive law that produces a Drucker–Prager-like yield condition in the limit of vanishing flow. A friction law for non-steady flows, incorporating rate-dependence and dilation, is derived and implemented within the constitutive law. We compare our SPH simulations with experimental data, demonstrating that they can capture both steady and non-steady dynamic flow behavior, notably including transient column collapse profiles. In conclusion, this technique may therefore be attractive for modeling the time-dependent evolution of natural and industrial flows.

42 ENGINEERING↗

Shear rate dependency on flowing granular biomass material

The commercialization of bioenergy has been significantly limited by various material handling issues due to the poor flowability of granular biomass materials. A good understanding of flow physics and robust constitutive models to predict flow behavior across multiple regimes are essential to address these issues. In this study, we investigated the multi-regime flow behavior of loblolly pine chips, a widely used bioenergy feedstock, through comprehensive inclined plane flow experiments and simulations. A quasi-static hypoplastic model and a cross-regime Drucker-Prager-µ(I) model were calibrated and validated against the physical experiments to investigate the quasi-static and dense flow behavior. The results show that for granular biomass, 1) plane flow (iso-thickness along the plane) exists within a smaller range of inclination while heap flow (varying thickness along the plane) exists within a broader range of inclination, as compared with conventional granular materials (e.g., glass beads); 2) the scaling law of granular biomass flowing on an inclined plane (Froude Number versus dimensionless thickness) forms a bi-linear trend with the turning point governed by the quasi-static and dense flow regimes; 3) the multi-regime DP-µ(I) model can capture the flow behavior in both regimes well at the cost of extra calibration. In conclusion, these findings advance the scientific understanding of the multi-regime flow behavior of granular biomass materials and shed light on formulating novel constitutive models to assist granular biomass handling in the bioenergy industry.

09 BIOMASS FUELS↗

Assessment of Models of Chemically Reacting Granular Flows

A report presents an assessment of a general mathematical model of dense, chemically reacting granular flows like those in fluidized beds used to pyrolize biomass. The model incorporates submodels that have been described in several NASA Tech Briefs articles, including "Generalized Mathematical Model of Pyrolysis of Biomass" (NPO-20068) NASA Tech Briefs, Vol. 22, No. 2 (February 1998), page 60; "Model of Pyrolysis of Biomass in a Fluidized-Bed Reactor" (NPO-20708), NASA Tech Briefs, Vol. 25, No. 6 (June 2001), page 59; and "Model of Fluidized Bed Containing Reacting Solids and Gases" (NPO- 30163), which appears elsewhere in this issue. The model was used to perform computational simulations in a test case of pyrolysis in a reactor containing sand and biomass (i.e., plant material) particles through which passes a flow of hot nitrogen. The boundary conditions and other parameters were selected for the test case to enable assessment of the validity of some assumptions incorporated into submodels of granular stresses, granular thermal conductivity, and heating of particles. The results of the simulation are interpreted as partly affirming the assumptions in some respects and indicating the need for refinements of the assumptions and the affected submodels in other respects.

Bellan, Josette↗

A study in three-dimensional chaotic dynamics: Granular flow and transport in a bi-axial spherical tumbler

We study three-dimensional (3D) chaotic dynamics through an analysis of transport in a granular flow in a half-full spherical tumbler rotated sequentially about two orthogonal axes (a bi-axial “blinking” tumbler). The flow is essentially quasi-two-dimensional in any vertical slice of the sphere during rotation about a single axis, and we provide an explicit exact solution to the model in this case. Hence, the cross-sectional flow can be represented by a twist map, allowing us to express the 3D flow as a linked twist map (LTM). We prove that if the rates of rotation about each axis are equal, then (in the absence of stochasticity) particle trajectories are restricted to two-dimensional (2D) surfaces consisting of a portion of a hemispherical shell closed by a “cap''; if the rotation rates are unequal, then particles can leave the surface they start on and traverse a volume of the tumbler. The period-one structures of the governing LTM are examined in detail: analytical expressions are provided for the location of period-one curves, their extent into the bulk of the granular material, and their dependence on the protocol parameters (rates and durations of rotations). Exploiting the restriction of trajectories to 2D surfaces in the case of equal rotation rates about the axes, a method is proposed for identifying and constructing 3D Kolmogorov--Arnold--Moser (KAM) tubes around the normally elliptic period-one curves. The invariant manifold structure arising from the normally hyperbolic period-one curves is also examined. When the motion is restricted to 2D surfaces, the structure of manifolds of the hyperbolic points in the bulk differs from that corresponding to hyperbolic points in the flowing layer. Each is reminiscent of a template provided by a non-integrable perturbation to a Hamiltonian system, though the governing LTM is not. This highlights the novel 3D chaotic behaviors observed in this model dynamical system.

42 ENGINEERING↗

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↗

BDEM (Discrete-element-simulator for high-solids granular flows) [SWR-22-72]

BDEM is a discrete element method based simulation tool developed specifically for modeling high-solids granular flows that include polydispersity, heat-transfer, moving boundaries and chemistry. Our solver provides facilities for simulating spherical/non-spherical particles with modified contact and friction models in complex dynamic geometries defined using level-sets or triangulated files. The solver is developed on top of NREL's open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Simulations at the scale of several millions to billion particles have been performed using this software on large scale computing resources. This software can be applied to non-reacting solids dominant flows in silos, hoppers and screw conveyors as well as in high temperature reacting systems such as screw kilns and auger reactors.

Sitaraman, Hariswaran↗

Shear dispersion in dense granular flows

We formulate and solve a model problem of dispersion of dense granular materials in rapid shear flow down an incline. The effective dispersivity of the depth-averaged concentration of the dispersing powder is shown to vary as the Péclet number squared, as in classical Taylor–Aris dispersion of molecular solutes. An extension to generic shear profiles is presented, and possible applications to industrial and geological granular flows are noted.

42 ENGINEERING↗

Experimental demonstration of a 5 kW th granular-flow reactor for solar thermochemical energy storage with aluminum-doped calcium manganite particles

A two-step cycle was discussed for solar thermochemical energy storage based on aluminum-doped calcium manganite reduction/oxidation reactions for direct integration into Air Brayton cycles. The two steps encompassed 1) the storage of concentrated solar direct irradiation via the thermal reduction of aluminum-doped calcium manganite and 2) the delivery of heat to an Air-Brayton cycle via re-oxidation of oxygen-deficient aluminum-doped calcium manganite. The re-oxidized aluminum-doped calcium manganite was fed back to the first step to complete the cycle. A 5 kWth solar thermochemical reactor operating under vacuum was fabricated and tested to examine the first cycle reduction step. Reactor operating conditions and high-flux solar simulator control were tuned for continuous reactor operation with particle temperatures > 1073 K. Continuous operation was achieved using intermittent, dense granular flows. A maximum absorption efficiency of 64.7% was shown, accounting for both sensible and chemical heat storage.

25 ENERGY STORAGE↗

Granular Flow Constitutive Model Add-ons For Abaqus

The software is a list of Abaqus User MATerial subroutines (VUMAT) for modeling granular flow physics. It includes (1) density dependent Mohr-Coulomb model, (2) density dependent Drager-Prager\Cap model, (3)Gudehus-Bauer hypoplastic model, and (4) critical state-based NorSand model.

Wencheng, Wencheng↗

SPH modeling of biomass granular flow: Theoretical implementation and experimental validation

The commercialization of biomass-derived energy is impeded by flowability challenges arising from the feeding and handling of granular biomass materials in full-scale biorefineries. To overcome these obstacles, a robust and accurate model to simulate the flow of granular biomass is indispensable. However, conventional mesh-based numerical codes are limited by inherent mesh distortion in simulating large deformation that commonly occurs in granular biomass handling. Here, in this study, we propose a graphics processing unit (GPU)-accelerated meshless Smoothed Particle Hydrodynamics (SPH) code to model the flow of granular biomass materials. A modified void ratio-based mass conversation, a hybrid particle-to-particle/surface frictional boundary treatment, and a hypoplastic constitutive model are implemented. Four numerical examples, an elastic block sliding on inclined planes, sand column collapse, Angle of Repose, and axial compression tests for pine chips, were simulated using the developed SPH code. The results demonstrate good agreement between numerical predictions and analytical and experimental data for all four examples, validating the SPH code and increasing confidence that it can be applied to simulate more complex granular biomass handling processes, such as hopper feeding or auger conveyance.

09 BIOMASS FUELS↗

Granular flow in pebble bed reactors: Dust generation and scaling

In experimental prototypes of pebble bed reactors, significant quantities of graphite dust have been observed due to rubbing between pebbles as they flow through the core. At the high temperatures and pressures in these reactors, little data is available to understand the frictional properties of the pebble surfaces, and as a result, the Paul Scherrer Institut (Switzerland) proposes a conceptual design of a scaled-down version of a pebble bed reactor to investigate this issue in detail. In this paper, simulations of granular flow in pebble bed reactors using the discrete-element method are presented. Simulations in the full geometry (using 440,000 pebbles) are compared to those in geometries scaled down by 3:1 and 6:1. The simulations show complex behavior due to discrete pebble packing effects, meaning that pebble flow and dust generation in a scaled-down facility may be significantly different. The differences between velocity profiles, packing geometry, and pebble wear at the different scales are discussed. The results can aid in the design of the prototypical facility to more accurately reproduce the flow in a full-size reactor. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Initial Validation of a Gas-Granular Flow Solver Using a Subscale, Reduced Pressure Plume Surface Interaction Ground Test

With NASA’s goal to land the next humans on the lunar surface in the next few years, it is vitally important to have a better understanding of the plume surface interaction (PSI) between the landing vehicles and the lunar regolith. The Fluid Dynamics Branch at NASA/MSFC has previously used the gas-granular flow solver Loci/GGFS to qualitatively predict crater formation due to PSI effects in a lunar (near vacuum) ambient environment. In this paper, validation of Loci/GGFS crater width and depth predictions in ambient near-lunar conditions are provided using experimental data collected at MSFC during the Physics-Focused Ground Test 1 (PFGT-1) campaign in 2022. To observe sensitivity to soil models, simulations were conducted with both monodisperse glass bead (MGB) and BP-1 lunar regolith simulant soil models in Loci/GGFS. Crater depth and width comparisons are made with PFGT-1 Run 56, which used BP-1 soil. The Loci/GGFS BP-1 soil model performed slightly better with a mean predicted crater depth within 10% of the experiment. Both soil models predicted crater width within 10%. Cratering occurred more quickly with the MGB soil model. Mesh and spatial order sensitivity are also examined for the BP-1 soil model.

Validation↗

Validation Assessment of Loci/GGFS Gas Granular Flow Solver Predictions of Ejecta from Physics Focused Ground Test

NASA is preparing to return humans to the Moon to establish a sustained Lunar presence through the Artemis program. One area of concern for Lunar landers is the plume surface interaction (PSI) environment that poses several risks during a propulsive landing. Understanding the PSI environment caused by the landers is therefore important to designing successful landing missions. Towards this end, a two-phase, gas granular flow solver, Loci/GGFS, has been developed to predict the cratering and ejecta physics expected during a Lunar landing. The focus of this paper is on a validation assessment of Loci/GGFS predictions of ejecta with monodisperse glass beads (MGB). The validation assessment is performed with respect to the Physics Focused Ground Test 1 (PFGT1) conducted at NASA Marshall Space Flight Center (MSFC) in 2021. It is found that Loci/GGFS predicts similar initial cratering and ejecta features, with predictions of average velocities that are within 25% of experimental values.

Plume Surface Interaction↗

Validation Assessment of Loci/GGFS Gas Granular Flow Solver Predictions of Ejecta From Physics Focused Ground Test

NASA is preparing to return humans to the Moon to establish a sustained Lunar presence through the Artemis program. One area of concern for Lunar landers is the plume surface interaction (PSI) environment that poses several risks during a propulsive landing. Understanding the PSI environment caused by the landers is therefore important to designing successful landing missions. Towards this end, a two-phase, gas granular flow solver, Loci/GGFS, has been developed to predict the cratering and ejecta physics expected during a Lunar landing. The focus of this paper is on a validation assessment of Loci/GGFS predictions of ejecta with monodisperse glass beads (MGB). The validation assessment is performed with respect to the Physics Focused Ground Test 1 (PFGT1) conducted at NASA Marshall Space Flight Center (MSFC) in 2021. It is found that Loci/GGFS predicts similar initial cratering and ejecta features, with predictions of average velocities that are within 25% of experimental values.

Plume Surface Interaction↗

Granular Flow in Novel Octet Shape–Based Lattice Frame Material

Moving packed-bed heat exchangers in concentrated solar power (CSP) plants involves heat transfer between heated falling particles and supercritical carbon dioxide. The overall effective thermal conductivity of the moving packed bed and particle-side channel contact resistances are still the bottlenecks in achieving the desirable thermal transport levels. To this end, a novel moving packed bed heat exchanger consisting of an Octet lattice packed between the walls of the particle-side channel is proposed in this study. Granular flow analysis in Octet lattice moving packed bed heat exchanger (OLHX) was conducted through experiments and discrete element method (DEM)-based numerical simulations. The experimental images clearly demonstrated stagnation regions upstream of lattice fibers, void regions downstream of the fiber junctions, and wavy-type unobstructed flow on the lateral sides of the fibers. DEM simulations were successful in capturing all these critical flow phenomena. Larger flow velocities were observed on the lateral sides of the fibers in the simulations. Also, when the particles in the silo were emptied, the final images showed an accumulation of particles on the inter-fiber as well as fiber–channel wall junctions. Moreover, the fiber connections resulted in some regions devoid of particle contact on the channel endwall, which means that these regions would suffer from poor thermal exchange. Finally, the overall mass flowrate increased with increasing porosity for a fixed particle diameter.

14 SOLAR ENERGY↗

Modified kinetic theory applied to the shear flows of granular materials

Here, granular materials are characterized by large collections of discrete particles, where the particle-particle interactions are significantly more important than the particle-fluid interactions. The current kinetic theory captures fairly accurately the granular flow behavior in the dilute case, when only binary interactions are significant, but is not accurate at all in the dense flow regime, where multi-particle interactions and contacts must be modeled. To improve the kinetic theory results for granular flows in the dense flow regime, we propose a Modified Kinetic Theory (MKT) model that utilizes the contact duration or cut-off time to account for the complex particle-particle interactions in the dense regime. The contact duration model, also called TC model, is originally proposed by Luding and McNamara to solve the inelastic collapse issue existing in the Inelastic Hard Sphere (IHS) model. This model defines a cut-off time t c such that dissipation is not counted if the time between two consecutive contacts is less than t c . As shown in their study, the use of a cut-off time t c can also reduce the dissipation during multi-particle contacts. In this paper we relate the TC model with the Discrete Element Method (DEM) by choosing the cut-off time t c to be the duration of contact calculated from the linear-spring-dashpot soft-sphere model of the DEM. We examine two types of granular flows: simple shear flow and the plane shear flow, and compare the results of the classical Kinetic Theory (KT) model, the present MKT model, and the DEM model. Here, we show that the MKT model entails a significant improvement over the KT model for simple shear flows at inertial regimes. With the MKT model the calculations are close to the DEM results at solid fractions as high as 0.57. Even for the plane shear flows, where shear rate and solid fraction are inhomogeneous, the results of the MKT model agree very well with the DEM results.

36 MATERIALS SCIENCE↗