Maximally dense random packings of spherocylinders

Evident structural anisotropies arising from near-zero particle asphericity in granular spherocylinder packings

With magnetic resonance imaging experiments, we study packings of granular spherocylinders with merely 2% asphericity. Evident structural anisotropies across all length scales are identified. Most interestingly, the global nematic order decreases with increasing packing fraction, while the local contact anisotropy shows an opposing trend. We attribute this counterintuitive phenomenon to a competition between gravity-driven ordering aided by frictional contacts and a geometric frustration effect at the marginally jammed state. It is also surprising to notice that such slight particle asphericity can trigger non-negligible correlations between contact-level and mesoscale structures, manifested in drastically different nonaffine structural rearrangements upon compaction from that of granular spheres. These observations can help improve statistical mechanical models for the orientational order transformation of nonspherical granular particle packings, which involves complex interplays between particle shape, frictional contacts, and external force field.

Dense Disordered Jammed Packings of Hard Spherocylinders with a Low Aspect Ratio: A Characterization of Their Structure

This work numerically investigates dense disordered (maximally random) jammed packings of hard spherocylinders of cylinder length L and diameter D by focusing on L/D ∈ [0,2]. It is within this interval that one expects that the packing fraction of these dense disordered jammed packings ϕMRJ hsc attains a maximum. This work confirms the form of the graph ϕMRJ hsc versus L/D: here, comparably to certain previous investigations, it is found that the maximal ϕMRJ hsc = 0.721 ± 0.001 occurs at L/D = 0.45 ± 0.05. Furthermore, this work meticulously characterizes the structure of these dense disordered jammed packings via the special pair-correlation function of the interparticle distance scaled by the contact distance and the ensuing analysis of the statistics of the hard spherocylinders in contact: here, distinctly from all previous investigations, it is found that the dense disordered jammed packings of hard spherocylinders with 0.45 ≲ L/D ≤ 2 are isostatic.

Optimal shapes of disk assembly in saturated random packings

Particle morphology is one of the most significant factors influencing the packing structures of granular materials. With certain targeted properties or optimization criteria, inverse packing problems have drawn extensive attention in terms of their adaptability to many material design tasks. An important question hard to answer is which particle shape, especially within given shape families, forms the densest (loosest) random packing? In this paper, we address this issue for the disk assembly model in two dimensions with an infinite variety of shapes, which are simulated in the random sequential adsorption process to suppress crystallization. Via a unique shape representation method, particle shapes are transformed into genotype sequences in the continuous shape space where we utilize the genetic algorithm as an efficient shape optimizer. Specifically, we consider three representative species of disk assembly, i.e., congruent tangent disks, incongruent tangent disks, and congruent overlapping disks, and carry out shape optimization on their packing densities in the saturated random state. We numerically search optimal shapes in the three species with a variable number of constituent disks which yield the maximal and minimal packing densities. We obtain an isosceles circulo-triangle and an unclosed ring for the maximal and minimal packing density in saturated random packings, respectively. The perfect sno-cone and isosceles circulo-triangle are also specifically investigated which give remarkably high packing densities of around 0.6, much denser than those of ellipses. This study is beneficial for guiding the design of particle shapes as well as the inverse design of granular materials.

Effect of sample height and particle elongation in the sedimentation of colloidal rods

We study theoretically the effect of a gravitational field on the equilibrium behaviour of a colloidal suspension of rods with different length-to-width aspect ratios. The bulk phases of the system are described with analytical equations of state. The gravitational field is then incorporated via sedimentation path theory, which assumes a local equilibrium condition at each altitude of the sample. The bulk phenomenology is significantly enriched by the presence of the gravitational field. In a suspension of elongated rods with five stable phases in bulk, the gravitational field stabilizes up to fifteen different stacking sequences. The sample height has a non-trivial effect on the stable stacking sequence. New layers of distinct bulk phases appear either at the top, at the bottom, or simultaneously at the top and the bottom when increasing the sample height at constant colloidal concentration. We also study sedimentation in a mass-polydisperse suspension in which all rods have the same shape but different buoyant masses.

Designing architectural materials: from granular form to functional granular material

Designed granular materials are a novel class of architectural material system. Following one of the key paradigms of designed matter, material form and material function are closely interrelated in these systems. In this context, the article aims to contribute a parametric particle design model as an interface for this interrelation. A granular material is understood as an aggregation of large numbers of individual particles between which only short-range repulsive contact forces are acting. Granular materials are highly pertinent material systems for architecture. Due to the fact that they can act both as a solid and a liquid, they can be recycled and reconfigured multiple times and are thus highly sustainable. Designed granular materials have the added potential that the function of the granular material can be calibrated through the definition of the particles' form. Research on the design of granular materials in architecture is nascent. In physics they have been explored mainly with respect to different particle shapes. However, no coherent parametric particle design model of designed particle shapes for granular material systems in architecture has yet been established which considers both fabrication constraints and simulation requirements. The parametric particle design model proposed in this article has been based on a design system which has been developed through feasibility tests and simulations conducted in research and teaching. Based on this design system the parametric particle design model is developed integrating both fabrication constraints for architecture-scale particle systems and the geometric requirements of established simulation methods for granular materials. Initially the design system and related feasibility tests are presented. The parametric particle design model resulting from that is then described in detail. Directions of further research are discussed especially with respect to the integration of the parametric particle design model in 'inverse' design methods.

Structural analysis of disordered dimer packings

Jammed disordered packings of non-spherical particles show significant variation in the packing density as a function of particle shape for a given packing protocol. Rotationally symmetric elongated shapes such as ellipsoids, spherocylinders, and dimers, e.g., pack significantly denser than spheres over a narrow range of aspect ratios, exhibiting a characteristic peak at aspect ratios of α_max ≈ 1.4-1.5. However, the structural features that underlie this non-monotonic behaviour in the packing density are unknown. Here, we study disordered packings of frictionless dimers in three dimensions generated by a gravitational pouring protocol in LAMMPS. Focusing on the characteristics of contacts as well as orientational and translational order metrics, we identify a number of structural features that accompany the formation of maximally dense packings as the dimer aspect ratio α is varied from the spherical limit. Our results highlight that dimer packings undergo significant structural changes as α increases up to α_max manifest in the reorganisation of the contact configurations between neighbouring dimers, increasing nematic order, and decreasing local translational order. Remarkably, for α > α_max our metrics remain largely unchanged, indicating that the peak in the packing density is related to the interplay of structural rearrangements for α < α_max and subsequent excluded volume effects with unchanged structure for α > α_max.

DEM study of the angle of repose and porosity distribution of cylindrical particles with different lengths

Effect of the length of cylindrical particle on repose angle and porosity of a pile was numerically studied using discrete element method. The variation of repose angle and porosity with coefficient of sliding and rolling friction were also discussed. The results shown that compared with sphere particle, the bottom size of cylindrical pile is smaller, while the height of cylinder pile is larger and the heap is steeper. With the increase of the length of cylinder, the contour line of the pile becomes steep, and the angle of repose increases. The repose angle shows a positive correlation with coefficient of sliding and rolling friction. The porosity increases with the increase of the length of cylinders. The trends of porosity are basically consist with that of repose angle, and with increase of friction coefficient, the average porosity increases.

Algorithms to generate saturated random sequential adsorption packings built of rounded polygons

We present the algorithm for generating strictly saturated random sequential adsorption packings built of rounded polygons. It can be used in studying various properties of such packings built of a wide variety of different shapes, and in modeling monolayers obtained during irreversible adsorption processes of complex molecules. Here, we apply the algorithm to study the densities of packings built of rounded regular polygons. Contrary to packings built of regular polygons, where the packing fraction grows with an increasing number of polygon sides, here the packing fraction reaches its maximum for packings built of rounded regular triangles. With a growing number of polygon sides and increasing rounding radius, the packing fractions tend to the limit given by a packing built of disks. However, they are still slightly higher, even for the rounded 25-gon, which is the highest-sided regular polygon studied here.

Start-up shear of spherocylinder packings: Effect of friction

We study the response to shear deformations of packings of long spherocylindrical particles that interact via frictional forces with friction coefficient μ. The packings are produced and deformed with the help of molecular dynamics simulations combined with minimization techniques performed on a GPU. We calculate the linear shear modulus g_{∞}, which is orders of magnitude larger than the modulus g_{0} in the corresponding frictionless system. The motion of the particles responsible for these large frictional forces is governed by and increases with the length ℓ of the spherocylinders. One consequence of this motion is that the shear modulus g_{∞} approaches a finite value in the limit ℓ→∞, even though the density of the packings vanishes, ρ∝ℓ^{-2}. By way of contrast, the frictionless modulus decreases to zero, g_{0}∼ℓ^{-2}, in accordance with the behavior of density. Increasing the strain beyond a value γ_{c}∼μ, the packing strain weakens from the large frictional to the smaller frictionless modulus when contacts saturate at the Coulomb inequality and start to slide. In this regime, sliding friction contributes a "yield stress" σ_{y}=g_{∞}γ_{c} and the stress behaves as σ=σ_{y}+g_{0}γ. The interplay between static and sliding friction gives rise to hysteresis in oscillatory shear simulations.

Carbon nanothreads enable remarkable enhancement in the thermal conductivity of polyethylene

Polymer nanocomposites with high thermal conductivity have been increasingly sought after in the electronic industry. Based on molecular dynamics simulations, this work assesses the thermal transport in polyethylene (PE) nanocomposites with the presence of a new one-dimensional nanofiller-a carbon nanothread (NTH). It is found that the axial thermal conductivity of PE nanocomposites increases linearly with the content of regularly aligned NTH fillers, while the aggregated pattern suppresses the enhancement effect. This phenomenon is explained by a stronger filler-filler interaction that reduces the intrinsic thermal conductivity of the NTH. Results show that the randomly dispersed NTHs can hardly promote heat transfer because effective heat transfer channels are lacking. Strikingly, surface functionalization has an adverse effect on the thermal conductivity due to the presence of additional voids. The presence of voids answers a long-standing open question that functionalization of the heat conductive filler only slightly improves the thermal conductivity of the polymer composite. Additionally, the transverse thermal conductivity degrades in the presence of the NTH and exhibits no clear correlation with the filler content or the distribution pattern. Overall, this study provides an in-depth understanding of the heat transfer within the polymer nanocomposites, which opens up possibilities for the preparation of highly conductive polymers.

Impact of the Drying Temperature and Grinding Technique on Biomass Grindability

The process of biomass compaction depends on many factors, related to material and process. One of the most important is the proper fragmentation of the raw material. In most cases, more fragmented raw material makes it easier to achieve the desired quality parameters of pellets or briquettes. While the chipping of biomass prefers moist materials, for grinding, the material needs to be dried. As drying temperature changes the properties of the material, these may affect the grinding process. The aim of this work was to determine the influence of the drying temperature of biomass raw material in the range of 60–140 °C on the biomass grindability. To only determine this effect, without the influence of moisture, grinding was carried out on the material in a dry state. The research was carried out on a mill with a knife and hammer grinding system, which is the most popular in the fragmentation of biomass. The analysis of particle size distribution and bulk density of the obtained material was carried out. The energy demand for the grinding process was determined and it was shown that drying temperature, grinding system, and mainly type of biomass affects the grindability.

Hard convex lens-shaped particles: Characterization of dense disordered packings

Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is "optimal" in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction ϕ_{MRJ}≃0.73 [G. Cinacchi and S. Torquato, Soft Matter 14, 8205 (2018)1744-683X10.1039/C8SM01519H]. This value is only a few percent lower than ϕ_{DKP}=0.76210⋯, the packing fraction of the corresponding densest-known crystalline (degenerate) packings [G. Cinacchi and S. Torquato, J. Chem. Phys. 143, 224506 (2015)JCPSA60021-960610.1063/1.4936938]. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method-based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows us to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated here can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of ten contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve.

Predicting maximally random jammed packing density of non-spherical hard particles via analytical continuation of fluid equation of state

We developed a formalism for accurately predicting the density of MRJ packing state of a wide spectrum of congruent non-spherical hard particles in 3D via analytical fluid EOS.

Random sequential adsorption of particles with tetrahedral symmetry

We study random sequential adsorption (RSA) of a class of solids that can be obtained from a cube by specific cutting of its vertices, in order to find out how the transition from tetrahedral to octahedral symmetry affects the densities of the resulting jammed packings. We find that in general solids of octahedral symmetry form less dense packing; however, the lowest density was obtained for the packing built of tetrahedra. The densest packing is formed by a solid close to a tetrahedron but with vertices and edges slightly cut. Its density is θ_{max}=0.41278±0.00059 and is higher than the mean packing fraction of spheres or cuboids but is lower than that for the densest RSA packings built of ellipsoids or spherocylinders. The density autocorrelation function of the studied packings is typical for random media and vanishes very quickly with distance.

Vibrational properties of two-dimensional dimer packings near the jamming transition

Jammed particulate systems composed of various shapes of particles undergo the jamming transition as they are compressed or decompressed. To date, sphere packings have been extensively studied in many previous works, where isostaticity at the transition and scaling laws with the pressure of various quantities, including the contact number and the vibrational density of states, have been established. Additionally, much attention has been paid to nonspherical packings, and particularly recent work has made progress in understanding ellipsoidal packings. In this work, we study the dimer packings in two dimensions, which have been much less understood than systems of spheres and ellipsoids. We first study the contact number of dimers near the jamming transition. It turns out that packings of dimers have "rotational rattlers," each of which still has a free rotational motion. After correcting this effect, we show that dimers become isostatic at the jamming, and the excess contact number obeys the same critical law and finite-size scaling law as those of spheres. We next study the vibrational properties of dimers near the transition. We find that the vibrational density of states of dimers exhibits two characteristic plateaus that are separated by a peak. The high-frequency plateau is dominated by the translational degree of freedom, while the low-frequency plateau is dominated by the rotational degree of freedom. We establish the critical scaling laws of the characteristic frequencies of the plateaus and the peak near the transition. In addition, we present detailed characterizations of the real space displacement fields of vibrational modes in the translational and rotational plateaus.

Packing Characteristics of Aggregate with Consideration of Particle size and Morphology

The packing characteristics of aggregates are very important for aggregate blend design, which is closely related to the performance of mixtures. In this study, an indoor packing experiment was designed to investigate the behaviors of single-size particle packing and two-size particle packing. The effects of particle composition, particle size and size ratio, particle morphology on packing characteristics were also evaluated. Two kinds of aggregates (crushed stone and gravel) with significant morphological differences were selected for the test. The angularity of the aggregates was quantitatively analyzed using the variance of mean curvature ( S C m 2 ) of particle surface in accordance with the 3-D scanning measurements. Based on the test results, the packing characteristics of aggregates were analyzed using the air void content (Va) and the packing function index (Ipf) proposed in this paper. It is shown that the analysis results of packed ideal spheres cannot be directly used to describe the packing characteristics of aggregates. Particle morphology has a significant impact on packing characteristics, especially on the Va. The Va of packed aggregates with poor angularity is significantly smaller than that with good angularity. Ipf can be used to quantitatively distinguish the packing function of particles. The test results show that the packing function of particles cannot be simply divided into the skeleton building and air voids filling. Generally, the particles in packed blend have both of these functions. The packing function of particles depends not only on the particle size, but also on the composition of particles with different size. When the size ratio and volume ratio are the same, the packing characteristics of the two-size particle blends will still change with the change of the particle size. The exploration of packing behaviors of single- size and two- size particle aggregates is helpful for analyzing the packing behaviors of blends with multi-size particles.

Hard convex lens-shaped particles: metastable, glassy and jammed states

We generate and study dense positionally and/or orientationally disordered, including jammed, monodisperse packings of hard convex lens-shaped particles (lenses). Relatively dense isotropic fluid configurations of lenses of various aspect ratios are slowly compressed via a Monte Carlo method based procedure. Under this compression protocol, while 'flat' lenses form a nematic fluid phase (where particles are positionally disordered but orientationally ordered) and 'globular' lenses form a plastic solid phase (where particles are positionally ordered but orientationally disordered), 'intermediate', neither 'flat' nor 'globular', lenses do not form either mesophase. In general, a crystal solid phase (where particles are both positionally and orientationally ordered) does not spontaneously form during lengthy numerical simulation runs. In correspondence to those volume fractions at which a transition to the crystal solid phase would occur in equilibrium, a 'downturn' is observed in the inverse compressibility factor versus volume fraction curve beyond which this curve behaves essentially linearly. This allows us to estimate the volume fraction at jamming of the dense non-crystalline packings so generated. These packings are nematic for 'flat' lenses and plastic for 'globular' lenses, while they are robustly isotropic for 'intermediate' lenses, as confirmed by the calculation of the τ order metric, among other quantities. The structure factors S(k) of the corresponding jammed states tend to zero as the wavenumber k goes to zero, indicating they are effectively hyperuniform (i.e., their infinite-wavelength density fluctuations are anomalously suppressed). Among all possible lens shapes, 'intermediate' lenses with aspect ratio around 2/3 are special because they are those that reach the highest volume fractions at jamming while being positionally and orientationally disordered and these volume fractions are as high as those reached by nematic jammed states of 'flat' lenses and plastic jammed states of 'globular' lenses. All of their attributes, taken together, make such 'intermediate' lens packings particularly good glass-forming materials.

Mesoscale Anisotropy in Porous Media Made of Clay Minerals. A Numerical Study Constrained by Experimental Data

The anisotropic properties of clay-rich porous media have significant impact on the directional dependence of fluids migration in environmental and engineering sciences. This anisotropy, linked to the preferential orientation of flat anisometric clay minerals particles, is studied here on the basis of the simulation of three-dimensional packings of non-interacting disks, using a sequential deposition algorithm under a gravitational field. Simulations show that the obtained porosities fall onto a single master curve when plotted against the anisotropy value. This finding is consistent with results from sedimentation experiments using polytetrafluoroethylene (PTFE) disks and subsequent extraction of particle anisotropy through X-ray microtomography. Further geometrical analyses of computed porous media highlight that both particle orientation and particle aggregation are responsible of the evolution of porosity as a function of anisotropy. Moreover, morphological analysis of the porous media using chord length measurements shows that the anisotropy of the pore and solid networks can be correlated with particle orientation. These results indicate that computed porous media, mimicking the organization of clay minerals, can be used to shed light on the anisotropic properties of fluid transfer in clay-based materials.

Evolutions of packing properties of perfect cylinders under densification and crystallization

Cylindrical particles are ubiquitous in nature and industry, and a cylinder is a representative shape of rod-like particles. However, the disordered packing results of cylinders in previous studies are quite inconsistent with each other. In this work, we obtain the MRJ (maximally random jammed) packings and the MDRPs (maximally dense random packings) of perfect cylinders with the aspect ratio (height/diameter) 0.2 ≤ w ≤ 6.0 using the ASC (adaptive shrinking cell) algorithm and the IMC (inverse Monte Carlo) method, respectively. The optimal aspect ratio corresponding to the maximal packing density is w = 0.9 in the MRJ state, while the value is w = 1.2 in the MDRP state. Then we investigate the evolutions of packing properties of perfect cylinders under densification and crystallization. We compare the different final packing states generated via the two methods with different compression rates and order constraints. In the densification procedure, we generate jammed and random packings of cylinders with various compression rates via the ASC and IMC method, respectively. When decreasing the compression rate, we find that the packing density increases but the optimal w remains the same in both methods. In the crystallization procedure, the order constraint in the IMC method is gradually released which means the degree of order in the packings is allowed to increase, and we find that the optimal w shifts from 1.2 to 0.9 while the packing density increases as well. Meanwhile, the random packings evolve into the jammed packings in the crystallization procedure which reflects the competition mechanism between the randomness and jamming. These results also indicate that the optimal w is solely related to the degree of order in the cylinder packings but not determined by the protocol or packing density. Furthermore, a uniform shape elongation effect on the random-packing densities of various shaped particles is found via a new proposed definition of the scaled aspect ratio. Finally, a rough linear relationship between the mean and standard deviation of the reduced Voronoi cell volumes is obtained only for the random packings. Our findings should lead to a better understanding toward the jammed and random packings and are helpful in guiding the granular material design.

Hypostatic jammed packings of frictionless nonspherical particles

We perform computational studies of static packings of a variety of nonspherical particles including circulo-lines, circulo-polygons, ellipses, asymmetric dimers, dumbbells, and others to determine which shapes form packings with fewer contacts than degrees of freedom (hypostatic packings) and which have equal numbers of contacts and degrees of freedom (isostatic packings), and to understand why hypostatic packings of nonspherical particles can be mechanically stable despite having fewer contacts than that predicted from naive constraint counting. To generate highly accurate force- and torque-balanced packings of circulo-lines and cir-polygons, we developed an interparticle potential that gives continuous forces and torques as a function of the particle coordinates. We show that the packing fraction and coordination number at jamming onset obey a masterlike form for all of the nonspherical particle packings we studied when plotted versus the particle asphericity A, which is proportional to the ratio of the squared perimeter to the area of the particle. Further, the eigenvalue spectra of the dynamical matrix for packings of different particle shapes collapse when plotted at the same A. For hypostatic packings of nonspherical particles, we verify that the number of "quartic" modes along which the potential energy increases as the fourth power of the perturbation amplitude matches the number of missing contacts relative to the isostatic value. We show that the fourth derivatives of the total potential energy in the directions of the quartic modes remain nonzero as the pressure of the packings is decreased to zero. In addition, we calculate the principal curvatures of the inequality constraints for each contact in circulo-line packings and identify specific types of contacts with inequality constraints that possess convex curvature. These contacts can constrain multiple degrees of freedom and allow hypostatic packings of nonspherical particles to be mechanically stable.

Dense crystalline packings of ellipsoids

An ellipsoid, the simplest nonspherical shape, has been extensively used as a model for elongated building blocks for a wide spectrum of molecular, colloidal, and granular systems. Yet the densest packing of congruent hard ellipsoids, which is intimately related to the high-density phase of many condensed matter systems, is still an open problem. We discover an unusual family of dense crystalline packings of self-dual ellipsoids (ratios of the semiaxes α:sqrt[α]:1), containing 24 particles with a quasi-square-triangular (SQ-TR) tiling arrangement in the fundamental cell. The associated packing density ϕ exceeds that of the densest known SM2 crystal [ A. Donev et al., Phys. Rev. Lett. 92, 255506 (2004)10.1103/PhysRevLett.92.255506] for aspect ratios α in (1.365, 1.5625), attaining a maximal ϕ≈0.75806... at α=93/64. We show that the SQ-TR phase derived from these dense packings is thermodynamically stable at high densities over the aforementioned α range and report a phase diagram for self-dual ellipsoids. The discovery of the SQ-TR crystal suggests organizing principles for nonspherical particles and self-assembly of colloidal systems.

Maximally dense random packings of cubes and cuboids via a novel inverse packing method

The packings of cubes and cuboids (i.e., "elongated" or "compressed" cubes) are ubiquitous in nature. The high symmetry and space-tiling nature of such particles make them easily packable in dense configurations with a high degree of orientational and translational order. In this paper, we devise a novel inverse packing method that enables one to generate dense hard-particle packings with a controllable degree of disorder/order quantified by predefined order metrics via stochastic Monte Carlo optimizations. We employ the inverse packing method to generate and investigate the maximally dense random packings (MDRPs) of hard cubes and cuboids with aspect ratio α, in which a series of newly introduced normalized local cubatic order parameters sensitive to the onset of any spatial order in packings of cubes and cuboids is minimized. The density of the MDRP of cubes is φ ≈ 0.637, which increases as the shape deviates from the cube limit (α = 1) and reaches the maximal values for cuboids with aspect ratios α = 0.7 or 1.5. These special α values associated with local density extrema are almost identical for those associated with the random packings of spherocylinders, spheroids and superellipsoids, suggesting a universal influence of shape elongation on random packing density. Our inverse packing method can be readily utilized to study the MDRPs of other hard particles and the normalized local cubatic order parameter introduced here is applicable to other shaped particles characterized by three principal axes.