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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗
DOE OSTI2025

Suppression of Composition g -modes in Chemically Equilibrating Warm Neutron Stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.

Neutron stars↗
DOE OSTI2023

Coupled chemo-mechanical modeling of point-defect diffusion in a crystal plasticity fast Fourier transform framework

Below the yield strength and at moderate-to-high homologous temperatures, the inelastic deformation of metals is mostly governed/rate-controlled by vacancy diffusion-mediated processes. As a function of grain size, stress, temperature and dislocation content, vacancies (or atoms) can adopt preferential diffusion paths across grain interiors, along grain boundaries, or towards and along dislocations, resulting in climb and self-climb. In the steady state and under constant load, grain boundary and grain bulk vacancy diffusion-mediated plasticity have been described in seminal works by Coble and by Nabarro and Herring, respectively. Yet, the interplay between all aforementioned potential diffusion pathways has not been comprehensively mapped. This work presents a thermodynamically-consistent full-field model integrated within a voxel-based elasto-viscoplastic fast Fourier transform framework, which considers the coupling between the diffusion-mediated plasticity mechanisms. In the proposed approach, the kinetics and kinematics of plastic deformation due to vacancy diffusion along grain boundaries and grain bulk, as well as the exchange between grain boundaries and bulk are described explicitly. A homogenization approach at the voxel level is further introduced to simultaneously consider bulk and grain boundary diffusion in a numerically efficient fashion. The new formulation predicts the expected strain rate dependencies and the scaling of the steady-state creep rate with respect to grain size, temperature, and stress. Finally, the model predicts the transition from grain bulk to grain boundary-dominated diffusion with reduction in grain size, a significant step towards capturing transitions in deformation behavior without any phenomenological or ad-hoc adjustments.

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

Liquefying Lignins: Determining Phase-Transition Temperatures in the Presence of Aqueous Organic Solvents

Under appropriate conditions of temperature and solvent composition, single-phase aqueous organic solvents are capable of liquefying lignin to form a lignin-rich liquid phase separate from the solvent phase, creating liquid–liquid equilibrium. The extent to which this phenomenon occurs when Kraft lignin is combined with the solvents methanol, ethanol, isopropanol, and acetone was determined, with solid–liquid (SL) to liquid–liquid (LL) phase-transition temperatures being measured over the range of applicable organic–water compositions. Lignins tested from a variety of biomass sources were all discovered to form this liquid–liquid equilibrium with aqueous ethanol solutions. Furthermore, the formation of such phase behavior between lignin and economical, green solvents creates new opportunities for the large-scale purification and fractionation of raw bulk (i.e., technical) lignins.

09 BIOMASS FUELS↗
DOE OSTI2020

Structure and dynamics of ethane confined in silica nanopores in the presence of CO 2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C 2 H 6 ) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C 2 H 6 and CO 2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C 2 H 6 , C 2 H 6 :CO 2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C 2 H 6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C 2 H 6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO 2 molecules are more strongly attracted to the pore surface compared to C 2 H 6 . The C 2 H 6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO 2 molecules play in enhancing C 2 H 6 diffusion (“molecular lubrication”) at high CO 2 :C 2 H 6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

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

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗
DOE OSTI2023

Influence of Single-Crystalline Morphology on the Electrochemical Behavior of High-Nickel Layered Oxide Cathodes

High-Nickel layered oxides are being pursued for their high specific capacities. With increasing nickel content, the layered oxides also suffer from structural and morphological instabilities, which degrade the electrochemical performance. Single-crystalline morphology presents an attractive solution to these problems, simultaneously reducing surface area available for parasitic reactions and preventing electrolyte penetration into the bulk with the removal of grain boundaries. Single-crystal cathodes have already shown promising results with NMC composition cathodes, but there are still gaps in the fundamental understanding of how single-crystalline morphology alters electrochemical behavior. We attempt to fill some of these gaps by studying in detail the electrochemical operation of polycrystalline and single-crystalline LiNiO 2 (LNO). Single-crystal LNO is prepared through a molten salt method yielding high-quality, distinct crystals with comparable lattice chemistry to polycrystalline LNO. The single-crystal LNO achieves far greater stability during long-term cycling, reaching 500 cycles with 82.5% capacity retention. Interestingly, the single-crystalline LNO also displays superior rate performance, delivering 157 mA h g −1 at 10C discharge rate. Investigation of the phase evolution behavior during cycling strongly suggests that the absence of grain boundaries in single-crystalline LNO is responsible for the superior performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Thermal-Strain-Enabled Enhanced Emission from UV Laser-Induced Defect Levels near the Surface of Multilayer MoS 2

Monolayer two-dimensional (2D) materials have been intensively studied while research on multilayers is still in its infancy. Here, we induce defects inside bulk MoS 2 through thermal annealing and near the surface of multilayer MoS 2 using 375 nm laser irradiation, and investigate their photoluminescence (PL) and fluorescence lifetime imaging (FLIM). Enhanced emission is limited within a certain MoS 2 thickness. The observed enhanced emission is evidenced by a threshold behavior in super-linear PL intensity increase, strong polarization effects, and increased lifetime of defect peak. The laser power threshold for enhanced emission is much smaller in defects near the surface than that inside the bulk of multilayer MoS 2 . The mechanical strain from a wrinkle of the sample further lowers the laser power threshold for enhanced emission. By exciting with a 639 nm laser that is close to the fundamental gap between the conduction band minimum and the valence band maximum, the lifetime of defect enhanced emission increased by 5 times. Furthermore, one of the competing indirect bandgap emissions disappears, and the defect emission peak dominates the PL spectrum in the wrinkle area with a strain. Furthermore, the discovered principle can be applied to future studies on the integration of enhanced emission and single photon emission involving selectively depopulating the conduction band of the host crystal to defect levels for quantum emitters.

2D materials↗
DOE OSTI2020

Unexpected Conformational Behavior of Poly(poly(ethylene glycol) methacrylate)-Poly(propylene carbonate)-Poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA) Amphiphilic Block Copolymers in Micellar Solution and at the Air-Water Interface

Self-assemblies (bulk micelles and Langmuir monolayers) of a new amphiphilic block copolymer, poly(poly(ethylene glycol) methacrylate)-poly(propylene carbonate)-poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA), were investigated in dilute aqueous solution and at the air-water interface. An important feature that distinguishes PPEGMA-PPC-PPEGMA from conventional linear PEG-based polymer surfactants is that in the hydrophilic PPEGMA block the PEG moieties exist as side chains attached to the poly(methacrylate) (PMA) backbone. This study was performed using three PPEGMA-PPC-PPEGMA samples having an identical PPC molecular weight (5.6 kDa) and different PPEGMA molecular weights (7.3, 2.8 and 2.1 kDa on either side) (named as “G7C6G7”, “G3C5G3”, and “G2C6G2”, respectively). Cryo-TEM images revealed that G2C6G2 micelles form large clusters, whereas G7C6G7 micelles exist as unimers; the behavior of G3C6G3 is intermediate. The average core surface distance between two adjacent G2C6G2 micelles in a cluster was estimated to be only about 3.7 nm, which suggests that in the micelle corona layer only the PEG side chains (but not the PMA backbone) exist as hydrated brush chains. A calculation showed that at the small micelle separation distance observed, the van der Waals attractive forces become significant, and are thus likely responsible for the cluster formation. Steady shear rheometry measurements confirmed that G2C6G2 micelles are indeed only weakly bound in clusters, and these micelle clusters are thus prone to breakup even under mild steady shear. The micelle breakup force estimated from the shear viscosity data is in reasonable agreement with an estimate obtained from the theoretical van der Waals interaction potential. Langmuir monolayers formed at the air-water interface by PPEGMA-PPC-PPEGMA were also characterized; surface pressure-area measurements were performed. Detailed features of the surface pressure-area isotherms further supported that only the PEG side chains are hydrated in the subphase, while the PMA backbone remains unhydrated and situated on the water surface.

Lee, Jaewon↗
DOE OSTI2025

Elucidating the Competitive Hydrodeoxygenation of Lignin-Derived Oxygenates over Bulk MoO 3 Catalyst through Kinetic Analysis

Ambient pressure hydrodeoxygenation (HDO) of lignin-derived oxygenates over molybdenum oxide-based catalysts is an effective strategy to produce chemicals that can be directly integrated into our existing petrochemical infrastructure. Complexities pertaining to the simultaneous kinetic and mechanistic analysis of HDO have limited research endeavors to single-compound systems. Although valuable insight into the catalytic reaction has been gained through this approach, it provides limited understanding of the competitive adsorption behavior manifest in a realistic multioxygenate reaction environment. To address this shortcoming, simultaneous gas-phase acetone and anisole HDO was performed at 330 °C and ≤1 bar H 2 partial pressure over bulk MoO 3 . Propene, propane, and benzene were the HDO products formed, showing a similar product distribution to the single-compound system. Selectivity to propene and propane was ∼14 times higher than benzene, even at three times higher anisole partial pressure compared to acetone. A negative anisole HDO (−0.97 ± 0.22) rate order with varying acetone partial pressure suggested a strong inhibition effect on anisole HDO by acetone. Conversely, with increasing anisole partial pressure, a rate order of −0.07 ± 0.12 was observed for acetone HDO, implying a weak impact of anisole cofeed on acetone HDO. A kinetic-driven approach was taken to estimate the relative adsorption constants of the oxygenates. Acetone exhibited a 6.4 times higher adsorption propensity on the HDO active site than anisole. Relative adsorption constants for phenolics increased with increasing basicity of the oxygenate but decreased for aliphatic molecules, suggesting a volcano-shaped relationship. The results suggest the possibility of an optimal electron density around the molecule’s oxygen atom to maximize the molecule’s adsorption strength.

acid sites↗
DOE OSTI2021

Characteristics of Non-Fullerene Acceptor-Based Organic Photovoltaic Active Layers Using X-ray Scattering and Solid-State NMR

With the recent development of non-fullerene acceptors, power conversion efficiencies of bulk-heterojunction organic solar cells have exceeded 18%. The morphology of the BHJ active layer, including packing, ordering, orientation, and phase behavior of the donor(s) and acceptor(s), plays critical roles in determining the device performance. We characterized the morphology of active layers consisting of mixtures of a PTB7/PTB7-Th (donor) and ITIC (acceptor) using grazing incidence wide-angle X-ray scattering, resonant soft X-ray scattering, and solid-state nuclear magnetic resonance (ssNMR) to correlate the morphology with device performance. PTB7-Th/ITIC shows a better device performance than that of PTB7/ITIC due to a smaller π-πstacking distance and smaller domain size of the phase separated morphology. One of PTB7/ITIC samples, processed from mixed solvents (chlorobenzene/benzene = 1:1 by volume), shows a possible partial miscibility and smaller domain size as revealed by ssNMR T 1ρ relaxation times, which is detrimental to device performance. ssNMR results showed that ITIC could crystallize into different forms depending on processing conditions, which may have implications on the manufacturing of devices using it as an ingredient, as well as its long-term stability in service.

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

First-Principles Characterization of Surface Phonons of Halide Perovskite CsPbI 3 and Their Role in Stabilization

The stability of halide perovskites has been a long-standing issue for their real-world application. Approaches to improve stability include nanostructuring, dimensionality reduction, and strain engineering, where surfaces play an important role in the formation of a stable structure. To understand the mechanism we compute the lattice dynamics of the surface of CsPbI 3 using density functional theory. We demonstrate, for the first time, that CsPbI 3 crystals exhibit surface phonons that are localized on the outermost layers of the slabs, and we perform a complete symmetry characterization including an identification of the Raman/IR active modes. These surface phonons are present in the optically active cubic phase but are absent in the optically inactive "yellow" phase. Furthermore, we show that the surface suppresses bulk instabilities by hardening soft modes of the bulk cubic phase, resulting in phase stabilization and quenching of dynamical disorder. Finally, this study is fundamental for understanding the structural behavior of halide perovskite materials with high surface area-to-volume ratios, and for guiding stabilization strategies.

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

To Substitute Rather Than Intercalate: Chimie douce Approach to Induce Ferromagnetism in Metastable Pt 0.8 M 0.2 Se 2 ( M = Cr, Co, Ni)

Two-dimensional (2D) magnetic materials with exotic magnetic properties have garnered significant interest due to their potential applications in spintronics and data storage technologies. However, the limited availability of intrinsic 2D magnetic materials has driven efforts to induce and manipulate magnetism in otherwise nonmagnetic 2D systems through approaches such as chemical intercalation, defect engineering, and substitutional doping. Herein, we present a facile, chimie douce method for incorporating 3d transition metals (Cr, Co, and Ni) into the nonmagnetic PtSe 2 sublattice. This synthetic approach enables control over layer thickness of Pt 1–x M x Se 2 (M = Cr, Co, Ni) nanosheets by varying the M identity and annealing conditions. Comprehensive scattering and spectroscopic characterizations confirm the successful and homogeneous substitution of M atoms at the Pt site, rather than intercalation, and reveal a strong correlation between nanosheet thickness and the identity of the substituting metal. High-temperature annealing of the nanosheets promotes an irreversible transformation toward the bulk phase, allowing for detailed characterization of structural and magnetic properties. A case study of Pt 0.8 Cr 0.2 Se 2 reveals that nanosheet thickness plays a critical role in modulating local magnetic interactions. While Cr atoms in the as-synthesized few-layers-thick nanosheets exhibit predominantly short-range antiferromagnetic interactions, the emergence of short-range ferromagnetic exchange is revealed in the bulk material. Detailed ac susceptibility and remanent magnetization measurements further demonstrate that bulk Pt 0.8 Cr 0.2 Se 2 adopts a frustrated magnetic ground state with clear signatures of ferromagnetic cluster-glass behavior. The systematic investigation presented herein establishes a clear and robust protocol for the synthesis and in-depth characterization of 2D transition-metal-substituted PtSe2 materials with varying layer thickness and paves a path toward their realization in spintronic and magnetic device applications.

crystallinity↗
DOE OSTI2023

Sequence-dependent self-assembly of supramolecular nanofibers in periodic dynamic block copolymers

Block copolymers (BCPs) can spontaneously self-assemble into various nanostructured morphologies which enable their diverse applications in selective membranes, polymer electrolytes, and optoelectronics. To expand the range of nanostructures accessible to block copolymers, we designed dynamic block copolymers (DBCPs) that combine the phase separation of traditional block copolymers with the supramolecular self-assembly of periodic dynamic polymers. We demonstrate that DBCPs synthesized with a periodic block sequence self-assemble into high aspect ratio supramolecular nanofibers with well-ordered PEG and PDMS domains, in contrast to those synthesized with a random sequence which do not form nanofibers but have disordered morphologies. The periodicity of the block sequence ensures regular placement of dynamic bonds along the polymer chain, which enables stacking of the hydrogen bonding units into ordered 1D supramolecular nanofibers and delays the onset of terminal flow by up to 60 °C compared to the random block sequence. The hierarchically assembled supramolecular nanofibers display complex mechanical and thermal phase behavior arising from the interplay between phase separation of the dissimilar polymer backbones and supramolecular interactions between dynamic bonds. Despite identical bulk composition, the periodic DBCPs demonstrate ionic conductivity values over two orders of magnitude higher than their random counterparts due to the formation of well-ordered, interconnected, high aspect ratio ion-transporting PEG domains. In conclusion, these results highlight the potential for DBCPs as an emerging material platform to achieve advanced, self-assembled nanostructured morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Lattice chain theories for dynamics of acoustic flexural phonons in nonpolar nanomaterials

Knowledge of lattice dynamics of acoustic flexural phonons is critically important in understanding the mechanical, electrical, and thermal properties of nanomaterials. However, the general nature of acoustic flexural phonons for the span of possible nanomaterials and the dependence of these on microscopic material properties remain unclear. Therefore, we develop general lattice chain theories to understand acoustic flexural vibrations in covalent and van der Waals lattices. Furthermore, the analytic theories suggest that the flexure branch has quadratic dispersion for low dimensional lattices in the long-wavelength limit, quite different from linear branches in bulk van der Waals bonded lattices like graphite. Surprisingly, interlayer van der Waals interactions in multilayer lattices do not determine the behavior of the flexural dispersion. The bonds of covalent interactions in multilayer and buckled monolayer lattices play different roles in determining the flexural dispersion with those in flat monolayer lattices, e.g., graphene. Further, based on the correspondence between microscopic and continuum dynamics, we provide a universal and effective approach to characterize the intrinsic bending rigidity of low dimensional nanomaterials using dispersion data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multiscale Shear Properties and Flow Performance of Milled Woody Biomass

One dominant challenge facing the development of biorefineries is achieving consistent system throughput with highly variant biomass feedstock quality and handling performance. Current handling unit operations are adapted from other sectors (primarily agriculture), where some simplifying assumptions about granular mechanics and flow performance do not translate well to a highly compressible and anisotropic material with nonlinear time- and stress-dependent properties. This work explores the shear and frictional properties of loblolly pine at multiple experimental test apparatus and particle scales to elucidate a property window that defines the shear behavior over a range of material attributes (particle size, size distribution, moisture content, etc.). In general, it was observed that the bulk internal friction and apparent cohesion depend strongly on both the stress state of the sample in granular shear testers and the overall particle size and distribution span. For equipment designed to characterize the quasi-static shear stress failure of bulk materials ranging from 50 to 1,000 ml in test volume, similar test results were observed for finely milled particles (50% passing size of 1.4 mm) with a narrow size distribution (span between 10 and 90% passing size of 0.9 mm), while stress chaining and over-torque issues persisted for the bench-scale test apparatus for larger particle sizes or widely dispersed sample sizes. Measurement of the anisotropic particle–particle friction ranged from coefficients of approximately 0.20 to 0.45 and resulted in significantly higher and more variable friction measurements for larger particle sizes and in perpendicular alignment orientations. To supplement these laboratory-scale properties, this work explores the flow of loblolly pine and Douglas fir through a pilot-scale wedge-shaped hopper and a screw feeder. For the gravity-driven hopper flow, the critical arching distance and mass discharge rate ranged from approximately 10 to 30 mm and 2 to 16 tons/hour, respectively, for both materials, where the arching distance depends strongly on the overall particle size and depends less on the hopper inclination angle. Comparatively, the auger feeder was found to be much more impacted by the size of the particles, where smaller particles had a more consistent and stable flow while consuming less power.

09 BIOMASS FUELS↗
DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗