Search NASA⌕ Search

SEARCH · Search NASA

Results for “granular flow 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗
DOE OSTI2024

Convolutional L2LFlows: generating accurate showers in highly granular calorimeters using convolutional normalizing flows

Abstract In the quest to build generative surrogate models as computationally efficient alternatives to rule-based simulations, the quality of the generated samples remains a crucial frontier. So far, normalizing flows have been among the models with the best fidelity. However, as the latent space in such models is required to have the same dimensionality as the data space, scaling up normalizing flows to high dimensional datasets is not straightforward. The prior L2LFlows approach successfully used a series of separate normalizing flows and sequence of conditioning steps to circumvent this problem. In this work, we extend L2LFlows to simulate showers with a 9-times larger profile in the lateral direction. To achieve this, we introduce convolutional layers and U-Net-type connections, move from masked autoregressive flows to coupling layers, and demonstrate the successful modelling of showers in the ILD Electromagnetic Calorimeter as well as Dataset 3 from the public CaloChallenge dataset.

Instruments & Instrumentation↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part I: Discrete-Particle Models

Biomass is a renewable and sustainable energy resource. Current design of biomass handling and feeding equipment leverage both experiments and numerical modeling. This paper reviews the state-of-the-art discrete element methods (DEM) for the flow of milled biomass (Part I), accompanied by a comprehensive review on continuum-based computational models (Part II). The present review on DEM is primarily focused on the features and suitability of various particle shape models for different types of milled biomass because particle shape is the predominant attribute controlling the flow behavior of complex-shaped granular material. The general strengths and weaknesses in the applicability of those models for the milled biomass modeling are summarized. In particular, comments are provided to balance the numerical model capabilities and the computational cost for the development of DEM models. To our best knowledge, this is the first-of-its-kind review on DEM specifically for biomass. Our study indicates that the current DEM models require further development, calibration, and validation based on a deep understanding of biomass particle contact mechanics and experimental data support before they can be reliably used for predictive simulations in handling and feeding systems.

42 ENGINEERING↗
NASA NTRSDate not supplied

Overview of the Predictive Simulation Capability Element of the Plume Surface Interaction Project

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the PSC Element of the PSI Project as well as providing descriptions of recent accomplishments and remaining work. The overall structure of the PSC Element is broken down into four areas of focus. The first area is the Prediction of Plume Flow in low pressure environments. The development approach taken is the augmentation of the existing production-mode computational fluid dynamics (CFD) tool Loci/Chem, with targeted extensions necessary to accurately model rarefied conditions found in both Martian and Lunar applications. Production readiness and validation of predictive capability are the major objectives of this task. The second area is the effect of mixed continuum/rarefied flow on crater development and ejecta sheets. A new CFD application, Loci/GGFS (Gas Granular Flow Solver), is being developed which implements an Eulerian/Eulerian two-phase model of gas- and soil-phases in order to simulate the soil erosion, crater formation, and soil ejecta transport in a fully coupled simulation. This task seeks to verify Loci/GGFS is production-ready as well as perform validation studies to determine the degree of predictive capability achieved by Loci/GGFS. The third area is focused on the details and extension of particle phase modeling of soil. In this task, Discrete Element Modeling (DEM) techniques are used to perform direct simulations of complex soil particles under the action of forcing similar to that to be cause by PSI. The simulation results are then used to construct closures to the Eulerian model of the soil phase used by Loci/GGFS. The fourth area is gas-particle interaction modeling. In this task, experiments are being conducted as well as detailed simulation results are being studied to further understand the complexities of gas-particle interactions in dilute, intermediate, and high soil volume fraction regimes. Improved models of particle drag and the particle turbulent kinetic energy (PTKE) resulting from the interaction of gas flows within particle clouds are the objective of this task.

Jeff West↗
DOE OSTI2022

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗
NASA NTRS1976

The effect of vertical distortion in the modeling of sedimentation phenomena - Martian crater wake streaks

Theoretical mass flow rate and the particle trajectory equations of motion of granular material in saltation are used to correlate experimental data for the rate of erosion in the wake regions of wind-tunnel-model Martian craters. Vertical geometric distortion is inherent in the simulation because of the equivalent roughness height characteristic of a turbulent boundary layer, which is affected by material in saltation. It is thus necessary to distort topographic model geometry in the vertical direction. A systematic similitude which is based on erosion rate and equivalent roughness in saltation is shown to correlate time-dependent model data as long as the model Reynolds number is higher than a critical value.

Iversen, J. D.↗
DOE OSTI2022

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗
NASA NTRS2002

Constitutive Relation in Transitional Granular Flows

To study the constitutive behavior of granular materials, the presence of gravity is detrimental. Although empirical relations have been obtained for engineering designs to control granular flows on Earth, it is not known how well these Earthbound relations can be used in another gravity field. Fundamental understanding must be derived to reliably design for granular flows in space exploration. There are two extremes of granular flows of which significant amount of knowledge is available. One deals with a dense and quasi-static situation where the deformation rate nearly vanishes. The other deals with dilute and rapidly fluctuating grain velocities where particle inertia dominates. This project, funded by the NASA Microgravity Fluid Physics Program, aims to study this transitional regime via physical experiments and computer simulations. A conceptual model has been established as described below. There are two natural time scales in a granular flow. One is the travel time between two consecutive collisions and the other is the duration of a collision contact. At a very low shear-rate, the shear-induced particle velocity is low. Hence the travel time between collisions is longer than the contact time between colliding particles. Binary collisions prevail. As the shear-rate increases, the traveling time between collisions reduces and the probability of multiple collisions goes up. These particle groups disperse shortly after and new groups form. When shear-rate is further increased, clusters grow in size due to an increasing chance for free particles to join before groups have the time to disperse. The maximum cluster size may depend on the global concentration and material properties. As the solid concentration approaches zero, the cluster size goes to one particle diameter. The maximum possible cluster size under any condition is the container size, provided that the shear flow is inside a container. The critical shear-rates that dictate the initiation of the multiple contacts, and the size and lifetime of the collision clusters, are functions of the concentration also. A 'regime' theory has been proposed by Babic et al. This theory suggested that both the solid concentration, C, and the non-dimensional shear-rate, B, are important in determining the regimes of the granular constitutive law. Additional information is included in the original extended abstract.

Shen, Hayley H.↗
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 OSTI2022

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗
DOE OSTI2020

Aluminum-doped calcium manganite particles for solar thermochemical energy storage: Reactor design, particle characterization, and heat and mass transfer modeling

A two-step cycle is discussed for solar thermochemical energy storage based on particulate aluminum-doped calcium manganite reduction/oxidation reactions for direct integration into Air Brayton cycles. The two steps encompass 1) the storage of concentrated solar irradiation within endothermic reduction of aluminum-doped calcium manganite and 2) the delivery of heat to an Air-Brayton cycle via exothermic re-oxidation of oxygen-deficient aluminum-doped calcium magnanite. A 5 kW th scale solar thermochemical reactor operating under vacuum was designed, modeled, and optimized to thermally reduce a continuous, gravity-driven flow of aluminum-doped calcium manganite particles. The granular flows were characterized in a tilt-flow rig, and particle image velocimetry was used to determine flow properties via frictional and velocity scaling relationships. Flow properties were integrated into a detailed heat and mass transfer model of the solar thermochemical reactor. A reactor design with 31° inclination angle, 230 g/min of particles, and 5.2 kW th radiative input from the high-flux solar simulator was discovered to produce an outlet flow temperature of 1158 K, with stoichiometric deviations of 0.076 and a storage efficiency of 0.628 while avoiding particle overheating and promoting longer particle residence times.

14 SOLAR ENERGY↗
NASA NTRS1989

The role of particle collisions in pneumatic transport

A model of dilute gas-solid flow in vertical risers is developed in which the particle phase is treated as a granular material, the balance equations for rapid granular flow are modified to incorporate the drag force from the gas, and boundary conditions, based on collisional exchanges of momentum and energy at the wall, are employed. In this model, it is assumed that the particle fluctuations are determined by inter-particle collisions only and that the turbulence of the gas is unaffected by the presence of the particles. The model is developed in the context of, but not limited to, steady, fully developed flow. A numerical solution of the resulting governing equations provides concentration profiles generally observed in dilute pneumatic flow, velocity profiles in good agreement with the measurements of Tsuji, et al. (1984), and an explanation for the enhancement of turbulence that they observed.

Mastorakos, E.↗
DOE OSTI2021

Corona splashing triggered by a loose monolayer of particles

In nature, high-speed raindrops often impact and spread on particulate surfaces (e.g., soil, plant leaves with spores or pollen). We study the dynamics of droplet impact on a loosely packed monolayer of particles by combining experimental and mathematical approaches. We find that the presence of mobile particles lowers the critical impact velocity at which the droplet exhibits corona splashing, as the particle area fraction is systematically increased. We rationalize this experimental observation by considering the jamming of frictional particles at the spreading rim. Elucidating the splashing transition of the drop on a particulate bed can lead to a better understanding of soil loss and erosion from falling raindrops

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2018

Gas-Granular Simulation Framework for Spacecraft Landing Plume-Surface Interaction and Debris Transport Analysis

The Gas-Granular Flow Solver (GGFS) multi-phase flow computational framework has been developed to enable simulations of particle flows complex extra-terrestrial regolith materials. Particle flows of interest include the damage of unprepared landing sites from rocket plume impingement on Moon, Mars, and asteroids. The flow solver implements an Eulerian-Eulerian two-fluid model with fluid representation of the gas phase and granular phase to avoid the need to model billions of particle interactions. The granular phase is modeled as an Eulerian fluid with constituent physics closure models derived from first-principle Discrete Element Model (DEM) particle interaction simulations that capture the complex, non-linear granular particle interaction effects. Granular phase constituent models have been developed and integrated that address the complex, non-linear granular material mechanics complexities resulting from both: the irregular, jagged particle shapes and poly-disperse mixture effects encountered in extra-terrestrial regolith, with lunar regolith as the extreme. The GGFS capabilities are being integrated into a proven NASA plume-surface interaction and debris transport simulation framework featuring the Loci/CHEM CFD program and Debris Transport Analysis (DTA) post-processing tools for applications in robotic and human Moon and Mars lander development. Integration of the three simulation tool components. Loci/CHEM, GGFS, and DTA, into a coordinated simulation framework will enable time-accurate spacecraft landing simulations that account for the alteration of the landing surface through plume-induced cratering and the resulting redirection of plume impingement flow and debris transport. Initial implementation of this simulation framework and application examples will be presented.

Liever, Peter A.↗
DOE OSTI2022

Modeling Needs for High Power Target

The next generation of high power targets will use more complex geometries, novel materials, and new concepts (like flowing granular materials); however, the current numerical approaches will not be sufficient to converge towards a reliable target design that satisfies the physical requirements. We will discuss what can be improved in the next 10 years in target modeling to support high power (MW class) targets.

43 PARTICLE ACCELERATORS↗