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.
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.
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.
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.
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.
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.
Self-Organized Stress Distributions in Polycrystalline Materials [Dissertation]
Understanding stress distributions in solid materials is complicated by the fact that most materials are polycrystalline in nature, with each crystal having an elastically anisotropic reaction to force. This study is to gain a better understanding on how external forces placed on a polycrystal are related to internal reactions within and between the grains. The hypothesis is that the stress distribution in porous to fully dense materials are self-organized based on strong contacts between and within the individual grains created by force chains. Force chains, commonly known in loaded granular materials, and could be the phenomenon that connect micro to macro deformation. Scale bridging measurements conducted through Raman spectroscopy, Atomic Force Microscopy, and Digital Image Correlation will be used to create stress maps, modulus maps, and elastic strain maps across a variety of geological and pharmaceutical polycrystals. When possible, the resulting maps will be compared to current homogenization schemes and a full field models. Finite element modeling will be used to assess whether the patterning seen in the experimental map is a reasonable approximation based on the orientation data of the samples used. A minimum of three publications is projected to be accomplished focusing each on a different method to experimentally test and analyze stress distributions.
STS-42 Space Shuttle mission report
The STS-42 Space Shuttle Program Mission Report contains a summary of the vehicle subsystem operations during the forty-fifth flight of the Space Shuttle Program and the fourteenth flight of the Orbiter vehicle Discovery (OV-103). In addition to the Discovery vehicle, the flight vehicle consisted of the following: an External Tank (ET) designated as ET-52 (LWT-45); three Space Shuttle main engines (SSME's), which were serial numbers 2026, 2022, and 2027 in positions 1, 2, and 3, respectively; and two Solid Rocket Boosters (SRB's) designated as BI-048. The lightweight redesigned Solid Rocket Motors (RSRM's) installed in each one of the SRB's were designated as 360L020A for the left SRM and 360Q020B for the right SRM. The primary objective of the STS-42 mission was to complete the objectives of the first International Microgravity Laboratory (IML-1). Secondary objectives were to perform all operations necessary to support the requirements of the following: Gelation of Sols: Applied Microgravity Research (GOSAMR); Student Experiment 81-09 (Convection in Zero Gravity); Student Experiment 83-02 (Capillary Rise of Liquid Through Granular Porous Media); the Investigation into Polymer Membrane Processing (IPMP); the Radiation Monitoring Equipment-3 (RME-3); and Get-Away Special (GAS) payloads carried on the GAS Beam Assembly.
Iodine Capture by Ag-Loaded Solid Sorbents Followed by Ag Recycling and Iodine Immobilization: An End-to-End Process
Here, this article demonstrates the complete process of I 2(g) capture using Ag-exchanged faujasite (AgX) and Ag-exchanged mordenite (AgZ) zeolite sorbents, the recovery of Ag for recyclability of iodine capture, and the immobilization of iodide (from the elution process) into an iodosodalite waste form, which is chemically and environmentally stable. The gaseous iodine, I 2(g) , was captured in the AgX via chemisorption by Ag and converted to AgI within the aluminosilicate zeolite frameworks. The elution reaction in Na 2 S resulted in the precipitation of Ag 2 S and dissolution of I – and Na + into the aqueous solution, where powdered sorbent showed higher conversion efficiency than the granular form, which is attributed to Na 2 S-solution diffusion limitations in the granules. The Ag 2 S-containing aluminosilicate precipitates were filtered out of the solution, and the filtrate (aqueous solution) was used to synthesize iodosodalite to immobilize iodide. The iodine capture using the recovered Ag 2 S-containing sorbents shows equivalent iodine loadings to the AgX starting material of Q e = 0.355 g/g for powdered material and 0.255 g/g for granular material, demonstrating the potential recycling of Ag for iodine capture as Ag 2 S-based sorbents.
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.
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.
A finite micro-rotation material point method for micropolar solid and fluid dynamics with three-dimensional evolving contacts and free surfaces
This paper introduces an explicit material point method designed specifically for simulating the micropolar continuum dynamics in the finite deformation and finite microrotation regime. The material point method enables us to simulate large deformation problems while circumventing the potential mesh distortion without remeshing. To eliminate rotational motion damping and loss of angular momentum during the projection, we introduce the mapping for microinertia and angular momentum between particles and grids through the affine particle-in-cell approach. The microrotation and the curvature at each particle are updated through zero-order forward integration of the microgyration and its spatial gradient. We show that the microinertia and the angular momentum are conserved during the projections between particles and grids in our formulation. We verify the formulation and implementation by comparing with the analytical dispersion relation of micropolar waves under the small strain and small microrotation, as well as the analytical soliton solution for solids undergoing large deformation and large microrotation. Additionally we also demonstrate the capacity of the proposed computational framework to handle a wide spectrum of simulations that exhibit size effects in the geometrical nonlinear regime through three representative numerical examples, i.e., a cantilever beam torsion problem, a fragment-impact penetration problem, and a micropolar fluid discharging problem.
Size distributions in two porous chondritic micrometeorites
Quantitative size measurements of granular units (GUs), and nm-sized minerals in these units, in two porous chondritic micrometeorites are investigated. The matrix of these micrometeorites consist of loosely packed, 0.1 micron-sized, GUs. These objects were a major component of the solar nebula dust that accreted into protoplanets. The matrix in micrometeorite W7010*A2 has a fractal dimension with a small coefficient that supports efficient sticking of carbon-rich GUs during accretion. The fractal nature of the matrix provides a way to calculate the density using the aggregate size. The resulting very low density for porous chondritic micrometeorites is 0.08-0.14 g/cu cm, which supports the view that they are the solid debris from unconsolidated solar system bodies. Chondritic GUs contain ultrafine olivines, pyroxenes, and sulfides, embedded in hydrocarbons and amorphous carbons. Nanocrystals in the micrometeorites W7010*A2 and U2015*B show log normal size distributions. The high incidence of disk-shaped grains, a changeover from disk-shaped to euhedral grains, the unevolved nature of the size distributions, and multiple populations for grains less than 127 nm in size, are consistent with continuous postaccretion nucleation and growth in amorphous GUs, including coarsening via Ostwald ripening.
Constitutive and Stability Behavior of Soils in Microgravity Environment
All aspects of soil stability, bearing capacity, slope stability, the supporting capacity of deep foundations, and penetration resistance depend on soil strength. The stress-deformation and stress-deformation-time behavior of soils are of importance in any problem where ground movements are of interest. In most engineering materials, the strength is derived from internal chemical and physico-chemical forces of interaction, which bond the atoms, molecules, and particles together. In soils, the constitutive relations are mainly derived from interparticle friction between particles and particle groups and dilatancy, and to a lesser extent from particle bonding by weak electrostatic, physico-chemical, and coulomb forces. For engineering purposes, soils are classified as cohesive (clays and silts; typical particle sizes range from 10 nm to 10 micrometers) and cohesionless (sand and gravel; typical particle sizes range from 10 micrometers to 75 mm). The mechanical or constitutive properties of cohesionless soils or granular materials are highly fabric-dependent, highly non-linear, and non-conservative with engineering properties primarily depending on the effects of gravity through self-weight and on the tractions or forces applied to the soil mass. Under moderate-to-high stress levels, the influence of gravity on the behavior of laboratory test specimens may not be pronounced and, therefore, the test results in terrestrial (1-g) environment may be sufficiently conclusive. However at low interparticle stresses, which can result either from low applied (confining) stresses or from excess pore fluid pressures developed within the soil mass without corresponding changes in the applied stresses, the presence of gravitational body forces acting on solid particles and interstitial fluids exerts a pronounced influence on movement of individual particles or particle groups. Such motions, in turn, cause changes in soil fabric which results in significant changes in the interparticle friction forces contributing to the soil's strength and deformation characteristics.
Evaluating the Adsorptive Capabilites of Chemsorb 1000 and Chemsorb 1425
The removal of trace contaminants from spacecraft cabin air is necessary for crew health and comfort during long duration space exploration missions. The air revitalization technologies used in these future exploration missions will evolve from current ISS ISS State-of-Art (SOA) and is being designed and tested by the Advanced Exploration Systems (AES) Program's Atmosphere Resource Recovery and Environmental Monitoring (ARREM) project. The ARREM project is working to mature optimum atmosphere revitalization and environmental monitoring system architectures to enable exploration beyond Lower Earth Orbit (LEO). The Air Revitalization Lab at KSC is one of six NASA field centers participating in the ARREM that specializes in adsorbent and catalyst characterization with simulated spacecraft gas streams using combinations of pressure, O2 partial pressure, CO2 partial pressure, and humidity that are representative of a range of anticipated cabin atmospheric conditions and loads. On board ISS, the Trace Contaminant Control Subassembly (TCCS) provides active control of trace contaminants from the cabin atmosphere utilizing physical adsorption, thermal catalytic oxidation, and chemical adsorption processes. High molecular weight contaminants and ammonia (NH3) are removed a granular activated carbon treated with approx. 10% by weight phosphoric acid (H3PO4) (B-S Type 3032 4×6 mesh), which is expendable and is periodically refurbished. The Type 3032 granular activated carbon bed is no longer commercially available and therefore it is important to characterize the efficiency and capacity of commercially available NH3 sorbents. This paper describes the characterization of two Molecular Products LTD activated carbons: Chemsorb 1000 and Chemsorb 1425. Untreated activated carbons (e.g. Chemsorb 1000) remove contaminants by physisorption, which concentrates the contaminant within the pores of the carbon while letting air to pass through the sorbent4. Low molecular weight or polar gases (e.g. HCl, SO2, formaldehyde, and NH3) are not removed by physisorption and typically require chemisorption for removal. Treated activated carbons (e.g. Chemsorb 1425) are impregnated with a a chemical agent (e.g. phosphoric acid) that reacts with those gases, converting them to solids or salts within the carbon and removes them from the air stream. This process occurs via neutralization or catalysis reactions and adsorption capacity is exhaustedwhen the available impregnated chemicals are consumed. Moisture affects removal performance since adsorption sites within the pores are filled with water. The performance of impregnated carbons may be enhanced by moisture content because the mechanisms of contaminant removal are chemical reactions that occur in reagents contained within the pores. The adsorptive capacity data (mol/kg) of Chemsorb 1000 and 1425 for gas mixtures (ethanol, acetone, toluene, acetaldehyde, dichloromethane, and xylene) was measured with 40% relative humidity at 23 deg C air temperature. The adsorptive capacity data (mol/kg) of Chemsorb 1425 was measured using NH3 gas streams.
Evidence for Methane Segregation at the Surface of Pluto
In May 1995, a set of spectrophotometric curves of the system Pluto-Charon was recorded with the UKIRT telescope equipped with the spectrometer CGS4. As for the previous observations, the spectra cover a part of the near infrared range, between 1.4 and 2.55 micrometers, but with a higher resolution of approximately 700. In both the 1992 and 1995 data, the existence of solid methane is confirmed by numerous absorption bands, and the carbon monoxide and the nitrogen ices are identified by their respective signatures at 2.35 and 2.15 um. The solid nitrogen seems to be the principal icy component and forms a matrix in which the CH4 and CO molecules are diluted. However a spectroscopic analysis of the 1995 observations indicates that pure methane may coexist with its diluted phase in N2. In order to derive the horizontal and vertical distribution of these different species and to obtain some quantitative information about their characteristics, we have modeled the spectrum of May 15 that corresponds to the maximum of Pluto's visible light curve. This was achieved by means of a radiative transfer algorithm dealing with compact and stratified media. Among the various representations we have tested to describe the surface of Pluto, only a geographical mixture of three distinct units explains all the significant structures of the analyzed spectrum. The first unit is a thin granular layer of pure CH4 covering a compact polycrystalline substratum of N2-CH4-CO, which are in a molecular mixture (concentrations of and CO of the order of 0.45%, 0.1-0.2% respectively). It covers about 70% of the observed area and corresponds to volatile deposits that are sublimating under solar illumination. The second unit is either (a) a single thick layer of pure granular methane or (b) a unit similar to the first unit but with the two components inverted (i.e. with CH4 forming a substratum and the N2-CH4-CO mixture a superficial layer of fine grains). Covering 20% of the surface, it represents some old surfaces that have been sublimated for a long time, and eventually recovered later by very small amounts of fresh deposits of the molecular mixture N2-CH4-CO. Finally, the third unit may result from the condensation of very fine grains of nearly pure N2. It covers the remainder of the surface (about 10%). All these results allow a better understanding of the processes of deposition, metamorphism, sublimation and transport affecting the different ices detected on Pluto during its climatic cycles.
Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering
Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.
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.