{"id":"9247d818-426a-439f-a809-ee0842edeb23","arxiv_id":"2502.07376","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A new analytical model computes acoustic microstreaming and shear stress on a viscoelastic particle near an oscillating bubble, predicting stresses much larger than standard wall-based estimates.","lead":"This paper builds an analytical model of the acoustic microstreaming and shear stress produced when an oscillating gas bubble sits near a viscoelastic particle, removing earlier restrictions on separation and viscous layer thickness. It predicts stress values far above the commonly used Nyborg wall formula and identifies new possible rupture sites on the particle surface.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The particle stress is computed from the Lagrangian streaming velocity imported from the two-bubble theory, without a derivation that this equals the physical Eulerian-mean stress; the 25.2 Pa vs 1.1 Pa comparison therefore rests on an unvalidated transfer.","rationale":"The paper has a plausible and internally coherent first-order solution, and the series-convergence checks (Eq. 3.1, Fig. 3) show that the truncation error in the stress sums is controlled. The first-order boundary conditions (Eqs. 2.22-2.27) are standard, and the viscoelastic particle model (Eq. 2.12) is a reasonable idealization. The main risk is the second-order step, because the entire quantitative claim depends on the applicability of the two-bubble streaming equations after replacing the right bubble by the particle. The reader identified this transfer as the weakest assumption. I agree that it is load-bearing, but I specify it more narrowly: the particle stress is computed from the Lagrangian streaming velocity, and it is not shown that this equals the physical stress in the liquid. Because the particle surface oscillates, the Stokes drift at the surface is not obviously zero, so imposing v_L = 0 and using v_L in the stress formula may not be equivalent to the Eulerian no-slip condition. This is a correctness risk rather than a known error; the authors may have a convention in which the Lagrangian streaming is the reported quantity, but the physical stress claim requires a derivation. The proposed check (dropping the Stokes-drift terms in Eq. 2.112) is a direct, low-cost way to see whether the distinction matters. If the 25.2 Pa value is stable under this change, the concern largely resolves; if not, the central comparison is unsupported and the verdict should be conditional pending a corrected stress calculation or a DNS validation. This does not justify rejection, because the first-order machinery and the qualitative streaming patterns may still be valuable, and the issue is addressable.","tokens_in":26613,"tokens_out":22082,"duration_ms":206123,"concrete_test":"Recompute the Fig. 6(b) shear-stress curve using the Eulerian streaming velocity instead of the Lagrangian one: in Eq. (2.112), remove the V_S^(2)-dependent Stokes-drift terms, re-impose the no-slip condition on the Eulerian streaming velocity at r2=R20, and compare the maximum stress with the quoted 25.2 Pa. If the maximum changes by more than roughly 20%, the head-to-head with Nyborg's 1.1 Pa is not robust. An independent finite-volume or lattice-Boltzmann simulation of the same bubble-particle configuration would settle the transfer definitively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim (25.2 Pa vs 1.1 Pa, Figs. 6 and 7) rests on Section 2.5, where the second-order streaming solution of Doinikov et al. (2022) is transferred from a two-bubble geometry to a bubble and a no-slip viscoelastic particle. The transfer is asserted rather than derived. In the transfer, the particle boundary condition is imposed on the Lagrangian streaming velocity v_L = v_E + v_S (Eq. 2.93), and the time-averaged shear stress is then computed from v_L (Eqs. 2.109-2.112), with the Stokes-drift contributions V_S^(2) explicitly included. For a Newtonian liquid, the time-averaged shear stress on a material surface is fixed by the actual (Eulerian) mean velocity field and the first-order displacement-stress correlation; it is not in general equal to the gradient of the Lagrangian mean velocity. The paper supplies no derivation showing that the Lagrangian-based expression is the physical stress, nor that the free-slip two-bubble streaming solution form remains valid when one scatterer is replaced by a no-slip solid whose surface is oscillating. If the Eulerian and Lagrangian formulations differ, the quantitative stress predictions and the comparison with Nyborg's formula have no support.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytical theory for acoustic microstreaming and time-averaged shear stress produced by an oscillating gas bubble near a viscoelastic particle. The first-order problem is solved by coupling linearized viscous-liquid equations to a linear viscoelastic-solid model, with boundary conditions that prescribe the bubble surface motion, zero tangential stress at the bubble, and continuity of velocity and stress at the particle. The linear scattering coefficients are obtained from a truncated system of equations. The second-order streaming is taken from a previous two-bubble theory by Doinikov et al. (2022), with the linear coefficients replaced and the particle boundary condition changed from zero tangential stress to zero tangential Lagrangian velocity. Numerical examples illustrate streaming patterns for different bubble modes and particle materials, and the shear stress on the particle surface is computed. For a radially oscillating bubble near a large particle, the predicted maximum stress is 25.2 Pa, versus 1.1 Pa from Nyborg's formula.","tokens_in":26865,"tokens_out":9517,"duration_ms":81138,"significance":"If valid, the theory would provide a general analytical framework for computing acoustic microstreaming and particle-surface shear stress without the usual restrictions on the ratios of radii to viscous penetration depth or separation. The prediction that Nyborg's wall-based formula can underestimate peak bubble-particle stress by more than an order of magnitude is practically relevant for sonoporation and microfluidic applications. Strengths of the paper include a carefully posed linear scattering problem with a clear truncation scheme and a reported convergence check, and the absence of fitted parameters in the numerical examples. However, the second-order results depend on an unvalidated transfer of the two-bubble streaming equations to the bubble-particle system and on an unjustified identification of the physical shear stress with a Lagrangian-velocity-based expression, so the quantitative claims are not yet fully supported.","major_comments":[{"comment":"The second-order streaming solution is transferred from the two-bubble theory of Doinikov et al. (2022) by simply substituting the linear scattering coefficients and changing the tangential boundary condition at the particle from zero stress to zero velocity. The authors do not derive the second-order problem for the bubble–particle geometry or show that the particular solutions and forcing terms used in the two-bubble derivation remain valid when one scatterer is a no-slip viscoelastic solid. Since the entire streaming field and all stress results in Section 3 rest on this transfer, the authors should either provide a self-contained derivation of the second-order equations for the bubble–particle system or demonstrate explicitly that the governing equations and boundary conditions reduce to the same functional form as in Doinikov et al. (2022).","section":"§2.5, Eqs. (2.93)–(2.98)"},{"comment":"The time-averaged shear stress is computed by inserting the Lagrangian streaming velocity v_L = v_E + v_S into the Newtonian viscous-stress formula. For a Newtonian liquid, the physical time-averaged stress on a material surface is determined by the Eulerian mean velocity and the first-order oscillatory displacement–stress correlation; it is not in general equal to the gradient of the Lagrangian mean velocity. The paper supplies no derivation showing that the Lagrangian-based expression equals the physical stress, and the no-slip condition is imposed on the Lagrangian tangential velocity (Eq. (2.93)) rather than on the Eulerian mean velocity. The central quantitative claim, including the 25.2 Pa versus 1.1 Pa comparison with Nyborg's formula, depends on this identification, so the authors must either justify it or recompute the stress from the correct Eulerian-mean velocity field.","section":"§2.6, Eqs. (2.109)–(2.112)"},{"comment":"The comparison between the present theory and Nyborg's formula is presented as evidence that wall-based estimates underestimate the stress by more than an order of magnitude. The comparison is only conclusive if the stress formula in (2.112) is verified, which is not established due to the issue in Major Comment 2. In addition, the paper does not discuss the range of validity of Nyborg's thin-boundary-layer approximation for the chosen parameters (gap 10 µm and viscous penetration depth of about 0.8 µm), which is needed to interpret the difference. The authors should provide a quantitative justification for the applicability of Nyborg's formula in this configuration.","section":"§3, Figs. 6 and 7"},{"comment":"The boundary condition at the particle is imposed on the Lagrangian streaming velocity at the equilibrium radius R20, but the physical no-slip condition applies to the instantaneous velocity and the mean-flow condition for an oscillating boundary is not self-evident. The authors do not derive the correct mean boundary condition for a viscoelastic particle whose surface is oscillating, nor do they show that the Lagrangian and Eulerian formulations coincide at the surface. This is a load-bearing part of the derivation and should be clarified.","section":"§2.5, Eqs. (2.93) and (2.97)"}],"minor_comments":[{"comment":"The statement that scaling results by the product of the mode amplitudes ab is sufficient to obtain results for amplitudes a and b is not generally correct for multimode excitation, because the streaming contains quadratic terms s_n s_m for all pairs, including self-interactions, so the total is not a simple product of the two amplitudes.","section":"§3, text after Eq. (3.1)"},{"comment":"The convergence of the truncated system (2.72) is demonstrated only for the stress in Section 3; the authors should also report the dependence of the streaming velocity field on the truncation order N, especially for the large-particle case in Figure 6 where N=35 is used.","section":"§2.4"},{"comment":"The claim that no restrictions are imposed on the ratios of radii to viscous penetration depth and separation distance is not supported by a convergence analysis in extreme parameter regimes; the numerical examples cover only a limited range of these ratios, so the claim should be either demonstrated or softened.","section":"Abstract and §1"},{"comment":"The comparison of stress for different particle materials would benefit from a table of the maximum stress values and their locations, since the stress curves are difficult to read in a single panel.","section":"§3, Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the first-order derivation appears sound. The main risk is the unverified transfer of the streaming equations from the two-bubble case and the treatment of the mean stress; both are fixable in principle but require substantial additional derivation or justification. I would also recommend that the authors clarify the relationship between their Lagrangian formulation and the physical stress on the particle surface, and that they report convergence of the velocity field as well as the stress."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline is this: the paper is not a repackaging. Solving the linear scattering problem for an oscillating bubble next to a finite viscoelastic particle, with no restrictions on radii, separation, or viscous penetration depth, is genuinely new. Earlier cell models used an infinite wall; the authors' own 2022 work treated two bubbles. The six boundary conditions, the viscoelastic constitutive law, the truncation scheme, and the convergence checks are all handled carefully. I also credit the absence of fitting: bubble mode amplitudes are inputs, material constants come from the literature, and the Nyborg comparison is an external benchmark. That is honest and reproducible in the analytical sense.\n\nThe soft spot is exactly where the stress-test note lands: Section 2.5. The paper takes the second-order streaming solution from the two-bubble theory and asserts that it carries over to a bubble and a no-slip particle, with only the constants re-solved and the tangential boundary condition changed. That is plausible—the governing second-order equations are the same—but it is not derived, and the solution class from two free-slip bubbles may not automatically cover an oscillating no-slip solid surface. The related issue is that Eq. (2.109) computes the time-averaged shear stress from the Lagrangian streaming velocity. The stress-test note is right that the physical stress on a material surface is not in general the gradient of the Lagrangian mean velocity. The paper neither cites nor confronts the Eulerian-versus-Lagrangian stress distinction. I would not call this fatal; the linear scattering part looks sound, and the streaming transfer is a reasonable hypothesis. But the 25.2 Pa vs 1.1 Pa claim rests on that hypothesis, so the quantitative headline is less secure than the abstract implies.\n\nA secondary weakness is the absence of any direct numerical simulation or experiment. The comparison against Nyborg's approximation is an analytical benchmark, and Nyborg's approximation is itself known to be rough at small separations. One DNS check for a representative case would substantially raise my confidence.\n\nThe citation pattern is legitimate. The heavy reliance on Doinikov et al. (2022) is a foundation, not a flaw, given that the paper is explicitly extending that machinery.\n\nWho is this for? The ultrasound-microbubble and microfluidics communities, especially people modeling sonoporation or interpreting shear-stress-driven membrane effects. The paper deserves a serious referee. I would send it to review and ask the authors to either derive the streaming transfer more carefully or benchmark it numerically, and to address the Lagrangian/Eulerian stress point. With those addressed, the paper would be a solid contribution.\n\nRecommendation: accept for peer review, with the expectation of a revision.","headline":"A serious analytical extension that solves a genuinely new boundary-value problem, with a concrete and striking stress prediction, but the load-bearing transfer of the two-bubble streaming equations to a bubble–particle system is asserted rather than derived.","tokens_in":27387,"tokens_out":2542,"would_cite":true,"duration_ms":28612,"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 analytical theory with no size or distance restrictions shows that an oscillating bubble exerts a peak shear stress of 25.2 Pa on a nearby viscoelastic particle, more than twenty times the 1.1 Pa given by the standard wall formula.","keywords":["acoustic microstreaming","gas bubble","viscoelastic particle","shear stress","shape modes","sonoporation","microfluidics"],"falsifier":"Measure the steady streaming velocity field around a 5-µm bubble oscillating at 500 kHz with 1-µm radial amplitude next to a 5-µm steel sphere in water with 20 µm centre-to-centre separation; the theory predicts a peak time-averaged shear stress of 25.2 Pa on the sphere, whereas the standard wall formula gives 1.1 Pa. A micro-PIV or particle-tracking experiment that yields a peak stress near 1.1 Pa would falsify the central claim.","tokens_in":26386,"feed_emoji":"🫧","tokens_out":8552,"duration_ms":68770,"temperature":0.7,"pith_summary":"This paper develops an analytical theory for the steady vortex flow, acoustic microstreaming, that an oscillating gas bubble drives around a nearby viscoelastic particle, and for the time-averaged shear stress this flow exerts on the particle surface. The theory is stated without restrictions on the bubble or particle radii relative to the viscous penetration depth or on the bubble–particle separation, so it covers configurations where the usual thin-boundary-layer approximation is not valid. Its central quantitative claim is that for 5 µm bubble and particle radii, 20 µm centre-to-centre separation, 500 kHz driving, and 1 µm oscillation amplitudes, the peak shear stress is 25.2 Pa, against 1.1 Pa from the standard wall-based formula. If the theory is right, shear-stress estimates for bubble-driven cell damage and sonoporation that treat the cell as a rigid wall can underestimate the peak stress by more than an order of magnitude.","feed_headline":"Bubble shears nearby particles 20x harder than standard estimate","feed_subtitle":"A no-restriction analytic model puts peak stress at 25.2 Pa, not the usual 1.1 Pa.","key_machinery":"The load-bearing object is the set of linear scattering coefficients $a_n^{(j)}$ and $b_n^{(j)}$ that encode how the liquid velocity field is re-radiated by the bubble and the particle, in expansions over axisymmetric spherical harmonics with scalar and vector potentials. These coefficients solve a truncated infinite system of boundary conditions: the normal velocity matches the bubble surface motion, the tangential stress vanishes on the bubble, and both velocity and stress are continuous at the particle. The second-order step reuses the authors' earlier two-bubble streaming derivation, replacing the second bubble with the particle and changing the boundary condition at the particle from zero tangential stress to zero tangential Lagrangian streaming velocity; this modified condition fixes a new set of constants that enter the final shear-stress formula on the particle surface.","core_discovery":"The authors claim that the shear stress on a viscoelastic particle near an oscillating bubble can be computed within an analytical framework that imposes no small-parameter assumptions linking the viscous penetration depth to the particle radius or to the bubble–particle gap. They solve the first-order linear scattering problem by expanding the liquid velocity in scalar and vector potentials and the particle displacement in viscoelastic wave solutions, then feed those coefficients into a second-order streaming calculation adapted from their earlier two-bubble theory. The only conceptual change at the particle is that the Lagrangian streaming velocity, not the tangential stress, is required to vanish at its surface. For a steel particle with $R_{10}=R_{20}=5$ µm, $d=20$ µm, $f=500$ kHz, and $s_0=s_1=1$ µm, the computed maximum shear stress is 25.2 Pa, compared with 1.1 Pa from the standard wall formula; at smaller gaps the theory gives 155.4 Pa and 534.5 Pa versus 7.4 Pa and 49.7 Pa.","pith_inferences":["Editorial extension: because the theory imposes no gap or penetration-depth restrictions, it could provide quantitative shear-stress inputs for sonoporation models where cells are positioned at a controlled standoff from a bubble; the paper itself stops at the stress calculation.","Editorial extension: the finding that soft particles shift the rupture site suggests a testable prediction for experiments with hydrogel or lipid-shell particles of varying stiffness: the spatial distribution of surface damage should shift from the bubble-facing pole to the equator as the Young's modulus drops.","Editorial extension: a direct experimental check could use the computed streaming patterns, not just the stress: the theory predicts four small counter-rotating vortices at the particle's bubble-facing side for a radially pulsating bubble, a feature absent in single-bubble streaming."],"forward_implications":["Bubble-induced shear stress on finite particles is not captured by the rigid-wall approximation; the maximum stress depends strongly on the finite particle geometry and can exceed wall-based estimates by more than an order of magnitude at typical separations.","The discrepancy shrinks as the bubble approaches the particle: at 5 µm and 1 µm surface-to-surface gaps the theory gives stresses about 21 and 11 times the wall-formula values, so the largest errors occur at larger gaps.","For stiff particles the most probable rupture site is the pole facing the bubble, but for erythrocyte-like soft particles the equator becomes a comparable rupture location, a result that follows only when the particle is treated as finite and deformable.","Increasing the particle's shear viscosity raises the maximum shear stress and the stress gradient at the pole, implying that more viscous cells are expected to be more susceptible to rupture under the same acoustic drive."],"supporting_citations":[{"why":"Supplies the second-order streaming equations and boundary-condition machinery that the bubble–particle theory adapts to a viscoelastic particle.","marker":"Doinikov et al. (2022)"},{"why":"The standard boundary-layer approximation whose time-averaged wall shear stress is the baseline that the paper's predictions are compared against.","marker":"Nyborg (1958)"},{"why":"Derives the rigid-wall stress formula used for the quantitative comparison; the paper's 25.2 Pa versus 1.1 Pa discrepancy is against this formula.","marker":"Doinikov & Bouakaz (2014)"},{"why":"Provides the single-bubble streaming solution for modes 0 and 1, used to identify that a lone radially pulsating bubble produces no streaming and to interpret the particle-induced flow.","marker":"Doinikov et al. (2019)"},{"why":"Gives the quadrupole streaming pattern for a translationally oscillating bubble, the reference pattern for one of the mode cases.","marker":"Longuet-Higgins (1998)"},{"why":"Supplies the linearized viscous-fluid equations and stress definitions on which the first-order and shear-stress calculations are built.","marker":"Landau & Lifshitz (1987)"}],"fun_headline_variants":["Bubble-induced shear on particles beats standard formula by 20x","Analytic model: bubble streaming stress up to 20x Nyborg's estimate","No-restriction theory reveals stronger bubble shear on particles","Bubble streaming stress on particles: 20x higher than classic formula","Shear from bubble streaming on particles up to 20x Nyborg's"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The theory stands or falls on the assumption that the second-order streaming equations derived for two interacting bubbles remain valid for a bubble and a viscoelastic particle once the linear scattering coefficients are replaced and the particle boundary condition is switched from zero tangential stress to zero tangential velocity; if this transfer is not legitimate, the quantitative stress predictions have no foundation.","fun_headline_variants_meta":{"raw":{"variants":["Bubble-induced shear on particles beats standard formula by 20x","Analytic model: bubble streaming stress up to 20x Nyborg's estimate","No-restriction theory reveals stronger bubble shear on particles","Bubble streaming stress on particles: 20x higher than classic formula","Shear from bubble streaming on particles up to 20x Nyborg's"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000419,"raw_usage":{"total_tokens":2147,"prompt_tokens":926,"completion_tokens":1221,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":542,"completion_tokens_details":{"reasoning_tokens":1127}},"tokens_in":542,"tokens_out":1221,"duration_ms":8257,"temperature":1.0,"reasoning_tokens":1127,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T12:56:48.327687+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the steady streaming velocity field around a 5-µm bubble oscillating at 500 kHz with 1-µm radial amplitude next to a 5-µm steel sphere in water with 20 µm centre-to-centre separation; the theory predicts a peak time-averaged shear stress of 25.2 Pa on the sphere, whereas the standard wall formula gives 1.1 Pa. A micro-PIV or particle-tracking experiment that yields a peak stress near 1.1 Pa would falsify the central claim.","supporting_citations":[{"cited_title":"1958 Acoustic streaming near a boundary","cited_arxiv_id":null,"evidence_quote":"The standard boundary-layer approximation whose time-averaged wall shear stress is the baseline that the paper's predictions are compared against."}],"review_version":1}