Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bulk behavior”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Let’s Stick Together: Interplay of Moisture, Particle Size, and Anatomical Fraction in the Flowability of Corn Stover Crumbles

Handling of agricultural waste biomass feedstocks (such as corn stover) is a persistent challenge in the production of biofuels and chemicals at integrated biorefineries. Inconsistent flowability and high feedstock variability creates equipment downtime and processing challenges. This work investigates the influence of anatomical fractionation and moisture content on the bulk solids handling behavior of corn stover crumbles. A combination of rotational shear testing, tribological measurements, and surface and sorption characterization techniques (inverse gas chromatography (IGC), dynamic vapor sorption (DVS), and electron microscopy) were used to evaluate the contributions of particle composition, size, and surface properties to flow behavior. Shear and frictional testing revealed that moisture content has the greatest effect on flowability, significantly increasing unconfined yield strength (ƒ c ) and reducing the flow function coefficient (F F C), regardless of anatomical fraction. While differences in surface energy were observed between fractions, particularly in polar contributions (γ AB ), these did not correspond to meaningful differences in flow behavior. Coarser 6 mm particles exhibited poorer flow performance than 4 mm crumbles, likely due to broader particle size distributions and increased particle interlocking. Taken together, these results suggest that anatomical fractionation provides minimal benefit from a bulk handling perspective. Moisture content and particle size heterogeneity dominate handling behavior, reinforcing the utility of unfractionated corn stover in processing environments. These insights can inform the design of more efficient and reliable feedstock processing systems for agricultural biomass in integrated biorefineries.

09 BIOMASS FUELS↗
DOE OSTI2024

Additive friction stir deposition induced stress ripples in aluminum alloy

In additive manufacturing, modification of the viscosity of softened or molten metal plays a significant role in the development of mechanical properties distribution. Similar to soft polymers, molten or softened metal usually exhibits vastly different viscoelastic properties. During scratching tests, viscoelastic polymers undergo spatiotemporal stress development and relaxation, causing stress ripples to form on their surfaces. Here, it is observed that friction stir-based solid-state manufacturing can yield similar stress ripples on softened aluminum alloys due to the linear motion of the deposition tool. Softened aluminum alloys exhibit viscoelastic and viscoplastic responses similar to polymeric material during the solid-state deposition process. In the final product, residual stress ripples appear on the surface of the metal alloy due to periodic residual stress relaxation, leading to oscillatory behavior from the bulk mechanical properties of the material. The orientation of the in-built stress pattern formed in the bulk of the additively manufactured metal alloy underneath the visually evident onion ring-like surface pattern formation is diametrically in the opposite direction. The experimentally observed results are consistent with analytical calculation and numerical simulations in terms of behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanoscale Superconducting States in the Fe-Based Filamentary Superconductor of Pr-Doped CaFe 2 As 2

The low-temperature scanning tunneling microscope and spectroscopy (STM/STS) are used to visualize superconducting states in the cleaved single crystal of 9% praseodymium-doped CaFe 2 As 2 (Pr-Ca122) with Tc ≈ 30 K. The spectroscopy shows strong spatial variations in the density of states (DOS), and the superconducting map constructed from spectroscopy discloses a localized superconducting phase, as small as a single unit cell. The comparison of the spectra taken at 4.2 K and 22 K (below vs. close to the bulk superconducting transition temperature) from the exact same area confirms the superconducting behavior. Nanoscale superconducting states have been found near Pr dopants, which can be identified using dI/dV conductance maps at +300 mV. There is no correlation of the local superconductivity to the surface reconstruction domain and surface defects, which reflects its intrinsic bulk behavior. We, therefore, suggest that the local strain of Pr dopants is competing with defects induced local magnetic moments; this competition is responsible for the local superconducting states observed in this Fe-based filamentary superconductor.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

From bulk to interface: electrochemical phenomena and mechanism studies in batteries via electrochemical quartz crystal microbalance

Understanding the bulk and interfacial behaviors during the operation of batteries (e.g., Li-ion, Na-ion, Li–O 2 batteries, etc.) is of great significance for the continuing improvement of the performance. Electrochemical quartz crystal microbalance (EQCM) is a powerful tool to this end, as it enables in situ investigation into various phenomena, including ion insertion/deinsertion within electrodes, solid nucleation from the electrolyte, interphasial formation/evolution and solid–liquid coordination. As such, EQCM analysis helps to decipher the underlying mechanisms both in the bulk and at the interface. This tutorial review will present the recent progress in mechanistic studies of batteries achieved by the EQCM technology. Here, the fundamentals and unique capability of EQCM are first discussed and compared with other techniques, and then the combination of EQCM with other in situ techniques is also covered. In addition, the recent studies utilizing EQCM technologies in revealing phenomena and mechanisms of various batteries are reviewed. Perspectives regarding the future application of EQCM in battery studies are given at the end.

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

Understanding the limitations and potential of micro tensile testing of tungsten and needs for crystal plasticity modelling

Tungsten is the leading plasma facing material candidate due to its exceptional properties. Understanding the response to neutron irradiation is crucial for the lifetime evaluation of tungsten. The limited space in nuclear reactors and the high levels of radioactivity of the specimens after irradiation are significant barriers to accurately measuring the mechanical properties after neutron irradiation. Testing of micro tensile specimens is one approach to reduce the total amount of irradiated material needed for a set of mechanical testing experiments. However, micro samples are not necessarily measuring the bulk material properties, as size effects produce higher measured mechanical properties than is observed in engineering size counterparts. Determining the minimal sample size for reliable bulk property measurement, combined with modeling efforts, is essential. We conducted room temperature tensile tests on tungsten micro tensile specimens fabricated with focused ion beam, plasma focused ion beam and femto-second laser ablation system to dimension of 2×2×7?µm³, 7×7×18?µm³, and 80×100×233?µm³ (width×thickness×gauge-length), respectively. Elevated tensile testing was performed up to 475°C on 5×5×18?µm³ specimens fabricated with focused ion beam. The smallest specimens exhibited a high degree of ductility and strength, whereas the largest specimens demonstrated behavior akin to bulk tungsten. To investigate the effects introduced by the micro specimen fabrication processes and to obtain the necessary bulk dislocation density for the crystal plasticity model, we employed X-ray diffraction depth profiling. This measurement was performed using different X-ray sources with varying penetration depths. Finally, the potential and limitations of micro mechanical tests for tungsten will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Development of Mechanistic Fission Gas Release and Swelling Models for UN Fuels in BISON

This report describes the work in NEAMS (Nuclear Energy Advanced Modeling and Simulation) to develop a mechanistic fission gas model for uranium nitride fuels in BISON. The existing Sifgrs (Simple integrated fission gas release and swelling) model tracks the average properties of two bubble populations in the bulk and at the grain boundaries. It was recognized that dislocations play a crucial role in the fission gas swelling of UN, and an irradiation-induced dislocation density model was needed, as well as a model describing how fission gas interacts with dislocations creating a third population of bubbles along dislocations. A mechanistic model that tracks dislocation bubbles was implemented in Sifgrs. This model was used to simulate fission gas swelling and release in UN. Lower-length-scale calculations and experimental observations from carbide fuel were leveraged to help populate the model with essential parameters. As a placeholder, an empirical function was formulated for the evolution of the dislocation network. Because of the difference in evolution of this network at different temperatures, the dislocation bubbles are able to capture behavior that the bulk intragranular bubbles cannot. It was found that the bulk bubbles dominate microscopic swelling at low temperatures, and the dislocation bubbles dominate at higher temperatures. Based on this model, the transition between the two bubble types is the main factor behind the breakaway swelling phenomenon in UN. In order to test the model, a couple of assessments were run. For the lower temperature JOYO pins, the model produces reasonable fission gas release and swelling values. For the higher temperature SP1 pin, the model dramatically underestimates the fission gas release. This issue is attributed to the gas being trapped inside dislocation bubbles, unable to escape to grain boundaries to cause fission gas release, and could be remedied by a mechanistic dislocation model that allows the dislocation density to decrease at very high temperatures. A preliminary mechanistic dislocation model was developed supported by first-principle calculations. These calculations provided valuable insight into how interstitial defects cluster in UN in the {110} orientation, and may eventually form dislocation loops, leading to the conclusion that dislocation loops may nucleate from these clusters. An estimate of the dislocation line energy in UN was obtained and will be improved in future work. The free-energy cluster dynamics code Centipede was used to track interstitial and vacancy absorption at dislocation loops and calculate their growth. In addition, improvements to Centipede were made including new convergence criteria for transient simulations and sink driving force updates.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Signature of magnetoelectric coupling driven finite momentum pairing in 3D ising superconductor

The finite momentum superconducting paring states (FMPs) represent a forefront of condensed matter physics. Here we report experimental evidence of FMP in a locally noncentrosymmetric bulk superconductor 4H b -TaS 2 . Using hard X-ray diffraction and angle-resolved photoemission spectroscopy, we reveal unusual 2D ferro-rotational charge density wave (CDW) and weak interlayer hopping in 4H b -TaS 2 . The superconducting upper critical field, H c2 , linearly increases via decreasing temperature, and well exceeds the Pauli limit, suggesting the dominant orbital pair-breaking mechanism. Remarkably, we observed evidence of field-induced superconductivity-to-superconductivity transition that breaks continuous rotational symmetry of the s-wave uniform pairing in the Bardeen-Cooper-Schrieffer theory down to the six-fold rotation symmetry. Ginzburg-Landau free energy analysis shows that magnetoelectric coupling, induced by 2D ferro-rotational CDW, stabilizes FMP that provides an explanation of the lowering rotation symmetry. Our results provide a new understanding of unconventional superconducting behaviors of the bulk quantum heterostructure 4H b -TaS 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗
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 OSTI2023

Fraction of Free-Base Nicotine in Simulated Vaping Aerosol Particles Determined by X-ray Spectroscopies

A new generation of electronic cigarettes is exacerbating the youth vaping epidemic by incorporating additives that increase the acidity of generated aerosols, which facilitate uptake of high nicotine levels. We need to better understand the chemical speciation of vaping aerosols to assess the impact of acidification. Here we used X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy to probe the acid-base equilibria of nicotine in hydrated vaping aerosols. Here we show that, unlike the behavior observed in bulk water, nicotine in the core of aqueous particles was partially protonated when the pH of the nebulized solution was 10.4, with a fraction of free-base nicotine (α FB ) of 0.34. Nicotine was further protonated by acidification with equimolar addition of benzoic acid (α FB = 0.17 at pH 6.2). By contrast, the degree of nicotine protonation at the particle surface was significantly lower, with 0.72 < α FB < 0.80 in the same pH range. The presence of propylene glycol and glycerol completely eliminated protonation of nicotine at the surface (α FB = 1) while not affecting significantly its acid-base equilibrium in the particle core. These results provide a better understanding of the role of acidifying additives in vaping aerosols, supporting public health policy interventions.

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

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2021

Comparison of Biotite Elastic Properties Recovered by Spherical Nanoindentations and Atomistic Simulations — Influence of Nano-Scale Defects in Phyllosilicates

Phyllosilicate minerals, due to their sheets structure and morphology, are known to cause anisotropy in bulk rock properties and make the bulk rock more compliant. Accurately characterizing the micromechanical behavior of phyllosilicate minerals from laboratory observations, which eventually translates to the bulk rock behavior, is still challenging due to their fine-grained nature. Recent advances in atomistic simulations open the possibility of theoretically investigating such mineral mechanical behavior. We compare the elastic properties of biotites recovered by spherical nanoindentation with those predicted from density functional theory (DFT) simulations to investigate to what extent theoretical predictions reproduce actual phyllosilicate properties. Spherical nanoindentation was conducted using schist rocks from Poorman Formation, South Dakota, USA, to recover continuous indentation stress-strain curves. Loading in the layer-normal orientation shows an average indentation modulus (M) of about 35 GPa, while loading in the layer-parallel orientation gives a higher average of about 95 GPa. To facilitate comparison, the elastic stiffness constants (cij) determined from DFT were converted to indentation modulus (M) using solutions proposed in this study. The majority of the nanoindentation modulus results are below the values inferred from the simulation results representing ideal defect-free minerals. We suggest that crystal defects present at the nano-scale, potentially ripplocations, are the dominant cause of the lower indentation modulus recovered from nanoindentation compared to those inferred from DFT simulations. Results highlight the importance of acknowledging the defects that exist down to the nano-scale as it modifies the mechanical properties of phyllosilicates compared to its pure defect-free form.

58 GEOSCIENCES↗
DOE OSTI2022

Suppressing electrolyte-lithium metal reactivity via Li + -desolvation in uniform nano-porous separator

Lithium reactivity with electrolytes leads to their continuous consumption and dendrite growth, which constitute major obstacles to harnessing the tremendous energy of lithium-metal anode in a reversible manner. Considerable attention has been focused on inhibiting dendrite via interface and electrolyte engineering, while admitting electrolyte-lithium metal reactivity as a thermodynamic inevitability. Here, we report the effective suppression of such reactivity through a nano-porous separator. Calculation assisted by diversified characterizations reveals that the separator partially desolvates Li + in confinement created by its uniform nanopores, and deactivates solvents for electrochemical reduction before Li 0 -deposition occurs. The consequence of such deactivation is realizing dendrite-free lithium-metal electrode, which even retaining its metallic lustre after long-term cycling in both Li-symmetric cell and high-voltage Li-metal battery with LiNi 0.6 Mn 0.2 Co 0.2 O 2 as cathode. The discovery that a nano-structured separator alters both bulk and interfacial behaviors of electrolytes points us toward a new direction to harness lithium-metal as the most promising anode.

25 ENERGY STORAGE↗
DOE OSTI2021

Excess-iron driven spin glass phase in Fe 1 + y Te 1 – x Se x

The iron-chalcogenide superconductor FeTe 1– x Se x displays a variety of exotic features distinct from iron pnictides. Although much effort has been devoted to understanding the interplay between magnetism and superconductivity near x = 0.5, the existence of a spin glass phase with short-range magnetic order in the doping range ( x ~ 0.1–0.3) has rarely been studied. Here, we use DC/AC magnetization and (quasi) elastic neutron scattering to confirm the spin-glass nature of the short-range magnetic order in a Fe 1.07 Te 0.8 Se 0.2 sample. The AC-frequency dependent spin-freezing temperature T f generates a frequency sensitivity Δ T f ( ω )/[ T f ( ω )Δlog 10 ω ] ≈ 0.028 and the description of the critical slowing down with τ = τ 0 ( T f / T SG – 1) - z v gives T SG ≈ 22 K and zv ≈ 10, comparable to that of a classical spin-glass system. We have also extended the frequency-dependent T f to the smaller time scale using energy-resolution-dependent neutron diffraction measurements, in which the T N of the short-range magnetic order increases systematically with increasing energy resolution. We find that by removing the excess iron through annealing in oxygen, the spin-freezing behavior disappears, and bulk superconductivity is realized. Thus, the excess Fe is the driving force for the formation of the spin-glass phase detrimental to bulk superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Photoemission and dynamical mean field theory study of electronic correlations in a t 2 g 5 metal SrRhO 3 thin film

Perovskite rhodates are characterized by intermediate strengths of both electronic correlation as well as spin-orbit coupling (SOC) and usually behave as moderately correlated metals. A recent publication on epitaxial SrRhO 3 thin films reported a bad-metallic behavior and suggested the occurrence of antiferromagnetism below 100 K. In this work, we have further studied this SrRhO 3 thin film by hard x-ray photoemission spectroscopy and found a very small density of states (DOS) at the Fermi level $(E_F)$, which is consistent with the previously reported bad-metallic behavior. However, the absence of DOS persists up to room temperature, which contradicts the explanation of antiferromagnetic transition at ~100 K. We also employed electronic structure calculations within the framework of density functional theory and dynamical mean-field theory. In contrast to the photoemission results, our calculations indicate metallic behaviors of both bulk SrRhO 3 and the SrRhO 3 thin film, and a stronger correlation effect was observed in the thin film than that in the bulk. The calculated uniform magnetic susceptibility is substantially larger in the thin film than that in the bulk. We also investigated the role of SOC and found only a moderate modulation of the band structure. Hence SOC is not expected to significantly affect the electronic correlation in SrRhO 3 . Extrinsic effects of finite-thickness effects beyond our calculations and localization effects may play important roles to induce the negligible spectral weight at $(E_F)$.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Orientation-dependent stabilization of MgCr 2 O 4 spinel thin films

AB 2 O 4 normal spinels with a magnetic B site can host a variety of magnetic and orbital frustrations leading to spin-liquid phases and field-induced phase transitions. Here, we report the epitaxial growth of (111)-oriented MgCr 2 O 4 thin films. By characterizing the structural and electronic properties of films grown along the (001) and (111) directions, the influence of growth orientation has been studied. Despite distinctly different growth modes observed during deposition, the comprehensive characterization reveals no measurable disorder in the cation distribution nor multivalency issue for Cr ions in either orientation. Contrary to a naive expectation, the (111) stabilized films exhibit a smoother surface and a higher degree of crystallinity than (001)-oriented films. The preference in growth orientation is explained within the framework of heteroepitaxial stabilization in connection to a significantly lower (111) surface energy. Furthermore, these findings open broad opportunities in the fabrication of two-dimensional kagome-triangular heterostructures with emergent magnetic behavior inaccessible in bulk crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

3D polycatenated architected materials

Architected materials derive their properties from the geometric arrangement of their internal structural elements. Their designs rely on continuous networks of members to control the global mechanical behavior of the bulk. Here, in this study, we introduce a class of materials that consist of discrete concatenated rings or cage particles interlocked in three-dimensional networks, forming polycatenated architected materials (PAMs). We propose a general design framework that translates arbitrary crystalline networks into particle concatenations and geometries. In response to small external loads, PAMs behave like non-Newtonian fluids, showing both shear-thinning and shear-thickening responses, which can be controlled by their catenation topologies. At larger strains, PAMs behave like lattices and foams, with a nonlinear stress-strain relation. At microscale, we demonstrate that PAMs can change their shapes in response to applied electrostatic charges. The distinctive properties of PAMs pave the path for developing stimuli-responsive materials, energy-absorbing systems, and morphing architectures.

36 MATERIALS SCIENCE↗
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↗