Compressibility, sound and model limits

✓ Solved (inherited)

✓ 已完成（继承）

Completion credit carries forward from your earlier accepted answers. Covered work: wall Mach number, acoustic crossing time, and incompressible-model limitations.

完成记录沿用此前已接受的作答。已涵盖内容：泡壁马赫数、声传播时间及不可压缩模型的局限。

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

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

7. Compressibility, sound, and model limits 7.1 Two different small parameters 7.2 Keller–Miksis with constant ambient pressure 7.3 Solving the implicit acceleration correctly 7.4 Gilmore's enthalpy formulation 7.5 What the far-field pressure means 7.6 Stop before you manufacture a result

7. Compressibility, sound, and model limits

可压缩性、声辐射与模型极限

7.1 Two different small parameters

两个不同的小参数

The wall Mach number $M=|\dot R|/c$ measures how fast the interface moves relative to sound. For harmonic forcing, $kR=\omega R/c$ compares bubble size with acoustic wavelength. A low wall speed does not automatically make a very high-frequency forcing field spatially uniform. Likewise, a bubble may spend most of its life at small $M$ but leave that regime during its final collapse. [S413, S409]

界面Mach数 $M=|\dot R|/c$ 衡量界面运动相对于声速的快慢。对简谐激励， $kR=\omega R/c$ 比较气泡尺寸与声波波长。界面速度较低，并不自动保证高频激励场在空间上均匀。同样，一个气泡可能在大部分寿命中保持较小 $M$ ，却在塌缩末期离开这一范围。[S413, S409]

A radial compressible equation can include outgoing acoustic energy while still treating the bubble as spherical. It does not thereby become a three-dimensional shock-resolving simulation. Model dimensionality, thermodynamic closure, and compressibility order are separate choices.

可压缩径向方程可以在保持球形假设的同时计入向外传播的声能，但这并不使它成为三维激波解析仿真。模型维度、热力学闭合和可压缩性近似阶次，是彼此不同的选择。

7.2 Keller–Miksis with constant ambient pressure

恒定环境压力下的Keller–Miksis方程

For the numerical laboratory we deliberately use a constant far-field pressure $p_0$ . One common first-order compressible form is

在数值实验中，我们有意采用恒定远场压力 $p_0$ 。一种常见的一阶可压缩形式为

\left(1-\frac{\dot R}{c}\right)R\ddot R +\frac32\left(1-\frac{\dot R}{3c}\right)\dot R^2 =\left(1+\frac{\dot R}{c}\right)\frac{p_L-p_0}{\rho} +\frac{R}{\rho c}\frac{dp_L}{dt}. \tag{7.1}

7.2 Keller–Miksis with constant ambient pressure / 恒定环境压力下的Keller–Miksis方程; Ṙ: radial wall velocity; c: sound speed / heat capacity; R: instantaneous wall radius; R̈: radial wall acceleration; p_L: liquid pressure next to wall; p_0: reference ambient pressure; ρ: liquid / layer / droplet density; t: physical time

Here $p_L=p_g+p_v-2\sigma/R-4\mu\dot R/R$ is the liquid-side interfacial pressure, not just gas pressure. As $c\to\infty$ , equation (7.1) reduces to the incompressible radial equation. This limiting check is implemented in the supplied code. [S413]

这里 $p_L=p_g+p_v-2\sigma/R-4\mu\dot R/R$ 是界面液侧压力，而不只是气体压力。当 $c\to\infty$ 时，式(7.1)退化为不可压缩径向方程。配套代码实际执行了这一极限检查。[S413]

Published acoustically forced forms depend on the convention for the imposed incident pressure, its evaluation time, and retained asymptotic terms. The original Keller–Miksis derivation explicitly contains an incident-field timing convention. Do not take a forcing term from one convention and a pressure derivative from another. The laboratory avoids that ambiguity by implementing free ring-down at constant $p_0$ ; extending it to a measured acoustic waveform requires a deliberate derivation. [S413]

文献中受声激励形式的差异，与入射压力的定义、求值时刻及保留的渐近项有关。Keller–Miksis原始推导明确包含入射场的时间约定。不能从一种约定中取激励项，再从另一种约定中取压力导数。配套实验通过实现恒定 $p_0$ 下的自由衰减振荡来避免这种混用；要扩展至实测声压波形，必须重新明确推导约定。[S413]

7.3 Solving the implicit acceleration correctly

正确求解隐式加速度

Let $U=\dot R$ . For constant $p_v$ , constant $\sigma$ , and the polytropic gas law, differentiation gives

令 $U=\dot R$ 。当 $p_v$ 和 $\sigma$ 恒定，且气体满足多方定律时，求导得到

\frac{dp_L}{dt}=-\frac{3\kappa p_gU}{R} +\frac{2\sigma U}{R^2}+\frac{4\mu U^2}{R^2} -\frac{4\mu}{R}\dot U. \tag{7.2}

7.3 Solving the implicit acceleration correctly / 正确求解隐式加速度; p_L: liquid pressure next to wall; t: physical time; κ: gas polytropic exponent / curvature; p_g: noncondensable gas pressure; U: velocity / velocity scale; R: instantaneous wall radius; σ: interface tension / solid stress; μ: dynamic viscosity / absorption coefficient; U̇: velocity / velocity scale

The derivative contains the unknown acceleration. Omitting the last term is not the same model. Define $D_p$ as the first three terms in (7.2); then the explicit acceleration used in the code is

该导数中包含未知加速度。漏掉最后一项，就不再是同一个模型。将式(7.2)前三项记作 $D_p$ ，配套代码采用的显式加速度为

\dot U=\frac{ (1+U/c)(p_L-p_0)/\rho+RD_p/(\rho c) -\tfrac32(1-U/(3c))U^2} {R(1-U/c)+4\mu/(\rho c)}. \tag{7.3}

7.3 Solving the implicit acceleration correctly / 正确求解隐式加速度; U̇: velocity / velocity scale; U: velocity / velocity scale; c: sound speed / heat capacity; p_L: liquid pressure next to wall; p_0: reference ambient pressure; ρ: liquid / layer / droplet density; R: instantaneous wall radius; D_p: bending rigidity / dose / residual derivative; μ: dynamic viscosity / absorption coefficient

This is a useful example of a general computational principle: when a constitutive quantity depends on velocity, its time derivative may contain acceleration. Symbolically rearrange the equation before passing it to an explicit right-hand-side function.

这说明了一个普遍计算原则：当某个本构量依赖速度时，其时间导数可能包含加速度。在将方程交给显式右端函数之前，应先完成符号整理。

Free ring-down from the same expanded initial radius under Rayleigh–Plesset and Keller–Miksis models. The difference includes radiation-related compressibility effects; the curves are teaching simulations. / 相同膨胀初态下Rayleigh–Plesset与Keller–Miksis模型的自由衰减响应。差异包含与辐射相关的可压缩效应；图中为教学仿真。 — Free ring-down from the same expanded initial radius under Rayleigh–Plesset and Keller–Miksis models. The difference includes radiation-related compressibility effects; the curves are teaching simulations.
相同膨胀初态下Rayleigh–Plesset与Keller–Miksis模型的自由衰减响应。差异包含与辐射相关的可压缩效应；图中为教学仿真。

7.4 Gilmore's enthalpy formulation

Gilmore的焓形式

For stronger compressibility, a liquid equation of state can replace the constant-density pressure difference with an enthalpy difference. Define

在更强可压缩条件下，可通过液体状态方程，用焓差替代恒密度压力差。定义

H=\int_{p_\infty}^{p_L}\frac{dp}{\rho(p)},\qquad C^2=\left.\frac{dp}{d\rho}\right|_{p_L}. \tag{7.4}

7.4 Gilmore's enthalpy formulation / Gilmore的焓形式; H: enthalpy / absorbed dose / liquid depth; p_L: liquid pressure next to wall; p: pressure at the indicated location; ρ: liquid / layer / droplet density

A standard Gilmore form is

一种标准Gilmore形式为

\left(1-\frac{U}{C}\right)R\dot U +\frac32\left(1-\frac{U}{3C}\right)U^2 =\left(1+\frac{U}{C}\right)H +\left(1-\frac{U}{C}\right)\frac{R}{C}\dot H. \tag{7.5}

7.4 Gilmore's enthalpy formulation / Gilmore的焓形式; U: velocity / velocity scale; C: gas-content constant / damping; R: instantaneous wall radius; U̇: velocity / velocity scale; H: enthalpy / absorbed dose / liquid depth; Ḣ: enthalpy / absorbed dose / liquid depth

For a Tait-type liquid relation $(p+B)/(p_\infty+B)=(\rho/\rho_\infty)^n$ , with constant reference ambient state,

对于Tait类液体关系 $(p+B)/(p_\infty+B)=(\rho/\rho_\infty)^n$ ，在参考远场状态恒定时，

H=\frac{n(p_\infty+B)}{(n-1)\rho_\infty} \left[\left(\frac{p_L+B}{p_\infty+B}\right)^{(n-1)/n}-1\right], \qquad C^2=c_\infty^2+(n-1)H. \tag{7.6}

7.4 Gilmore's enthalpy formulation / Gilmore的焓形式; H: enthalpy / absorbed dose / liquid depth; n: mode / bubble number density; p∞: ambient liquid pressure; B: kinetic pressure coefficient; ρ∞: liquid / layer / droplet density; p_L: liquid pressure next to wall; C: gas-content constant / damping; c∞: sound speed / heat capacity

The parameters $B,n$ and the valid pressure range must be appropriate to the liquid. The extensive laser-cavitation analysis in your collection uses an extended Gilmore framework and carefully treats initialization and energy partition. The existence of a more sophisticated radial model does not remove uncertainties in deposited energy, gas content, or phase change. [S409]

参数 $B,n$ 及其适用压力范围必须与具体液体相符。你提供的长篇激光空化分析采用扩展Gilmore框架，并仔细处理初始化与能量分配。更复杂的径向模型并不会自动消除沉积能量、气体含量或相变过程中的不确定性。[S409]

7.5 What the far-field pressure means

远场压力究竟表示什么

For a compact spherical source in a weakly nonlinear or linear far-field propagation approximation, the monopole pressure is

在紧致球形声源以及弱非线性或线性远场传播近似下，单极子声压为

p'(r,t)\simeq\frac{\rho}{4\pi r}\ddot V(t-r/c) =\frac{\rho}{r}\left[R^2\ddot R+2R\dot R^2\right]_{t-r/c}. \tag{7.7}

7.5 What the far-field pressure means / 远场压力究竟表示什么; p': pressure at the indicated location; r: liquid observation radius / embryo radius; t: physical time; ρ: liquid / layer / droplet density; π: circle geometry constant; V̈: bubble or observation volume; c: sound speed / heat capacity; R: instantaneous wall radius; R̈: radial wall acceleration; Ṙ: radial wall velocity

The observation distance $r$ is measured from the bubble center and must be in a suitable radiation region. A local wall-impact pressure cannot be obtained by putting a tiny $r$ into this far-field formula. Nor does equation (7.7) resolve shock steepening, attenuation, reflection, sensor averaging, or nonspherical jets.

观测距离 $r$ 从气泡中心量起，且必须位于适当的辐射区域。不能通过在远场公式中代入极小 $r$ 来获得局部壁面冲击压力。式(7.7)也不解析激波陡化、衰减、反射、传感器平均效应或非球形射流。

For a freely outgoing, approximately linear spherical wave, an acoustic-energy estimate is

对于自由向外传播、近似线性的球面波，声能估计为

E_{\mathrm{ac}}\simeq\frac{4\pi r^2}{\rho c}\int p'(r,t)^2\,dt. \tag{7.8}

7.5 What the far-field pressure means / 远场压力究竟表示什么; E_ac: energy / elastic modulus as defined; π: circle geometry constant; r: liquid observation radius / embryo radius; ρ: liquid / layer / droplet density; c: sound speed / heat capacity; p': pressure at the indicated location; t: physical time

Use this only with its propagation assumptions. An arbitrary near-field pressure trace may include reactive rather than radiated energy. Your laser-bubble source [S409] is useful precisely because it distinguishes liquid compression, acoustic emission, and other energy channels rather than equating them all to $pV$ .

使用该式时必须保留其传播假设。任意近场压力历程可能包含反应性储能，而不是纯辐射能量。你提供的激光气泡文献[S409]的重要价值，正是在于区分液体压缩、声发射及其他能量通道，而不是把它们都等同于 $pV$ 。

Approximate monopole pressure calculated from the teaching Keller–Miksis trajectory at 10 mm. This is not a resolved collapse shock or wall-impact prediction. / 由教学Keller–Miksis轨迹计算的10毫米处近似单极子声压。它不是解析得到的塌缩激波，也不是壁面冲击预测。 — Approximate monopole pressure calculated from the teaching Keller–Miksis trajectory at 10 mm. This is not a resolved collapse shock or wall-impact prediction.
由教学Keller–Miksis轨迹计算的10毫米处近似单极子声压。它不是解析得到的塌缩激波，也不是壁面冲击预测。

7.6 Stop before you manufacture a result

在制造伪结果之前停止

A solver can integrate beyond the validity of its model. The supplied laboratory uses a conservative teaching guard at $|U|/c=0.10$ and a finite minimum radius. The value 0.10 is a chosen warning threshold, not a universal accuracy boundary for Keller–Miksis. A stopped calculation is preferable to reporting a spurious pressure spike generated by an invalid extrapolation.

求解器可能在模型适用范围之外继续积分。配套数值实验在 $|U|/c=0.10$ 处设置保守教学保护，并设置有限最小半径。0.10是人为选择的警戒阈值，不是Keller–Miksis模型普适的精度分界线。与其报告无效外推产生的虚假压力尖峰，不如及时停止计算。

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.

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

A 10-µm bubble has wall speed 150 m/s and a 1-MHz signal component. For c = 1482 m/s compute Mw and kR. Are these interchangeable checks?
10 µm 气泡壁面速率为 150 m/s，信号含 1 MHz 分量。取 c = 1482 m/s，计算 Mw 与 kR。这两项检查可以互换吗？
Reference answer / 参考答案
$M_w=|\dot R|/c,\qquad kR=\omega R/c$
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.
图示标出物理体系，图例列出公式中的符号。若某字母有多种作用，应以邻近推导中的定义为准。
Mw = 0.1012, kR = 0.0424. Mw compares wall speed with sound speed; kR compares source size with wavelength. Small kR does not rescue a fast collapse; the book’s 0.10 stop rule is a conservative teaching guard, not a universal accuracy threshold.
Mw = 0.1012，kR = 0.0424。Mw 比较壁面速率和声速；kR 比较源尺寸和波长。kR 小不能挽救快速塌缩；教材 0.10 停止规则是保守教学保护条件，不是普适精度阈值。
Starting from pL = pg + pv − 2σ/R − 4μṘ/R, isolate the hidden acceleration in Keller–Miksis. At R = 10 µm, U = −10 m/s, calculate its denominator.
从 pL = pg + pv − 2σ/R − 4μṘ/R 出发，分离 Keller–Miksis 中隐藏的加速度。在 R = 10 µm、U = −10 m/s 时计算其分母。
Reference answer / 参考答案
$\dot p_L=D_p-4\mu\ddot R/R,\quad D_p=-3\kappa p_g\dot R/R+2\sigma\dot R/R^2+4\mu\dot R^2/R^2,\quad \mathcal D=R(1-\dot R/c)+4\mu/(\rho c)$
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.
图示标出物理体系，图例列出公式中的符号。若某字母有多种作用，应以邻近推导中的定义为准。
Using μ = 0.001, ρ = 998, c = 1482, the denominator is 1.00702×10⁻⁵ m, including 2.705×10⁻⁹ m from viscosity. Move the −4μR̈/R derivative term to the left before dividing; leaving it in a supposedly known RHS is circular.
使用 μ = 0.001、ρ = 998、c = 1482，分母为 1.00702×10⁻⁵ m，其中黏性贡献 2.705×10⁻⁹ m。除法前应将导数中的 −4μR̈/R 移至左端；把它留在所谓已知右端会造成循环定义。
For volume acceleration V̈ = 0.001 m³/s², estimate linear monopole pressure at r = 1 mm and 2 mm. Why is it not a jet-impact pressure?
体积加速度 V̈ = 0.001 m³/s²，在 r = 1 mm 和 2 mm 估算线性单极声压。它为何不是射流冲击压力？
Reference answer / 参考答案
$p'(r,t)=\rho\ddot V(t-r/c)/(4\pi r)$
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.
图示标出物理体系，图例列出公式中的符号。若某字母有多种作用，应以邻近推导中的定义为准。
At appropriately retarded times the amplitudes are 79.42 Pa and 39.71 Pa. The 1/r relation is a far-field linear radiation result; local impact, shocks, boundary reflection and compressible near flow require other models.
在相应迟滞时刻，幅值为 79.42 Pa 和 39.71 Pa。1/r 关系是远场线性辐射结果；局部冲击、激波、边界反射及可压缩近场需要其他模型。

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.

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