{"id":"a8a66efe-bb92-43ca-9b4d-7ad53c172953","arxiv_id":"1908.01560","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A 2D elastic cantilever beam in shear flow is simulated with the immersed boundary method, including a new porous-beam extension that lets fluid flow through the structure.","lead":"This paper simulates how a flexible cantilever beam bends in a channel flow using the immersed boundary method, and extends the method to let fluid pass through a porous beam. It could be useful for modeling biofilms, cilia, and similar flexible biological structures in flow.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The EI–σ_b equivalence is assumed, not derived; without an independent calibration the Euler–Bernoulli scaling claim remains an empirical power-law fit.","rationale":"The most load-bearing concern is the unproven equivalence between the discrete spring network and a continuum Euler–Bernoulli beam, because the paper's headline result (reproduction of EB scaling) and its inference that EI ∝ σ_b both rest on it. The authors themselves list this as future work in Section 5, which is an explicit admission that the central comparison is incomplete. An independent static calibration of EI from the spring network would settle whether the FSI scaling data genuinely reflect beam-like elasticity or are coincidence. This does not impugn the numerical work: the simulations appear carefully set up, the smoothing of corner effects is a sensible remedy, and the observed power laws are plausible. But the central claim is conditional on this calibration, exactly as the reader concluded. I agree with the reader's weakest assumption and verdict.","tokens_in":14308,"tokens_out":10097,"duration_ms":110282,"concrete_test":"Run a static, fluid-free calibration: take the same triangular spring network (same mesh, same σ_b), clamp the base, apply a uniform horizontal load to all nodes, solve the spring equilibrium, and measure the tip deflection. Use the EB cantilever formula d = q0 H_b^4/(8 EI) to extract EI. Repeat for σ_b = 70, 140, 280, 560, 1120, 2240 and for H_b = 0.0056, 0.0077, 0.014. If the extracted EI is not proportional to σ_b with the same proportionality across H_b values, the paper's interpretation of the FSI scaling data is unsupported. Also recompute the FSI power-law fits using only cases with steady-state d/H_b < 0.2 and check whether the slopes move closer to −1 and 4.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim is that the IB spring-network model reproduces Euler–Bernoulli deflection scaling, d ∝ σ_b^{−1} and d ∝ H_b^4. This interpretation requires that the effective flexural rigidity EI of the triangulated spring network is proportional to σ_b and independent of H_b. The authors do not derive this equivalence; Section 5 explicitly defers it to future work. Without such a calibration, the fitted exponents (−0.8995 and 3.8094) are only empirical power laws, and the agreement with beam theory (to within 10% and 5%) is not a quantitative validation. The issue is compounded because the fits include cases with d/H_b up to 0.65, outside the small-deflection regime where EB theory applies; the deviations of the exponents from −1 and 4 may be systematic artifacts of including nonlinear deflections. A concrete check is to calibrate EI by applying a known static load to the isolated spring network (no fluid) and fitting the EB formula, then verify that EI ∝ σ_b and is independent of H_b. If that proportionality fails, the scaling comparison in Sections 3.2–3.3 does not support the EB interpretation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14572,"tokens_out":11348,"duration_ms":116245,"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":[{"comment":"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.","section":"Section 3.2, Section 5"},{"comment":"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.","section":"Section 3.2, Figure 4"},{"comment":"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.","section":"Section 3.1, Section 3.3, Figure 5"},{"comment":"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.","section":"Section 4.1, Eqs. (14)–(15)"}],"minor_comments":[{"comment":"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.","section":"Section 3.4, table"},{"comment":"The phrase 'In this thesis' should read 'In this paper.'","section":"Section 4.1, paragraph"},{"comment":"'Filled or round' should read 'filleted or rounded' in the sentence about smoothing the corners.","section":"Section 5, conclusions"},{"comment":"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.","section":"Section 2.1"},{"comment":"The 'normalized RMS error' is not defined; please state the normalization used and the number of data points in the least-squares fits.","section":"Section 3.2, Section 3.3"},{"comment":"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.","section":"Figure 2b"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a revised version of a 2020 preprint. The work is within the scope of Applied Mathematics and Computation. The primary concerns are the uncalibrated EI-σ_b relation and the unvalidated porous model; both are addressable with additional simulations. I do not see any citation concerns, although the authors could engage more explicitly with the slender-body load analysis of Pozrikidis when discussing the load distribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, honest IB simulation paper with one genuinely new piece (a 2D porous solid extension of the IB method) and one central interpretation that is not as secure as the abstract suggests.\n\nWhat it does well: the base IB model is standard, but the parametric study is careful. Varying σ_b, H_b, and shear rate yields power laws close to Euler–Bernoulli expectations, and the corner-effect work is practical and clearly presented—rounded fillets remove non-physical kinks with less than 1% change in tip deflection. The porous extension in Section 4 is the real novelty: instead of the 1D membrane with pores normal to the fiber (Kim & Peskin, Stockie), the authors allow a 2D porous region with flow driven by the local pressure gradient via Darcy's law. That is a genuine generalization, and the simulations show plausible behavior: at high permeability the beam deflects less because fluid passes through. The paper is also refreshingly candid about its limitations.\n\nWhere it is soft: the stress-test note is right. The claim that the spring-network beam reproduces EB scaling rests on the assumption that effective EI is proportional to σ_b and independent of H_b. That equivalence is not derived; the authors defer it to future work in Section 5. So the fitted exponents are empirical power laws, not a validation of EB theory. Worse, the fits include deflections up to 65% of the beam length—outside the small-deflection regime where EB applies—so the agreement could be partly fortuitous. A simple static-load calibration of the isolated spring network would settle this, and it is missing. Also, no code/data or grid-convergence study is included, making the quantitative claims hard to check. For the porous model, the authors do not explain how ∇p is evaluated on the Lagrangian points, and there is no benchmark or experimental comparison—the model is a physical postulate. These are addressable issues, not fatal flaws.\n\nWho it is for: people working on IB simulations of flexible or porous structures, especially in biofilm or cilia contexts. It is a useful contribution despite the caveats.\n\nRecommendation: I would send it to a serious referee. The porous extension is new and the paper is honest, but it needs revision: add the EI calibration, exclude or flag points outside the linear regime, and provide some verification of the porous implementation.","headline":"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.","tokens_in":15050,"tokens_out":2513,"would_cite":false,"duration_ms":26058,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["immersed boundary method","fluid-structure interaction","cantilever beam","Euler-Bernoulli beam theory","porosity","Darcy's law","shear flow","biofilm mechanics"],"falsifier":"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.","tokens_in":14109,"feed_emoji":"🌊","tokens_out":5736,"duration_ms":51957,"temperature":0.7,"pith_summary":"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%.","feed_headline":"Simulated cantilever bending tracks beam theory","feed_subtitle":"A spring-network model reproduces beam theory's 1/stiffness and length^4 laws, then adds porosity.","key_machinery":"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.","core_discovery":"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%.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the spring-network model of biofilms that is used for the beam's elastic forces (Eq. 9).","marker":"[1]"},{"why":"Provides the Euler-Bernoulli beam theory and the tip-deflection formulas against which simulations are compared.","marker":"[2]"},{"why":"The immersed-boundary formulation (spreading, interpolation, force density) follows Peskin's method summarized here.","marker":"[13]"},{"why":"Introduces porosity in the IB framework via a normal slip velocity, the starting point for the present porous extension.","marker":"[18]"},{"why":"Incorporates porosity directly via Darcy's law for a 1D membrane, which the paper extends to a 2D solid region.","marker":"[19]"},{"why":"Generates the uniform triangular mesh (DistMesh) used to discretize the beam and define the spring network.","marker":"[22]"}],"fun_headline_variants":["Fluid-sheared beam simulation matches Euler-Bernoulli scaling","Beam bend simulation confirms 1/stiffness and length^4 laws","Porosity in simulated beams: flow alters deflection in surprising ways","Fillets smooth corner artifacts in immersed-boundary beam model","Simulated cantilever bending: theory matched, then extended to porous"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Fluid-sheared beam simulation matches Euler-Bernoulli scaling","Beam bend simulation confirms 1/stiffness and length^4 laws","Porosity in simulated beams: flow alters deflection in surprising ways","Fillets smooth corner artifacts in immersed-boundary beam model","Simulated cantilever bending: theory matched, then extended to porous"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000805,"raw_usage":{"total_tokens":3500,"prompt_tokens":874,"completion_tokens":2626,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":490,"completion_tokens_details":{"reasoning_tokens":2535}},"tokens_in":490,"tokens_out":2626,"duration_ms":22747,"temperature":1.0,"reasoning_tokens":2535,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:09:02.740797+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Alpkvist, I","cited_arxiv_id":null,"evidence_quote":"Supplies the spring-network model of biofilms that is used for the beam's elastic forces (Eq. 9)."},{"cited_title":"Timoshenko, D","cited_arxiv_id":null,"evidence_quote":"Provides the Euler-Bernoulli beam theory and the tip-deflection formulas against which simulations are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The immersed-boundary formulation (spreading, interpolation, force density) follows Peskin's method summarized here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces porosity in the IB framework via a normal slip velocity, the starting point for the present porous extension."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Incorporates porosity directly via Darcy's law for a 1D membrane, which the paper extends to a 2D solid region."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Generates the uniform triangular mesh (DistMesh) used to discretize the beam and define the spring network."}],"review_version":1}