Search NASA⌕ Search

SEARCH · Search NASA

Results for “bulk 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Direct observations of passages of the distant neutral line (80-140 RE) following substorm onsets - ISEE-3

Specific events recorded by ISEE-3 passage through the earth magnetotail at 80-140 earth radii are discussed. The data were taken during March 20-28, 1983, when Kp ranged from +6 to -40, and included both plasma and magnetic field signatures. Magnetic field polar angle reversals from N-S flows were detected, along with alteration in the bulk plasma flow from tailward to stagnation. Tailward flow was associated with negative field values, while stagnation was mainly present with positive field values. The tail plasma at 100 radii exhibited changes shortly after a substorm event which featured particle ejection at 6.6 earth radii. The plasma sheet swept by the events could have a 25-50 radii length scale.

Baker, D. N.↗
DOE OSTI2021

Flow behavior characterization of biomass Feedstocks

With the advancement of renewable energy, the processing and handling of biomass feedstocks has drawn enormous research interest. The key to economically viable and operationally reliable biomass handling processes involves a complete understanding of the flow behavior of different feedstocks as functions of the inherent critical material properties environmental conditions, and critical processing parameters. This review introduces the basic principles and methods in bulk solids flow and highlights recent efforts to understand the flow behavior of biomass feedstocks. The effects of inherent properties, environmental parameters, and handling equipment on the flow behavior of biomass feedstocks are presented. Particle size, shape, moisture, consolidation stress, temperature, storage time at rest, chemical treatment, and the material and geometry of the handling equipment all substantially influence biomass flowability. The limitations of current studies and the potential methods to overcome them are presented. Finally, an outlook for future research effort is provided, with proposed opportunities for improving the characterization of biomass flow.

Bulk solids testing↗
NASA NTRS1979

Magnetic fields in flowing magnetotail plasmas and their significance for magnetic reconnection

A statistical study is made of the magnetic field components that occur in the presence of high-velocity flows. Results are reported for a statistical analysis of 324 intervals of intense bulk plasma flows in the plasma sheet of the geomagnetic tail, measured by the Imp 6 spacecraft. A pronounced tendency for southward fields to be associated with tailward flow more often than with earthward flow yields support for the merging concept, though not in its simplest form. Histograms of 15.36-s averages of vector magnetic field measurements during plasma flows reveal no significant changes in the distribution of beta-y between earthward or tailward plasma flow directions.

Caan, M. N.↗
NASA NTRS1991

Gasdynamic modeling of the Venus magnetotail

A gasdynamic, convected magnetic field model of the solar wind interaction with Venus is used to model the steady state Venus magnetotail. The flow obstacle surface is approximated as a tangential discontinuity. An initial obstacle shape is defined by balancing a hydrostatic equilibrium approximation for the internal plasma pressure with an external flow pressure approximation. These approximations produce a cylindrical obstacle in the distant tail. A refined obstacle shape that attempts to balance this internal pressure with the calculated external flow pressure tapers inward toward the tail axis downstream of the terminator. The bulk plasma flow and magnetic field properties compare well with experimental observations. The model predicts central magnetotail oxygen plasma number densities of about 0.2/cu cm and temperatures on the order of 10 to the 6th K flowing tailward at speeds as low as 200 m/s.

Moore, K. R.↗
NASA NTRS1977

Storage and feeding of coal

Reliable feeding of coal from storage bins to process requires the knowledge of the behavior of coal during flow. The study of the flow of bulk solids was undertaken in the 1950's and led to the development of flow ability testing equipment and of the Mass Flow concept of design for reliable flow. The theory has since been expanded to two-phase, solids-gas system, and has found world wide application in the design of storage and feeding systems.

Jenike, A. W.↗
NASA NTRS2019

The Lunar Geophysical Network Mission

Overarching Principles: Must be better than Apollo (coverage, duration, instrument performance); Learn from the Apollo experience. Lunar Geophysical Network (LGN) New Frontiers (NF)-class mission, as part of the NF-5 call. “This mission consists of several identical landers distributed across the lunar surface, each carrying geophysical instrumentation. The primary science objectives are to characterize the Moon’s internal structure, seismic activity, global heat flow budget, bulk composition, & magnetic field.” Global distribution of multiple stations. Each station should contain a seismometer, heat flow probe, electromagnetic sounder, laser retroreflector (lunar nearside). Each station must be long-lived (e.g., approximately10 years)to allow other stations (from other countries?) to be integrated with the anchor nodes to form the International Lunar Network. Why LGN? Planetary Science: Moon represents an end-member in planetary evolution (large small body, small rocky planet); Primary planetary differentiation preserved; Key to understanding terrestrial planet initial differentiation. Lunar Science: Heat flow probes yield crustal heat budget estimates; Combined with EMS (ElectroMagnetic Sounding), the temperature profile of the deep interior can be modeled along with mineralogy; Seismic and LLR (Lunar Laser Ranging) data also yield structure and compositional information of the lunar interior; High fidelity data from LGN would enhance the usefulness of the GRAIL (Gravity Recovery and Interior Laboratory) and SELENE (Selenological and Engineering Explorer) gravity data. Human Exploration: LGN must be established prior to renewed human lunar activity - we do not know the exact locations or causes of the shallow moonquakes (SMQs) - the largest magnitude seismic events recorded by Apollo (1 event per year of magnitude greater than or equal to 5); Establishing surface infrastructure near SMQ epicenters must be avoided.

Neal, C. R.↗
NASA NTRS1991

Recent developments in oscillatory thermocapillary flows

Recent research aimed at understanding the transition of thermocapillary flows from steady to oscillatory states is described. Particular attention is given to conditions for the onset of the oscillations and how they are influenced by changes in configuration such as geometry and thermal signature. It is shown that, for high Prandtl number fluids, the oscillations are caused by the coupling among the imposed thermal signature, the surface flow, and the surface deformations. The surface deformations cause a time lag between the surface and bulk-return flows and alter the driving force (surface-temperature gradient).

Ostrach, Simon↗
NASA NTRS1988

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The influence of bulk turbulence and boundary layer thickness on the secondary flow development in a square, 90 degree turning duct was investigated. A three-dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity at six cross-planes in the duct. The results from this investigation, with entrance boundary layer thickness of 20 percent, were compared with the thin boundary layer results documented in NASA CR-174811. The axial velocity profiles, cross-flow velocities, and turbulence intensities were compared and evaluated with regard to the influence of bulk turbulence intensity and boundary layer thickness, and the influence was significant. The results of this investigation expand the 90 degree curved duct experimental data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, R. A.↗
DOE OSTI2023

Investigation of Technologies to Improve Condenser Heat Transfer and Performance in a Relevant Coal-Fired Power Plant

Improvements in thermal power-generating plant performance is correlated directly to societal benefits including lower cost of reduced fuel consumption, resulting in lower cost of electricity for the consumer and reduced carbon emissions to the atmosphere. Warm steam exhausted from low-pressure steam turbines is condensed to liquid water on the exterior of thin-walled metal condenser tubes with cooling water passing through the tube interior. The condensation of steam creates a vacuum that supports turbine rotation and the concurrent generation of electricity. This vacuum is optimized when heat transfer across the wall of condenser tube is maximized. Common hindrances to heat transfer include foulants in cooling water that may form and adhere to the interior of condenser tubes, including mineral scale, microbiological films, and particulate deposition. Flowing cooling water may also include a laminar layer at the interior metal surfaces that travels more slowly than bulk water flow, serving to impede heat transfer. On the tube exterior, condensing steam forms an insulating layer of water that flows down the tube and reduces the effectiveness of cooling. Both the interior and exterior barriers to optimal heat transfer may be alleviated to some extent by surface treatments. On the tube interior, hydrophobic coatings may be applied that can reduced the adherence of foulants and of the laminar flow layer to the tube surface. On the tube exterior, hydrophobic coatings or mechanical treatments can be applied that may result in the termination of droplet growth and the departure of droplets from the surface rather than coalescence into a continuous layer of flowing water. Fourteen surface treatments were applied to condenser tubes in this study, including eight interior coatings and six exterior treatments, five of which were coatings and one a microstructural texture. Heat transfer measurement equipment simulating conditions in the condenser of an operating power plant was used to determine heat transfer coefficients by measuring sufficient flow, temperature, water chemistry and other data. Several of the tubes with interior surface treatments showed improvement in heat transfer coefficients compared with a plain (uncoated) tube, and several of the tubes with exterior surface treatments also showed enhanced heat transfer coefficients.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1984

Lean limit phenomena

The concept of flammability limits in the presence of flame interaction, and the existence of negative flame speeds are discussed. Downstream interaction between two counterflow premixed flames of different stoichiometries are experimentally studied. Various flame configurations are observed and quantified; these include the binary system of two lean or rich flames, the triplet system of a lean and a rich flame separated by a diffusion flame, and single diffusion flames with some degree of premixedness. Extinction limits are determined for methane/air and butane/air mixtures over the entire range of mixture concentrations. The results show that the extent of flame interaction depends on the separation distance between the flames which are functions of the mixtures' concentrations, the stretch rate, and the effective Lewis numbers (Le). In particular, in a positively-stretched flow field Le 1 ( 1) mixtures tend to interact strongly (weakly), while the converse holds for flames in a negatively-stretched flow. Also established was the existence of negative flames whose propagation velocity is in the same general direction as that of the bulk convective flow, being supported by diffusion alone. Their existence demonstrates the tendency of flames to resist extinction, and further emphasizes the possibility of very lean or rich mixtures to undergo combustion.

Law, C. K.↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2020

Measurement of flow properties coupled to experimental and numerical analyses of dense, granular flows for solar thermal energy storage

Granular flows of sintered bauxite proppants were researched along an inclined plane for solar thermal energy storage applications. Granular flow properties needed to drive numeric granular models were measured for improved numerical model predictions for Carbobead CP 50/140, 40/100, and 30/60 particles. Particle shape and size distributions were determined by coupling optical microscopy to an in-house image processing algorithm. The impulse excitation technique was used to measure elastic and shear moduli, and compute Poisson’s ratio. The coefficient of static sliding friction was measured using the slip-stick method, and the static rolling friction was determined from measured shear on particles positioned between two hot-pressed plates. The coefficient of restitution was measured by dropping particles on a surface and determining the kinetic energy before and after impact with the surface using high resolution particle tracking velocimetry. Particle size did not significantly impact the coefficients of restitution and static rolling friction, however, particle shape distribution resulted in a large variation in measurements. An inclined flow experiment was performed to characterize granular flows of Carbobead CP 30/60 particles using particle image velocimetry. Numerical models of the experiment using discrete element method were generated with the measured mechanical properties as inputs for comparison with experimental results. A constant directional torque rolling friction model best predicted bulk granular flow behavior. Furthermore, good agreement between the model and experiment was achieved at ambient, steady state conditions, with average velocity differences <10%.

14 SOLAR ENERGY↗
DOE OSTI2020

Laminar drag reduction in microchannels with liquid infused textured surfaces

Liquid infused surfaces (LIS) achieve non-wetting properties through asperity structures containing pockets of a lubricating liquid rather than air. Although studies have demonstrated several potential applications of LIS, the effect of liquid infused texture on the flow of liquid in microchannels has not been systematically addressed. The present work focuses on the effect of liquid-infused texture geometry and the infused liquid properties on the flow characteristics of a bulk fluid in rectangular and round microchannel geometries. Considering four different textured geometries of transverse ridges, longitudinal ridges, in-lined square posts and staggered squares, a detailed parametric study is presented to assess the resulting drag reduction and friction factor, expressed in terms of the Poiseuille number, on the flow in microchannels. The effect of geometrical parameters is investigated for different asperity solid fraction (Φ s ), asperity height (h s ) and channel thickness (H) or tube diameter (D), while the effect of the infused fluid is studied in terms of a ratio of its viscosity relative to the viscosity of the flowing bulk fluid. Based on the results of the parametric study an extension of the conventional Moody diagram for non-wetting surfaces is developed. Finally, a design of experiments study is conducted to relate the drag reduction and Poiseuille number to the various governing parameters, using which conditions that maximize drag reduction are identified.

42 ENGINEERING↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

A different view of plasma flow inside P/Halley

The Giotto spacecraft carried two different instruments - the JPA and the IMS - for the observation of hot ions in the coma of P/Halley. Although there are many similarities in the time and distance profiles of the plasma flow parameters (bulk velocity, number density, and temperature) computed from the two data sets, there are also some significant differences, especially at cometocentric distances less than 500,000 km. The principal discrepancies between the JPA results presented by Formisano et al. (1990) and the IMS observations are: (1) the IMS did not detect the levelling off of the speed and temperature profiles that Formisano et al. interpreted as flow stabilization; (2) the IMS detected differential north-south flow between the solar wind and cometary ions for only a brief interval when the magnetic field was oriented nearly southward, whereas Formisano et al. reported more extensive differential north-south flow that was independent of the direction of the field; (3) the JPA ion densities were factors of 2 to 4 higher than the IMS ion densities which, in turn, were an order of magnitude greater than theoretical values.

Neugebauer, M.↗
NASA NTRSDate not supplied

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗
DOE OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part II: Continuum-Mechanics Models

The design of efficient material-handling systems for milled lignocellulosic biomass is challenging due to their complex particle morphologies and frictional interactions. Computational modeling, including the discrete element method (DEM) and continuum-based finite-element/volume methods, may offer scientific insight and predictive capabilities for the flow of milled biomass in hoppers and feeders. Herein, this article (Part II) presents a review of current state-of-the-art continuum models for the flow of milled biomass, whereas DEM models are reviewed in a companion article (Part I). Advances of numerical methods to solve the global governing equations are discussed first, followed by a comprehensive review of constitutive models for granular materials, including Drucker–Prager, hypoplastic, Cambridge-type, inertial-rheology, and nonlocal granular fluidity models. Specifically, we provide in-depth discussion on the suitability of those models for milled lignocellulosic biomass materials in terms of nonlinear elasticity, dependence of flow strength on pressure, density and shear rate, and compaction (dilation) associated with hardening (softening). Furthermore, our study shows that, despite the recent advances in continuum granular flow modeling, the most suitable constitutive models still need further development to account for material parametrization, multiflow regimes, and multiscale behavior before they can be reliably used to optimize the design and operation of biomass handling systems.

09 BIOMASS FUELS↗