Search NASA⌕ Search

SEARCH · Search NASA

Results for “granular media”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

Using limited neural networks to assess relative mechanistic influence on shock heating in granular solids

The rapid compaction of granular media results in localized heating that can induce chemical reactions, phase transformations, and melting. However, there are numerous mechanisms in play that can be dependent on a variety of microstructural features. Machine learning techniques such as neural networks offer a ubiquitous method to develop models for physical processes. Limiting what kind of microstructural information is used as an input and assessing normalized changes in network error, the relative importance of different mechanisms can be inferred. Here we utilize binned, initial density information as network inputs to predict local shock heating in a granular high explosive trained from large-scale molecular dynamics simulations. In this study, the spatial extent of the density field used in the network is altered to assess the importance and relevant length scales of the physical mechanisms in play, where different microstructural features result in different predictive capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Closed-form existence conditions for bandgap resonances in a finite periodic chain under general boundary conditions

Bragg scattering in periodic media generates bandgaps, frequency bands where waves attenuate rather than propagate. Yet, a finite periodic structure may exhibit resonance frequencies within these bandgaps. This is caused by boundary effects introduced by the truncation of the nominal infinite medium. Previous studies of discrete systems determined existence conditions for bandgap resonances, although the focus has been limited to mainly periodic chains with free–free boundaries. In this paper, we present closed-form existence conditions for bandgap resonances in discrete diatomic chains with general boundary conditions (free–free, free–fixed, fixed–free, or fixed–fixed), odd or even chain parity (contrasting or identical masses at the ends), and the possibility of attaching a unique component (mass and/or spring) at one or both ends. The derived conditions are consistent with those theoretically presented or experimentally observed in prior studies of structures that can be modeled as linear discrete diatomic chains with free–free boundary conditions. An intriguing case is a free–free chain with even parity and an arbitrary additional mass at one end of the chain. Introducing such an arbitrary mass underscores a transition among a set of distinct existence conditions, depending on the type of chain boundaries and parity. The proposed analysis is applicable to linear periodic chains in the form of lumped-parameter models, examined across the frequency spectrum, as well as continuous granular media models, or similar configurations, examined in the low-frequency regime.

Bastawrous, Mary V.↗
DOE OSTI2020

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗
DOE OSTI2025

The role of fast and slow dynamics in nonlinear resonant ultrasound spectroscopy of consolidated granular materials

Abstract Elastic nonlinearity observed in consolidated granular media can be attributed to the combination of slow and fast effects, which give rise to hysteresis and relaxation of both modulus and damping after the sample is perturbed. A consequence is a high level of complexity in the measurements of the sample linear and nonlinear elastic parameters. The results of experiments are dependent on the experimental protocol that is adopted to measure the relevant quantities and it is hard to quantify parameters with accuracy and repeatability. Here we focus on examining Nonlinear Resonant Ultrasound Spectroscopy, showing experimentally the role of slow dynamics in the process and quantifying/discussing its influence on the quantification of nonlinearity. We also propose a model to describe the process, which shows that different contributions to nonlinearity (e.g., classical and hysteretic) could be due to physical features (defects) relaxing with different relaxation times.

Science & Technology - Other Topics↗
DOE OSTI2020

A massively parallel and scalable multi-CPU material point method

Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid (G2P2G) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array (AoSoA) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 10243 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Wang, Xinlei↗
DOE OSTI2020

Adhesion-force micro-scale study of desiccating granular material

Experiments on five-, four-, three- and two-wet-hydrophilic-grain clusters were performed to investigate evolution of adhesion of granular media during drying on the micro-scale. The experiments show that the adhesion-force of a cluster initially grows at most to three times the original value before decreasing to zero by the end of evaporation. The adhesion-force is composed of capillary pressure force acting over the liquid/solid contact surface area, and surface tension forces acting over the three-phase contact perimeter length. This is in contrast with most macro-scale phenomenological models, in which the only desaturation process variables affecting strength are suction and saturation. Both the contact surface area and contact perimeter length are reduced to zero upon complete liquid evaporation. The morphology of an evaporating water body evolves through slow flow controlled by evaporation rate, interrupted by various modes of fast air entry, which are non-equilibrium jumps of liquid/gas interfaces (Haines jumps). The instabilities involve large adhesion force discontinuities and substantial water mass reconfiguration with water flow in an extremely short time, which makes the process transient. The reconfigurations can reduce the original multi-grain water clusters to four-, three- and two-grain clusters by way of three different instability modes: of thin-sheet instability, or meniscus snap-through instability, depending on the sign of the Gauss curvature of the liquid surface, or finally, for two-grain bridges only, a liquid wire pinch-off. For larger meso-scale assemblies, however, the global adhesion-force evolution is little affected by the jumps. The air entries are potential sites for drying cracks. The (approximately) calculated capillary pressure for two- and three-grain clusters, in no cases is seen to reach high values, predicted from water retention curves.

Engineering↗
DOE OSTI2025

Compact Counter Flow Fluidized Bed Heat Exchanger (Final Report)

Generation 3 concentrating solar power systems which utilize solid granular media as the heat transfer and storage media hold promise in reducing the levelized cost of electricity relative to generation 2 systems which utilize molten salt. The particles can be heated to temperatures exceeding 800 °C which enables the use of high efficiency supercritical CO 2 Brayton power cycles. A key challenge associate with the use of solid media is its low bulk thermal conductivity which limits the ability to transfer heat from the particles to a working fluid, such as supercritical CO 2 . Tradition shell and tube and shell and plate heat exchangers feature gravity driven particle flow which limits the maximum particle side heat transfer coefficient in these systems to approximately 400 W/m 2 K. A fluidized bed heat exchanger does not rely solely on gravity for particle flow. Instead, enhanced mixing and horizontal conveyance is achieved through fluidization induced through air injection at the bottom of a bed of particles. A tube bundle within the bed conveys the supercritical CO 2 which is heated via the turbulent particle motion at the tube wall.

14 SOLAR ENERGY↗
DOE OSTI2022

Thermal protection barrier

A thermal protection barrier including a base layer having a high melting temperature granular media having free standing granules disposed within the base layer and a distributor grid layer having a plurality of spaced-apart distributors composed of a high melting temperature material disposed adjacent to the base layer. The spaced-apart distributors are dimensioned and configured to provide load bearing support to the thermal protection barrier prior to a breach of the thermal protection barrier and to divide and disperse a superheated flowing mass to the free standing granules of the base layer upon occurrence of a breach. In certain embodiments, the thermal protection barrier is positioned adjacent to a nuclear system such that, upon occurrence of a loss of containment accident, the distributor grid layer is positioned and configured to divide and disperse a superheated flowing nuclear mass to the free standing granules.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

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

Electrochemical reduction of per- and polyfluorinated alkyl substances (PFAS): is it possible? Applying experimental and quantum mechanical insights from the reductive defluorination literature

Remediation of per- and polyfluorinated alkyl substances (PFAS) in global water systems is a critical human and environmental health challenge facing society. PFAS consumption is associated with a litany of adverse health effects, and our knowledge of these dangers is still evolving. Current techniques to remove PFAS from water include adsorption to media (e.g. granular activated carbon, ion-exchange resin), nanofiltration, and reverse osmosis. However, these processes create a concentrated PFAS residual that requires further management. Destructive techniques are therefore needed to detoxify these residuals. Oxidative techniques have garnered the most attention (e.g. supercritical water oxidation, electrochemical oxidation) but are energy intensive and potentially form toxic by-products. As an alternative, several groups have researched advanced reduction processes that form aqueous electrons, but these processes are still chemical and energy intensive (e.g. ultraviolet/SO 3 2− , electron beam). This concise review therefore focuses on whether electrochemical reduction — a chemical-free, modular process — could be technically feasible for PFAS destruction.

King, Jacob F.↗
DOE OSTI2021

Hybrid-energy approach enabled by heat storage and oxy-combustion to generate electricity with near-zero or negative CO 2 emissions

We assess a hybrid-energy approach that modifies a steam-turbine power plant to use renewable energy sources (electricity, plus options for geothermal and solar heat), plus fossil fuel (natural gas and coal) and/or waste biomass (e.g., Douglas fir woodchips). Heat storage allows heat to be created during periods of excess energy supply and for that heat to be converted to electricity when demanded. Excess electricity, such as from variable renewable energy (VRE), is used to generate oxygen for oxy-combustion furnaces that create very-hot, high-purity CO 2 that heats granular rock beds in insulated vessels. Cool CO 2 leaving the beds is dried, sent to compressors powered by excess VRE electricity, before being sent by pipeline to geologic CO 2 storage. Very-hot CO 2 transfers high-grade heat from storage to the power plant in a closed loop that returns medium-grade heat back to storage, allowing low- and medium-grade renewable-heat sources to be stacked beneath combustion heat, with all heat sources being converted to electricity at the same (high) thermal efficiency. Reliable, on-demand power may be generated with near-zero CO 2 emissions with fossil fuel and with negative CO 2 emissions with waste biomass. Our analyses show that with fossil fuel, up to 35% of gross power can be derived from renewable sources, while for waste biomass, it can be entirely derived from renewable sources. Because our approach has the potential to ensure that grids have a continuous supply of clean energy and because electricity is only generated once, when demanded, it could serve as an efficient alternative to bulk energy storage.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗
DOE OSTI2023

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated 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 concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗
DOE OSTI2024

In‐operando thermal transport characterization of moving particle bed heat exchanger

The project aimed to develop modulated photothermal radiometry (MPR) for conducting in-operando measurements of thermal transport in moving particle bed heat exchangers (MPBE) located at Sandia National Laboratory (SNL). The plan was to initially conduct lab-scale measurements on a 1-meter-long particle channel at UCSD, followed by in- operando measurements at SNL on a 20 kWth MPBE prototype built for this project, and if time permitting, potential final deployment on SNL’s 1 MWth G3P3 MPBE. The first two goals were successfully achieved but due to changes with commissioning timeline of the G3P3 MPBE, in-operando measurements on the G3P3 MPBE could not be conducted. The measurement campaigns have resulted in some promising results and the MPR tool has been handed over to trained SNL personnel for future use. MPR measurements on UCSD’s 1-meter-long particle channel showed a peculiar length dependence of the particle bed thermal conductivity and near-wall thermal resistance. The bed thermal conductivity increased, and the near-wall thermal resistance decreased, with increasing distance from the channel inlet, both eventually reaching a plateau. This was found to be due to hydrostatic pressure screening by the channel walls, known as the “Janssen effect”, and is the first experimental study on its role in granular heat transfer. A 20 kWth small-scale MPBE prototype, the first of its kind with 3 mm particle channels and 2 mm sCO2 plates, was designed and manufactured by Vacuum Process Engineering (VPE) for this project. This MPBE prototype was installed in SNL’s HEX test stand and the MPR tool was integrated into the test stand for in-operando measurements. Several modifications have been made to the MPR setup for facile installation and easy operation at SNL. Measurements with Carbo HSP40/70 particles were conducted with various particle and sCO2 temperatures and flow rates. Further, apart from regular MPBE operation, fault scenarios such as plate warping, particle stagnation, and choking were also detected and MPR showed excellent sensitivity to the detection of such conditions. Through the course of this project, MPR has evolved from a sensitive laboratory technique to a rugged industrial tool with various applications in high temperature measurements and diagnostics.

14 SOLAR ENERGY↗
DOE OSTI2022

Life Expectancy of Evaporating Capillary Bridges Predicted by Tertiary Creep Modeling

The evaporation of capillary bridges is experimentally investigated at the microscale through a three-grain capillary cluster. This setting provides the minimum viable description of Haines jumps during evaporation, that is, capillary instabilities stemming from air entry into a saturated granular material. The displacement profile of a meniscus is obtained via digital image correlation for different grain materials, geometries, and separations. While it is well known that Haines jumps are triggered at the pore throat, we find that these instabilities are of three types depending on the separation. We also provide a temporal characterization of Haines jumps; we find that they are accurately described, as tertiary creep instabilities, by Voight’s relation, similarly to landslides and volcanic eruptions. This finding extends the description of capillary instabilities beyond their onset predicted by Laplace equilibrium. Our contribution also paves the way for a microscopically-informed description of desiccation cracks, of which Haines jumps are the precursors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗