Jamming is the transition from a fluid‐like state to a solid‐like state of a packing system. Recent studies have shown that jamming transition depends upon many factors: particle shape, friction/cohesion between particles, particle size dispersity, the stress of the packing, etc. This study aims to contribute to this growing area of research by exploring the jamming density of soil with strong dispersity. In analogous to Gibbs excess energy, we introduce excess volume‐potentials for each species. We then proposed a mathematical model to quantitatively compute the jamming density based on the second law of equilibrium in thermodynamics. This approach is validated using experimental results on glass beads and on silty sand. It is hoped that this study will provide to a deeper understanding of the link between jamming density, packing dispersity and the second law of thermodynamics.
This content will become publicly available on January 1, 2025
- NSF-PAR ID:
- 10528003
- Editor(s):
- Bellwied, R; Geurts, F; Rapp, R; Ratti, C; Timmins, A; Vitev, I
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- EPJ Web of Conferences
- Volume:
- 296
- ISSN:
- 2100-014X
- Page Range / eLocation ID:
- 14014
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
To explore the curvature dependence of solid–fluid interfacial thermodynamics, we calculate, using Grand Canonical Monte Carlo simulation, the surface free energy for a 2d hard-disk fluid confined in a circular hard container of radius R as a function of the bulk packing fraction η and wall curvature C̄=−1/R. (The curvature is negative because the surface is concave.) Combining this with our previous data [Martin et al., J. Phys. Chem. B 124, 7938–7947 (2020)] for the positive curvature case (a hard-disk fluid at a circular wall, C̄=+1/R), we obtain a complete picture of surface thermodynamics in this system over the full range of positive and negative wall curvatures. Our results show that γ is linear in C̄ with a slope that is the same for both positive and negative wall curvatures, with deviations seen only at high negative curvatures (strong confinement) and high density. This observation indicates that the surface thermodynamics of this system is consistent with the predictions of so-called morphometric thermodynamics at both positive and negative curvatures. In addition, we show that classical density functional theory and a generalized scaled particle theory can be constructed that give excellent agreement with the simulation data over most of the range of curvatures and densities. For extremely high curvatures, where only one or two disks can occupy the container at maximum packing, it is possible to calculate γ exactly. In this limit, the simulations and density functional theory calculations are in remarkable agreement with the exact results.more » « less
-
more » « less
Abstract The 3D reference interaction site model (3D‐RISM) of molecular solvation is a powerful tool for computing the equilibrium thermodynamics and density distributions of solvents, such as water and co‐ions, around solute molecules. However, 3D‐RISM solutions can be expensive to calculate, especially for proteins and other large molecules where calculating the potential energy between solute and solvent requires more than half the computation time. To address this problem, we have developed and implemented treecode summation for long‐range interactions and analytically corrected cut‐offs for short‐range interactions to accelerate the potential energy and long‐range asymptotics calculations in non‐periodic 3D‐RISM in the AmberTools molecular modeling suite. For the largest single protein considered in this work, tubulin, the total computation time was reduced by a factor of 4. In addition, parallel calculations with these new methods scale almost linearly and the iterative solver remains the largest impediment to parallel scaling. To demonstrate the utility of our approach for large systems, we used 3D‐RISM to calculate the solvation thermodynamics and density distribution of 7‐ring microtubule, consisting of 910 tubulin dimers, over 1.2 million atoms.
-
more » « less
Conservation of total, kinetic, and thermal energy in the atmosphere is revisited, and the derived thermal energy balance is examined with observations. Total energy conservation (TEC) provides a constraint for the sum of kinetic, thermal, and potential energy changes. In response to air thermal expansion/compression, air density variation leads to vertical density fluxes and potential energy changes, which in turn impact the thermal energy balance as well as the kinetic energy balance due to the constraint of TEC. As vertical density fluxes can propagate through a large vertical domain to where local thermal expansion/compression becomes negligibly small, interactions between kinetic and thermal energy changes in determining atmospheric motions and thermodynamic structures can occur when local diabatic heating/cooling becomes small. The contribution of vertical density fluxes to the kinetic energy balance is sometimes considered but that to the thermal energy balance is traditionally missed. Misinterpretation between air thermal expansion/compression and incompressibility for air volume changes with pressure under a constant temperature would lead to overlooking important impacts of thermal expansion/compression on air motions and atmospheric thermodynamics. Atmospheric boundary layer observations qualitatively confirm the contribution of potential energy changes associated with vertical density fluxes in the thermal energy balance for explaining temporal variations of air temperature.
-
null (Ed.)Here we revisit the thermodynamics of the Kitaev quantum spin liquid realized on the honeycomb lattice. We address two main questions: First, we investigate whether there are observable thermodynamic signatures of the topological Majorana boundary modes of the Kitaev honeycomb model. We argue that for the time-reversal invariant case the residual low-temperature entropy is the primary thermodynamic signature of these Majorana edge modes. Using large-scale Monte Carlo simulations, we verify that this residual entropy is present in the full Kitaev model. When time-reversal symmetry is broken, the Majorana edge modes are potentially observable in more direct thermodynamic measurements such as the specific heat, though only at temperatures well below the bulk gap. Second, we study the energetics, and the corresponding thermodynamic signatures, of the flux excitations in the Kitaev model. Specifically, we study the flux interactions on both cylinder and torus geometries numerically and quantify their impact on the thermodynamics of the Kitaev spin liquid by using a polynomial fit for the average flux energy as a function of flux density and extrapolating it to the thermodynamic limit. By comparing this model to Monte Carlo simulations, we find that flux interactions have a significant quantitative impact on the shape and the position of the low-temperature peak in the specific heat.more » « less
