pith. sign in

arxiv: 2604.15517 · v1 · submitted 2026-04-16 · ❄️ cond-mat.soft · physics.flu-dyn

Formation of cylindrical shells via sphere packing from fluidized beds

Vin\'icius Pereira da Silva Oliveira , Danilo da Silva Borges , Erick de Moraes Franklin , Jorge Manuel Peixinho This is my paper

Pith reviewed 2026-05-10 09:23 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.flu-dyn
keywords fluidized bedssphere packingcylindrical shellshexagonal latticepolydispersity effectsparticle frictionwall settlingcontact forces
0
0 comments X

The pith

In narrow pipes, fluidized spheres spontaneously settle into stable hexagonal cylindrical shells along the walls.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates how particles in a fluidized bed inside a narrow vertical pipe, with pipe-to-particle diameter ratios of 4.3 and 4.7, transition from a mixed state to one where they settle preferentially on the wall. This settling reduces bed height and fluctuations, producing either disordered glass-like or ordered crystal-like shells that match a hexagonal lattice when unwrapped. Shell stability depends on low polydispersity and sufficient particle friction, with contact force calculations showing particle-particle interactions as the main support. A sympathetic reader would care because the finding identifies a spontaneous packing route in confined flows that could influence transport, filtration, or deposition processes.

Core claim

Starting from a steady fluidized state, particles settle on the pipe wall to form a cylindrical shell. When unwrapped to a plane, the shell shows a hexagonal lattice whose defect density rises with increasing polydispersity. Formation is hindered by polydispersity, with a critical threshold above which crystal-like order becomes unstable. In bidisperse cases the ordered shell requires high enough particle-particle friction. Contact forces within chains and at the wall confirm that particle-particle forces dominate over particle-wall forces in sustaining the shell.

What carries the argument

Spontaneous wall settling of spheres into a hexagonal lattice shell sustained by interparticle contact chains in a narrow-pipe fluidized bed.

If this is right

  • Shells form only within limited ranges of flow velocity and bed height.
  • Hexagonal order is lost above a critical polydispersity level.
  • Bidisperse mixtures form ordered shells only when friction is sufficiently high.
  • Particle-particle contact forces are the primary mechanism holding the shell in place.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same settling process might be tuned by changing pipe diameter or flow rate to control deposit thickness in industrial pipelines.
  • Polydispersity thresholds could serve as a design parameter for avoiding unwanted wall buildup in granular transport systems.
  • Unwrapping analysis of the shell suggests the packing is essentially a 2D hexagonal lattice rolled into a cylinder, which may generalize to other curved confinements.

Load-bearing premise

The numerical model with the chosen diameter ratios and interaction rules faithfully reproduces the real settling dynamics and shell stability seen in experiments.

What would settle it

Perform experiments at the same D/d ratios, flow velocities, and polydispersity values and observe either no wall accumulation or a lattice defect density that does not rise with polydispersity.

read the original abstract

The results of a numerical investigation of fluidized beds of spherical particles in a narrow vertical cylindrical pipe, with particular attention to the spontaneous settling along the wall, are reported. Starting from a steady fluidized state, the particles fluctuate because of fluid-particle, particle-particle, and particle-wall interactions. The particles are heavier than the fluid, with diameters d yielding ratios of pipe to particle diameters D/d=4.3 and 4.7. For given ranges of flow velocities and bed sizes, particles settle on the wall, with a decrease in the bed height and particle fluctuations. Either a glass- or crystal-like shell forms along the pipe wall, in qualitative agreement with previous experiments. The polydispersity and the particle-particle friction are varied to test the stability of the particulate shell formation. The shell structure is analyzed by unwrapping it in a plane and locating all particles and their contact points, and we find that it exhibits a hexagonal lattice with a defects density that increases with polydispersity. The shell formation is hindered by polydispersity, and there exists a critical point for polydispersity above which a crystal-like shell is unstable. In a particular case of bidisperse beds, the crystal-like shell only appears when the particle-particle friction is high enough. Finally, we compute the contact forces within particle-particle chains and in particle-wall contacts, which sustain the cylindrical shell, highlighting the dominant role of particle-particle forces.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports a numerical investigation of fluidized beds of spherical particles in narrow vertical cylindrical pipes (D/d=4.3 and 4.7). Starting from a steady fluidized state, particles settle along the pipe wall with reduced bed height and fluctuations, forming either glass-like or crystal-like shells. The shell is analyzed via unwrapping to a plane, revealing a hexagonal lattice whose defect density increases with polydispersity. Shell formation is hindered by polydispersity, with a critical polydispersity threshold above which crystal-like shells are unstable; in bidisperse cases, sufficient particle-particle friction is required for crystal-like order. Contact-force computations show particle-particle forces dominate over particle-wall contacts in sustaining the shells. Results are presented in qualitative agreement with prior experiments.

Significance. If the numerical model faithfully reproduces the relevant force balances, the work provides mechanistic insight into spontaneous wall ordering in confined granular-fluid systems through direct simulation with independently varied parameters (flow velocity, bed size, polydispersity, friction). The unwrapping analysis yielding explicit hexagonal-lattice identification and defect-density trends, together with the force decomposition, supplies concrete, falsifiable observations that could guide targeted experiments or industrial control of particle deposition in pipes. The parameter sweeps and emphasis on particle-particle chains constitute a strength for an exploratory numerical study.

major comments (2)
  1. [Abstract / Results (unwrapping analysis)] Abstract and results description of shell stability: the claim that 'there exists a critical point for polydispersity above which a crystal-like shell is unstable' is load-bearing for the central conclusion yet is supported only by qualitative trends; no specific polydispersity values, defect-density plots, or statistical measures (e.g., order-parameter thresholds or error bars across realizations) are referenced to define or locate the transition.
  2. [Methods / Results (force analysis)] Methods and validation: the central claims rest on the numerical scheme (CFD-DEM or equivalent) accurately capturing fluid drag, lubrication, and contact forces at D/d=4.3–4.7, but only 'qualitative agreement with previous experiments' is stated. Absence of quantitative benchmarks—bed-height reduction, fluctuation amplitudes, contact-force distributions, or defect densities versus measured values—undermines that the observed hexagonal shells and critical polydispersity are physical rather than artifacts of the drag law or friction calibration.
minor comments (2)
  1. [Abstract] The phrase 'glass- or crystal-like' is used repeatedly without an accompanying order parameter or quantitative metric (e.g., bond-orientational order or coordination-number distribution) to distinguish the two regimes.
  2. [Introduction / Methods] Notation for the pipe-to-particle ratio is given as D/d=4.3 and 4.7; a brief statement of how these specific values were chosen relative to experimental setups would improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight opportunities to strengthen the presentation of our results on polydispersity effects and the validation of the numerical approach. We address each major comment below and outline the revisions we will make.

read point-by-point responses

  1. Referee: [Abstract / Results (unwrapping analysis)] Abstract and results description of shell stability: the claim that 'there exists a critical point for polydispersity above which a crystal-like shell is unstable' is load-bearing for the central conclusion yet is supported only by qualitative trends; no specific polydispersity values, defect-density plots, or statistical measures (e.g., order-parameter thresholds or error bars across realizations) are referenced to define or locate the transition.

    Authors: We agree that the critical polydispersity threshold is central and would benefit from more explicit quantification. The manuscript already varies polydispersity and reports increasing defect density in the unwrapped hexagonal lattice, but the transition is described qualitatively. In the revision we will add the specific polydispersity values examined, reference the defect-density trends with data from multiple independent realizations, and introduce a simple order-parameter threshold (e.g., based on the fraction of six-fold coordinated particles) to locate the point above which long-range hexagonal order is lost. These additions will be placed in both the abstract and the results section. revision: yes

  2. Referee: [Methods / Results (force analysis)] Methods and validation: the central claims rest on the numerical scheme (CFD-DEM or equivalent) accurately capturing fluid drag, lubrication, and contact forces at D/d=4.3–4.7, but only 'qualitative agreement with previous experiments' is stated. Absence of quantitative benchmarks—bed-height reduction, fluctuation amplitudes, contact-force distributions, or defect densities versus measured values—undermines that the observed hexagonal shells and critical polydispersity are physical rather than artifacts of the drag law or friction calibration.

    Authors: We acknowledge that stronger validation would increase confidence. The study employs a standard CFD-DEM formulation with well-documented drag and contact models; the force decomposition (particle-particle versus particle-wall) is an internal diagnostic that directly supports the claim that particle-particle chains dominate shell stability. Quantitative experimental data for defect densities and contact-force statistics at these exact narrow-pipe ratios are not available in the cited literature, limiting direct benchmarking. In the revised manuscript we will expand the methods section with additional model-validation tests (e.g., bed-expansion curves against literature correlations for similar D/d) and a brief sensitivity analysis to drag-law parameters. We will also note the exploratory character of the work while emphasizing that the observed trends with polydispersity and friction are robust within the simulation framework. revision: partial

Circularity Check

0 steps flagged

No circularity detected in simulation-derived results

full rationale

The manuscript reports outcomes from direct numerical simulations (fluidized-bed CFD-DEM type) in which pipe-to-particle diameter ratio, flow velocity, bed height, polydispersity, and friction coefficient are set as independent inputs. Shell formation, hexagonal ordering, defect density, and force-chain analysis are observed as emergent consequences of the integrated particle-fluid and particle-particle dynamics; no equation, fit, or self-citation is invoked to derive these structures from themselves. The central claims therefore remain non-circular.

Axiom & Free-Parameter Ledger

4 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the (unspecified) numerical scheme for modeling fluidized-bed interactions; free parameters include chosen D/d ratios and flow velocities tuned to produce steady states and shell formation.

free parameters (4)
  • pipe-to-particle diameter ratio D/d
    Set to 4.3 and 4.7 to enable narrow-pipe shell formation.
  • flow velocities
    Varied to maintain steady fluidized state before settling.
  • polydispersity level
    Varied to identify stability threshold.
  • particle-particle friction coefficient
    Varied in bidisperse cases to test shell formation.
axioms (1)
  • domain assumption The numerical method accurately reproduces fluid-particle, particle-particle, and particle-wall interactions for the studied conditions.
    Invoked to justify the observed settling and force balance.

pith-pipeline@v0.9.0 · 5569 in / 1323 out tokens · 50811 ms · 2026-05-10T09:23:29.897069+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages · 1 internal anchor

  1. [1]

    Cambridge University Press, Cambridge, UK (2013)

    Andreotti, B., Forterre, Y., Pouliquen, O.: Granular Flow: Between Fluid and Solid. Cambridge University Press, Cambridge, UK (2013)

    work page 2013

  2. [2]

    C´ u˜ nez, F.D., Franklin, E.M.: Crystallization and jamming in narrow fluidized beds. Phys. Fluids32(8), 083303 (2020) https://doi.org/ 10.1063/5.0015410

    work page doi:10.1063/5.0015410 2020

  3. [3]

    Royall, C.P., Charbonneau, P., Dijkstra, M., Russo, J., Smallenburg, F., Speck, T., Valeriani, C.: Colloidal hard spheres: Tri- umphs, challenges, and mysteries. Rev. Mod. Phys.96(4), 045003 (2024) https://doi.org/ 10.1103/RevModPhys.96.045003

    work page doi:10.1103/revmodphys.96.045003 2024

  4. [4]

    Pickett, G.T., Gross, M., Okuyama, H.: Spontaneous chirality in simple systems. Phys. Rev. Lett.85(17), 3652 (2000) https: //doi.org/10.1103/PhysRevLett.85.3652

    work page doi:10.1103/physrevlett.85.3652 2000

  5. [5]

    Mughal, A., Chan, H.K., Weaire, D.: Phyl- lotactic description of hard sphere packing in cylindrical channels. Phys. Rev. Lett. 106(11), 115704 (2011) https://doi.org/10. 1103/PhysRevLett.106.115704

    work page 2011

  6. [6]

    Mughal, A., Chan, H.K., Weaire, D., Hutzler, S.: Dense packings of spheres in cylinders: Simulations. Phys. Rev. E (5), 051305 (2012) https://doi.org/10.1103/ PhysRevE.85.051305

    work page 2012

  7. [7]

    Soft Matter12(9), 2505–2514 (2016) https://doi.org/10.1039/ C5SM02875B

    Fu, L., Steinhardt, W., Zhao, H., Socolar, J.E.S., Charbonneau, P.: Hard sphere pack- ings within cylinders. Soft Matter12(9), 2505–2514 (2016) https://doi.org/10.1039/ C5SM02875B

    work page 2016

  8. [8]

    Lee, T., Gizynski, K., Grzybowski, B.A.: Non-equilibrium self-assembly of monocomponent and multicomponent tubular structures in rotating fluids. Adv. Mater.29(47), 1704274 (2017) https://doi.org/10.1002/adma.201704274

    work page doi:10.1002/adma.201704274 2017

  9. [9]

    Soft Matter21(11), 2049–2058 (2025) https://doi.org/10.1039/ D4SM01359J

    Kurban, E., Vescovi, D., Berzi, D.: Crys- tallization in load-controlled shearing flows of monosized spheres. Soft Matter21(11), 2049–2058 (2025) https://doi.org/10.1039/ D4SM01359J

    work page 2049

  10. [10]

    Soft Matter13(18), 3296–3306 (2017) https://doi.org/10.1039/ C7SM00316A

    Fu, L., Bian, C., Shields, C.W., Cruz, D.F., L´ opez, G.P., Charbonneau, P.: Assembly of hard spheres in a cylinder: a computational and experimental study. Soft Matter13(18), 3296–3306 (2017) https://doi.org/10.1039/ C7SM00316A

    work page 2017

  11. [11]

    Langmuir20(5), 2033–2035 (2004) https://doi.org/10.1021/la0358015

    Moon, J.H., Kim, S., Yi, G.-R., Lee, Y.-H., Yang, S.-M.: Fabrication of ordered macro- porous cylinders by colloidal templating in microcapillaries. Langmuir20(5), 2033–2035 (2004) https://doi.org/10.1021/la0358015

    work page doi:10.1021/la0358015 2033

  12. [12]

    Oliveira, V.P.S., Borges, D.S., Franklin, E.M.: Crystallization and refluidization in very-narrow fluidized beds. Phys. Fluids 35(9), 093306 (2023) https://doi.org/10. 1063/5.0163555

    work page 2023

  13. [13]

    Beverloo, W.A., Leniger, H.A., Velde, J.: The flow of granular solids through orifices. Chem. Eng. Sci.15(3-4), 260–269 (1961) https:// doi.org/10.1016/0009-2509(61)85030-6

    work page doi:10.1016/0009-2509(61)85030-6 1961

  14. [14]

    Zuriguel, I., Garcimart´ ın, A., Maza, D., Pug- naloni, L.A., Pastor, J.: Jamming during the discharge of granular matter from a silo. Phys. Rev. E71(5), 051303 (2005) https: //doi.org/10.1103/PhysRevE.71.051303

    work page doi:10.1103/physreve.71.051303 2005

  15. [15]

    arXiv preprint 10 arXiv:2505.16469 (2025) https://doi.org/10

    Zhou, Y., Li, M., Wang, Y., Guan, Y., Liu, Y., Zou, Z.: Clogging-unclogging tran- sition in 2d vertical pipe. arXiv preprint 10 arXiv:2505.16469 (2025) https://doi.org/10. 48550/arXiv.2505.16469

    work page arXiv 2025

  16. [16]

    Pournin, L., Weber, M., Tsukahara, M., Fer- rez, J.-A., Ramaioli, M., Liebling, T.M.: Three-dimensional distinct element simula- tion of spherocylinder crystallization. Granul. Matter7, 119–126 (2005) https://doi.org/10. 1007/s10035-004-0188-4

    work page 2005

  17. [17]

    Pournin, L., Ramaioli, M., Folly, P., Liebling, T.M.: About the influence of friction and polydispersity on the jamming behav- ior of bead assemblies. Eur. Phys. J. E 23, 229–235 (2007) https://doi.org/10.1140/ epje/i2007-10176-5

    work page 2007

  18. [18]

    Zhang, S., Lin, P., Wang, M., Wan, J.- f., Peng, Y., Yang, L., Hou, M.: Flow- induced surface crystallization of granular particles in cylindrical confinement. Sci. Rep. 11(1), 13227 (2021) https://doi.org/10.1038/ s41598-021-92136-9

    work page 2021

  19. [19]

    Kloss, C., Goniva, C., Hager, A., Amberger, S., Pirker, S.: Models, algorithms and vali- dation for opensource DEM and CFD–DEM. Prog. Comput. Fluid Dyn.12(2-3), 140–152 (2012) https://doi.org/10.1504/PCFD.2012. 047457

    work page doi:10.1504/pcfd.2012 2012

  20. [20]

    Par- ticuology10(5), 582–591 (2012) https://doi

    Goniva, C., Kloss, C., Deen, N.G., Kuipers, J.A., Pirker, S.: Influence of rolling friction on single spout fluidized bed simulation. Par- ticuology10(5), 582–591 (2012) https://doi. org/10.1016/j.partic.2012.05.002

    work page doi:10.1016/j.partic.2012.05.002 2012

  21. [21]

    Powder Technology71(3), 239–250 (1992) https:// doi.org/10.1016/0032-5910(92)88030-L

    Tsuji, Y., Tanaka, T., Ishida, T.: Lagrangian numerical simulation of plug flow of cohe- sionless particles in a horizontal pipe. Powder Technology71(3), 239–250 (1992) https:// doi.org/10.1016/0032-5910(92)88030-L

    work page doi:10.1016/0032-5910(92)88030-l 1992

  22. [22]

    Discrete particle simulation of two-dimensional fluidized bed,

    Tsuji, Y., Kawaguchi, T., Tanaka, T.: Discrete particle simulation of two- dimensional fluidized bed. Powder Technology77(1), 79–87 (1993) https: //doi.org/10.1016/0032-5910(93)85010-7

    work page doi:10.1016/0032-5910(93)85010-7 1993

  23. [23]

    Liu, D., Liu, X., Fu, X., G., W.: Quantifi- cation of the bed load effects on turbulent open-channel flows. J. Geophys. Res. Earth Surf.121, 767–789 (2016) https://doi.org/ 10.1002/2015JF003723

    work page doi:10.1002/2015jf003723 2016

  24. [24]

    Masset, R

    C´ u˜ nez, F.D., Franklin, E.: Mimicking layer inversion in solid-liquid fluidized beds in narrow tubes. Powder Technology364, 994–1008 (2020) https://doi.org/10.1016/j. powtec.2019.09.089

    work page doi:10.1016/j 2020

  25. [25]

    Dickerson

    S. Oliveira, V.P., Borges, D.S., Franklin, E.M., Peixinho, J.: Role of Friction on the Formation of Confined Granular Structures (2025). https://doi.org/10.48550/arXiv. 2505.14382

    work page internal anchor Pith review doi:10.48550/arxiv 2025

  26. [26]

    Zuriguel, I., Gray, J., Peixinho, J., Mullin, T.: Pattern selection by a gran- ular wave in a rotating drum. Phys. Rev. E73(6), 061302 (2006) https: //doi.org/10.1103/PhysRevE.73.061302

    work page doi:10.1103/physreve.73.061302 2006

  27. [27]

    Gray, J.M.N.T.: Particle segregation in dense granular flows. Ann. Rev. Fluid Mech. 50(1), 407–433 (2018) https://doi.org/10. 1146/annurev-fluid-122316-045201

    work page 2018

  28. [28]

    Vo, T.-T., Vu, T.L., Mutabaruka, P.: Effects of size polydispersity on segregation of spher- ical particles in rotating drum. Eur. Phys. J. E44(6), 86 (2021) https://doi.org/10.1140/ epje/s10189-021-00091-0

    work page 2021

  29. [29]

    Academic Press, San Diego, California (1994)

    Gidaspow, D.: Multiphase Flow and Fluidiza- tion: Continuum and Kinetic Theory Descrip- tions, 1st edn. Academic Press, San Diego, California (1994)

    work page 1994

  30. [30]

    Lomholt, S., Maxey, M.R.: Force-coupling method for particulate two-phase flow: Stokes flow. J. Comput. Phys.184(2), 381–405 (2003) https://doi.org/10.1016/ S0021-9991(02)00021-9

    work page 2003

  31. [31]

    Abbas, M., Climent, E., Simonin, O.: Fully coupled simulations of non-colloidal monodis- perse sheared suspensions. Chem. Eng. Res. Des.85(6), 778–791 (2007) https://doi.org/ 10.1205/cherd06114 11

    work page doi:10.1205/cherd06114 2007