pith. sign in

arxiv: 2502.00863 · v2 · submitted 2025-02-02 · ⚛️ physics.flu-dyn

Nonspherical oscillations of an encapsulated magnetic microbubble

Pith reviewed 2026-05-23 04:18 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords encapsulated magnetic microbubblesnonspherical oscillationsmembrane theorystability regionMaxwell stressaxisymmetric deformationspolymer coatingmagnetic particles
0
0 comments X

The pith

The applied magnetic field does not influence the stability region for encapsulated magnetic microbubbles.

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

This paper develops a model for nonspherical oscillations of encapsulated magnetic microbubbles using membrane theory for thin weakly magnetic membranes. The model assumes generalized plane stress where only the applied magnetic field contributes to Maxwell stress. It computes the pressure-frequency stability region for axisymmetric linear oscillations and finds that the second mode dominates. Computational results indicate that varying the applied magnetic field does not change the stability region, although material properties and magnetic susceptibility do affect oscillation amplitudes. A reader would care if designing magnetic bubbles for applications where stability under magnetic influence is key.

Core claim

The study focuses on axisymmetric deformations under symmetrically arranged magnetic coils and restricts non-spherical oscillations to the linear regime. The pressure-frequency stability region is computationally determined and its variation with material properties and applied magnetic field is analyzed. The natural frequency of each mode is estimated using boundary layer approximation. Time-series analysis of the second mode amplitude reveals a significant oscillation amplitude relative to the bubble radius. Estimation indicates that the interface magnetic susceptibility and initial bubble radius enhance the amplitude of second-mode oscillations. Computational findings suggest that the 0.1

What carries the argument

Membrane theory for thin weakly magnetic membranes under generalized plane stress, with Maxwell stress contributed only by the applied magnetic field.

Load-bearing premise

The membrane is under generalized plane stress and only the applied magnetic field contributes to the Maxwell stress.

What would settle it

Observing whether the pressure-frequency stability boundaries change when an external magnetic field is applied in experiments with real encapsulated magnetic microbubbles.

Figures

Figures reproduced from arXiv: 2502.00863 by Arun Krishna B. J., Ganesh Tamadapu.

Figure 1
Figure 1. Figure 1: (a) Deformed and (b) undeformed configurations of a bubble surface. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Two coils placed symmetrically about e1−e2 plane carrying currents in opposite direction as har = H2 r 2a (3 cos2 θ − 1), haθ = H2 r 2a (−3 sin θ cos θ). (2.32) In-spherical co-ordinates the magnetic forces become, F mag r = (µoχ∇ha · ha)r = µoχ  har ∂har ∂r + haθ r ∂har ∂θ − h 2 aθ r  = µoχ H2 2 r 4a 2 (3 cos2 θ + 1), F mag θ = (µoχ∇ha · ha)θ = µoχ  har ∂haθ ∂r + haθ r ∂haθ ∂θ + harhaθ r  = µoχ H2 2 r… view at source ↗
Figure 3
Figure 3. Figure 3: Convergence of stability diagram with number of modes and results of linear [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Variation of natural frequency with mode numbers plot [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Variation of natural frequency with mode numbers for shear modulus [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Plots for an undeformed bubble radius R0 = 10 µm with different values of surface encapsulation shear modulus c1. (a) Variation of natural frequency with mode numbers. (b) Critical excitation pressure (ϵp) versus driving frequency (ωd) curves representing the stability boundaries. occurs because, in the absence of initial disturbances in the shape modes, the a2 mode is initially dominant. The amplitude of … view at source ↗
Figure 7
Figure 7. Figure 7: Variation of radial mode amplitude R(t) and second mode amplitude a2(t) with time for a bubble of initial radius R0 = 10 µm at various values of encapsulation susceptibility. The simulations are conducted at ωd = 1 MHz and ϵp = 0.5. contains c1 and R¨ contains −ϵp. So, lower pressure amplitude ϵp is enough to cause exponential blow up at higher c1. 5.2. Time-series of mode shapes Radial and shape oscillati… view at source ↗
Figure 8
Figure 8. Figure 8: Variation of radial mode amplitude R(t) and second mode amplitude a2(t) with time for a bubble encapsulation of susceptibility χ = 1 at different values of undeformed bubble radius R0. The simulations are conducted at ωd = 1 MHz and ϵp = 0.5. 5.2.1. Variation with magnetic susceptibility Simulations are conducted with the following parameters: R0 = 10 µm, c1 = 0.1 N m−1 , c2 = 0.12 N m−1 , N1 = 1000 turns,… view at source ↗
read the original abstract

This paper presents a model for nonspherical oscillations of encapsulated bubbles coated with a polymer infused with magnetic particles, developed using membrane theory for thin weakly magnetic membranes. According to this theory, only the applied magnetic field significantly contributes to the Maxwell stress and membrane is under generalized plane stress. The study focuses on axisymmetric deformations of bubbles under symmetrically arranged magnetic coils. Non-spherical oscillations of the bubble are restricted to the linear regime, with the second mode dominating within the pressure range of the stability region. The pressure-frequency stability region is computationally determined, and its variation with different material properties and applied magnetic field is analyzed. The natural frequency of each mode is estimated using boundary layer approximation. Time-series analysis of the second mode amplitude reveals a significant oscillation amplitude relative to the bubble radius. Estimation using the model indicates that the interface magnetic susceptibility and initial bubble radius enhance the amplitude of second-mode oscillations. Computational findings suggest that the applied magnetic field does not influence the stability region for exponential stability.

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 / 1 minor

Summary. The manuscript develops a model for nonspherical linear oscillations of encapsulated magnetic microbubbles coated with polymer infused with magnetic particles, using membrane theory under generalized plane stress where only the applied magnetic field contributes to Maxwell stress. It restricts analysis to axisymmetric deformations and the linear regime (second mode dominant), computationally determines the pressure-frequency stability region for exponential stability via numerical continuation, estimates natural frequencies with a boundary-layer approximation, analyzes time series of second-mode amplitude, and reports that the applied magnetic field does not influence the stability region while interface magnetic susceptibility and initial bubble radius enhance oscillation amplitudes.

Significance. If the modeling assumptions and numerical results hold, the finding of magnetic-field independence for the stability region could be relevant for applications of magnetic microbubbles in fluid dynamics. The computational determination of stability boundaries and boundary-layer frequency estimates are positive technical elements, but the absence of shown derivations, error estimates, or validation against nonlinear simulations limits the strength of the contribution.

major comments (2)
  1. [Abstract] Abstract (computational findings paragraph): the conclusion that the applied magnetic field does not influence the stability region for exponential stability rests directly on the modeling choice that only the applied field enters the Maxwell stress under generalized plane stress; without a comparison to the full Maxwell stress (including particle-induced terms) or an alternative stress state, the eigenvalue problem for stability could shift and the reported independence is not shown to be robust.
  2. [Abstract] Abstract: central claims on the stability region, natural frequencies, and amplitude scaling lack any shown derivation steps, error estimates, or validation against full nonlinear simulations, so the accuracy of the numerical continuation and boundary-layer results cannot be assessed from the provided information.
minor comments (1)
  1. [Abstract] Abstract: grammatical phrasing 'and membrane is under generalized plane stress' should be 'and the membrane is under generalized plane stress'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the opportunity to respond to the referee's report. We address the major comments point by point below, providing clarifications and indicating revisions where appropriate.

read point-by-point responses
  1. Referee: [Abstract] Abstract (computational findings paragraph): the conclusion that the applied magnetic field does not influence the stability region for exponential stability rests directly on the modeling choice that only the applied field enters the Maxwell stress under generalized plane stress; without a comparison to the full Maxwell stress (including particle-induced terms) or an alternative stress state, the eigenvalue problem for stability could shift and the reported independence is not shown to be robust.

    Authors: We thank the referee for highlighting this important point. The independence result is indeed tied to our modeling assumption, as stated in the manuscript, that under generalized plane stress for a thin weakly magnetic membrane, only the applied field contributes significantly to the Maxwell stress. This is based on the physical setup with symmetrically arranged coils and the membrane properties. To strengthen the presentation, we will revise the abstract and add a dedicated discussion section explaining the rationale for this approximation and acknowledging that a full Maxwell stress treatment (including induced fields from the particles) might alter the stability boundaries. This will make the scope of the claim clearer without changing the core analysis. revision: partial

  2. Referee: [Abstract] Abstract: central claims on the stability region, natural frequencies, and amplitude scaling lack any shown derivation steps, error estimates, or validation against full nonlinear simulations, so the accuracy of the numerical continuation and boundary-layer results cannot be assessed from the provided information.

    Authors: The derivations for the linear membrane model, the eigenvalue problem for stability via numerical continuation, and the boundary-layer approximation for natural frequencies are detailed in the methods and results sections of the full manuscript. However, we agree that the abstract could better indicate these. In revision, we will modify the abstract to briefly reference the numerical methods used and include a statement on the linear regime limitations. Error estimates for the boundary-layer approximation can be added in the text. Full validation against nonlinear simulations is beyond the scope of this linear analysis paper but will be noted as future work. revision: partial

Circularity Check

0 steps flagged

No circularity: stability region obtained by direct numerical continuation under explicit modeling assumptions.

full rationale

The paper states its core modeling choices (membrane theory under generalized plane stress with only applied field contributing to Maxwell stress) upfront and then computes the pressure-frequency stability region via numerical continuation of the linearized system. Natural frequencies are estimated separately via boundary-layer approximation. No quoted step shows a 'prediction' reducing to a fitted parameter by construction, no self-citation chain is load-bearing for the independence claim, and no ansatz is smuggled via prior work. The reported lack of magnetic-field influence on the stability region is therefore a computational output of the chosen equations rather than a definitional tautology.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on standard thin-shell assumptions plus two domain-specific choices whose independent support is not shown in the abstract.

free parameters (2)
  • interface magnetic susceptibility
    Fitted or chosen to match material; directly scales reported amplitude.
  • initial bubble radius
    Enters amplitude scaling; treated as a free geometric parameter.
axioms (2)
  • domain assumption Only the applied magnetic field contributes to Maxwell stress; internal fields neglected.
    Stated in the membrane-theory paragraph of the abstract.
  • domain assumption Deformations remain in the linear regime.
    Explicitly restricts analysis to linear oscillations.

pith-pipeline@v0.9.0 · 5697 in / 1393 out tokens · 42833 ms · 2026-05-23T04:18:24.412882+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

30 extracted references · 30 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    Acta Mechanica 191 (1), 1--19

    Barham, M , Steigmann, DJ , McElfresh, M & Rudd, RE 2007 Finite deformation of a pressurized magnetoelastic membrane in a stationary dipole field . Acta Mechanica 191 (1), 1--19

  4. [4]

    International Journal of Non-Linear Mechanics 47 (2), 185--196

    Barham, M , Steigmann, DJ & White, Dan 2012 Magnetoelasticity of highly deformable thin films: theory and simulation . International Journal of Non-Linear Mechanics 47 (2), 185--196

  5. [5]

    Journal of Applied Physics 66 (12), 5691--5701

    Boridy, Elie 1989 Magnetic fields generated by axially symmetric systems . Journal of Applied Physics 66 (12), 5691--5701

  6. [6]

    Theranostics 8 (2), 341

    Chertok, Beata & Langer, Robert 2018 Circulating magnetic microbubbles for localized real-time control of drug delivery by ultrasonography-guided magnetic targeting and ultrasound . Theranostics 8 (2), 341

  7. [7]

    The Journal of the Acoustical Society of America 155 (4), 2445--2459

    Dash, Nehal & Tamadapu, Ganesh 2024 Nonspherical oscillations of an encapsulated microbubble with interface energy under the acoustic field . The Journal of the Acoustical Society of America 155 (4), 2445--2459

  8. [8]

    Applied Acoustics 216 , 109792

    Du, Zhiwei , Pan, Ruiqi , Li, Fan , Mo, Runyang & Wang, Chenghui 2024 Local control of magnetic microbubble behavior by magnetic field and pulsed ultrasound . Applied Acoustics 216 , 109792

  9. [9]

    Journal of Fluid Mechanics 857 , 681--703

    Gu \'e dra, Matthieu & Inserra, Claude 2018 Bubble shape oscillations of finite amplitude . Journal of Fluid Mechanics 857 , 681--703

  10. [10]

    Physical Review E 94 (5), 053115

    Gu \'e dra, Matthieu , Inserra, Claude , Mauger, Cyril & Gilles, Bruno 2016 Experimental evidence of nonlinear mode coupling between spherical and nonspherical oscillations of microbubbles . Physical Review E 94 (5), 053115

  11. [11]

    Physics of Fluids 11 (6), 1309--1317

    Hao, Yue & Prosperetti, Andrea 1999 The effect of viscosity on the spherical stability of oscillating gas bubbles . Physics of Fluids 11 (6), 1309--1317

  12. [12]

    Journal of Fluid Mechanics 844 , 567--596

    Liu, Yunqiao , He, Dongdong , Gong, Xiaobo & Huang, Huaxiong 2018 Deformation of an encapsulated bubble in steady and oscillatory electric fields . Journal of Fluid Mechanics 844 , 567--596

  13. [13]

    Journal of Fluid Mechanics 691 , 315--340

    Liu, Yunqiao , Sugiyama, Kazuyasu , Takagi, Shu & Matsumoto, Yoichiro 2012 Surface instability of an encapsulated bubble induced by an ultrasonic pressure wave . Journal of Fluid Mechanics 691 , 315--340

  14. [14]

    Ultrasound in Medicine & Biology 38 (5), 864--875

    Mulvana, Helen , Eckersley, Robert J , Tang, Meng-Xing , Pankhurst, Quentin & Stride, Eleanor 2012 Theoretical and experimental characterisation of magnetic microbubbles . Ultrasound in Medicine & Biology 38 (5), 864--875

  15. [15]

    Interface Focus 5 (5), 20150001

    Owen, Joshua , Rademeyer, Paul , Chung, Daniel , Cheng, Qian , Holroyd, David , Coussios, Constantin , Friend, Peter , Pankhurst, Quentin A & Stride, Eleanor 2015 Magnetic targeting of microbubbles against physiologically relevant flow conditions . Interface Focus 5 (5), 20150001

  16. [16]

    Journal of Fluid Mechanics 440 , 269--291

    Pozrikidis, C 2001 Effect of membrane bending stiffness on the deformation of capsules in simple shear flow . Journal of Fluid Mechanics 440 , 269--291

  17. [17]

    Quarterly of Applied Mathematics 34 (4), 339--352

    Prosperetti, Andrea 1977 Viscous effects on perturbed spherical flows . Quarterly of Applied Mathematics 34 (4), 339--352

  18. [18]

    Royal Society Open Science 3 (8), 160063

    Sciallero, Claudia , Balbi, Luca , Paradossi, Gaio & Trucco, Andrea 2016 Magnetic resonance and ultrasound contrast imaging of polymer-shelled microbubbles loaded with iron oxide nanoparticles . Royal Society Open Science 3 (8), 160063

  19. [19]

    Physics of Fluids 18 (7), 072104

    Shaw, SJ 2006 Translation and oscillation of a bubble under axisymmetric deformation . Physics of Fluids 18 (7), 072104

  20. [20]

    Physics of Fluids 21 (2), 022104

    Shaw, SJ 2009 The stability of a bubble in a weakly viscous liquid subject to an acoustic traveling wave . Physics of Fluids 21 (2), 022104

  21. [21]

    Physics of Fluids 29 (12)

    Shaw, SJ 2017 Nonspherical sub-millimeter gas bubble oscillations: Parametric forcing and nonlinear shape mode coupling . Physics of Fluids 29 (12)

  22. [22]

    Journal of Engineering Mathematics 65 , 291--310

    Shaw, SJ , Spelt, PDM & Matar, OK 2009 Electrically induced bubble deformation, translation and collapse . Journal of Engineering Mathematics 65 , 291--310

  23. [23]

    Nonlinear Dynamics and Systems Theory 11 (2), 183--198

    Slane, Jean & Tragesser, Steven 2011 Analysis of periodic nonautonomous inhomogeneous systems . Nonlinear Dynamics and Systems Theory 11 (2), 183--198

  24. [24]

    International Journal of Non-Linear Mechanics 39 (7), 1193--1216

    Steigmann, DJ 2004 Equilibrium theory for magnetic elastomers and magnetoelastic membranes . International Journal of Non-Linear Mechanics 39 (7), 1193--1216

  25. [25]

    Ultrasound in Medicine & Biology 35 (5), 861--868

    Stride, Eleanor , Porter, Colin , Prieto, Ana Garcia & Pankhurst, Quentin 2009 Enhancement of microbubble mediated gene delivery by simultaneous exposure to ultrasonic and magnetic fields . Ultrasound in Medicine & Biology 35 (5), 861--868

  26. [26]

    Physics of Fluids 23 (1)

    Tsiglifis, Kostas & Pelekasis, Nikos A 2011 Parametric stability and dynamic buckling of an encapsulated microbubble subject to acoustic disturbances . Physics of Fluids 23 (1)

  27. [27]

    Chinese Physics B 31 (3), 034302

    Zhao, Lixia , Shi, Huimin , Bello, Isaac , Hu, Jing , Wang, Chenghui & Mo, Runyang 2022 Nonlinear oscillation characteristics of magnetic microbubbles under acoustic and magnetic fields . Chinese Physics B 31 (3), 034302

  28. [28]

    @esa (Ref

    \@ifclassloaded aguplus natbib The aguplus class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later remove the command natbib from the document \@ifclassloaded nlinproc natbib The nlinproc class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later r...

  29. [29]

    @stdbsttrue NAT@ctr \@lbibitem[ NAT@ctr ] \@lbibitem[#1]#2 \@extra@b@citeb \@ifundefined br@#2\@extra@b@citeb \@namedef br@#2 \@nameuse br@#2\@extra@b@citeb \@ifundefined b@#2\@extra@b@citeb @num @parse #2 [ @natanchorstart #2\@extra@b@citeb \@biblabel @num @natanchorend] @ifcmd#1(@)(@)\@nil #2 @lbibitem\@undefined @lbibitem\@lbibitem \@lbibitem[#1]#2 @lb...

  30. [30]

    " @ +G |!ϯ w x2 -]qT[C8-⨶[ hڦ (B`0B@m C5үh_ @Dڙ jCJ A58Anƨi 6PO⧶ :l S |B<mJ zw =5*Zq @0@z O H x 1BL

    @open @close @open @close and [1] URL: #1 \@ifundefined chapter * \@mkboth \@ifundefined NAT@sectionbib * \@mkboth * \@mkboth\@gobbletwo \@ifclassloaded amsart * \@ifclassloaded amsbook * \@ifundefined bib@heading @heading NAT@ctr thebibliography [1] @ \@biblabel NAT@ctr \@bibsetup #1 NAT@ctr 0 @openbib .11em \@plus.33em \@minus.07em 4000 4000 `\.=1000 \@...