pith. sign in

arxiv: 2510.12571 · v2 · submitted 2025-10-14 · ⚛️ physics.flu-dyn

Low Reynolds number flow in a packed bed of rotated bars

Wojciech Sadowski , Christin Velten , Maximilian Br\"ommer , Hakan Demir , Kerstin H\"ulz , Francesca di Mare , Katharina Z\"ahringer , Viktor Scherer This is my paper

Pith reviewed 2026-05-18 07:35 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords packed bedrotated barslow Reynolds numberparticle image velocimetryparticle-resolved simulationvoid spacesfreeboard jets
0
0 comments X

The pith

The flow inside a packed bed of rotated square bars is largely independent of Reynolds number and determined by the geometry of the void spaces.

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

This paper investigates gas flow through an experimental packed bed reactor made of square bars arranged in layers, each rotated by 30 degrees. Measurements using Particle Image Velocimetry at Reynolds numbers of 100 and 200 reveal that the internal flow patterns remain consistent across these values. Two particle-resolved numerical methods confirm the experimental data, showing that the complex void spaces shape the flow more than the flow speed does. This matters for designing reactors where predictable flow inside the packing is needed without strong dependence on operating conditions.

Core claim

The flow inside the bed is largely independent from the Reynolds number and seems to be determined by the geometry of the void spaces. The flow in the freeboard is dominated by the presence of slowly dissipating jets downstream of the bed, which are characterized by unsteady oscillations at the higher Reynolds number. The numerical results obtained with both boundary-conforming meshing and blocked-off methods are in good agreement with the measurements.

What carries the argument

The geometry of the void spaces formed by layers of square bars each rotated by 30 degrees, which dictates the internal flow patterns.

Load-bearing premise

That the numerical artifacts in the two simulation strategies do not systematically affect the agreement with experimental data especially in the freeboard.

What would settle it

If new PIV measurements at Re=50 or Re=300 inside the bed show velocity fields that differ substantially from those at 100 and 200 in a way not attributable to the void geometry alone.

Figures

Figures reproduced from arXiv: 2510.12571 by Christin Velten, Francesca di Mare, Hakan Demir, Katharina Z\"ahringer, Kerstin H\"ulz, Maximilian Br\"ommer, Viktor Scherer, Wojciech Sadowski.

Figure 1
Figure 1. Figure 1: a) Schematic of the assembly process of the packed bed using 30° rotation between modules. The red and blue lines mark the two measurement positions (P1 and P3) in different layers. b) The cut￾out illustrating the structure of void spaces in the bed. c) An optical module. d) The PIV setup including the packed bed and the measuring equipment. 2.1 Experimental setup To allow optical access inside the packed … view at source ↗
Figure 2
Figure 2. Figure 2: Averaged velocity fields at P1 inside layer 17, visualized using streamlines. Color denotes the vertical component 𝑤/ 𝑤 . Plotted data has been gathered using PIV and simulation approaches: bound￾ary-conforming and blocked-off methods, at Re𝑝 = 100 (top) and 200 (bottom) [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Averaged velocity fields at P3 inside layer 17, visualized using streamlines. Color denotes the vertical component 𝑤/ 𝑤 . Plotted data has been gathered using PIV and simulation approaches: bound￾ary-conforming and blocked-off methods, at Re𝑝 = 100 (top) and 200 (bottom). The flow inside the packed bed is visualized in detail in the Figs. 2 and 3, illustrating the velocity field at positions P1 and P3, res… view at source ↗
Figure 4
Figure 4. Figure 4: Averaged vertical (𝑤) and horizontal (𝑢𝜉) velocity components inside the bed for Re𝑝 = 100, plotted along the lines at vertical coordinates 𝑧/𝐵 = 16.25, 16.5, 16.75 (annotated by gray numbers) in planes P1 (left) and P3 (right) in the 17th layer [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Averaged vertical (𝑤) and horizontal (𝑢𝜉) velocity components inside the bed for Re𝑝 = 200, plotted along the lines at vertical coordinates 𝑧/𝐵 = 16.25, 16.5, 16.75 (annotated by gray numbers) in planes P1 (left) and P3 (right) in the 17th layer. The differences in the flow field near the side walls, are further presented quantitatively in Figs. 4 and 5 corresponding to Re𝑝 = 100 and 200, respectively. Bot… view at source ↗
Figure 6
Figure 6. Figure 6: An overview of the flow in the freeboard: velocity fields from the simulations employing boundary-conforming mesh for Re𝑝 = 100 (a) and 200 (b) at the planes at 𝑧/𝐵 = 19, 21 and 23. Colormap denotes the vertical velocity component 𝑤/ 𝑤 and the red lines visualize the streamlines of the flow field in a selection of jets. Green lines indicate sampling lines at the positions P1 and P3. The freeboard flow fiel… view at source ↗
Figure 7
Figure 7. Figure 7: Velocity field in the freeboard for Re𝑝 = 100 and 200 at P1 (left) and P3 (right), visualized using streamlines. Colormap denotes the vertical component 𝑤/ 𝑤 . Plotted data has been gathered using PIV and simulation approaches: boundary-conforming and blocked-off methods. At the sides of the P1 plane at Re𝑝 = 100, slow-moving gas is drawn from the outlet into the freeboard and below the bed surface, locate… view at source ↗
Figure 8
Figure 8. Figure 8: Averaged vertical (𝑤) and horizontal (𝑢𝜉) velocity components in the freeboard for Re𝑝 = 100, plotted along the lines at vertical coordinates 𝑧/𝐵 = 19, 21, 23 (annotated by gray numbers) in planes P1 (left) and P3 (right) [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Averaged vertical (𝑤) and horizontal (𝑢𝜉) velocity components in the freeboard for Re𝑝 = 200, plotted along the lines at vertical coordinates 𝑧/𝐵 = 19, 21, 23 (annotated by gray numbers) in planes P1 (left) and P3 (right). Examining first the results for lower Re𝑝 , near the bed surface (at 𝑧/𝐵 = 19), the boundary-conforming method reproduces the measured vertical velocity 𝑤 and shapes of the jets very acc… view at source ↗
read the original abstract

The present study focuses on the gas flow through an experiment-scale modular packed bed reactor consisting of square bars, arranged in layers. Each layer is rotated by $30^\circ$ resulting in a complex shape of the void spaces between the bars. Particle Image Velocimetry measurement results inside and on top of the studied system are presented for particle-based Reynolds numbers of 100 and 200, and used as validation data for two sets of particle-resolved numerical simulations, using boundary conforming meshing strategy and treating the solid boundaries via the blocked-off method. The flow inside the bed is largely independent from the Reynolds number and seems to be determined by the geometry of the void spaces. The flow in the freeboard is dominated by the presence of slowly dissipating jets downstream of the bed, which are characterized by unsteady oscillations at the higher Reynolds number. The numerical results obtained with both methods are in good agreement with the measurements, both inside and above the bed. However, stronger deviations between the results can be observed in the freeboard and can be traced to numerical properties of the current simulation approaches.

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

1 major / 2 minor

Summary. The manuscript presents PIV measurements of gas flow through a modular packed bed of square bars rotated by 30° per layer at particle-based Reynolds numbers of 100 and 200. These data serve as validation for two particle-resolved simulation approaches (boundary-conforming meshing and blocked-off method). The central claim is that internal-bed flow is largely independent of Reynolds number and controlled by void-space geometry, whereas freeboard flow is dominated by slowly dissipating jets that exhibit unsteady oscillations at the higher Re; both numerical methods agree well with PIV inside the bed but show larger deviations above it.

Significance. If the Re-independence result holds, the work would simplify modeling of low-Re packed-bed reactors by emphasizing geometric control over inertial effects. The combination of direct PIV validation with two independent simulation strategies supplies a useful benchmark for complex void geometries; the explicit acknowledgment of freeboard discrepancies is a strength.

major comments (1)
  1. [Abstract and Results] Abstract and Results (velocity-field comparisons at Re=100 vs. Re=200): the claim that internal flow is 'largely independent from the Reynolds number' rests on visual agreement between only two relatively close values lying in the same inertial regime. No quantitative similarity metric (e.g., normalized L2 difference, correlation coefficient, or integrated kinetic-energy difference between the two Re cases) is reported, so the assertion that geometry of the void spaces is the sole determinant remains qualitative and could be undermined by gradual Re dependence outside this narrow window.
minor comments (2)
  1. [Freeboard flow discussion] The description of jet unsteadiness in the freeboard would be strengthened by reporting a specific frequency or Strouhal number extracted from the time-resolved data or simulations.
  2. [Figures] Figure captions and axis labels should explicitly state the normalization used for velocity (e.g., inlet velocity or superficial velocity) to facilitate direct comparison between PIV and the two simulation sets.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comment point by point below and will incorporate revisions to strengthen the presentation.

read point-by-point responses

  1. Referee: [Abstract and Results] Abstract and Results (velocity-field comparisons at Re=100 vs. Re=200): the claim that internal flow is 'largely independent from the Reynolds number' rests on visual agreement between only two relatively close values lying in the same inertial regime. No quantitative similarity metric (e.g., normalized L2 difference, correlation coefficient, or integrated kinetic-energy difference between the two Re cases) is reported, so the assertion that geometry of the void spaces is the sole determinant remains qualitative and could be undermined by gradual Re dependence outside this narrow window.

    Authors: We agree that the independence statement would benefit from quantitative support rather than relying solely on visual inspection. In the revised manuscript we will add explicit metrics for the internal-bed region, including the normalized L2 difference between the time-averaged velocity fields at Re=100 and Re=200 together with the spatial correlation coefficient. These quantities will be reported in the Results section alongside the existing figures. Within the narrow inertial window examined, the similarity is consistent with geometric control of the void-scale flow, but we accept that the added metrics will make the claim more objective. revision: yes

standing simulated objections not resolved
  • We cannot address possible gradual Reynolds-number dependence at values well outside the Re=100–200 range without performing additional experiments or simulations at other particle-based Reynolds numbers.

Circularity Check

0 steps flagged

No circularity: empirical PIV validation and dual numerical methods provide independent benchmarks for Re-independence claim

full rationale

The paper's central observation—that internal flow is largely independent of Reynolds number and governed by void-space geometry—is drawn directly from comparing PIV measurements at Re=100 and Re=200 against two distinct particle-resolved simulation strategies (boundary-conforming meshing and blocked-off method). These measurements serve as external experimental data rather than fitted inputs, and the simulations are validated against them instead of being tuned to reproduce a target result. No self-definitional steps, fitted parameters renamed as predictions, or load-bearing self-citations appear in the derivation chain; the agreement between methods and data confirms numerical fidelity without reducing the physical claim to an internal tautology. The analysis remains self-contained against the reported benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

This is an experimental and numerical validation study relying on standard low-Re fluid assumptions rather than new fitted parameters or postulated entities.

axioms (1)
  • domain assumption Gas flow treated as incompressible at the studied low Reynolds numbers
    Standard modeling choice for low-Re internal flows implied by the abstract.

pith-pipeline@v0.9.0 · 5748 in / 1184 out tokens · 41926 ms · 2026-05-18T07:35:58.092932+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

12 extracted references · 12 canonical work pages

  1. [1]

    L. R. Rosseau, V. Middelkoop, H. A. Willemsen, I. Roghair, M. van Sint Annaland, Front. chem. eng. 2022, 4, 834547. DOI: 10.3389/fceng.2022.834547

    work page doi:10.3389/fceng.2022.834547 2022

  2. [2]

    Moghaddam, E

    E. Moghaddam, E. Foumeny, A. Stankiewicz, and J. Padding, Chem. Eng. Sci. 2019, 1, 100006. DOI: 10.1016/j.cesx.2019.100006

    work page doi:10.1016/j.cesx.2019.100006 2019

  3. [3]

    Jurtz, M

    N. Jurtz, M. Kraume, and G. D. Wehinger, Rev. Chem. Eng. 2019, 35 (2), 139–190. DOI: 10.1515/revce-2017-0059

    work page doi:10.1515/revce-2017-0059 2019

  4. [4]

    Sadowski, M

    W. Sadowski, M. Sayyari, F. di Mare, C. Velten, K. Zähringer, Phys. Fluids 2024, 36 (2), 023330. DOI: 10.1063/5.0188247

    work page doi:10.1063/5.0188247 2024

  5. [5]

    Patankar, Numerical Heat Transfer and Fluid Flow , 1st ed., CRC Press, Boca Raton 1980

    S. Patankar, Numerical Heat Transfer and Fluid Flow , 1st ed., CRC Press, Boca Raton 1980. Chemical Engineering & Technology 12

    work page 1980

  6. [6]

    Gorges, M

    C. Gorges, M. Brömmer, C. Velten, S. Wirtz, E. I. Mahiques, V. Scherer, K. Zähringer B. van Wachem, Particuology 2024, 86, 1–12. DOI: 10.1016/j.partic.2023.04.006

    work page doi:10.1016/j.partic.2023.04.006 2024

  7. [7]

    Brömmer, M

    M. Brömmer, M. Scharnowski, E. Illana Mahiques, S. Wirtz, V. Scherer, Particuology 2024, 85, pp. 89–101. DOI: 10.1016/j.partic.2023.03.022

    work page doi:10.1016/j.partic.2023.03.022 2024

  8. [8]

    Illana Mahiques, M

    E. Illana Mahiques, M. Brömmer, S. Wirtz, V. Scherer, Chem. Eng. Technol. 2023, 46, 1362-1372. DOI: 10.1002/ceat.202200622

    work page doi:10.1002/ceat.202200622 2023

  9. [9]

    Lovreglio, S

    P. Lovreglio, S. Das, K. A. Buist, E. A. J. F. Peters, L. Pel, and J. A. M. Kuipers, AIChE J. 2018, 64 (5), pp. 1896–1907. DOI: 10.1002/aic.16127

    work page doi:10.1002/aic.16127 2018

  10. [10]

    Khayamyan, T

    S. Khayamyan, T. S. Lundström, P. Gren, H. Lycksam, J. G. I. Hellström, Transp. Porous Media 2017, 117 (1), 45–67. DOI: 10.1007/s11242-017-0819-y

    work page doi:10.1007/s11242-017-0819-y 2017

  11. [11]

    Velten, K

    C. Velten, K. Zähringer, Transp. Porous Media 2023, 150, 307–326. DOI: 10.1007/s11242-023- 02010-7

    work page doi:10.1007/s11242-023- 2023

  12. [12]

    Velten, K

    C. Velten, K. Hülz, K. Zähringer, submitted to Particuology 2025

    work page 2025