REVIEW 4 major objections 4 minor 5 references
Acoustic microstreaming and shear stress produced by the interaction of an oscillating gas bubble with a viscoelastic particle
T0 review · 4 major / 4 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [§2.5, Eqs. (2.93)–(2.98)] 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).
- [§2.6, Eqs. (2.109)–(2.112)] 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.
- [§3, Figs. 6 and 7] 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.
- [§2.5, Eqs. (2.93) and (2.97)] 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.
minor comments (4)
- [§3, text after Eq. (3.1)] 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.
- [§2.4] 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.
- [Abstract and §1] 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.
- [§3, Fig. 4] 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.
Circularity Check
No circularity: the claimed stress values are computed outputs from prescribed mode amplitudes and a prior independent streaming derivation; self-citation is not self-fulfilling.
full rationale
No load-bearing step in this paper reduces to its own inputs. The first-order scattering coefficients are obtained by solving a linear system with prescribed bubble mode amplitudes s_n and material constants; the streaming constants are re-calculated under a changed boundary condition (zero tangential Lagrangian velocity at the particle surface, Eq. 2.93), and the shear stress is then a computed output (Eqs. 2.109-2.112). No parameter is fitted to the 25.2 Pa result, and the Nyborg-based formula (3.5) is an independent external benchmark rather than an input to the model. The reliance on Doinikov et al. (2022) for the second-order streaming equations is a substantive prior derivation by overlapping authors, but the present paper does not define its target result in terms of that citation; it transfers and adapts the framework with a new boundary condition. Whether that transfer is physically valid is a correctness and rigor concern, not a circularity concern. The self-citation is load-bearing in the sense of providing the streaming formalism, but it is not unverified in a way that makes the present stress prediction equivalent to the paper's own assumptions. Accordingly, the appropriate finding is a low score reflecting normal reliance on prior work, with no circular step identified.
Assumptions & free parameters
free parameters (1)
- Bubble mode amplitudes s_n =
s0=s1=s2=1 um in numerical examples (500 kHz for modes 0,1 and 250 kHz for mode 2)
assumptions (5)
- domain assumption Liquid is viscous, incompressible, Newtonian; first-order flow obeys linearized Navier-Stokes equations (2.2)-(2.3).
- domain assumption Bubble and particle are spherical at rest; motion is axisymmetric; gas viscosity is negligible so tangential stress vanishes at bubble surface (2.23).
- domain assumption Particle is a homogeneous isotropic linear viscoelastic solid governed by equation (2.12) with constant Lame coefficients and viscosities.
- ad hoc to paper The two-bubble streaming framework of Doinikov et al. (2022) applies to bubble-particle systems with only the tangential boundary condition replaced by zero Lagrangian tangential velocity at the particle surface.
- standard math Time averages are taken over the oscillation period and the Lagrangian streaming velocity is the sum of Eulerian streaming and Stokes drift.
Cite this review
Pith. "Pith review of Acoustic microstreaming and shear stress produced by the interaction of an oscillating gas bubble with a viscoelastic particle." pith.science (2026). https://pith.science/paper/RKFTOR4I
@misc{pith2026250207376,
author = {Pith},
title = {Pith review of: Acoustic microstreaming and shear stress produced by the interaction of an oscillating gas bubble with a viscoelastic particle},
year = {2026},
howpublished = {\url{https://pith.science/paper/RKFTOR4I}},
note = {Machine review of arXiv:2502.07376}
}
read the original abstract
An analytical theory is developed that describes acoustic microstreaming produced by the interaction of an oscillating gas bubble with a viscoelastic particle. The bubble is assumed to undergo axisymmetric oscillation modes, which can include radial oscillation, translation and shape modes. The oscillations of the particle are excited by the oscillations of the bubble. No restrictions are imposed on the ratio of the bubble and the particle radii to the viscous penetration depth and the separation distance, as well as on the ratio of the viscous penetration depth to the separation distance. Capabilities of the developed theory are illustrated by computational examples. The shear stress produced by the acoustic microstreaming on the particle surface is calculated. It is shown that this stress is much higher than the stress predicted by the formula of Nyborg, which is commonly used to evaluate the time-averaged shear stress produced by a bubble on a rigid wall.
Figures
Reference graph
Works this paper leans on
-
[1]
Abramowitz, M. & Stegun, I. A. 1972 Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables. Dover. Deng, C.X. , Sieling, F., Pan, H. & Cui, J. 2004 Ultrasound-induced cell membrane porosity . Ultrasound Med. Biol. 30, 519–526. 33 Doinikov, A.A. & Bouakaz, A. 2010 Theoretical investigation of shear stress generated by a contrast m...
work page 1972
-
[152]
Tachibana, K., Uchida, T., Ogawa, K., Yamashita, N. & Tamura, K. 1999 Induction of cell - membrane porosity by ultrasound. Lancet 353,
work page 1999
-
[157]
1958 Acoustic streaming near a boundary
Nyborg, W.L. 1958 Acoustic streaming near a boundary. J. Acoust. Soc. Am. 30, 329–339. Nyborg, W.L. 1965 Acoustic streaming. In Physical Acoustics (ed. W.P. Mason), vol. IIB, pp. 266–331. Academic. Rooney, A. 1972 Shear as a mechanism for sonically induced biological effects. J. Acoust. Soc. Am. 52, 1718–1724. Saeidi, D., Saghafian, M., Javanmard, S.H. & ...
work page 1958
-
[222]
34 Nilsson, J., Evander, M., Hammarström, B. & Laurell, T. 2009 Review of cell and particle trapping in microfluidic systems. Anal. Chim. Acta 649, 141 –
work page 2009
-
[1409]
2007 Effects of shear stress on endothelial cells: Possible relevance for ultrasound applications
VanBavel, E. 2007 Effects of shear stress on endothelial cells: Possible relevance for ultrasound applications. Prog. Biophys. Mol. Biol. 93, 374–383. van Wamel, A., Kooiman, K., Harteveld, M., Emmer, M., ten Cate, F.J., Versluis, M. & de Jong, N. 2006 Vibrating microbubbles poking individual cells: Drug transfer into cells via sonoporation. J. Controlled...
work page 2007
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.