pith. sign in

arxiv: 2511.05157 · v2 · pith:PC46WKQZnew · submitted 2025-11-07 · ⚛️ physics.flu-dyn

A multiple-scales framework for branched channel filters

T. Fastnedge , C. J. W. Breward , I. M. Griffiths This is my paper

Pith reviewed 2026-05-18 00:33 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords multiple scalesbranched channel filterseffective boundary conditionricochet separationparticle trajectoriesStokes numberhigh-Reynolds-number flowleakage boundary condition
0
0 comments X

The pith

Multiple-scales analysis produces an effective leakage boundary condition for flow over branched channel filters.

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

The paper develops a multiple-scales method for high-Reynolds-number laminar flow past a surface containing many discrete branched channels. It averages the local flow disturbances caused by the branches to obtain a simplified leakage boundary condition. This condition then supplies closed-form expressions for the velocity and pressure fields at distances away from the surface. The same framework tracks spherical particles under a force-balance model with wall bounces to determine how the fraction entering the branches depends on the Stokes number. The resulting model therefore links device efficiency directly to geometric and operating parameters without repeated full-scale simulations.

Core claim

We use multiple-scales techniques to derive an effective leakage boundary condition, which smooths out localised effects in the flow velocity and pressure that arise due to the discrete branched channels, and then use this boundary condition to explicitly determine the flow away from the boundary. We find that our explicit solution compares well with an analogous numerical solution containing a discrete set of branched channels. We further consider the behaviour of individual spherical particles in the device, with their trajectories determined via a simple force balance model with a wall-bounce condition, and explore the dependence of the fraction of particles that flow into the branched渠道s

What carries the argument

Effective leakage boundary condition obtained by multiple-scales averaging of discrete branched channels.

If this is right

  • Velocity and pressure fields are available in closed form outside a thin layer near the filter surface.
  • The smoothed boundary condition reproduces the far-field behavior of fully resolved discrete-channel computations.
  • Particle entry into branches is governed by Stokes number through the integrated trajectory equations.
  • Overall filter efficiency becomes a direct function of branch spacing, flow speed, and particle properties.

Where Pith is reading between the lines

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

  • The same averaging procedure could be applied to optimize branch angles or spacings for higher ricochet efficiency.
  • Extensions to unsteady or lower-Reynolds-number regimes would require checking whether the scale separation still permits a simple leakage condition.
  • The explicit flow solution offers a fast surrogate model for exploring parameter space before committing to detailed simulations.

Load-bearing premise

The separation between branches is much larger than the thickness of the viscous boundary layer.

What would settle it

A full numerical simulation of the discrete branched geometry that produces flow fields or particle-capture fractions differing substantially from the explicit multiple-scales solution.

Figures

Figures reproduced from arXiv: 2511.05157 by C. J. W. Breward, I. M. Griffiths, T. Fastnedge.

Figure 1
Figure 1. Figure 1: Layout schematic of a branched channel filter preceding a dead-end filter. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: 2-dimensional repeatable T-junction domain, Ωˆ , given by a main channel compartment with 𝑁 perpendicular branched channels on the bottom wall. Inlet and outlets are indicated by dashed black lines, the T-junction spacing is indicated by dashed red lines and boundary walls are denoted by 𝜕Ωˆ 𝑤, in solid black lines. The domain design parameters are indicated as ℎ1, ℎ2, 𝐿, 𝐿1, 𝐿2 and 𝑁. We apply an outlet p… view at source ↗
Figure 3
Figure 3. Figure 3: Reduced dimensionless geometry, with point-sinks replacing each branched [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Outer flow domain with boundary conditions, including the effective boundary [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Conformal map of the semi-infinite half strip inner region to the positive [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Numerical solution for the magnitude of the flow velocity, [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Numerical solutions, zoomed into individual branched channels, for (a) the [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Numerical solution for the pressure, 𝑝, corresponding to figure 6, solved via the Navier–Stokes equations (2.13)–(2.14). Here, Pout = 0.4, Re = 1000, 𝜖 = 0.04, 𝛿 = 0.1, 𝜆 = 0.1 and 𝛾 = 0.5. The colour range for the pressure is restricted to [0.397, 0.401], and white otherwise, to highlight the asymptotic observation of constant pressure to leading-order in 𝜖 over the main channel. The branched channels are… view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of the leading-order asymptotic composite solution ( [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Comparison of the leading-order asymptotic composite solution ( [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Leading-order asymptotic composite solution ( [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of the leading-order asymptotic composite solution ( [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Numerical solutions, zoomed into individual branched channels, for (a) the [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Comparison of the leading-order asymptotic composite solution in red dashed [PITH_FULL_IMAGE:figures/full_fig_p019_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Comparison of the (a) flux through each branched channel, recorded at [PITH_FULL_IMAGE:figures/full_fig_p020_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Angled geometry with (a) constant branch width [PITH_FULL_IMAGE:figures/full_fig_p020_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Comparison of the total flux through the discrete branches, [PITH_FULL_IMAGE:figures/full_fig_p021_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Examples of possible bounces for St = 1, with the flow field given by the outer solution 𝒖 (𝑜) 0 (𝑥, 𝑦) (blue), the composite solution 𝒖𝑐 (𝑥, 𝑦) (red) and the full numerical flow (black) from figure 6. We vary the initial position, 𝑦0 to show the minimal difference between trajectories when St = O (1), for particles initialise both inside and outside of the boundary layer indicated by black dashed lines. … view at source ↗
Figure 19
Figure 19. Figure 19: Examples of possible bounces for St = 0.04, with the flow field given by the outer solution 𝒖 (𝑜) 0 (𝑥, 𝑦) (blue), the composite solution 𝒖𝑐 (𝑥, 𝑦) (red) and the full numerical flow (black) from figure 6. We vary the initial position, 𝑦0 to show the minimal difference between trajectories when St = O (𝜖), for particles initialise both inside and outside of the boundary layer indicated by black dashed line… view at source ↗
Figure 20
Figure 20. Figure 20: Example trajectory in the composite solution for the flow ( [PITH_FULL_IMAGE:figures/full_fig_p025_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Phase diagram for St = ∞, with the limiting trajectory (red dashed) having inlet position 𝒙 𝑝 = (0, 𝑦0) and velocity 𝒖(0, 𝑦0), for whether a particle will bounce. 𝛿𝜖/2 either side of a point-sink, we remove the particle from the simulation and assume that the particle will have gone down a branched channel (see figure 20 for an example particle removal via a branched channel). To calculate the proportion … view at source ↗
Figure 22
Figure 22. Figure 22: The 𝑥-intercept, calculated using (10.4), for a ballistic particle given an inlet position 𝑦0, in the case where St = ∞. We show the position of the branched channels by black horizontal lines, extended along to the 𝑦0 values for which the 𝑥-intercept falls at the centre of the branched channel. Note that 𝑦0 ≠ 0 and so the first 𝑥-intercept is after the first two branched channels. that a large proportion… view at source ↗
Figure 23
Figure 23. Figure 23: The proportion of particles, K, leaving through the main channel compared to the total at the inlet, calculated using (10.1), plotted as black points between (a) St ∈ [0, 0.4] with a spacing of 0.01 and (b) St ∈ [1, 10] with a spacing of 1. Here, we have taken 𝜖 = 0.04, 𝛿 = 0.1, 𝛾 = 0.5 and 𝑄 = 1/3. The limiting values for K are shown by red dashed lines. An explanation of the behaviour near St = 0.26 is … view at source ↗
Figure 24
Figure 24. Figure 24: The proportion of particles, R, leaving through the main channel rather than through the branched channels, calculated using (10.5), and plotted as black points between (a) St ∈ [0, 0.4] with spacing 0.01 and (b) St ∈ [1, 10] with spacing 1. We plot the corresponding values of R for figure 23. Here, we have taken 𝜖 = 0.04, 𝛿 = 0.1, 𝛾 = 0.5 and 𝑄 = 1/3. The limiting values for R are shown by red dashed lin… view at source ↗
Figure 25
Figure 25. Figure 25: The proportion of particles, R, leaving through the main channel rather than through the branched channels, calculated using (10.5), and plotted as black points on a larger range of St, on a log scale, over figure 24. Here, we have taken 𝜖 = 0.04, 𝛿 = 0.1, 𝛾 = 0.5 and 𝑄 = 1/3. The limiting value at St = ∞ for R is shown by a red dashed line. showing that it is possible to divert a reasonable fraction of t… view at source ↗
Figure 26
Figure 26. Figure 26: Flow solution structure. The full high-Reynolds-number structure is given by [PITH_FULL_IMAGE:figures/full_fig_p031_26.png] view at source ↗
Figure 27
Figure 27. Figure 27: Numerical solution (in solid lines) and asymptotic parabolic prediction (in red [PITH_FULL_IMAGE:figures/full_fig_p032_27.png] view at source ↗
Figure 28
Figure 28. Figure 28: The inlet positions 𝑦0 of particles that exit the device through a branched channel versus the Stokes number, St. For each St, we release 𝑖 = 19, 999 particles at slightly perturbed initial points around 𝑖 equispaced points, and run 60 separate simulations, plotting each point once. Here, Pout = 0.4, Re = 1000, 𝜖 = 0.04, 𝛿 = 0.1, 𝜆 = 0.1 and 𝛾 = 0.5. (a) (b) 0.22 0.23 0.24 0.25 0.26 0.27 0.00 0.01 0.02 0.… view at source ↗
Figure 29
Figure 29. Figure 29: Zoomed in versions of figures 23 and 24 around St ∈ [0.22, 0.27]. that result in particles passing through the branched channels. We see that the large deviation in K and R at St = 0.26 is aligned with the steepest change in the 𝑦0 values in figure 28. We note that there are a finite number rapidly changing 𝑦0 regions. In figure 29, we zoom in on the behaviour of K and R for St near 0.26. We see that, wit… view at source ↗
read the original abstract

Fibres shed from our clothes during a washing machine cycle constitute around 35% of the primary microplastics in our oceans. Current conventional dead-end washing machine filters clog relatively quickly and require frequent cleaning. We consider a new concept, ricochet separation, inspired by the feeding process of manta rays, to reduce the cleaning frequency. In such a device, some fluid is diverted through branched channels whilst particles ricochet off the wall structure, forcing them back into the main flow and then into the dead-end filter. In this paper, we consider a simple branched-channel filter structure beneath a high-Reynolds-number laminar flow, in the case where the branch separation is much larger than the thickness of the viscous boundary layer. We use multiple-scales techniques to derive an effective leakage boundary condition, which smooths out localised effects in the flow velocity and pressure that arise due to the discrete branched channels, and then use this boundary condition to explicitly determine the flow away from the boundary. We find that our explicit solution compares well with an analogous numerical solution containing a discrete set of branched channels. We further consider the behaviour of individual spherical particles in the device, with their trajectories determined via a simple force balance model with a wall-bounce condition. We explore the dependence of the fraction of particles that flow into the branched channels on the particle's Stokes number. The resulting combined model is able to predict the relationship between the efficiency of a ricochet filter device and the design and operating parameters, avoiding the need to conduct extensive numerically challenging simulations.

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 paper develops a multiple-scales analysis for high-Reynolds-number laminar flow over a branched-channel filter where branch spacing greatly exceeds the viscous boundary-layer thickness. It derives an effective leakage boundary condition that homogenizes discrete-branch effects, obtains an explicit outer-flow solution, compares this solution to discrete-branch numerical simulations, and then uses a simple force-balance particle model with wall bounces to predict the fraction of particles entering the branches as a function of Stokes number, thereby relating device efficiency to design and operating parameters.

Significance. If the scale-separation assumption holds and the comparison is quantitatively sound, the work supplies a computationally inexpensive, parameter-free route to efficiency predictions that avoids full discrete simulations. The derivation from the governing equations under an explicit asymptotic assumption, together with the direct link to particle trajectories, constitutes a practical strength for microplastic-filter design.

major comments (2)
  1. [§4] §4 (numerical comparison): the claim that the explicit solution 'compares well' with the discrete-branch numerical solution is not accompanied by quantitative error measures (e.g., L2 or L∞ norms of velocity or pressure differences, or tabulated maximum relative errors). Without such metrics the support for the accuracy of the homogenized boundary condition remains qualitative and limits assessment of the approximation's robustness.
  2. [§5] §5 (particle trajectories): the dependence of branched-channel capture fraction on Stokes number is presented, yet the manuscript does not report the sensitivity of these curves to the wall-bounce restitution coefficient or to the precise form of the drag law; because the efficiency prediction rests on these trajectories, the lack of such checks weakens the quantitative reliability of the final design relations.
minor comments (2)
  1. [Figures 3–5] Figure captions and axis labels in the flow-field comparisons would benefit from explicit indication of which curves are analytic and which are numerical, together with a statement of the Reynolds number and branch-spacing ratio used.
  2. [Introduction] The introduction would be strengthened by a brief reference to prior multiple-scales treatments of perforated or ribbed walls to clarify the novelty of the leakage-condition derivation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive review and recommendation of minor revision. We address each major comment below and will incorporate the suggested quantitative support in the revised manuscript.

read point-by-point responses

  1. Referee: [§4] §4 (numerical comparison): the claim that the explicit solution 'compares well' with the discrete-branch numerical solution is not accompanied by quantitative error measures (e.g., L2 or L∞ norms of velocity or pressure differences, or tabulated maximum relative errors). Without such metrics the support for the accuracy of the homogenized boundary condition remains qualitative and limits assessment of the approximation's robustness.

    Authors: We agree that quantitative error measures would strengthen the comparison. In the revised manuscript we will compute and report the L2 and L∞ norms of the velocity and pressure differences between the explicit multiple-scales solution and the discrete-branch numerical simulation. We will also tabulate maximum relative errors at representative locations along the boundary to provide a more rigorous quantification of the homogenized boundary condition's accuracy. revision: yes

  2. Referee: [§5] §5 (particle trajectories): the dependence of branched-channel capture fraction on Stokes number is presented, yet the manuscript does not report the sensitivity of these curves to the wall-bounce restitution coefficient or to the precise form of the drag law; because the efficiency prediction rests on these trajectories, the lack of such checks weakens the quantitative reliability of the final design relations.

    Authors: We acknowledge that the particle-capture results depend on modeling choices for bounces and drag. In the revised manuscript we will add a sensitivity study showing the capture-fraction curves for restitution coefficients between 0.7 and 1.0 and will compare the baseline Stokes-drag results with those obtained when an added-mass term is included. These additional checks will be presented to demonstrate that the reported Stokes-number dependence remains robust. revision: yes

Circularity Check

0 steps flagged

Multiple-scales derivation of effective leakage BC is self-contained asymptotic analysis

full rationale

The paper applies standard multiple-scales homogenization to the high-Re laminar flow equations under the explicit assumption that branch separation greatly exceeds viscous boundary-layer thickness. This produces an effective leakage boundary condition directly from the governing PDEs and scale-separation ansatz; the outer flow solution and subsequent particle trajectories follow from that BC without any parameter fitting to target data or renaming of known results. No self-citation chain, uniqueness theorem, or fitted-input-called-prediction appears in the derivation. The numerical comparison is performed inside the same asymptotic regime and serves as validation rather than input. The derivation chain therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard continuum fluid assumptions and an explicit scale-separation hypothesis; no new entities are postulated and no parameters are fitted to the target efficiency data.

axioms (2)
  • domain assumption High-Reynolds-number laminar flow
    Invoked to justify the flow regime and the thin viscous boundary layer relative to branch spacing.
  • domain assumption Branch separation much larger than viscous boundary layer thickness
    Required for the validity of the multiple-scales reduction to an effective boundary condition.

pith-pipeline@v0.9.0 · 5586 in / 1357 out tokens · 44670 ms · 2026-05-18T00:33:56.009051+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

3 extracted references · 3 canonical work pages

  1. [1]

    R., YANG, H

    ACHARYA, S., SHAIDA, S. R., YANG, H. & NOUREDDINE, A. 2021 Microfibers from synthetic textiles as a major source of microplastics in the environment: A review. Text. Res. J. 91 (17-18), 2136–2156. AKARSU, C., KUMBUR, H. & KIDEYS, A. E. 2021 Removal of microplastics from wastewater through electrocoagulation-electroflotation and membrane filtration process...

    work page 2021

  2. [2]

    V., STROTHER, J

    DIVI, R. V., STROTHER, J. A. & PAIG-TRAN, E. W. M. 2018 Manta rays feed using ricochet separation, a novel nonclogging filtration mechanism. Sci. Adv. 4 (9). DRIS, R., IMHOF, H. I., SANCHEZ, W. & GASPERI, J. 2015 Beyond the ocean: contamination of freshwater ecosystems with (micro-) plastic particles. Environ. Chem. 12 (5), 539–550. ENTEN, A. C., LEIPNER,...

    work page 2018

  3. [3]

    2023 Suspension feeders as biological models to develop biomimetic filter modules and reduce microplastic emissions

    HAMANN, L. 2023 Suspension feeders as biological models to develop biomimetic filter modules and reduce microplastic emissions. PhD thesis, Universitäts-und Landesbibliothek Bonn. HAZLEHURST, A., TIFFIN, L., SUMNER, M. & TAYLOR, M. 2023 Quantification of microfibre release from textiles during domestic laundering. Environ. Sci. Pollut. Res. 30 (15), 43932...

    work page 2023