Pith. sign in

REVIEW 4 major objections 6 minor 28 references

Immersed boundary simulations of fluid shear-induced deformation of a cantilever beam

T0 review · 4 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read An immersed-boundary spring network reproduces Euler-Bernoulli cantilever deflection scaling, and adding porosity through a Darcy-type slip velocity reduces deformation at high permeability.

desk verdict A competent IB cantilever-beam study whose porous 2D extension is genuinely new, but whose Euler–Bernoulli scaling agreement is an empirical fit pending a proper EI–σ_b calibration. read the letter →

arxiv 1908.01560 v2 pith:SAYUZDVI submitted 2019-08-05 physics.flu-dyn physics.comp-ph

classification physics.flu-dynphysics.comp-ph
keywords immersedboundarymethodfluid-structureinteractioncantileverbeamEuler-BernoullitheoryporosityDarcy'slawshearflowbiofilmmechanics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper asks whether a simple immersed-boundary spring network can reproduce the classical Euler-Bernoulli beam law when a cantilever is bent by a shear flow. It simulates a two-dimensional elastic beam in a Couette-type channel, driving the flow with a moving top wall, and measures steady-state tip deflection as beam stiffness, length, and shear speed vary. The simulated deflections follow the predicted power laws: inversely proportional to spring stiffness and proportional to about the fourth power of beam length. The paper then extends the fiber equation with a Darcy-type slip velocity to model a porous beam, finding that permeability below $10^{-8}$ $cm^{2}$ leaves the beam effectively solid, while high permeability lets fluid pass through and reduces deflection. It also shows that rounding the beam's corners removes non-physical shape irregularities without changing the deflection by more than 1%.

What carries the argument

The engine of the model is the immersed-boundary force density: a triangular spring network on the beam (Eq. 9) whose link stiffness sigma_b converts local stretch into an elastic force spread onto the fluid grid, plus stiff tether springs for the fixed wall and the beam-wall attachment. The beam's motion is set by the fiber evolution equation (Eq. 4), which is modified for porosity by subtracting a Darcy slip velocity u_p = -(K/mu) grad p (Eqs. 14-15). The Euler-Bernoulli equation (Eq. 13) supplies the theoretical deflection exponents (d proportional to 1/EI and d proportional to $H_b^{4}$) against which the measured deflections are compared.

What would settle it

Measure the effective flexural rigidity of the same triangulated spring network by applying a known static end load in the numerical solver (or by computing the strain energy of a bent configuration) and check whether EI/sigma_b is constant across mesh sizes and stiffness values. If that ratio varies, or if the deflection exponents drift away from -1 and 4 under mesh refinement, the claimed correspondence with Euler-Bernoulli theory would break down.

Watch

Extended reading notes

Core claim

The central claim is that the immersed-boundary spring-network model, with a uniform stiffness sigma_b on each link, produces cantilever deflection that scales as d proportional to $sigma_b^{{-0.8995}}$ and d proportional to $H_b^{{3.8094}}$, in close agreement with the Euler-Bernoulli predictions d proportional to 1/EI and d proportional to $H_b^{4}$; the paper interprets this as evidence that the effective flexural rigidity of the triangulated structure is proportional to sigma_b. For the porous extension, replacing the no-slip fiber equation with the Darcy slip velocity u_p = -(K/mu) grad p gives a deflection that is unchanged for K less than about $10^{-8}$, slightly increased for intermediate permeability, and strongly decreased for high permeability as flow passes through the beam. The paper further shows that replacing sharp corners by circular fillets eliminates the non-physical kinks and protrusions at the free and fixed ends, with a relative change in tip deflection under 1%.

Load-bearing premise

The results rest on the assumption that the uniform spring network with stiffness sigma_b behaves as a continuum elastic beam whose flexural rigidity EI is a fixed constant multiple of sigma_b; the paper explicitly leaves the derivation of EI proportional to sigma_b to future work.

Editorial extensions

If this is right

  • Within the reported parameter range, the spring network can serve as a faithful stand-in for a thin elastic cantilever in shear flow, so subsequent studies of biologically motivated beams can trust the deflection scaling without resolving the continuum solid.
  • The porosity model gives a handle on how much flow passes through a deformable structure: below a permeability threshold the beam is effectively solid, and above it the same shear produces less deformation.
  • Smoothing the beam's corners is a cheap fix for the non-physical end distortions, with a tip-deflection penalty below 1%.
  • The three-regime response to shear velocity (linear, transitional, saturated) can be used to predict when a beam ceases to deform further as flow strength grows.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The observed exponents (-0.8995 and 3.8094) deviate from the ideal integers by a few percent; this may reflect finite-thickness or discrete-network corrections that could be tested by refining the triangle mesh and checking whether the exponents converge to exactly -1 and 4.
  • The non-monotonic deflection in the mid-range of permeability suggests two competing effects—enhanced horizontal fluid transport that increases the shear load versus reduced blockage—and could be mapped to a dimensionless Darcy number for a collapse of the data.
  • If the EI proportional to sigma_b relation can be derived analytically, the method would become a parameter-free bridge between a discrete spring network and continuum beam theory, allowing the same computational setup to be used for arbitrary beam shapes without recalibration.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. This paper develops a two-dimensional immersed boundary (IB) model of a cantilever beam attached to a wall inside a Couette shear flow. The beam is discretized as a triangulated network of springs with stiffness σb, and the fluid is described by the incompressible Navier-Stokes equations with IB forcing. The authors measure the steady tip deflection d as σb, beam length Hb, and top-wall velocity utop are varied; report power-law fits d ∝ σ_b^{-0.8995} and d ∝ H_b^{3.8094}; compare these with Euler-Bernoulli predictions; investigate corner irregularities and propose a smoothed 'fillet' shape; and extend the IB formulation to porous beams via Darcy slip (Eqs. 14–15), reporting the dependence of d on permeability K. The central stated goal is qualitative comparison with linear beam theory and a demonstration that the spring-network IB model yields physically plausible deflection behavior.

Significance. The paper addresses a relevant problem in fluid-structure interaction and bio-inspired flows. Its strengths are the clear presentation of the IB spring-network discretization, the explicit treatment of corner artifacts, and the extension of porous IB boundaries to a 2D solid region. The power-law fits are reported with normalized RMS errors and the flow visualizations are informative. However, the quantitative link between the spring-network stiffness and a continuum flexural rigidity is not established, the scaling exponents are fitted over a range that includes large deflections, and the porous model is an unvalidated postulate. The paper's contribution is therefore best viewed as a computational parameter study rather than a fully validated quantitative model. If the missing calibration and model validation are provided, the paper would have clear value for biofilm and cilia modeling.

major comments (4)
  1. [Section 3.2, Section 5] The paper interprets d ∝ σ_b^{-0.8995} as evidence that the effective flexural rigidity EI is proportional to σ_b. This inference requires the Euler-Bernoulli relation d ∝ 1/EI to hold for the spring network, but the equivalence between the spring-network parameters and EI is never established; Section 5 explicitly defers this to future work. Without an independent calibration, e.g., applying a known static load to the isolated spring network and fitting EI from the Euler-Bernoulli solution, the exponent agreement is an empirical power-law fit rather than a validation of continuum beam theory. This point is load-bearing for the paper's central comparison.
  2. [Section 3.2, Figure 4] The power-law fit in Figure 4b includes configurations with d/H_b as large as 0.65, which is outside the small-deflection assumption underlying Eq. (13). Because the large-deflection points are systematically included, the fitted exponent -0.8995 may deviate from -1 due to geometric nonlinearity rather than to any intrinsic relation between EI and σ_b. The authors should report fits restricted to the linear regime (e.g., d/H_b < 0.1) and show that the exponent is stable.
  3. [Section 3.1, Section 3.3, Figure 5] The Euler-Bernoulli comparisons in §3.1 use a constant load and a load that decreases linearly from a maximum at the wall. In the simulated linear shear flow with no-slip at the bottom wall, the undisturbed load on a vertical beam is zero at the wall and grows with height; for a load q(y) ∝ y, the tip deflection is d = 11 q' H_b^5/(120 EI), scaling as H_b^5 rather than H_b^4 at fixed shear rate. The measured H_b^{3.8094} is therefore not directly comparable to the m=4 benchmark unless the actual load distribution is computed and its dependence on H_b is accounted for.
  4. [Section 4.1, Eqs. (14)–(15)] The porous-beam model is introduced as a postulate. Equation (15) states that the fiber moves with velocity -u_p + ∫ u δ, with u_p = -(K/μ)∇p, but no derivation, limiting-case check, or comparison against existing porous IB models (e.g., the 1D membrane limit of Kim and Peskin [18] or Stockie [19]) is provided. The physical interpretation of the negative sign is asserted but not tested. Because the permeability results in Section 4 are presented as predictions of the model, the authors should validate the porous formulation in a simple geometry against a known solution or against the non-porous limit.
minor comments (6)
  1. [Section 3.4, table] For utop=0.003, the listed ratio ‖u‖/utop is 0.17, which is not close to 0.1; the statement that the ratio is approximately 0.1 for all cases should be corrected or rephrased.
  2. [Section 4.1, paragraph] The phrase 'In this thesis' should read 'In this paper.'
  3. [Section 5, conclusions] 'Filled or round' should read 'filleted or rounded' in the sentence about smoothing the corners.
  4. [Section 2.1] The units of fIB are written as [g/s^2], but as a body force in the Navier-Stokes equations it should have units of force per unit volume, e.g., [g/(cm^2 s^2)]; please check the units of fIB and FIB for consistency.
  5. [Section 3.2, Section 3.3] The 'normalized RMS error' is not defined; please state the normalization used and the number of data points in the least-squares fits.
  6. [Figure 2b] The caption says the fillets have diameter equal to the beam width, but the text in Section 3.5 refers to circular arcs; specify the radius or arc geometry explicitly.

Circularity Check

1 steps flagged · score 2.0 of 10

Minor built-in σ_b scaling presented as material-property inference; otherwise simulations are self-contained against externally stated beam theory.

  1. other [Section 3.2, Eq. (9)]
    "Assuming that the linear theory holds here (thin beam, small deflection) our computational results suggest that the effective flexural rigidity of our triangulated spring link structure is directly proportional to the spring stiffness, or in other words that EI ∝ σb."

    Equation (9) defines the beam IB force as a linear spring law: F_b^ℓ = Σ_m σ_b I_{ℓm} (d_{ℓm}(0) − d_{ℓm}) d_{ℓm}/d_{ℓm}. For a fixed external fluid load, any equilibrium displacement therefore scales inversely with σ_b by construction. The observed d ∝ σ_b^{−0.8995} and the resulting suggestion 'EI ∝ σb' thus restate the model's linear-spring input under the name of an effective flexural rigidity, rather than independently testing the Euler-Bernoulli relation d ∝ 1/EI. The paper itself acknowledges that this relationship is not derived and defers it to future work (Section 5). The H_b scaling and porous-beam results are not affected by this step.

full rationale

The paper is a numerical study whose central outputs are steady-state tip deflections computed from an IB spring-network model. The Euler-Bernoulli formulas in Section 3.1 are quoted from an external textbook and serve as an independent benchmark. The σ_b test is the only near-circular step: because Eq. (9) is a linear spring law with stiffness σ_b, equilibrium deflection under a roughly fixed fluid load scales as 1/σ_b by definition; the inference EI ∝ σb is therefore an interpretation of a built-in model property, not an independent validation. Section 5 explicitly lists EI ∝ σb as future work, so the paper does not conceal the gap. The H_b-scaling study (Section 3.3) is genuinely independent: mesh geometry and load distribution change with beam length, and the fitted exponent 3.8094 is compared with the externally stated EB value 4. The porous extension (Section 4) is a stated Darcy-type postulate with qualitative consequences that are not used to define permeability. Self-citations [20,21] supply numerical algorithm details but are not load-bearing for the scaling claims. Overall, the derivation chain is self-contained and the flagged step is minor.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

The model rests on standard fluid mechanics, a spring-network representation of the beam whose continuum limit is unproven, and an ad hoc Darcy-style porous extension (eq. 15) that is the paper's own postulate. The fitted power-law exponents are derived from the simulations, not imposed as inputs.

free parameters (2)
  • Power-law exponent for d vs sigma_b = -0.8995
    Least-squares slope on log-log plot of steady-state tip deflection versus beam spring stiffness (Fig. 4b); used to infer EI proportional to sigma_b.
  • Power-law exponent for d vs Hb = 3.8094
    Least-squares slope on log-log plot of deflection versus beam length (Fig. 5); compared with theoretical exponent 4.
assumptions (6)
  • standard math Incompressible Newtonian fluid satisfies Navier-Stokes equations (eqs. 1-2).
    Standard governing equations for the fluid; accepted as the physical model for water in the channel.
  • domain assumption The immersed boundary method uses periodic boundary conditions in x and y and stiff tether springs to hold the walls (Section 2.2.1).
    The periodic domain and spring tether representation are modeling choices that affect the computed flow and deflection.
  • domain assumption The beam is represented by a triangulated network of linear springs with uniform stiffness sigma_b, based on the Alpkvist-Klapper biofilm model (eq. 9).
    This network model is assumed to behave like an elastic solid; the mapping to Young's modulus or EI is not established.
  • domain assumption Euler-Bernoulli beam theory with constant or linearly varying load applies to the small-deflection cases used for comparison (Section 3.1).
    The theory assumes thin beam and small deflection; some simulations exceed these limits, with deflections up to 65% of the beam length.
  • domain assumption Darcy's law with a single scalar permeability K relates the porous slip velocity to the local pressure gradient (eq. 14).
    Adopted from Kim-Peskin and Stockie; extends 1D normal-pore models to a 2D region with arbitrary flow direction.
  • ad hoc to paper The porous fiber evolves according to dX/dt = -up + integral(u delta) (eq. 15), with the structure moving opposite to the Darcy slip velocity.
    This is the paper's own modeling postulate for porosity in a solid 2D region; it is not derived and may not conserve mass.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Immersed boundary simulations of fluid shear-induced deformation of a cantilever beam." pith.science (2026). https://pith.science/paper/SAYUZDVI

@misc{pith2026190801560,
  author       = {Pith},
  title        = {Pith review of: Immersed boundary simulations of fluid shear-induced deformation of a cantilever beam},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SAYUZDVI}},
  note         = {Machine review of arXiv:1908.01560}
}
read the original abstract

We derive a mathematical model and the corresponding computational scheme to study deflection of a two-dimensional elastic cantilever beam immersed in a channel, where one end of the beam is fixed to the channel wall. The immersed boundary method has been employed to simulate numerically the fluid-structure interaction problem. We investigate how variations in physical and numerical parameters change the effective material properties of the elastic beam and compare the results qualitatively with linear beam theory. We also pay careful attention to "corner effects" -- irregularities in beam shape near the free and fixed ends -- and show how this can be remedied by smoothing out the corners with a "fillet" or rounded shape. Finally, we extend the immersed boundary formulation to include porosity in the beam and investigate the effect that the resultant porous flow has on beam deflection.

Figures

Figures reproduced from arXiv: 1908.01560 by the authors.

Figure 1
Figure 1. A general elastic boundary Γ = S3 i=1 Γ i that consists of several possibly disconnected components is immersed within a doubly-periodic rectangular fluid domain Ω. • The motion of the fluid is governed by incompressible Navier-Stokes equations. The variables involved are: u [cm/s](fluid velocity), p [g/cm s2 ](pressure), ρ [g/cm3 ] (density), µ [g/cm s](dynamic viscosity), fIB (q, t) [g/s 2 ] ( IB forcing term whic… view at source ↗
Figure 2
Figure 2. Unstructured triangular mesh for the two beam shap [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Flow streamlines (left) and pressure contours (ri [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: Deflection of the solid cantilever beam as [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Plot of Hb and final deflection d for the non-porous cantilever. Parameters: σb = 560 g/cm s2 , utop = 0.02 cm/s. Computed results are shown using diamond points (⋄) on top of which is superimposed a solid line corresponding to the linear least squares fit (normalized …
Figure 6
Figure 6. Figure 6: Deflected beam for several values of utop, showing streamlines (left) and contours of velocity magnitude (right). Other parameter values are Hb = 0.0077 cm, σb = 560 g/cm s2 . 10 [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: Maximum tip deflection as a function of utop for the solid beam, with parameters Hb = 0.0077 cm and σb = 560 g/cm s2 . The results of these simulations, as well as for a number of other flow speeds, are summarized in [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]
Figure 8
Figure 8. Figure 8: Zoomed view of the large deformations near the corn [PITH_FULL_IMAGE:figures/full_fig_p012_8.png]
Figure 9
Figure 9. Figure 9: Zoomed view of the two ends of the smoothed beam when [PITH_FULL_IMAGE:figures/full_fig_p012_9.png]
Figure 10
Figure 10. Figure 10: Maximum tip deflection for the rectangular and smo [PITH_FULL_IMAGE:figures/full_fig_p013_10.png]
Figure 11
Figure 11. Figure 11: Streamline plots at the point of maximum deflectio [PITH_FULL_IMAGE:figures/full_fig_p015_11.png]
Figure 12
Figure 12. Figure 12: Final deflection (d) versus permeability (K) for the smoothed porous beam. surrounding viscous, incompressible fluid flow. The choice of cantilever was motivated by our efforts to model biofilm structures in the near future. Our study included careful validation of th…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [18]

    Y. Kim, C. S. Peskin, 2–d parachute simulation by the imm ersed boundary method, SIAM Journal on Scientific Computing 28 (6) (2006) 2294–2312

  2. [19]

    J. M. Stockie, Modelling and simulation of porous immer sed boundaries, Computers & Structures 87 (11) (2009) 701–7 09

  3. [1]

    Alpkvist, I

    E. Alpkvist, I. Klapper, Description of mechanical resp onse including detachment using a novel particle model of biofilm/flow interaction, W ater Science and Technology 55 (8 -9) (2007) 265–273

  4. [2]

    Timoshenko, D

    S. Timoshenko, D. Young, Elements of strength of materia ls (1968)

  5. [3]

    Kamrin, C

    K. Kamrin, C. H. Rycroft, J.-C. Nave, Reference map techn ique for finite-strain elasticity and fluid–solid interacti on, Journal of the Mechanics and Physics of Solids 60 (11) (2012) 1952–1969

  6. [4]

    Pozrikidis, Shear flow over cylindrical rods attached to a substrate, Journal of Fluids and Structures 26 (3) (2010 ) 393–405

    C. Pozrikidis, Shear flow over cylindrical rods attached to a substrate, Journal of Fluids and Structures 26 (3) (2010 ) 393–405

  7. [5]

    Pozrikidis, Shear flow past slender elastic rods attac hed to a plane, International Journal of Solids and Structur es 48 (1) (2011) 137–143

    C. Pozrikidis, Shear flow past slender elastic rods attac hed to a plane, International Journal of Solids and Structur es 48 (1) (2011) 137–143

  8. [6]

    Farjoun, D

    Y. Farjoun, D. G. Schaeffer, A thin cantilever beam in a flow , in: AIP Conference Proceedings, Vol. 1389, 2011, pp. 1446–1449. 16

Show all 28 references
  1. [7]

    L. P. Li, K. Schulgasser, G. Cederbaum, Theory of poroela stic beams with axial diffusion, Journal of the Mechanics and Physics of Solids 43 (12) (1995) 2023–2042

  2. [8]

    A. Goza, T. Colonius, A strongly-coupled immersed-boun dary formulation for thin elastic structures, Journal of Co mpu- tational Physics 136 (2017) 401–411

  3. [9]

    X. Yang, C. W ang, A nonlinear mathematical model for larg e deflection of incompressible saturated poroelastic beams , Applied Mathematics and Mechanics (English Edition) 28 (12 ) (2007) 1587–1595

  4. [10]

    L. P. Li, K. Schulgasser, G. Cederbaum, Large deflection analysis of poroelastic beams, International Journal of No nlinear Mechanics 33 (1) (1998) 1–14

  5. [11]

    X. Yang, Q. W en, Dynamic and quasi-static bending of sat urated poroelastic timoshenko cantilever beam, Applied Ma th- ematics and Mechanics (English Edition) 31 (8) (2010) 995–1 008

  6. [12]

    L. J. Fauci, R. Dillon, Biofluidmechanics of reproducti on, Annu. Rev. Fluid Mech. 38 (2006) 371–394

  7. [13]

    C. S. Peskin, The immersed boundary method, Acta numeri ca 11 (2002) 479–517

  8. [14]

    C. S. Peskin, Flow patterns around heart valves: a numer ical method, Journal of computational physics 10 (2) (1972) 252–271

  9. [15]

    Dillon, L

    R. Dillon, L. Fauci, A microscale model of bacterial and biofilm dynamics in porous media, Biotechnology and Bioengi - neering 68 (5) (2000) 536–547

  10. [16]

    R. H. Dillon, L. J. Fauci, C. Omoto, X. Yang, Fluid dynami c models of flagellar and ciliary beating, Annals of the New York Academy of Sciences 1101 (1) (2007) 494–505

  11. [17]

    E. A. Schwartz, M. L. Leonard, R. Bizios, S. S. Bowser, An alysis and modeling of the primary cilium bending response t o fluid shear, American Journal of Physiology-Renal Physiolo gy 272 (1) (1997) F132–F138

  12. [20]

    Ghosh, The immersed boundary method for simulating g ravitational settling and fluid shear-induced deformation of elastic structures, Ph.D

    S. Ghosh, The immersed boundary method for simulating g ravitational settling and fluid shear-induced deformation of elastic structures, Ph.D. thesis, Science: Department of M athematics (2013)

  13. [21]

    Ghosh, J

    S. Ghosh, J. M. Stockie, Numerical simulations of parti cle sedimentation using the immersed boundary method, Comm u- nications in Computational Physics 18 (02) (2015) 380–416

  14. [22]

    P. O. Persson, G. Strang, A simple mesh generator in matl ab, SIAM review 46 (2) (2004) 329–345

  15. [23]

    Z. J. Huang, J. M. Tarbell, Numerical simulation of mass transfer in porous media of blood vessel walls, American Jou rnal of Physiology-Heart and Circulatory Physiology 273 (1) (19 97) H464–H477

  16. [24]

    M. C. M. van Loosdrecht, J. J. Heijnen, H. Eberl, J. Kreft , C. Picioreanu, Mathematical modelling of biofilm structur es, Antonie van Leeuwenhoek 81 (2002) 245–256

  17. [25]

    V. T. Nguyen, E. Morgenroth, H. J. Eberl, A mesoscale mod el for hydrodynamics in biofilms that takes microscopic flow effects into account, W ater science and technology 52 (7) (20 05) 167–172

  18. [26]

    Thullner, P

    M. Thullner, P. Baveye, Computational pore network mod eling of the influence of biofilm permeability on bioclogging in porous media, Biotechnology and Bioengineering 99 (6) (200 8) 1337–1351

  19. [27]

    Davit, H

    Y. Davit, H. Byrne, J. Osborne, J. Pitt-Francis, D. Gava ghan, M. Quintard, Hydrodynamic dispersion within porous biofilms, Physical Review E 87 (1) (2013) 012718

  20. [28]

    A. T. Layton, Modeling water transport across elastic b oundaries using an explicit jump method, SIAM Journal on Scientific Computing 28 (6) (2006) 2189–2207. 17

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.