Bubble oscillations as vibration theory

✓ Solved

✓ 已完成

Your submitted answer was checked and accepted. Confirmed work: equilibrium gas pressure, small-amplitude radial oscillations, static response, viscous and total damping, and shape stability limits.

本章提交的作答已核对并接受。已确认内容：平衡气体压力、小振幅径向振荡、静态响应、黏性阻尼与总阻尼，以及形状稳定性的局限。

Source: supplied generic textbook, Chapter 6. The complete bilingual chapter follows, with its derivations and qualifications. Textbook source (Markdown) · Original Markdown.

来源：所给通用教材第 6。下方完整收录双语章节及推导和限定。教材原文（Markdown） · 原始 Markdown。

6. Bubble oscillations as vibration theory 6.1 Linearization around equilibrium 6.2 Harmonic response and phase 6.3 Three damping mechanisms 6.4 What becomes nonlinear first

6. Bubble oscillations as vibration theory

用振动理论理解气泡振荡

6.1 Linearization around equilibrium

平衡态附近的线性化

Let $R=R_0(1+x)$ with $|x|\ll1$ , and write $p_\infty=p_0+p_a(t)$ . Expand the gas law as $p_g\simeq p_{g0}(1-3\kappa x)$ and $2\sigma/R\simeq(2\sigma/R_0)(1-x)$ . Products of small perturbations are second order. The radial equation reduces to

令 $R=R_0(1+x)$ ，其中 $|x|\ll1$ ，并设 $p_\infty=p_0+p_a(t)$ 。展开气体定律 $p_g\simeq p_{g0}(1-3\kappa x)$ 和 $2\sigma/R\simeq(2\sigma/R_0)(1-x)$ ，忽略小扰动乘积等二阶项，径向方程化为

\ddot x+2\beta_\mu\dot x+\omega_0^2x =-\frac{p_a(t)}{\rho R_0^2}, \tag{6.1}

6.1 Linearization around equilibrium / 平衡态附近的线性化; ẍ: position / dimensionless radius; β_μ: damping / logistic coefficient; ẋ: position / dimensionless radius; ω_0: linear radial angular frequency; x: position / dimensionless radius; p_a: pressure at the indicated location; t: physical time; ρ: liquid / layer / droplet density; R_0: reference equilibrium radius

\omega_0^2=\frac{3\kappa p_{g0}-2\sigma/R_0}{\rho R_0^2} =\frac{3\kappa(p_0-p_v)+(3\kappa-1)2\sigma/R_0}{\rho R_0^2}, \qquad \beta_\mu=\frac{2\mu}{\rho R_0^2}. \tag{6.2}

6.1 Linearization around equilibrium / 平衡态附近的线性化; ω_0: linear radial angular frequency; κ: gas polytropic exponent / curvature; p_g0: reference gas pressure; σ: interface tension / solid stress; R_0: reference equilibrium radius; ρ: liquid / layer / droplet density; p_0: reference ambient pressure; p_v: vapor partial pressure; β_μ: damping / logistic coefficient; μ: dynamic viscosity / absorption coefficient

This is the familiar damped oscillator, but its stiffness is supplied mainly by gas compression and capillarity, and its inertia comes from surrounding liquid. A positive acoustic pressure compresses the bubble, which explains the minus sign in the forcing term. The restoring coefficient must be positive for stable small oscillations.

这就是熟悉的阻尼振子，但其刚度主要来自气体压缩和毛细作用，惯性则来自周围液体。正声压会压缩气泡，因此激励项前有负号。要产生稳定小振动，回复系数必须为正。

If surface tension and vapor pressure are negligible compared with $p_0$ , then

当表面张力和蒸气压相对于 $p_0$ 可以忽略时，

f_0\simeq\frac{1}{2\pi R_0}\sqrt{\frac{3\kappa p_0}{\rho}}. \tag{6.3}

6.1 Linearization around equilibrium / 平衡态附近的线性化; f_0: frequency / predictor / false detections; π: circle geometry constant; R_0: reference equilibrium radius; κ: gas polytropic exponent / curvature; p_0: reference ambient pressure; ρ: liquid / layer / droplet density

This is the familiar Minnaert-type scaling. The inverse-radius dependence is more important than memorizing a numerical coefficient for one liquid and one ambient pressure. For micron-scale bubbles, ignoring surface tension can produce a noticeable error.

这就是常见的Minnaert类尺度关系。与其记忆某一种液体、某一个环境压力下的数值系数，不如理解频率与半径成反比的规律。对于微米级气泡，忽略表面张力可能造成明显误差。

Worked example. For the code's teaching values and $R_0=10\ \mu\mathrm m$ , $p_{g0}=113386\ \mathrm{Pa}$ , $f_0=342.37\ \mathrm{kHz}$ , and $T_0=2.92084\ \mu\mathrm s$ . The viscous decay rate is $\beta_\mu=2.0040\times10^4\ \mathrm{s^{-1}}$ . These values are obtained by direct substitution into (6.2).

计算示例。采用配套代码中的教学参数，且 $R_0=10\ \mu\mathrm m$ ，可得 $p_{g0}=113386\ \mathrm{Pa}$ 、 $f_0=342.37\ \mathrm{kHz}$ 和 $T_0=2.92084\ \mu\mathrm s$ 。黏性衰减率为 $\beta_\mu=2.0040\times10^4\ \mathrm{s^{-1}}$ 。这些数值由式(6.2)直接代入计算得到。

6.2 Harmonic response and phase

简谐响应与相位

For $p_a=P_a\cos\omega t$ and total linear damping rate $\beta$ , the steady amplitude is

对于 $p_a=P_a\cos\omega t$ 和总线性阻尼率 $\beta$ ，稳态振幅为

|\hat x|=\frac{P_a/(\rho R_0^2)} {\sqrt{(\omega_0^2-\omega^2)^2+(2\beta\omega)^2}}. \tag{6.4}

6.2 Harmonic response and phase / 简谐响应与相位; x^: position / dimensionless radius; P_a: event probability; ρ: liquid / layer / droplet density; R_0: reference equilibrium radius; ω_0: linear radial angular frequency; ω: angular frequency; β: damping / logistic coefficient

The phase changes across resonance. In a pressure-driven bubble, remember the sign reversal between pressure and expansive force when interpreting phase. The frequency of maximum displacement response is not exactly $\omega_0$ when damping is appreciable. Nor is a nonlinear bubble's resonance curve guaranteed to match the linear one. [S413]

相位会跨越共振区发生变化。对于压力驱动气泡，解释相位时必须记住压力与膨胀驱动力之间的反号关系。当阻尼不可忽略时，最大位移响应对应的频率不再严格等于 $\omega_0$ 。非线性气泡的共振曲线也不一定与线性模型一致。[S413]

Normalized linear compliance for three illustrative damping ratios. These curves are oscillator calculations, not measured bubble spectra. / 三个示意阻尼比下的归一化线性柔度曲线。这些是振子计算，而非实测气泡频谱。 — Normalized linear compliance for three illustrative damping ratios. These curves are oscillator calculations, not measured bubble spectra.
三个示意阻尼比下的归一化线性柔度曲线。这些是振子计算，而非实测气泡频谱。

6.3 Three damping mechanisms

三种阻尼机制

Viscous dissipation is explicit in the incompressible radial equation. Thermal damping arises because pressure and volume are not related by a perfectly reversible instantaneous law when heat exchange is finite-rate. Acoustic radiation carries energy into the far field and is absent from an incompressible model. Near resonance, a leading small-bubble radiation estimate is $\beta_{\mathrm{rad}}\simeq\omega_0^2R_0/(2c)$ ; the detailed frequency-dependent impedance is more general. [S413, S409]

黏性耗散直接出现在不可压缩径向方程中。有限速率热交换使压力和体积不再服从完全可逆、瞬时的关系，从而产生热阻尼。声辐射将能量输运至远场，而不可压缩模型不包含这一机制。在共振附近，小气泡的一个首阶声辐射估计为 $\beta_{\mathrm{rad}}\simeq\omega_0^2R_0/(2c)$ ；更一般的描述需要考虑随频率变化的阻抗。[S413, S409]

For the 10-micrometer example, this radiation estimate is about $1.56\times10^4\ \mathrm{s^{-1}}$ , comparable in order of magnitude to the viscous rate. Changing $\kappa$ to fit a decaying envelope is not a clean way to model radiation loss: it changes stiffness and thermodynamic assumptions at the same time.

对于10微米示例，这个声辐射估计约为 $1.56\times10^4\ \mathrm{s^{-1}}$ ，与黏性衰减率处于相近量级。通过修改 $\kappa$ 来拟合衰减包络，并不是表达声辐射损失的合理方式，因为这样会同时改变刚度和热力学假设。

6.4 What becomes nonlinear first

哪些效应首先变得非线性

Large oscillations sample the nonlinear gas pressure law, the $\dot R^2$ inertia term, and a changing effective mass. Expansion can be relatively slow and collapse very sharp. Harmonics, subharmonics, multiple responses, and sensitivity to initial conditions can arise in appropriate forcing regimes, but they are not inevitable in every experiment. Keller and Miksis show why radiation matters for large-amplitude oscillations rather than acting as an optional cosmetic correction. [S413]

大振幅振荡会显现气体压力定律的非线性、 $\dot R^2$ 惯性项，以及变化的有效质量。膨胀可能相对缓慢，而塌缩十分尖锐。在适当激励条件下，可能出现谐波、次谐波、多重响应和初值敏感性，但并不是每个实验都会如此。Keller和Miksis的研究说明，声辐射会影响大振幅振荡，而不是可有可无的装饰性修正。[S413]

Checkpoint. A Fourier spectrum containing high frequencies does not by itself prove a shock wave. A sharp but continuous nonlinear waveform also contains harmonics. Shock identification requires a compatible propagation model or suitably resolved measurements.

自检。频谱中出现高频成分，并不能单独证明存在激波。尖锐但连续的非线性波形也会含有谐波。识别激波需要相容的传播模型，或具有足够分辨率的测量。

Integrated calculation and model check

综合计算与模型检查

Use the declared parameter set in each part. The values are teaching scenarios unless specifically identified as inherited project inputs. Each reference panel follows its question.

使用每小题明确给出的参数。除明确标为沿用项目输入外，数值均为教学情景。各参考解答紧随题目。

For R0 = 10 µm and the book water properties, κ = 1.4, calculate pg0, ω0, f0 and the slow response to an extra ambient 1000 Pa.
R0 = 10 µm，使用教材水参数、κ = 1.4，计算 pg0、ω0、f0 及环境压力缓慢增加 1000 Pa 后的响应。
Reference answer / 参考答案
$p_{g0}=p_0-p_v+2\sigma/R_0,\quad \omega_0^2=(3\kappa p_{g0}-2\sigma/R_0)/(\rho R_0^2),\quad x_{\mathrm{static}}=-\delta p/(3\kappa p_{g0}-2\sigma/R_0)$
The diagram locates the physical system; its key lists the symbols in the formula. Read the definition in the surrounding derivation when a letter has different roles.
图示标出物理体系，图例列出公式中的符号。若某字母有多种作用，应以邻近推导中的定义为准。
pg0 = 113,386 Pa, ω0 = 2.151×10⁶ rad/s, f0 = 342.37 kHz. xstatic = −0.0021653, so ΔR = −21.65 nm. This stiffness includes capillarity; Minnaert’s gas-only limit discards it.
pg0 = 113,386 Pa，ω0 = 2.151×10⁶ rad/s，f0 = 342.37 kHz。xstatic = −0.0021653，ΔR = −21.65 nm。此刚度包含毛细作用；仅含气体的 Minnaert 极限舍弃该项。
Take total β = 4.0×10⁴ s⁻¹ and forcing amplitude pa = 100 Pa at ω = ω0. Calculate relative oscillation amplitude and interpret doubling viscosity.
设总 β = 4.0×10⁴ s⁻¹，在 ω = ω0 下驱动幅值 pa = 100 Pa，计算相对振幅并解释黏度加倍。
Reference answer / 参考答案
$|x|=\frac{p_a/(\rho R_0^2)}{\sqrt{(\omega_0^2-\omega^2)^2+(2\beta\omega)^2}},\qquad \beta_\mu=2\mu/(\rho R_0^2)$
The diagram locates the physical system; its key lists the symbols in the formula. Read the definition in the surrounding derivation when a letter has different roles.
图示标出物理体系，图例列出公式中的符号。若某字母有多种作用，应以邻近推导中的定义为准。
|x| = 0.00582, a 58.2-nm radius amplitude, compatible with this small-amplitude approximation. βμ = 2.004×10⁴ s⁻¹. Doubling μ doubles this damping contribution; it does not directly double the restoring stiffness or necessarily total β.
|x| = 0.00582，对应半径振幅 58.2 nm，符合本例小振幅近似。βμ = 2.004×10⁴ s⁻¹。黏度加倍使该阻尼分量加倍，不会直接使回复刚度加倍，也不一定使总 β 加倍。
Explain why the radial frequency does not determine shape stability, and why a resonance peak alone cannot identify every damping mechanism.
解释径向频率为何不能确定形状稳定性，以及仅凭共振峰为何不能辨识所有阻尼机制。
Reference answer / 参考答案
$\beta=\beta_\mu+\beta_{\mathrm{thermal}}+\beta_{\mathrm{radiation}}$
The diagram locates the physical system; its key lists the symbols in the formula. Read the definition in the surrounding derivation when a letter has different roles.
图示标出物理体系，图例列出公式中的符号。若某字母有多种作用，应以邻近推导中的定义为准。
A fitted linewidth gives an effective sum unless independent information separates the terms. Shape modes require their own perturbation equations. Large response can change gas stiffness, detune resonance, and destabilize spherical shape; the linear radial oscillator is then insufficient.
若无独立信息分离各项，拟合线宽只能得到有效总和。形状模态需要独立扰动方程。大响应可改变气体刚度、使共振失谐并破坏球形稳定性，此时线性径向振子不足。

Brief conceptual synthesis

简短概念综合

After using the reference panels, give a short explanation of the physical mechanism, the strongest assumption, and one observation that could falsify this chapter’s model. Preparation of later chapters does not require waiting for that reply.

核对参考解答后，简述物理机制、最强假设，以及一个可证伪本章模型的观测。后续章节编写无需等待该回答。