Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “granular solids”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Application of Interparticle Spacing Model to Maximize Filler Content in Cementitious Pastes

One of the most impactful ways to reduce embodied carbon of concrete in the near term is by partially replacing Portland cement with supplementary cementitious materials and/or fillers. This article describes an approach to reduce cement content in concrete through the development of high filler, low water (HFLW) cementitious pastes. Pastes with cement content as low as 50% (% weight of fines) were designed by applying models that maximize the packing density of the granular system and that consider the characteristics of the materials to calculate interparticle spacing (IPS), thereby allowing to obtain pastes with higher solids concentration and lower water demand. Strong correlations were found between IPS, the consistency index of the Herschel-Bulkley rheological model, and cement hydration kinetic parameters of binary and ternary pastes containing a Portland limestone cement and up to two different types of ground limestone. The approach showed to be feasible for the design of HFLW concrete paste fraction.

Antunes da Silva, Denise↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of processing parameters on the morphologies of complex sesquioxide thin films

Controlling and predicting the morphology of lanthanide sesquioxides in thin film form is vital to their use in current applications. In the present study, single and codeposited Sm2O3, Er2O3, and Lu2O3 thin films were grown on yttria-stabilized zirconia (8%) substrates by radio frequency magnetron sputtering at room temperature and 500 °C. The effect of two different substrate temperatures and altering the oxide cation on the structural and morphological properties of the films was analyzed. The thin films were characterized by profilometry, scanning electron microscopy, transmission electron microscopy, and x-ray diffraction. The single-component Lu2O3 and Sm2O3 films obtained were of the cubic phase, and the Er2O3 was a mix of cubic and monoclinic phases. It was observed for both the Er2O3 and Lu2O3 films that increasing the substrate temperature to 500 °C resulted in larger grained polycrystalline films. In contrast, large grained polycrystalline films were obtained at both room temperature and 500 °C for Sm2O3 and uneven granularity increased as temperature increased. Codeposition of Lu2O3 and Sm2O3, and Lu2O3 and Er2O3 resulted in a cubic bixbyite phase (the C phase of the lanthanide sesquioxide) solid solution. It was observed that the structure and morphology of the films can be controlled by manipulating deposition parameters. Both substrate temperature and altering the oxide cation contributed to changes in crystallinity and grain structure, which can modify the chemical and physical properties of the films for their applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring pressure-dependent inelastic deformation and failure in bonded granular composites: An energetic materials perspective

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material’s mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. Further, to construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer’s glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A hybrid machine-learning approach for analysis of methane hydrate formation dynamics in porous media with synchrotron CT imaging

Fast multi-phase processes in methane hydrate bearing samples pose a challenge for quantitative micro-computed tomography study and experiment steering due to complex tomographic data analysis involving time-consuming segmentation procedures. This is because of the sample's multi-scale structure, which changes over time, low contrast between solid and fluid materials, and the large amount of data acquired during dynamic processes. Here, a hybrid approach is proposed for the automatic segmentation of tomographic data from time-resolved imaging of methane gas-hydrate formation in sandy granular media, which includes a deep-learning 3D U-Net model. To prepare a training dataset for the 3D U-Net, a technique to automate data labeling based on sample-specific information about the mineral matrix immobility and occasional fluid movement in pores is proposed. Automatic segmentation allowed for studying properties of the hydrate growth in pores, as well as dynamic processes such as incremental flow and redistribution of pore brine. Results of the quantitative analysis showed that for typical gas-hydrate stability parameters (100 bar methane pressure, 7°C temperature) the rate of formation is slow (less than 1% per hour), after which the surface area of contact between brine and gas increases, resulting in faster formation (2.5% per hour). Hydrate growth reaches the saturation point after 11 h of the experiment. Finally, the efficacy of the proposed segmentation scheme in on-the-fly automatic data analysis and experiment steering with zooming to regions of interest is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Understanding the trade-offs of national municipal solid waste estimation methods for circular economy policy

Policies embracing circular economy concepts have taken hold in national legislation around the world. As the number of governments and organizations adopting circular economy policies increases, so does the need for accurate and timely measurement of material resource flows. Since many countries do not have access to centrally reported municipal solid waste (MSW) data, estimation and modeling are critical in evaluating circular economy policy effectiveness. The purpose of this paper is to examine three modeling approaches estimating national MSW data in the United States, including industry-based material flow analysis, waste-extended input-output modeling, and aggregated regional waste reporting. We establish five criteria to guide the analysis through the context of policy monitoring (data quality, flow totality, update frequency, sensitivity to disruption, and product granularity) and use these criteria to analyze and score each model. We then use a literature search to identify five, internationally-implemented options for circular economy policy and determine the data and modeling components that are most helpful in evaluating policy effectiveness. Finally, we provide a crosswalk of the model scores and policy needs to inform the suitability of model selection by policy type. We found that data quality and update frequency are identified as critical components for evaluating circular economy policies within the models evaluated, and can both be fulfilled by aggregated regional waste reporting. Flow totality, sensitivity to disruption, and product granularity requirements vary by both model and policy types. While none of the evaluated models satisfy the combination of requirements for any of the five policies, industry-based material flow analysis offers flow totality for extended producer responsibility, landfill bans, and recycling rate target policies that typically require it. Here, the waste-extended input-output model can provide disruption sensitivity and product granularity as needed for policies like minimum recycled content and market restrictions. Policy developers in areas where strong centralized data collection is not an option should design policy action(s) with modeling tradeoffs in mind, including the potential hybridization of modeling approaches that may provide the most accurate national MSW estimates.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
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 OSTI2022

Effect of Local Structural Distortions on Antiferroelectric–Ferroelectric Phase Transition in Dilute Solid Solutions of K x Na 1– x NbO 3

The fundamental principles that govern antiferroelectric (AFE)–ferroelectric (FE) transitions are not well understood for many solid solutions of perovskite compounds. For example, crystal chemical considerations based on the average Goldschmidt tolerance factor or ionic polarizability do not precisely predict the boundary between the AFE and FE phases in dilute solid solutions of alkali niobates, such as K x Na 1– x NbO 3 ( x ≤ 0.02). Here, based on detailed structural analysis from neutron total scattering experiments, we provide insights about how the relative local distortions around the A- and B-sites of the ABO 3 perovskite structure affect the AFE/FE order of the average crystallographic phases in K x Na 1– x NbO 3 . Specifically, we show that a higher (lower) ratio of B-site-centered distortions over A-site-centered distortions drives transition toward a long-range FE (AFE) phase, which is based on a competition between the long-range polarizing field of the Nb–O dipoles and the disordering effect of local distortions around the A-site. Our study provides a predictive tool for designing complex solid-solution perovskites with tunable (anti)ferroelectric polarization properties, which can be of interest for various energy-related applications such as high-density energy storage and solid-state cooling.

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

Tuning Magnetic Symmetry and Properties in the Olivine Series Li 1–x Fe x Mn 1–x PO 4 through Selective Delithiation

In this work, we present the first in-depth study of the nuclear and magnetic structures of the fully lithiated LiFe x Mn 1–x PO 4 solid solutions as well as the partially delithiated Li 1– xFe x Mn 1–x PO 4 series using neutron powder diffraction. By using a weak oxidizing solution that only oxidizes Fe 2+ , we find that lithium can be removed from these olivines in a stepwise manner. Representational analysis was performed on the observed magnetic reflections and found that the lithiated solid solutions order in the magnetic space group Pnma', which is identical to triphylite (LiFePO 4 ). We report that, even when removing 20% of the lithium from the hcp lattice, the magnetic structure changes and multiple magnetic contributions are needed to fully encompass the observed magnetic order. We find that the partially delithiated x = 0.2 and 0.3 have a magnetic structure P2 1 '/n11, and when x surpasses 0.3, the magnetic structure changes to P2 1 2 1 2 1 . We give a detailed report on the effect of delithiation on magnetic order, bond distances/angles, and the super- and supersuperexchange interactions in the material. The symmetries of the magnetic order of these solid solutions are then discussed with regard to their potential as ferrotoroidic candidates.

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

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗
DOE OSTI2022

Effect of grain dissolution on sloping ground

The static and dynamic stability of natural or constructed slopes can be affected by dissolution or dissolution-like phenomena. Their underlying mechanisms, however, remain unclear. New experimental results and discrete element simulations provide particle-level and macroscale information on the consequences of mineral dissolution on slope behavior. At the microscale, load-carrying grain arches develop around dissolving particles, the porosity increases, and contact force chains evolve to form a honeycomb topology. At the macroscale, while vertical settlements are the prevailing deformation pattern, lateral granular movements that create mass wasting are prominent in sloping ground, even under the quasi-static granular loss. Horizontal grain displacement is maximum at the surface and decreases linearly with the distance from the slope surface to become zero at the bottom boundaries, much like vertical granular displacement along the depth. Sediments with smaller friction angles and steeper slopes experience greater displacement, both vertically and horizontally. Slopes become flatter after dissolution, with the reduction in slope angle directly related to the loss in ground elevation, ΔH/H o . Yet, because of the porous fabric that results from dissolution, vertical shortening is less than the upper bound, estimated from the loss in the solid mass fraction, ΔH/H o ≈SF. Under water-saturated conditions, the post-dissolution fabric may lead to sudden undrained shear and slope slide.

42 ENGINEERING↗
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 OSTI2022

Enhancement of Two Heterogeneous Reactions in the Production of the Medical Isotope 188 W

We report the efficiency of heterogeneous reactions is influenced by kinetic conditions with an emphasis on reactants in different phases to have effective interactions. Optimal interactions of reactants rely on the rational design of the equipment and engineering parameters for the reaction process, especially for the process on a large scale. 188 W production at the Oak Ridge National Laboratory is such an example. In processing the irradiated 186 W metal target for the 188 W product, two of the heterogeneous reactions are essential ones: (1) the gas–solid reaction of O 2 + W to form WO 3 and (2) the gas–liquid reaction of OsO 4 being absorbed by the NaOH sorbent. Better understanding of the two reactions has led to revisions to the design of the quartz reaction vessel (gas–solid reaction) and the scrubbing system (gas–liquid reaction). The optimized quartz vessel resulted in less heating hours, whereas a more efficient OsO 4 scrubbing system ensured that all volatile OsO 4 was qualitatively captured.

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

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2022

Thermal Modeling and Simulation of the Packed-bed Thermal Energy Storage Combined with INL Thermal Energy Distribution System

Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.

25 ENERGY STORAGE↗