Core route · step 16/21 · source G16

主路线 · 第 16/21 步 · 来源 G16

Thirty original worked problems

三十道原始详解习题

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

来源:所给通用教材第 16。下方完整收录双语章节及推导和限定。教材原文(Markdown) · 原始 Markdown。

16. Problems with worked solutions16.1 Laplace pressure at three scales16.2 Where the liquid kinetic energy resides16.3 Track a liquid particle16.4 Diagnose a Bernoulli error16.5 An isothermal Blake threshold16.6 Probability depends on observation volume16.7 Long-time inertial growth speed16.8 Collapse-time scaling16.9 Speed at half maximum radius16.10 When the incompressible collapse becomes suspect16.11 Static compliance of a microbubble16.12 Viscosity versus stiffness16.13 The isothermal gas potential16.14 The hidden acceleration in Keller–Miksis16.15 Radial frequency versus shape frequency16.16 A uniform pressure impulse16.17 Equal impulse, zero residual oscillation16.18 A necessary adhesion-energy scale16.19 Droplet expansion by mass conservation16.20 Geometrical spreading of acoustic energy16.21 A false temperature prediction16.22 A mesh-convergence calculation16.23 A nonidentifiable gas parameterization16.24 Why more frames do not automatically fix acceleration16.25 Source-mechanism classification16.26 A terminology audit16.27 Detect an overconstrained interface16.28 Avoid double-counting fluid inertia16.29 Design a discriminating experiment16.30 Build a falsifiable capstone study

16. Problems with worked solutions

习题与详细解答

Attempt each problem before reading its solution. Numerical questions use the teaching parameter set unless otherwise stated. The problems progress from direct balances to model criticism; the last category is essential for research rather than optional commentary.

请先独立作答,再阅读解答。除非另有说明,数值题均采用本书教学参数。题目从直接平衡逐步过渡到模型批判;对科研而言,最后一类能力是必要内容,而不是可有可无的补充。

16.1 Laplace pressure at three scales

三个尺度下的Laplace压力

Problem. Find the capillary pressure of clean spherical bubbles with radii 1, 10, and 100 micrometers for σ=0.072 N m−1\sigma=0.072\ \mathrm{N\,m^{-1}}. At which scales is it small compared with approximately 100 kPa ambient pressure?

题目。当σ=0.072 N m−1\sigma=0.072\ \mathrm{N\,m^{-1}}时,求半径分别为1、10和100微米的洁净球形气泡毛细压力。在哪些尺度下,它相对于约100 kPa环境压力较小?

Reference answer · original problem 1 / 原题参考答案
pσ=2σ/Rp_\sigma=2\sigma/R
16.1 Laplace pressure at three scales / 三个尺度下的Laplace压力; p_σ: pressure at the indicated location; σ: interface tension / solid stress; R: instantaneous wall radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Equation (1.1) gives 144, 14.4, and 1.44 kPa. It is dominant or comparable at 1 micrometer, a noticeable correction at 10 micrometers, and a small correction at 100 micrometers. “Small bubble” is therefore not a purely geometric description; it changes the pressure balance.

解答。由式(1.1),结果为144、14.4和1.44 kPa。在1微米处,它占主导或至少与环境压力相当;在10微米处,它是明显修正;在100微米处,它较小。因此,“小气泡”不只是几何描述,还意味着压力平衡发生变化。

16.2 Where the liquid kinetic energy resides

液体动能分布在哪里

Problem. For spherical incompressible motion, what fraction of the total liquid kinetic energy lies outside r=aRr=aR, with a>1a>1?

题目。对于球对称不可压缩运动,半径r=aRr=aR之外的液体动能占总液体动能多少,其中a>1a>1?

Reference answer · original problem 2 / 原题参考答案
T=12ρ∫R∞u24πr2dr,u=R2R˙/r2T=\tfrac12\rho\int_R^\infty u^2 4\pi r^2dr,\quad u=R^2\dot R/r^2
16.2 Where the liquid kinetic energy resides / 液体动能分布在哪里; T: temperature / kinetic energy as defined; ρ: liquid / layer / droplet density; u: liquid or solid particle velocity; π: circle geometry constant; r: liquid observation radius / embryo radius; R: instantaneous wall radius; Ṙ: radial wall velocity

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Since u2r2∝r−2u^2r^2\propto r^{-2}, the energy outside a lower radius rbr_b is proportional to 1/rb1/r_b. Therefore Tr>aR/T=1/aT_{r>aR}/T=1/a. Half the energy lies outside 2R2R. A computational boundary only a few bubble radii away can alter the inertia substantially unless an appropriate exterior representation is used.

解答。由于u2r2∝r−2u^2r^2\propto r^{-2},从下限rbr_b向外积分的动能与1/rb1/r_b成正比,因此Tr>aR/T=1/aT_{r>aR}/T=1/a。一半动能位于2R2R之外。因此,如果计算边界只距气泡几个半径,且未采用适当外域处理,就可能显著改变惯性。

16.3 Track a liquid particle

跟踪液体质点

Problem. Derive the relation between a liquid particle's radial position r(t)r(t) and the bubble radius.

题目。推导液体质点径向位置r(t)r(t)与气泡半径之间的关系。

Reference answer · original problem 3 / 原题参考答案
r˙=u(r,t)=R2R˙/r2\dot r=u(r,t)=R^2\dot R/r^2
16.3 Track a liquid particle / 跟踪液体质点; ṙ: liquid observation radius / embryo radius; u: liquid or solid particle velocity; r: liquid observation radius / embryo radius; t: physical time; R: instantaneous wall radius; Ṙ: radial wall velocity

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Along a particle trajectory, dr/dt=R2R˙/r2dr/dt=R^2\dot R/r^2. Multiply by 3r23r^2 to obtain d(r3)/dt=d(R3)/dtd(r^3)/dt=d(R^3)/dt. Thus r(t)3−R(t)3=r(t0)3−R(t0)3r(t)^3-R(t)^3=r(t_0)^3-R(t_0)^3. The relation expresses conservation of the liquid volume between the bubble and that material spherical surface.

解答。沿质点轨迹有dr/dt=R2R˙/r2dr/dt=R^2\dot R/r^2,乘以3r23r^2得d(r3)/dt=d(R3)/dtd(r^3)/dt=d(R^3)/dt,因此r(t)3−R(t)3=r(t0)3−R(t0)3r(t)^3-R(t)^3=r(t_0)^3-R(t_0)^3。该关系表达了气泡与该物质球面之间液体体积守恒。

16.4 Diagnose a Bernoulli error

诊断Bernoulli推导错误

Problem. A student differentiates ϕ(R(t),t)=−RR˙\phi(R(t),t)=-R\dot R and substitutes −R˙2−RR¨-\dot R^2-R\ddot R as ∂tϕ\partial_t\phi. What is missing?

题目。某同学对ϕ(R(t),t)=−RR˙\phi(R(t),t)=-R\dot R求导,并把−R˙2−RR¨-\dot R^2-R\ddot R代作∂tϕ\partial_t\phi。遗漏了什么?

Reference answer · original problem 4 / 原题参考答案
ϕ=−R2R˙/r,dϕ(R(t),t)/dt=∂tϕ+R˙∂rϕ\phi=-R^2\dot R/r,\quad d\phi(R(t),t)/dt=\partial_t\phi+\dot R\partial_r\phi
16.4 Diagnose a Bernoulli error / 诊断Bernoulli推导错误; φ: velocity potential / mode shape; R: instantaneous wall radius; Ṙ: radial wall velocity; r: liquid observation radius / embryo radius; t: physical time

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. The chain rule gives dϕ(R(t),t)/dt=∂tϕ+R˙∂rϕd\phi(R(t),t)/dt=\partial_t\phi+\dot R\partial_r\phi. Since ∂rϕ∣R=R˙\partial_r\phi|_R=\dot R, subtract R˙2\dot R^2 from the total derivative. The correct partial derivative is −2R˙2−RR¨-2\dot R^2-R\ddot R, which produces the 3/23/2 inertia coefficient after adding u2/2u^2/2.

解答。链式法则给出dϕ(R(t),t)/dt=∂tϕ+R˙∂rϕd\phi(R(t),t)/dt=\partial_t\phi+\dot R\partial_r\phi。由于∂rϕ∣R=R˙\partial_r\phi|_R=\dot R,应从全导数中再减去R˙2\dot R^2。正确偏导为−2R˙2−RR¨-2\dot R^2-R\ddot R,加上u2/2u^2/2后即可得到3/23/2惯性系数。

16.5 An isothermal Blake threshold

等温Blake阈值

Problem. Find the turning point for a reference gas nucleus of radius R0=5 μmR_0=5\ \mu\mathrm m at p0=101325 Pap_0=101325\ \mathrm{Pa}, with pv=2339 Pap_v=2339\ \mathrm{Pa} and κ=1\kappa=1.

题目。对于p0=101325 Pap_0=101325\ \mathrm{Pa}下参考半径R0=5 μmR_0=5\ \mu\mathrm m的气核,取pv=2339 Pap_v=2339\ \mathrm{Pa}、κ=1\kappa=1,求平衡曲线转折点。

Reference answer · original problem 5 / 原题参考答案
pg0=p0−pv+2σ/R0,RB=3pg0R03/(2σ),pBlake=pv−4σ/(3RB)p_{g0}=p_0-p_v+2\sigma/R_0,\quad R_B=\sqrt{3p_{g0}R_0^3/(2\sigma)},\quad p_{\mathrm{Blake}}=p_v-4\sigma/(3R_B)
16.5 An isothermal Blake threshold / 等温Blake阈值; p_g0: reference gas pressure; p_0: reference ambient pressure; p_v: vapor partial pressure; σ: interface tension / solid stress; R_0: reference equilibrium radius; R_B: bubble wall radius; p_Blake: critical ambient pressure

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. First calculate pg0=p0−pv+2σ/R0=127786 Pap_{g0}=p_0-p_v+2\sigma/R_0=127786\ \mathrm{Pa}. Equation (2.5) gives RB=18.242 μmR_B=18.242\ \mu\mathrm m and pBlake=−2923.5 Pap_{\mathrm{Blake}}=-2923.5\ \mathrm{Pa} absolute. This negative absolute pressure is a tensile metastable-liquid value within the model, not a negative gas pressure and not a universal measured inception threshold.

解答。先算得pg0=p0−pv+2σ/R0=127786 Pap_{g0}=p_0-p_v+2\sigma/R_0=127786\ \mathrm{Pa}。由式(2.5),得到RB=18.242 μmR_B=18.242\ \mu\mathrm m和绝对压力pBlake=−2923.5 Pap_{\mathrm{Blake}}=-2923.5\ \mathrm{Pa}。这个负绝对压力表示模型中亚稳液体的拉伸状态,不是负气体压力,也不是普适实测初生阈值。

16.6 Probability depends on observation volume

概率依赖观测体积

Problem. A spatially uniform, constant-rate nucleation process has 50% probability of at least one event in volume VV during time τ\tau. What is the probability in volume 2V2V during the same time?

题目。一个空间均匀、速率恒定的成核过程,在体积VV、时间τ\tau内至少出现一次事件的概率为50%。同一时间内,在体积2V2V中的概率是多少?

Reference answer · original problem 6 / 原题参考答案
P=1−e−JVτP=1-e^{-JV\tau}
16.6 Probability depends on observation volume / 概率依赖观测体积; P: event probability; J: nucleation rate / areal impulse; V: bubble or observation volume; τ: duration / dimensionless time / dummy time

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. P=1−e−JVτ=0.5P=1-e^{-JV\tau}=0.5 implies JVτ=ln⁡2JV\tau=\ln2. Doubling volume gives P2=1−e−2ln⁡2=0.75P_2=1-e^{-2\ln2}=0.75. A threshold defined by a fixed event probability therefore shifts when observation volume changes, even if the microscopic rate law does not.

解答。由P=1−e−JVτ=0.5P=1-e^{-JV\tau}=0.5得JVτ=ln⁡2JV\tau=\ln2。体积加倍后,P2=1−e−2ln⁡2=0.75P_2=1-e^{-2\ln2}=0.75。因此,即使微观速率定律不变,以固定事件概率定义的阈值也会随观测体积变化。

16.7 Long-time inertial growth speed

长时间惯性增长速度

Problem. Under constant positive internal-minus-external driving pressure Δp\Delta p, neglect all losses and gas-law changes. What constant-speed asymptote is compatible with the radial equation?

题目。在恒定正内外驱动压差Δp\Delta p下,忽略全部损失和气体定律变化,什么恒速渐近值与径向方程相容?

Reference answer · original problem 7 / 原题参考答案
RR¨+32R˙2=Δp/ρR\ddot R+\tfrac32\dot R^2=\Delta p/\rho
16.7 Long-time inertial growth speed / 长时间惯性增长速度; R: instantaneous wall radius; R̈: radial wall acceleration; Ṙ: radial wall velocity; Δ: change / pressure difference; p: pressure at the indicated location; ρ: liquid / layer / droplet density

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Set R¨→0\ddot R\to0 in RR¨+3R˙2/2=Δp/ρR\ddot R+3\dot R^2/2=\Delta p/\rho. Then U∞=2Δp/(3ρ)U_\infty=\sqrt{2\Delta p/(3\rho)}. For Δp=98986 Pa\Delta p=98986\ \mathrm{Pa} and ρ=998 kg m−3\rho=998\ \mathrm{kg\,m^{-3}}, this is 8.132 m s−18.132\ \mathrm{m\,s^{-1}}. The result assumes the driving pressure remains available; it does not apply to a fixed gas mass whose pressure falls as the bubble grows.

解答。在RR¨+3R˙2/2=Δp/ρR\ddot R+3\dot R^2/2=\Delta p/\rho中令R¨→0\ddot R\to0,得到U∞=2Δp/(3ρ)U_\infty=\sqrt{2\Delta p/(3\rho)}。取Δp=98986 Pa\Delta p=98986\ \mathrm{Pa}和ρ=998 kg m−3\rho=998\ \mathrm{kg\,m^{-3}},可得8.132 m s−18.132\ \mathrm{m\,s^{-1}}。该结果假定驱动压差始终存在,不适用于压力随膨胀下降的固定质量气体。

16.8 Collapse-time scaling

塌缩时间尺度

Problem. Starting from the 100-micrometer example, determine the collapse time after doubling maximum radius, or after quadrupling driving pressure while keeping radius fixed.

题目。以100微米示例为基准,分别求最大半径加倍,以及半径不变而驱动压差增至四倍时的塌缩时间。

Reference answer · original problem 8 / 原题参考答案
tc=0.914681Rmax⁡ρ/Δpt_c=0.914681R_{\max}\sqrt{\rho/\Delta p}
16.8 Collapse-time scaling / 塌缩时间尺度; t_c: collapse duration; R_max: maximum expansion radius; ρ: liquid / layer / droplet density; Δ: change / pressure difference; p: pressure at the indicated location

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Since tc∝Rmax⁡Δp−1/2t_c\propto R_{\max}\Delta p^{-1/2}, the first change gives 18.369 μs18.369\ \mu\mathrm s, and the second gives 4.592 μs4.592\ \mu\mathrm s. Under independent small uncertainties, stc2/tc2≈(sR/R)2+14(sρ/ρ)2+14(sΔp/Δp)2s_{t_c}^2/t_c^2\approx(s_R/R)^2+\tfrac14(s_\rho/\rho)^2+\tfrac14(s_{\Delta p}/\Delta p)^2. Correlated uncertainties require covariance terms.

解答。由tc∝Rmax⁡Δp−1/2t_c\propto R_{\max}\Delta p^{-1/2},第一种变化得到18.369 μs18.369\ \mu\mathrm s,第二种得到4.592 μs4.592\ \mu\mathrm s。对于独立小不确定度,有stc2/tc2≈(sR/R)2+14(sρ/ρ)2+14(sΔp/Δp)2s_{t_c}^2/t_c^2\approx(s_R/R)^2+\tfrac14(s_\rho/\rho)^2+\tfrac14(s_{\Delta p}/\Delta p)^2。若误差相关,还需加入协方差项。

16.9 Speed at half maximum radius

半最大半径处的速度

Problem. Find the inward wall speed when an ideal Rayleigh cavity has contracted to R=Rmax⁡/2R=R_{\max}/2 under the teaching pressure difference.

题目。在教学压差下,理想Rayleigh空腔收缩至R=Rmax⁡/2R=R_{\max}/2时,求界面向内速度。

Reference answer · original problem 9 / 原题参考答案
R˙2=2Δp3ρ[(Rmax⁡/R)3−1]\dot R^2=\frac{2\Delta p}{3\rho}[(R_{\max}/R)^3-1]
16.9 Speed at half maximum radius / 半最大半径处的速度; Ṙ: radial wall velocity; Δ: change / pressure difference; p: pressure at the indicated location; ρ: liquid / layer / droplet density; R_max: maximum expansion radius; R: instantaneous wall radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Equation (5.2) gives U=−(2Δp/3ρ)(8−1)=−21.514 m s−1U=-\sqrt{(2\Delta p/3\rho)(8-1)}=-21.514\ \mathrm{m\,s^{-1}}. The value is independent of Rmax⁡R_{\max} in this scale-free idealization at a fixed radius ratio. Collapse time, however, remains proportional to Rmax⁡R_{\max}.

解答。式(5.2)给出U=−(2Δp/3ρ)(8−1)=−21.514 m s−1U=-\sqrt{(2\Delta p/3\rho)(8-1)}=-21.514\ \mathrm{m\,s^{-1}}。在这一无额外长度尺度的理想化中,固定半径比处的速度与Rmax⁡R_{\max}无关,但塌缩时间仍与Rmax⁡R_{\max}成正比。

16.10 When the incompressible collapse becomes suspect

不可压缩塌缩何时值得怀疑

Problem. Use the Rayleigh speed to find R/Rmax⁡R/R_{\max} when ∣U∣/c=0.10|U|/c=0.10 for the teaching values.

题目。利用Rayleigh速度关系,求教学参数下∣U∣/c=0.10|U|/c=0.10时的R/Rmax⁡R/R_{\max}。

Reference answer · original problem 10 / 原题参考答案
Mw=∣R˙∣/c,R/Rmax⁡=[1+3ρc2Mw2/(2Δp)]−1/3M_w=|\dot R|/c,\quad R/R_{\max}=[1+3\rho c^2M_w^2/(2\Delta p)]^{-1/3}
16.10 When the incompressible collapse becomes suspect / 不可压缩塌缩何时值得怀疑; M_w: generalized mass / molar mass; Ṙ: radial wall velocity; c: sound speed / heat capacity; R: instantaneous wall radius; R_max: maximum expansion radius; ρ: liquid / layer / droplet density; Δ: change / pressure difference; p: pressure at the indicated location

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Rearranging (5.2) gives R/Rmax⁡=[1+3ρc2M2/(2Δp)]−1/3=0.1443R/R_{\max}=[1+3\rho c^2M^2/(2\Delta p)]^{-1/3}=0.1443. Therefore a numerical collapse continued to 0.05Rmax⁡0.05R_{\max} is already outside this conservative physical warning level. It can still be used as a mathematical benchmark of the incompressible equation, but not as a credible physical late-collapse trajectory.

解答。整理式(5.2)得R/Rmax⁡=[1+3ρc2M2/(2Δp)]−1/3=0.1443R/R_{\max}=[1+3\rho c^2M^2/(2\Delta p)]^{-1/3}=0.1443。因此,将塌缩继续积分至0.05Rmax⁡0.05R_{\max},已经超出这一保守物理警戒范围。它仍可作为不可压缩方程的数学基准,但不能作为可信的真实塌缩末期轨迹。

16.11 Static compliance of a microbubble

微气泡的静态柔度

Problem. Estimate the relative radius change of the 10-micrometer reference bubble under a slowly applied additional pressure of 1000 Pa, using the code's polytropic closure.

题目。在配套代码的多方闭合下,估计10微米参考气泡在缓慢增加1000 Pa压力后的相对半径变化。

Reference answer · original problem 11 / 原题参考答案
ω02=(3κpg0−2σ/R0)/(ρR02),x=−δp/(ρR02ω02)\omega_0^2=(3\kappa p_{g0}-2\sigma/R_0)/(\rho R_0^2),\quad x=-\delta p/(\rho R_0^2\omega_0^2)
16.11 Static compliance of a microbubble / 微气泡的静态柔度; ω_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; x: position / dimensionless radius; δ: separation / perturbation / penetration depth; p: pressure at the indicated location

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Set time derivatives to zero in (6.1): x=−Pa/[3κpg0−2σ/R0]=−0.0021653x=-P_a/[3\kappa p_{g0}-2\sigma/R_0]=-0.0021653. The radius decreases by approximately 21.65 nm21.65\ \mathrm{nm}. A truly slow physical experiment may be closer to isothermal than the code's κ=1.4\kappa=1.4; the calculation is conditional on the stated closure.

解答。在式(6.1)中令时间导数为零,得x=−Pa/[3κpg0−2σ/R0]=−0.0021653x=-P_a/[3\kappa p_{g0}-2\sigma/R_0]=-0.0021653,半径约减少21.65 nm21.65\ \mathrm{nm}。真实缓慢实验可能比代码中的κ=1.4\kappa=1.4更接近等温;此计算以所指定闭合为条件。

16.12 Viscosity versus stiffness

黏度与刚度

Problem. If viscosity doubles while all other parameters stay fixed, what changes in the linear model?

题目。如果黏度加倍而其他参数不变,线性模型中哪些量发生变化?

Reference answer · original problem 12 / 原题参考答案
βμ=2μ/(ρR02)\beta_\mu=2\mu/(\rho R_0^2)
16.12 Viscosity versus stiffness / 黏度与刚度; β_μ: damping / logistic coefficient; μ: dynamic viscosity / absorption coefficient; ρ: liquid / layer / droplet density; R_0: reference equilibrium radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. βμ\beta_\mu doubles, while the undamped ω0\omega_0 in (6.2) is unchanged. The damped free frequency ω02−β2\sqrt{\omega_0^2-\beta^2} and forced-response peak can still change because they depend on damping. Thus “viscosity does not change frequency” is only correct when referring specifically to the undamped stiffness-to-mass frequency.

解答。βμ\beta_\mu加倍,而式(6.2)中的无阻尼ω0\omega_0不变。但有阻尼自由频率ω02−β2\sqrt{\omega_0^2-\beta^2}和受迫响应峰值位置仍可能改变,因为它们依赖阻尼。因此,“黏度不改变频率”仅在专指由刚度与质量决定的无阻尼频率时成立。

16.13 The isothermal gas potential

等温气体势

Problem. Derive the mechanical potential associated with pgV=pg0V0p_gV=p_{g0}V_0.

题目。推导与pgV=pg0V0p_gV=p_{g0}V_0对应的力学势。

Reference answer · original problem 13 / 原题参考答案
Wg=∫−pg dV,pg=pg0V0/VW_g=\int -p_g\,dV,\quad p_g=p_{g0}V_0/V
16.13 The isothermal gas potential / 等温气体势; W_g: work / potential energy; p_g: noncondensable gas pressure; V: bubble or observation volume; p_g0: reference gas pressure; V_0: bubble or observation volume

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. A reversible potential satisfies dWg=−pg dVdW_g=-p_g\,dV. Integrating gives Wg=−pg0V0ln⁡(V/V0)+constantW_g=-p_{g0}V_0\ln(V/V_0)+\mathrm{constant}. For an isothermal ideal gas, internal energy is constant; this logarithmic expression is therefore not the gas internal energy. It is the mechanical potential appropriate to the prescribed pressure law and thermal reservoir.

解答。可逆势满足dWg=−pg dVdW_g=-p_g\,dV,积分得Wg=−pg0V0ln⁡(V/V0)+constantW_g=-p_{g0}V_0\ln(V/V_0)+\mathrm{constant}。等温理想气体的内能恒定,因此这个对数表达式不是气体内能,而是与所指定压力定律及热库相对应的力学势。

16.14 The hidden acceleration in Keller–Miksis

Keller–Miksis中的隐藏加速度

Problem. Why does viscosity appear in the denominator of (7.3)?

题目。为什么黏度会出现在式(7.3)的分母中?

Reference answer · original problem 14 / 原题参考答案
p˙L=Dp−4μR¨/R\dot p_L=D_p-4\mu\ddot R/R
16.14 The hidden acceleration in Keller–Miksis / Keller–Miksis中的隐藏加速度; ṗ_L: pressure at the indicated location; D_p: bending rigidity / dose / residual derivative; μ: dynamic viscosity / absorption coefficient; R̈: radial wall acceleration; R: instantaneous wall radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. The liquid-side pressure contains −4μU/R-4\mu U/R. Differentiation contributes −4μU˙/R+4μU2/R2-4\mu\dot U/R+4\mu U^2/R^2. Multiplication by R/(ρc)R/(\rho c) in (7.1) gives −4μU˙/(ρc)-4\mu\dot U/(\rho c) on the right. Moving it to the left adds 4μ/(ρc)4\mu/(\rho c) to the coefficient of U˙\dot U. The effect is algebraically required, not an empirical correction.

解答。液侧压力包含−4μU/R-4\mu U/R,求导产生−4μU˙/R+4μU2/R2-4\mu\dot U/R+4\mu U^2/R^2。在式(7.1)中乘以R/(ρc)R/(\rho c)后,右端出现−4μU˙/(ρc)-4\mu\dot U/(\rho c)。将其移到左端,就使U˙\dot U的系数增加4μ/(ρc)4\mu/(\rho c)。这一项由代数必然产生,并非经验修正。

16.15 Radial frequency versus shape frequency

径向频率与形状频率

Problem. Find the ℓ=2\ell=2 static shape frequency for a 10-micrometer bubble and compare it with the radial frequency.

题目。求10微米气泡的ℓ=2\ell=2静态形状频率,并与径向频率比较。

Reference answer · original problem 15 / 原题参考答案
ωl2=(l−1)(l+1)(l+2)σ/(ρR03)\omega_l^2=(l-1)(l+1)(l+2)\sigma/(\rho R_0^3)
16.15 Radial frequency versus shape frequency / 径向频率与形状频率; ω_l: angular frequency; l: shape harmonic order; σ: interface tension / solid stress; ρ: liquid / layer / droplet density; R_0: reference equilibrium radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Equation (8.1) gives f2=12σ/(ρR3)/(2π)=148.09 kHzf_2=\sqrt{12\sigma/(\rho R^3)}/(2\pi)=148.09\ \mathrm{kHz}, compared with 342.37 kHz342.37\ \mathrm{kHz} for the radial mode under the teaching gas closure. Different deformation families have different restoring mechanisms. A single measured oscillation frequency should not be assigned to radial motion without checking the shape.

解答。由式(8.1),f2=12σ/(ρR3)/(2π)=148.09 kHzf_2=\sqrt{12\sigma/(\rho R^3)}/(2\pi)=148.09\ \mathrm{kHz},而教学气体闭合下径向频率为342.37 kHz342.37\ \mathrm{kHz}。不同变形类型具有不同回复机制。若未检查形状,就不能把某个实测振荡频率直接认定为径向振动。

16.16 A uniform pressure impulse

均匀压力冲量

Problem. Does a spatially uniform scalar pressure impulse create a velocity jump in an unbounded incompressible liquid?

题目。空间均匀的标量压力冲量,是否会在无界不可压缩液体中产生速度跃变?

Reference answer · original problem 16 / 原题参考答案
Ip=Δpτ,ΔU=Ip/(ρL)I_p=\Delta p\tau,\quad \Delta U=I_p/(\rho L)
16.16 A uniform pressure impulse / 均匀压力冲量; I_p: pressure impulse per unit area; Δ: change / pressure difference; p: pressure at the indicated location; τ: duration / dimensionless time / dummy time; U: velocity / velocity scale; ρ: liquid / layer / droplet density; L: length / latent heat / load disk radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Not through (9.1), because ∇Π=0\nabla\Pi=0. Boundaries, interfaces, or nonuniform loading are needed to generate the relevant pressure gradients. This does not say that uniform pressure on one side of a plate produces no force; that is a different boundary-value problem with a pressure difference across the structure.

解答。按式(9.1)不会,因为∇Π=0\nabla\Pi=0。需要边界、界面或非均匀加载形成相关压力梯度。这并不是说板一侧的均匀压力不产生力;后者是结构两侧存在压差的另一类边值问题。

16.17 Equal impulse, zero residual oscillation

相同冲量却无残余振动

Problem. For the undamped modal problem, what happens after a rectangular force pulse with ωnτp=2π\omega_n\tau_p=2\pi?

题目。对于无阻尼模态问题,若矩形力脉冲满足ωnτp=2π\omega_n\tau_p=2\pi,脉冲结束后会怎样?

Reference answer · original problem 17 / 原题参考答案
qfree=∣In∣Mnωn∣2sin⁡(ωnτp/2)ωnτp∣q_{\mathrm{free}}=\frac{|I_n|}{M_n\omega_n}\left|\frac{2\sin(\omega_n\tau_p/2)}{\omega_n\tau_p}\right|
16.17 Equal impulse, zero residual oscillation / 相同冲量却无残余振动; q_free: heat flux / mode displacement; I_n: pressure or modal impulse; M_n: generalized mass / molar mass; ω_n: structural modal angular frequency; τ_p: duration / dimensionless time / dummy time

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Equation (12.3) gives zero post-pulse free amplitude, because sin⁡π=0\sin\pi=0. The mode moved during the pulse; it did not remain stationary. The force first supplied and later recovered energy over one full oscillation. This is a counterexample to predicting structural response from total impulse alone when pulse duration is not negligible.

解答。式(12.3)给出脉冲后自由振幅为零,因为sin⁡π=0\sin\pi=0。模态在脉冲期间发生过运动,并不是始终静止。外力在完整振动周期内先输入能量,再回收能量。这说明当脉冲时长不可忽略时,不能只凭总冲量预测结构响应。

16.18 A necessary adhesion-energy scale

黏附能的必要尺度

Problem. Use mA=0.020 kg m−2m_A=0.020\ \mathrm{kg\,m^{-2}} and Gc=0.10 J m−2G_c=0.10\ \mathrm{J\,m^{-2}} to estimate the idealized velocity and pressure impulse required by the areal kinetic-energy comparison.

题目。取mA=0.020 kg m−2m_A=0.020\ \mathrm{kg\,m^{-2}}和Gc=0.10 J m−2G_c=0.10\ \mathrm{J\,m^{-2}},用面动能比较估计理想化速度与压力冲量。

Reference answer · original problem 18 / 原题参考答案
12mAv2=Gc,Ip=mAv\tfrac12m_Av^2=G_c,\quad I_p=m_Av
16.18 A necessary adhesion-energy scale / 黏附能的必要尺度; m_A: moving mass per unit area; v: velocity / squared wall speed as defined; G_c: fracture work per release area; I_p: pressure impulse per unit area

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. v=2Gc/mA=3.162 m s−1v=\sqrt{2G_c/m_A}=3.162\ \mathrm{m\,s^{-1}}, and Ip=mAv=0.06325 Pa sI_p=m_Av=0.06325\ \mathrm{Pa\,s}. These are not sufficient detachment thresholds. The estimate assumes all relevant kinetic energy becomes separation work and ignores bending, added liquid inertia, flaws, mode mix, and incomplete load coverage.

解答。v=2Gc/mA=3.162 m s−1v=\sqrt{2G_c/m_A}=3.162\ \mathrm{m\,s^{-1}},且Ip=mAv=0.06325 Pa sI_p=m_Av=0.06325\ \mathrm{Pa\,s}。这些数值不是充分脱离阈值。该估计假定相关动能全部转为分离功,并忽略弯曲、液体附加惯性、缺陷、混合模态和载荷覆盖不完整等因素。

16.19 Droplet expansion by mass conservation

由质量守恒估计液滴膨胀

Problem. For a hypothetical pure volatile liquid with ρd=1600 kg m−3\rho_d=1600\ \mathrm{kg\,m^{-3}}, M=0.300 kg mol−1M=0.300\ \mathrm{kg\,mol^{-1}}, Tb=300 KT_b=300\ \mathrm K, and pPF=105 Pap_{\mathrm{PF}}=10^5\ \mathrm{Pa}, estimate Rb/a0R_b/a_0 under complete ideal vaporization.

题目。对于一种假想纯挥发性液体,取ρd=1600 kg m−3\rho_d=1600\ \mathrm{kg\,m^{-3}}、M=0.300 kg mol−1M=0.300\ \mathrm{kg\,mol^{-1}}、Tb=300 KT_b=300\ \mathrm K和pPF=105 Pap_{\mathrm{PF}}=10^5\ \mathrm{Pa},在完全理想汽化假设下估计Rb/a0R_b/a_0。

Reference answer · original problem 19 / 原题参考答案
Rb/a0=[ρdRuT/(MpPF)]1/3R_b/a_0=[\rho_dR_uT/(Mp_{PF})]^{1/3}
16.19 Droplet expansion by mass conservation / 由质量守恒估计液滴膨胀; R_b: bubble wall radius; a_0: droplet radius / film span / acceleration; ρ_d: liquid / layer / droplet density; R_u: bubble wall radius; T: temperature / kinetic energy as defined; M: generalized mass / molar mass; p_PF: pressure at the indicated location

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Equation (11.1) gives Rb/a0=5.105R_b/a_0=5.105, corresponding to a volume ratio of about 133.0. This calculation provides no evidence that sufficient heat exists to complete vaporization. The parameters deliberately do not identify a measured PFC formulation.

解答。式(11.1)给出Rb/a0=5.105R_b/a_0=5.105,对应体积比约133.0。这个计算并不能证明体系具有足够热量完成汽化。所列参数有意不对应任何已测量的特定全氟碳配方。

16.20 Geometrical spreading of acoustic energy

声能的几何扩散

Problem. In the linear outgoing-wave approximation, if observation distance doubles, how do pressure amplitude and the integrated spherical acoustic-energy estimate change?

题目。在线性出射波近似下,观测距离加倍时,压力振幅和球面声能积分估计如何变化?

Reference answer · original problem 20 / 原题参考答案
Ea=4πr2ρc∫p′2dtE_a=\frac{4\pi r^2}{\rho c}\int p'^2dt
16.20 Geometrical spreading of acoustic energy / 声能的几何扩散; E_a: 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

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Pressure amplitude halves because p′∝1/rp'\propto1/r. The sphere area increases by four, while p′2p'^2 decreases by four, so (7.8) remains unchanged in the absence of absorption. This reasoning fails if near-field terms, reflections, nonlinear dissipation, or sensor averaging dominate.

解答。由于p′∝1/rp'\propto1/r,压力振幅减半。球面积增至四倍,而p′2p'^2减至四分之一,因此无吸收时式(7.8)保持不变。如果近场项、反射、非线性耗散或传感器平均效应占主导,这个推理就会失效。

16.21 A false temperature prediction

错误的温度预测

Problem. A student extrapolates T/T0=(R0/R)3(κ−1)T/T_0=(R_0/R)^{3(\kappa-1)} to arbitrarily small RR and reports the result as an experimental collapse temperature. Diagnose the mistake.

题目。某同学将T/T0=(R0/R)3(κ−1)T/T_0=(R_0/R)^{3(\kappa-1)}外推至任意小RR,并把结果报告为实验塌缩温度。错误在哪里?

Reference answer · original problem 21 / 原题参考答案
Tg/Tg0=(R0/R)3(κ−1)T_g/T_{g0}=(R_0/R)^{3(\kappa-1)}
16.21 A false temperature prediction / 错误的温度预测; T_g: temperature / kinetic energy as defined; T_g0: temperature / kinetic energy as defined; R_0: reference equilibrium radius; R: instantaneous wall radius; κ: gas polytropic exponent / curvature

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. The expression presumes a fixed-mass ideal gas and a specified polytropic process. At severe collapse, gas nonideality, phase change, thermal gradients, compressibility, chemistry, and nonsphericity may invalidate those premises. Even a valid model output is not an experimental measurement. The correct report labels it a conditional estimate and gives the model's validity limits.

解答。该表达式预设固定质量理想气体与指定多方过程。剧烈塌缩时,气体非理想性、相变、温度梯度、可压缩性、化学过程和非球形性都可能使这些前提失效。即使模型输出有效,它也不是实验测量。正确报告应将其标为条件性估计,并说明模型适用范围。

16.22 A mesh-convergence calculation

网格收敛计算

Problem. A scalar output is 1.00, 0.92, and 0.90 on meshes h,h/2,h/4h,h/2,h/4. Assuming a clean asymptotic error Qh=Q∗+ChpQ_h=Q_*+Ch^p, estimate pp and Q∗Q_*.

题目。某标量输出在网格h,h/2,h/4h,h/2,h/4上分别为1.00、0.92和0.90。假设具有清晰渐近误差Qh=Q∗+ChpQ_h=Q_*+Ch^p,估计pp和Q∗Q_*。

Reference answer · original problem 22 / 原题参考答案
p=log⁡2[(Qh−Qh/2)/(Qh/2−Qh/4)],Q∗=Qh/4+(Qh/4−Qh/2)/(2p−1)p=\log_2[(Q_h-Q_{h/2})/(Q_{h/2}-Q_{h/4})],\quad Q_*=Q_{h/4}+(Q_{h/4}-Q_{h/2})/(2^p-1)
16.22 A mesh-convergence calculation / 网格收敛计算; p: pressure at the indicated location; Q_h: heat or generalized load as defined; Q_h/2: heat or generalized load as defined; Q_h/4: heat or generalized load as defined; Q: heat or generalized load as defined

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. 2p=(1.00−0.92)/(0.92−0.90)=42^p=(1.00-0.92)/(0.92-0.90)=4, so p=2p=2. Richardson extrapolation gives Q∗=0.90+(0.90−0.92)/(4−1)=0.89333Q_*=0.90+(0.90-0.92)/(4-1)=0.89333. The inference is conditional on the asymptotic model. Interface topology changes, stochastic nucleation, and unresolved shock peaks can make three-grid extrapolation misleading.

解答。2p=(1.00−0.92)/(0.92−0.90)=42^p=(1.00-0.92)/(0.92-0.90)=4,因此p=2p=2。Richardson外推得Q∗=0.90+(0.90−0.92)/(4−1)=0.89333Q_*=0.90+(0.90-0.92)/(4-1)=0.89333。这一推断依赖渐近误差模型;界面拓扑变化、随机成核和未解析激波峰值,都可能使三网格外推产生误导。

16.23 A nonidentifiable gas parameterization

不可辨识的气体参数化

Problem. Suppose the gas law is written pg=C/R3κp_g=C/R^{3\kappa} with fixed κ\kappa. Can a radius trajectory independently determine pg0p_{g0} and R0R_0 when only their product C=pg0R03κC=p_{g0}R_0^{3\kappa} enters the model?

题目。气体定律写为pg=C/R3κp_g=C/R^{3\kappa},且κ\kappa固定。如果模型中只出现乘积C=pg0R03κC=p_{g0}R_0^{3\kappa},半径轨迹能否独立确定pg0p_{g0}与R0R_0?

Reference answer · original problem 23 / 原题参考答案
pg=C/R3κ,C=pg0R03κp_g=C/R^{3\kappa},\quad C=p_{g0}R_0^{3\kappa}
16.23 A nonidentifiable gas parameterization / 不可辨识的气体参数化; p_g: noncondensable gas pressure; C: gas-content constant / damping; R: instantaneous wall radius; κ: gas polytropic exponent / curvature; p_g0: reference gas pressure; R_0: reference equilibrium radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. No. Infinitely many pairs yield the same CC and hence the same pressure law. An additional equilibrium relation or independent measurement is needed. An optimizer returning two precise-looking numbers does not remove this structural nonidentifiability.

解答。不能。无穷多组参数都能给出相同CC,从而产生相同压力定律。需要附加平衡关系或独立测量。优化器即使返回两个看似精确的数值,也不能消除这种结构性不可辨识。

16.24 Why more frames do not automatically fix acceleration

为什么更多图像帧不自动解决加速度问题

Problem. If independent radius noise remains fixed while sampling interval halves, how does the noise scale in a naive second difference?

题目。如果独立半径噪声保持不变,而采样间隔减半,简单二阶差分中的噪声尺度如何变化?

Reference answer · original problem 24 / 原题参考答案
R¨i≈(Ri+1−2Ri+Ri−1)/Δt2\ddot R_i\approx(R_{i+1}-2R_i+R_{i-1})/\Delta t^2
16.24 Why more frames do not automatically fix acceleration / 为什么更多图像帧不自动解决加速度问题; R̈_i: bubble wall radius; R_i+1: bubble wall radius; R_i: bubble wall radius; R_i-1: bubble wall radius; Δ: change / pressure difference; t: physical time

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. It grows by a factor of four because the denominator is proportional to Δt2\Delta t^2. Higher temporal resolution is valuable, but acceleration estimation needs an appropriate observation model, smoothing with quantified bias, or direct fitting of a dynamical model. Do not confuse temporal resolution with derivative precision.

解答。由于分母与Δt2\Delta t^2成正比,噪声尺度增大至四倍。更高时间分辨率仍然有价值,但加速度估计需要适当观测模型、偏差可量化的平滑方法,或直接拟合动力学模型。不要将时间分辨率与导数精度混为一谈。

16.25 Source-mechanism classification

源机制分类

Problem. Classify [S390], [S72], [S343], and [S340] by their principal route to mechanical loading.

题目。按照主要机械加载路径,对[S390]、[S72]、[S343]和[S340]分类。

Reference answer · original problem 25 / 原题参考答案
Elaser⟶pacoustic(t)⟶R(t)E_{\mathrm{laser}}\longrightarrow p_{\mathrm{acoustic}}(t)\longrightarrow R(t)
16.25 Source-mechanism classification / 源机制分类; E_laser: energy / elastic modulus as defined; p_acoustic: pressure at the indicated location; t: physical time; R: instantaneous wall radius

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. [S390]: optoacoustic transmitter, focused ultrasound, nucleation, and liquid jetting. [S72]: photothermal phase change inside a hydrogel-composite stamp, surface bulging, and adhesion modification. [S343]: blister motion coupled to cavitation in a liquid-transfer setting. [S340]: stress-wave-mediated thin-film delamination. The common presence of a laser does not make the four source models interchangeable.

解答。[S390]:光声发射器、聚焦超声、成核与液体喷射。[S72]:水凝胶复合印章内部光热相变、表面鼓起与黏附改变。[S343]:液体转移环境中鼓包运动与空化耦合。[S340]:应力波介导的薄膜分层。共同使用激光,并不使四种源模型可以相互替代。

16.26 A terminology audit

术语审查

Problem. Why is it wrong to cite [S371] directly as evidence for a vapor-cavitation nucleation threshold?

题目。为什么不能直接引用[S371]作为蒸气空化成核阈值的证据?

Reference answer · original problem 26 / 原题参考答案
pPFVb=mPFRuT/Mp_{PF}V_b=m_{PF}R_uT/M
16.26 A terminology audit / 术语审查; p_PF: pressure at the indicated location; V_b: bubble or observation volume; m_PF: mass or areal mass; R_u: bubble wall radius; T: temperature / kinetic energy as defined; M: generalized mass / molar mass

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Its nucleation process forms liquid perfluorocarbon droplets through solvent exchange. The new phase and free-energy problem differ from vapor embryos forming in a liquid. A relevant title word is not enough; phase identity and process must match the claim.

解答。该文的成核过程通过溶剂交换形成液态全氟碳液滴。其新生相与自由能问题,不同于液体中形成蒸气胚核的过程。标题中存在相同关键词并不足够,物相与过程必须和论点匹配。

16.27 Detect an overconstrained interface

识别过度约束界面

Problem. In a fluid simulation, both a bubble radius history and an unrelated bubble pressure history are prescribed on the same interface. Why is this dangerous?

题目。在流体仿真中,同一界面同时指定气泡半径历程和一个无关的气泡压力历程,为什么危险?

Reference answer · original problem 27 / 原题参考答案
pB,required=p∞+ρ(RR¨+32R˙2)+2σ/R+4μR˙/Rp_{B,\mathrm{required}}=p_\infty+\rho(R\ddot R+\tfrac32\dot R^2)+2\sigma/R+4\mu\dot R/R
16.27 Detect an overconstrained interface / 识别过度约束界面; p_B,required: pressure at the indicated location; p∞: ambient liquid pressure; ρ: liquid / layer / droplet density; R: instantaneous wall radius; R̈: radial wall acceleration; Ṙ: radial wall velocity; σ: interface tension / solid stress; μ: dynamic viscosity / absorption coefficient

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. The motion and traction are linked by the fluid equations. Prescribing both independently usually overconstrains the boundary or injects energy inconsistent with either physical source. Prescribe kinematics and calculate reaction pressure, or prescribe a thermodynamic pressure law and solve for motion, unless a carefully formulated coupling supplies a consistent relation.

解答。运动与牵引通过流体方程相联系。独立指定二者通常会使边界过度约束,或注入与任何真实源都不一致的能量。除非具有严格相容的耦合关系,否则应指定运动学并计算反作用压力,或指定热力学压力定律并求解运动。

16.28 Avoid double-counting fluid inertia

避免重复计算流体惯性

Problem. A fully coupled fluid–structure model already resolves the water next to a film. A researcher also adds a hydrodynamic added mass to the film density. What should be checked?

题目。完整流固耦合模型已经解析薄膜旁的水域,研究者又给薄膜密度加上水动力附加质量。应该检查什么?

Reference answer · original problem 28 / 原题参考答案
(Ms+Ma)q¨+Cq˙+Kq=Q(t)(M_s+M_a)\ddot q+C\dot q+Kq=Q(t)
16.28 Avoid double-counting fluid inertia / 避免重复计算流体惯性; M_s: generalized mass / molar mass; M_a: generalized mass / molar mass; q̈: heat flux / mode displacement; C: gas-content constant / damping; q̇: heat flux / mode displacement; K: stiffness / collapse constant; q: heat flux / mode displacement; Q: heat or generalized load as defined; t: physical time

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Determine whether the added mass represents liquid already included in the solved fluid domain. If so, it is double-counted. Added-mass corrections belong to reduced models or omitted fluid regions, with a clearly defined partition. A coupled model should recover a known wet-versus-dry structural benchmark without arbitrary duplicate inertia.

解答。应判断附加质量是否表示已包含在求解流体域中的液体。如果是,就发生了重复计算。附加质量修正应当用于降阶模型或被省略的流体区域,并明确划分。耦合模型应能在没有任意重复惯性的情况下复现已知干态与浸液结构基准。

16.29 Design a discriminating experiment

设计具有区分能力的实验

Problem. Two models fit maximum bubble radius equally well, but one assumes constant vapor pressure and the other finite-rate thermal transport. Suggest a more discriminating measurement.

题目。两个模型同样准确地拟合最大气泡半径,但一个假设蒸气压恒定,另一个包含有限速率热输运。建议一种更有区分能力的测量。

Reference answer · original problem 29 / 原题参考答案
ta=ℓ/c,tth=ℓ2/αt_a=\ell/c,\quad t_{\mathrm{th}}=\ell^2/\alpha
16.29 Design a discriminating experiment / 设计具有区分能力的实验; t_a: acoustic crossing time; ℓ: absorption / cohesive-zone length; c: sound speed / heat capacity; t_th: physical time; α: thermal diffusivity

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Compare full growth and collapse histories, rebound radius and period, and their variation with independently controlled ambient temperature or pulse duration. A single extra peak is less informative than a controlled family of trajectories. Fit shared parameters on one subset and test on held-out conditions; do not refit every trace independently and call the result prediction.

解答。可比较完整增长与塌缩历程、回弹半径和周期,以及它们随独立控制的环境温度或脉冲时长的变化。相比只增加一个峰值,受控的一组轨迹更有信息。应在部分数据上拟合共享参数,再在留出条件上检验;不要对每条曲线分别重新拟合后仍称之为预测。

16.30 Build a falsifiable capstone study

构建可证伪综合研究

Problem. Propose a minimal study of the hypothesis that matching bubble-loading time to a film mode improves useful payload motion.

题目。针对“使气泡加载时间匹配薄膜模态可改善有效负载运动”的假设,提出最小研究方案。

Reference answer · original problem 30 / 原题参考答案
Θ=ωntb,Qn=∫AΔpϕndA\Theta=\omega_nt_b,\quad Q_n=\int_A\Delta p\phi_n dA
16.30 Build a falsifiable capstone study / 构建可证伪综合研究; Θ: loading-to-mode timing ratio; ω_n: structural modal angular frequency; t_b: physical time; Q_n: heat or generalized load as defined; Δ: change / pressure difference; p: pressure at the indicated location; φ_n: velocity potential / mode shape; A: area / volume-acceleration coefficient

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.

图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

Solution. Fix geometry, support, payload mass, and a declared source normalization such as fixed volume excursion or fixed work. Sweep loading duration; measure or calculate mode-projected force, payload velocity, bending energy, and losses. Verify the uncoupled oscillator and fluid limit first. Compare at least one held-out duration and one altered support stiffness. A peak in useful motion may support the hypothesis, but only after showing that it is not caused by a hidden change in source energy or by numerical resonance.

解答。固定几何、支承、负载质量,并声明源归一化方式,例如固定体积变化或固定功。扫描加载时长,测量或计算模态投影力、负载速度、弯曲能和损失。首先验证未耦合振子与流体极限,然后至少检验一个留出时长和一个改变支承刚度的条件。有效运动出现峰值可能支持假设,但必须先证明它不是由隐含源能量变化或数值共振引起的。

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.

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

  1. Combine the original problems: for Rmax = 100 µm, find ideal collapse time, half-radius speed and Mw at that radius. Separate a benchmark result from a pressure maximum.

    综合原习题:Rmax = 100 µm 时求理想塌缩时间、半半径速度及该处 Mw,区分基准结果与压力最大值。

    Reference answer / 参考答案
    tc=0.914681Rmax⁡ρ/Δp,R˙=−2Δp[(Rmax⁡/R)3−1]/(3ρ),Mw=∣R˙∣/ct_c=0.914681R_{\max}\sqrt{\rho/\Delta p},\quad \dot R=-\sqrt{2\Delta p[(R_{\max}/R)^3-1]/(3\rho)},\quad M_w=|\dot R|/c
    Thirty original worked problems · original governing laws; t_c: collapse duration; R_max: maximum expansion radius; ρ: liquid / layer / droplet density; Δ: change / pressure difference; p: pressure at the indicated location; Ṙ: radial wall velocity; R: instantaneous wall radius; M_w: generalized mass / molar mass; c: sound speed / heat capacity

    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.

    图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

    tc = 9.184 µs, Ṙ = −21.514 m/s, Mw = 0.01452 at half radius. These verify a finite-state ideal radial solution. They do not supply a finite singular-limit maximum or a boundary jet-impact pressure.

    半半径处 tc = 9.184 µs、Ṙ = −21.514 m/s、Mw = 0.01452。这些核查有限状态的理想径向解,不给出奇异极限有限最大值或边界射流冲击压力。

  2. Connect this calculation to a film with mA = 0.020 kg/m² and Gc = 0.10 J/m². State a necessary impulse scale and at least three additional conditions for actual release.

    将该计算与 mA = 0.020 kg/m²、Gc = 0.10 J/m² 的膜连接,给出必要冲量尺度及实际释放至少三个额外条件。

    Reference answer / 参考答案
    Ip,nom=2mAGc,Qn=∫AΔpϕn dAI_{p,\mathrm{nom}}=\sqrt{2m_AG_c},\qquad Q_n=\int_A\Delta p\phi_n\,dA
    Thirty original worked problems · original governing laws; I_p,nom: pressure or modal impulse; m_A: moving mass per unit area; G_c: fracture work per release area; Q_n: heat or generalized load as defined; Δ: change / pressure difference; p: pressure at the indicated location; φ_n: velocity potential / mode shape; A: area / volume-acceleration coefficient

    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.

    图示标出物理体系,图例列出公式中的符号。若某字母有多种作用,应以邻近推导中的定义为准。

    Ip,nom = 0.06325 Pa s. Need physically consistent local traction, spatial and temporal modal participation, a calibrated rate/mixed-mode cohesive law, and a propagation/complete-separation criterion. Read all thirty original questions below; this integrated check supplements them.

    Ip,nom = 0.06325 Pa s。需要物理自洽局部牵引、空间与时间模态参与、经标定的速率与混合模态黏聚律,以及传播和完全分离判据。下方全部三十道原题仍需阅读;本综合检查是其补充。

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.

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