Search NASA⌕ Search

SEARCH · Search NASA

Results for “granular matter”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Protocol-dependent frictional granular jamming simulations: cyclical, compression, and expansion

Granular matter takes many paths to pack in natural and industrial processes. The path influences the packing microstructure, particularly for frictional grains. We perform discrete element modeling simulations of different paths to construct packings of frictional spheres. Specifically, we explore four stress-controlled protocols implementing packing expansions and compressions in various combinations thereof. We characterize the eventual packed states through their dependence of the packing fraction and coordination number on packing pressure, identifying non-monotonicities with pressure that correlate with the fraction of frictional contacts. These stress-controlled, bulk-like particle simulations access very low-pressure packings, namely, the marginally stable limit, and demonstrate the strong protocol dependence of frictional granular matter.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational Investigation of Fuel Dispersal Phenomena during Large-Break Loss of Coolant Accident in Light-Water Reactors

In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Soft algebras for leaf amplitudes

Celestial MHV amplitudes are comprised of non-distributional leaf amplitudes associated to an AdS3 leaf of a foliation of flat spacetime. It is shown here that the leaf amplitudes are governed by the same infinite-dimensional soft ‘S-algebra’ as their celestial counterparts. Moreover, taking the soft limit of the smooth three-point MHV leaf amplitude yields a nondegenerate minus-minus two-point leaf amplitude. The two- and three-point MHV leaf amplitudes are used to compute the plus-minus-minus leaf operator product coefficients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Soft gravitons in the BFSS matrix model

BFSS proposed that asymptotically flat M-theory is dual to a large N limit of the matrix quantum mechanics describing N nonrelativistic D0-branes. Recent insights on the soft symmetries of any quantum theory of gravity in asymptotically flat space are applied to the BFSS matrix model. It is shown that soft gravitons are realized by submatrices whose rank is held fixed in the large N M-theory limit, rather than the usual linear scaling with N for hard gravitons. The soft expansion is identified with the large N expansion and the soft theorem becomes a universal formula for the quantum mechanical scattering of such submatrix excitations. This formula is shown to be the Ward identity of large type IIA U(1) RR asymptotic gauge symmetry in the matrix model, whose asymptotic boundaries are at future and past timelike infinity.

1/N Expansion↗
DOE OSTI2023

Goldstone bosons on celestial sphere and conformal soft theorems

In this paper, we study celestial amplitudes of Goldstone bosons and conformal soft theorems. Motivated by the success of soft bootstrap in momentum space and the important role of the soft limit behavior of tree-level amplitudes, our goal is to extend some of the methods to the celestial sphere. The crucial ingredient of the calculation is the Mellin transformation, which transforms four-dimensional scattering amplitudes to correlation functions of primary operators in the celestial CFT. The soft behavior of the amplitude is then translated to the singularities of the correlator. Only for amplitudes in “UV completed theories” (with sufficiently good high energy behavior) the Mellin integration can be properly performed. In all other cases, the celestial amplitude is only defined in a distributional sense with delta functions. We provide many examples of celestial amplitudes in UV-completed models, including linear sigma models and Z-theory, which is a certain completion of the SU(N) non-linear sigma model. We also comment on the BCFW-like and soft recursion relations for celestial amplitudes and the extension of soft bootstrap ideas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A rheological model for loose sands with insights from DEM

A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.

Elasticity↗
DOE OSTI2023

Detecting low-mass perturbers in cluster lenses using curved arc bases

ABSTRACT Strong gravitationally lensed arcs produced by galaxy clusters have been observationally detected for several decades now. These strong lensing constraints provided high-fidelity mass models for cluster lenses that include substructure down to $10^{9{-}10}\, \mathrm{M}_\odot$. Optimizing lens models, where the cluster mass distribution is modelled by a smooth component and subhaloes associated with the locations of individual cluster galaxies, has enabled deriving the subhalo mass function, providing important constraints on the nature and granularity of dark matter. In this work, we explore and present a novel method to detect and measure individual perturbers (subhaloes, line-of-sight haloes, and wandering supermassive black holes) by exploiting their proximity to highly distorted lensed arcs in galaxy clusters, and by modelling the local lensing distortions with curved arc bases. This method offers the possibility of detecting individual low-mass perturber subhaloes in clusters and haloes along the line of sight down to a mass resolution of $10^8\, \mathrm{M}_\odot$. We quantify our sensitivity to low-mass perturbers ($M\sim 10^{7{-}9}\, \mathrm{M}_\odot$) in clusters (M ∼ 1014–15 M⊙), by creating realistic mock data. Using three lensed images of a background galaxy in the cluster SMACS J0723, taken by JWST, we study the retrieval of the properties of potential perturbers with masses $M=10^{7{-}9}\, \mathrm{M}_\odot$. From the derived posterior probability distributions for the perturber, we constrain its concentration, redshift, and ellipticity. By allowing us to probe lower mass substructures, the use of curved arc bases can lead to powerful constraints on the nature of dark matter as discrimination between dark matter models appears on smaller scales.

Astronomy & Astrophysics↗
DOE OSTI2022

Quantifying local rearrangements in three-dimensional granular materials: Rearrangement measures, correlations, and relationship to stresses

Quantifying the ways in which local particle rearrangements contribute to macroscopic plasticity is one of the fundamental pursuits of granular mechanics and soft matter physics. Here we examine local rearrangements that occur naturally during the deformation of three samples of 3D granular materials subjected to distinct boundary conditions by employing in situ x-ray measurements of particle-resolved structure and stress. We focus on five distinct rearrangement measures, their statistics, interrelationships, contributions to macroscopic deformation, repeatability, and dependence on local structure and stress. Our most significant findings are that local rearrangements (1) are correlated on a scale of three to four particle diameters, (2) exhibit volumetric strain-shear strain and nonaffine displacement-rotation coupling, (3) exhibit correlations that suggest either rearrangement repeatability or that rearrangements span multiple steps of incremental sample strain, and (4) show little dependence on local stress but correlate with quantities describing local structure, such as porosity. Our results are presented in the context of relevant plasticity theories and are consistent with recent findings suggesting that local structure may play at least as important of a role as local stress in determining the nature of local rearrangements.

58 GEOSCIENCES↗
DOE OSTI2024

Qualification and Quantification of Porosity at the Top of the Fuel Pins in Metallic Fuels Using Image Processing

Approximately 130,000 metal fuel pins were irradiated in the Experimental Breeder Reactor II (EBR-II) during its 30 years of operation to develop and characterize existing and prospective fuels. For many of the metal fuel irradiation experiments, neutron radiography imaging was performed to characterize fuel behavior, such as fuel axial expansion. While several fuel expansion results obtained from neutron radiography imaging have been published, the analysis of neutron radiography for the purpose of describing statistical properties of porous matter formed on top of the fuel pins, also referred to as fluff in previous publications, is significantly less represented in the literature with just a single paper so far. This study aims to validate and augment results reported in previous publications using automated image processing. The paper describes the statistical properties of the porous matter in terms of nine parameters derived from radiography images and correlates those parameters with such fuel properties as composition, expansion, temperature, and burnup. The reported results are based on 1097 fuel pins of eight different fuel compositions. For three major fuel types, U-10Zr, U-8Pu-10Zr, and U-19Pu-10Zr, a clear negative correlation is found between the Pu content and five parameters describing the amount of porous matter generated. The parameters describing granularity properties, however, showed either negative correlation or nonlinear dependency from fuel composition. The parameters describing the amount showed a positive correlation with fuel axial expansion, while granularity parameters showed a negative correlation with axial expansion. The dependency on cladding temperature was found to be weak. A positive correlation is demonstrated for volume parameters and fuel burnup. In general, reported results confirm and validate findings published in previous studies using a much larger number of pins and automated processing techniques, which easily lend themselves to reproducibility, thus avoiding subjective bias.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying microbial functional guilds performing cryptic organotrophic and lithotrophic redox cycles in anaerobic granular biofilms

Granular biofilms used in anaerobic digester systems contain diverse microbial populations that interact to hydrolyze organic matter and produce methane within controlled environments. Prior research investigated the feasibility of utilizing granular biofilms obtained from an anaerobic digester to remove nitrate without the addition of exogenous electron donors. These granules possessed a unique structure of alternating light and dark iron sulfide and pyrite rich layers that potentially served as both an electron source and sink, linking carbon, nitrogen, sulfur, and iron cycles. To characterize the functional roles of diverse microbial populations enriched within these layered biofilms, we analyzed metagenomes obtained from three different granules. Comparisons between the functional gene content of forty metagenome assembled genomes (MAGs) identified phylogenetically cohesive functional guilds. Each of these functional MAG clusters was assigned to specific steps in anaerobic digestion (hydrolysis, acidogenesis, acetogenesis, and methanogenesis) and anaerobic respiration (denitrification and sulfate reduction). Comparisons with metagenomes derived from a variety of natural and engineered ecosystems confirmed that the enriched denitrifying bacteria were similar to populations typically found in wetlands and biological nitrogen removal systems. Analysis of read alignments to individual genes within the forty MAGs identified conserved genomic features that were representative of the functions that distinguished functional guilds. Overall, this research illustrates the utility of functional based classification of microorganisms for characterizing ecosystem functions and highlights the potential application of engineered ecosystems to serve as experimental models for complex natural ecosystems.

Ecosystem engineering↗
DOE OSTI2020

Granular packings with sliding, rolling, and twisting friction

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

Here, we examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coupled Dynamical Phase Transitions in Driven Disk Packings

Under the influence of oscillatory shear, a monolayer of frictional granular disks exhibits two dynamical phase transitions: a transition from an initially disordered state to an ordered crystalline state and a dynamic active-absorbing phase transition. Although there is no reason a priori for these to be at the same critical point, they are. Also, the transitions may also be characterized by the disk trajectories, which are nontrivial loops breaking time-reversal invariance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Empty perovskites as Coulomb floppy networks: Entropic elasticity and negative thermal expansion

Floppy networks (FNs) provide valuable insight into the origin of anomalous mechanical and thermal properties in soft matter systems, from polymers, rubber, and biomolecules to glasses and granular materials. Here, we use the same FN concept to construct a quantitative microscopic theory of empty perovskites, a family of crystals with ReO 3 structure, which exhibit a number of unusual properties. One remarkable example is ScF 3 , which shows a near-zero-temperature structural instability and large negative thermal expansion (NTE). Furthermore, we trace these effects to an FN-like crystalline architecture formed by strong nearest-neighbor bonds, which is stabilized by net electrostatic repulsion that plays a role similar to osmotic pressure in polymeric gels. NTE in these crystalline solids, which we conceptualize as Coulomb floppy networks, emerges from the tension effect of Coulomb repulsion combined with the FN's entropic elasticity and has the same physical origin as in gels and rubber. Our theory provides an accurate, quantitative description of phonons, thermal expansion, compressibility, and structural phase diagram, all in excellent agreement with experiments. The entropic stabilization of critical soft modes, which play only a secondary role in NTE, explains the observed phase diagram. Significant entropic elasticity resolves the puzzle of a marked, ≈50% discrepancy between the experimentally observed bulk modulus and ab initio calculations. The Coulomb FN approach is potentially applicable to other important materials with markedly covalent bonds, from perovskite oxides to iron chalcogenides, whose anomalous vibrational and structural properties are still poorly understood.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗