Estimating Porosity of Agglomerated Products Using Optimized Sphere Packing

All types of homogeneous sphere packing and interpenetrating sphere packings and layers were derived that correspond to point configurations of the 15 trivariant hexagonal lattice complexes. The respective sphere packings are assigned to 147 types. In total, sphere packings of 170 types can be realized with hexagonal symmetry. 103 types of sphere packing refer exclusively to trivariant hexagonal lattice complexes. For 23 of these types, the corresponding sphere packings can be generated only in hexagonal lattice complexes with less than three degrees of freedom or with trigonal or lower symmetry. In addition, seven types of interpenetrating sphere packings and two types of interpenetrating sphere layers were found. Interpenetrating 4.8(2) nets of spheres with 120 degrees angles between the nets were assumed to be not possible, so far. The sphere packings belonging to 85 of the 170 hexagonal types can be split up into parallel layers of spheres with mutual contact and can be characterized by symbols derived from those for the Shubnikov nets. The sphere packings of 135 hexagonal types may be subdivided into rod-like subsets of spheres with mutual contact. Such rods may be described by rolling up a plane net. Only 23 types of sphere packing cannot be symbolized on the basis of layers or rods of spheres with mutual contact. Examples are given for crystal structures that can be described by means of sphere packings.

In a series of recent breakthroughs, Allen-Zhu and Orecchia [Allen-Zhu/Orecchia, STOC 2015; Allen-Zhu/Orecchia, SODA 2015] leveraged insights from the linear coupling method [Allen-Zhu/Oreccia, arXiv 2014], which is a first-order optimization scheme, to provide improved algorithms for packing and covering linear programs. The result in [Allen-Zhu/Orecchia, STOC 2015] is particularly interesting, as the algorithm for packing LP achieves both width-independence and Nesterov-like acceleration, which was not known to be possible before. Somewhat surprisingly, however, while the dependence of the convergence rate on the error parameter epsilon for packing problems was improved to O(1/epsilon), which corresponds to what accelerated gradient methods are designed to achieve, the dependence for covering problems was only improved to O(1/epsilon^{1.5}), and even that required a different more complicated algorithm, rather than from Nesterov-like acceleration. Given the primal-dual connection between packing and covering problems and since previous algorithms for these very related problems have led to the same epsilon dependence, this discrepancy is surprising, and it leaves open the question of the exact role that the linear coupling is playing in coordinating the complementary gradient and mirror descent step of the algorithm. In this paper, we clarify these issues, illustrating that the linear coupling method can lead to improved O(1/epsilon) dependence for both packing and covering problems in a unified manner, i.e., with the same algorithm and almost identical analysis. Our main technical result is a novel dimension lifting method that reduces the coordinate-wise diameters of the feasible region for covering LPs, which is the key structural property to enable the same Nesterov-like acceleration as in the case of packing LPs. The technique is of independent interest and that may be useful in applying the accelerated linear coupling method to other combinatorial problems.

A simulation study of the effects of dynamic variables on the packing of spheres

We present a model of (modified) gravity on spacetimes with fractal structure based on packing of spheres, which are (Euclidean) variants of the packed swiss cheese cosmology models. As the action functional for gravity we consider the spectral action of noncommutative geometry, and we compute its expansion on a space obtained as an Apollonian packing of three-dimensional spheres inside a four-dimensional ball. Using information from the zeta function of the Dirac operator of the spectral triple, we compute the leading terms in the asymptotic expansion of the spectral action. They consist of a zeta regularization of the divergent sum of the leading terms of the spectral actions of the individual spheres in the packing. This accounts for the contribution of points 1 and 3 in the dimension spectrum (as in the case of a 3-sphere). There is an additional term coming from the residue at the additional point in the real dimension spectrum that corresponds to the packing constant, as well as a series of fluctuations coming from log-periodic oscillations, created by the points of the dimension spectrum that are off the real line. These terms detect the fractality of the residue set of the sphere packing. We show that the presence of fractality influences the shape of the slow-roll potential for inflation, obtained from the spectral action. We also discuss the effect of truncating the fractal structure at a certain scale related to the energy scale in the spectral action.

Algorithm for Random Close Packing of Spheres with Periodic Boundary Conditions

Based on a theory of porous solids previously developed by the author, the elasticity of a hexagonal close packing of equal spheres is treated. The packing is anisotropic and because of the weight of the spheres, also inhomogeneous. The velocities of propagation of elastic waves have been calculated for evacuated interspaces and for interspaces filled with a liquid or gas. In the case of evacuated or air-filled interspaces, the wave rays and travel times have been computed. The packing which has been treated may be of use as a model for a dry or wet loose material such as gravel or sand. Though the model is very simplified, the results obtained show some typical effects such as anisotropy, inhomogeneity, and a 90 degrees angle of emergence.

The permeability of packs of spheres is important in a wide range of physical scenarios. Here, we create numerically generated random periodic domains of spheres that are polydisperse in size and use lattice-Boltzmann simulations of fluid flow to determine the permeability of the pore phase interstitial to the spheres. We control the polydispersivity of the sphere size distribution and the porosity across the full range from high porosity to a close packing of spheres. We find that all results scale with a Stokes permeability adapted for polydisperse sphere sizes. We show that our determination of the permeability of random distributions of spheres is well approximated by models for cubic arrays of spheres at porosities greater than ∼0.38, without any fitting parameters. Below this value, the Kozeny-Carman relationship provides a good approximation for dense, closely packed sphere packs across all polydispersivity.

An effective optics-electrochemistry approach to random packing density of non-equiaxed ellipsoids

Constant gradient PFG sequence and automated cumulant analysis for quantifying dispersion in flow through porous media

Packing of Spheres: Packing of Equal Spheres

On hydraulic permeability of random packs of monodisperse spheres: Direct flow simulations versus correlations

Packings of granular materials are complex systems consisting of large sets of particles interacting via contact forces. Their internal structure is interesting for several theoretical and practical reasons, especially when the model system consists in a large amount (up to 105) of identical spheres. We herein present a method to process three-dimensional water density maps recorded in wet granular packings of mm-size spheres by magnetic resonance imaging (MRI). Packings of spheres with highly mono-dispersed diameter are considered and the implementation of an ad hoc reconstruction algorithm tailored for this feature allows for the determination of the position of each single sphere with an unprecedented precision (with respect to the scale of the system) while ensuring that all spheres are identified and no non-existing sphere is introduced in the reconstructed packing. The reconstruction of a 0.5 L sample containing about 2 × 104 spheres is presented to demonstrate the robustness of the method.

We studied the dynamic response of a two-dimensional square packing of uncompressed stainless steel spheres excited by impulsive loadings. We developed a new experimental measurement technique, employing miniature tri-axial accelerometers, to determine the stress wave properties in the array resulting from both an in-plane and out-of-plane impact. Results from our numerical simulations, based on a discrete particle model, were in good agreement with the experimental results. We observed that the impulsive excitations were resolved into solitary waves traveling only through initially excited chains. The observed solitary waves were determined to have similar (Hertzian) properties to the extensively studied solitary waves supported by an uncompressed, uniform, one-dimensional chain of spheres. The highly directional response of this system could be used as a basis to design granular crystals with predetermined wave propagation paths capable of mitigating stress wave energy.

Simulation of Drainage and Imbibition in a Random Packing of Equal Spheres

Comparison of NMR simulations of porous media derived from analytical and voxelized representations