超声近壁微气泡溃灭过程详解 [page::2][page::3]

• 超声作用引起气泡内外压差导致体积震荡及非对称溃灭。

- 一次溃灭过程产生最大射流速度约250 m/s，射流穿透气泡且气泡分裂为多个碎小气泡。

• 二次溃灭阶段射流速度更高，达到350 m/s，碎小气泡群继续膨胀压缩，溃灭过程复杂，伴随多处漩涡形成。

最大射流速度影响因素分析 [page::4][page::5]

| 初始半径 ($\upmu$m) | 频率 (kHz) | 无量纲近壁距离 $\gamma$ | 最大射流速度 (m/s) |
|---------------------|------------|--------------------------|---------------------|
| 60 | 30 | 1.1 - 1.6 | 190 - 270 |
| 80 | 30 | 1.1 - 1.6 | 210 - 380 |
| 100 | 30 | 1.1 - 1.6 | 200 - 410 |

• 气泡与壁面距离越近，最大射流速度越高，近壁距离$\gamma$与最大射流速度成正比。

- 超声频率升高，最大射流速度降低，因高频缩短溃灭周期，能量储存不足。

• 气泡初始半径越大，最大射流速度越高，因气泡原始势能更大。

数值模拟方法与模型验证 [page::1][page::2]

• 采用有限体积法结合VOF模型模拟二维轴对称近壁微气泡运动。

- 通过网格无关性测试及与Rayleigh-Plesset方程解析解比对，验证模型准确性。

研究报告详尽分析 — 超声对近壁微气泡溃灭过程的影响

---

一、元数据与报告概览

• 报告标题：《超声对近壁微气泡溃灭过程的影响》

- 作者：王舰航，陈韬厚，包福兵†，王月兵

• 机构：中国计量大学计量测试工程学院，杭州 310018

- 发表时间及刊物信息：2020年，中国计量大学相关学报，文献号1000-310X(2020)03-0329-07，DOI: 10.11684/j.issn.1000-310X.2020.03.002

• 主题：研究超声波作用下，近壁（靠近刚性边界面）微气泡的空化溃灭过程，通过数值模拟探讨气泡动力学和射流速度规律，为靶向给药和基因转染等医学应用提供理论依据。

核心论点：

• 气泡溃灭最大射流速度与近壁距离无量纲参数 $\gamma$、超声频率及气泡初始半径均存明显关系：

- 射流速度与 $\gamma$（1.1 ~ 1.6范围内）呈正比
- 射流速度与超声频率（10 ~ 60 kHz范围）成正比
- 射流速度与气泡初始半径（50 ~ 100 μm范围）成反比

• 气泡的二次溃灭阶段最大射流速度高于一次溃灭，暗示在过程整体影响中二次溃灭尤为重要。

- 本文采用有限体积法结合VOF模型进行二维轴对称数值模拟，仿真实现超声驱动下气泡体积、压力、温度及射流速度的时空演化。

主要意图传递超声参数及边界条件对微泡动力学影响规律，支撑临床靶向给药/基因治疗中微泡超声空化机制的基础研究。

---

二、逐节深度解读

0. 摘要与引言部分（第0-1页）

• 摘要明确指出研究目标是利用数值模拟探讨超声对近壁微气泡溃灭动力学的影响。当前技术瓶颈是基于Rayleigh-Plesset方程的球泡理论不适用于近壁微气泡溃灭的复杂不对称过程，故采用VOF模型结合有限体积法分析。
• 引言部分（第1页）回顾了基因治疗中超声靶向微泡破坏技术的进展和优势，指出其利用超声产生声孔效应促进基因转染的效率和安全性，强调微气泡超声空化在药物靶向输运上的重要性。[page::0,1]
• 引用大量文献论述近壁气泡的实验和数值模拟研究现状和不足，特别是强调了目前多聚焦在一次溃灭上，而二次溃灭研究较少，填补二次溃灭动态过程的研究空缺。
• 理论模型部分明确解释VOF模型的基本假设和方程体系，包括体积分数方程、连续性和动量方程，确定气相为可压缩流体，液相为不可压缩流体，模型适用于本文计算。

1. 数值计算方法（第1-2页）

• 利用二维轴对称矩形计算区（500μm × 250μm）进行模拟，设定刚性无滑移壁面作为边界条件，图1示意了计算域及气泡相对位置，定义无量纲近壁距离参数为 $\gamma = L / R0$，$L$为气泡中心与壁面的距离，$R0$为气泡初始半径。
• 边界条件设置中，超声场模拟为边界压力随时间的周期性变化，模拟压力为 $p = pl + pd \cos(2\pi f t)$，其中 $pl$为液体静压（默认为标准大气压），$pd$为声压振幅，$f$为超声频率。
• 初始气泡压力由表面张力及液相压力决定，采用Rayleigh理论近似气泡溃灭时间计算气泡所经历的周期。
• 网格无关性验证详尽，采用5套不同大小的结构化网格（从约1.7万到33万单元），结果显示当网格超过约23万单元时，气泡体积变化曲线收敛，且模拟结果与Rayleigh-Plesset方程解析解高度拟合，验证了数值模型的可靠性和稳定性。[page::1,2]

2. 结果与讨论（第2-5页）

2.1 超声作用下的微气泡近壁溃灭过程（第2-4页）

• 报告通过图3（图像展示）详细描绘了超声驱动下近壁微气泡溃灭的形态及流场演变：

- 初始阶段气泡逐渐被压缩，远壁端压力达到峰值，气泡出现凹陷与射流，射流速度可超过250 m/s，且产生旋涡结构。
- 气泡被射流穿透，分裂为两个较小气泡，碎小气泡沿中轴线远离壁面移动。
- 后续出现二次溃灭，分裂气泡再次被压缩，高速液流沿壁面产生复杂的漩涡，射流速度甚至超过一次溃灭（最高约350 m/s），反映了二次溃灭更强的动力效应。
- 末期，原大气泡完全崩解，变为多个小气泡群，流场复杂且呈现多重涡旋结构，标志着溃灭过程结束。

• 图4展现了气泡内部压力、体积及最大射流速度随时间变化的定量趋势：

- 一次溃灭时，压力峰值约10 MPa，射流速度峰值约250 m/s，温度约850 K。
- 之后气泡反弹膨胀，压力和速度下降。
- 二次溃灭时，压力与温度峰值较一次溃灭低，但射流速度达到更高峰值（约350 m/s），进一步证明二次溃灭动能释放更为剧烈。

这一动态过程分析为理解近壁微气泡溃灭的非对称射流及其力学效应提供了系统且微观的数值支持。[page::2,3,4]

2.2 最大射流速度影响因素分析（第4-5页）

报告重点分析了影响最大射流速度的三大关键参数：

1. 近壁距离 $\gamma$ 的影响（图5）：

- 随着 $\gamma$ 从1.1增加到1.6，最大射流速度明显下降。
- 解释为气泡离壁面越远，壁面对气泡非对称压缩的影响减弱，气泡溃灭过程趋于均匀，射流动能减少。
- 边界限制越明显，速度越大。

1. 超声频率的影响（图6）：

- 频率越高，气泡溃灭的最大射流速度越低。
- 理由是高频率缩短了溃灭周期，气泡缺乏足够时间积累能量，从而减少释放的动能，符合现有文献的研究趋势。

1. 气泡初始半径的影响（图7）：

- 气泡半径增大时，最大射流速度增加。
- 这一点与其堆积的势能成正比解释相符，初始半径越大，潜在能越大，因而释放时射流动能更高。

这三方面的定量关系揭示了超声参数与微泡溃灭动力学的关键控制因素，对于基于超声微气泡的精准医疗技术应用具备指导意义。[page::4,5]

3. 结论（第5页）

• 数值模拟清晰揭示了超声对近壁微气泡溃灭过程的多维度影响机制，完整揭示了溃灭过程中气泡形态、流场动量变化以及动力学特征。

- 最大射流速度呈现与无量纲近壁距离、超声频率和气泡初始半径的正/反比关系，为超声靶向给药、基因转染提供理论支撑。

• 二次溃灭阶段的射流速度及动能更加突出，强调完全考虑整个溃灭过程的重要性。

---

三、图表深度解读

图1 计算示意图（第1页）

• 描绘了二维轴对称计算域尺寸（500μm × 250μm），气泡的位置、高度与壁面的关系。

- 壁面为无滑移刚性边界，气泡中心距离壁面距离为 $L$。

• 表达了无量纲参数 $\gamma = L / R_0$ 用于后续参数化分析。

此图为数值模型几何设定提供直观基础。[page::1]

---

图2 不同网格数下气泡体积随时间变化（第2页）

• 展示不同网格数量对气泡体积变化曲线的影响，曲线与Rayleigh-Plesset方程近似解析解进行对比。

- 结果显示气泡达到最小体积前，网格影响不明显；溃灭后气泡回弹时，细化网格使体积预测更准确；

• 网格超过约23万时曲线趋于稳定，验证了网格无关性和计算准确性。

这一图表确保了报告数值模拟的可靠性，是数值方案合理的关键验证。[page::2]

---

图3 超声近壁微气泡溃灭过程（第3页）

• 8张子图结合流线（左半部分）与压力云图（右半部分），揭示气泡溃灭不同时间截面的流场与压力分布。

- 气泡由初始收缩到射流穿透，再到气泡分裂与碎小气泡扩散，动态清晰展现。

• 最大射流速度峰值出现在射流穿透初期（$t\approx 3.6 \rm\mu s$），达到约250 m/s，压力峰值在穿透点最大。

- 后续二次溃灭产生更多涡旋和复杂流动，碎小泡群形式多样。

此图集整体支持后续量化结果与物理分析，形象直观地揭示复杂的近壁微泡空化机制。[page::3]

---

图4 二次溃灭过程各参数随时间变化（第4页）

• 三个叠加子图描绘二次溃灭过程中气泡体积比例、内部压力和最大射流速度变化。

- 一次溃灭对应首个压力和速度高峰，二次溃灭出现第二个速度峰值且高于第一次，压力峰较低。

• 气泡体积呈振荡周期性变化，反映溃灭后气泡的膨胀与压缩交替状态。

- 最大速度峰值约350 m/s显著提升二次溃灭对系统整体影响的重视。

数据定量阐述了二次溃灭过程是关键的动力学阶段，强调不可忽视。[page::4]

---

图5 最大射流速度随不同壁面距离变化（第4页）

• 三条曲线代表不同气泡初始半径，x轴为无量纲近壁距离$\gamma$，y轴为最大射流速度。

- 整体趋势为随着$\gamma$增加（气泡远离壁面），最大射流速度递减。

• 射流速度随着气泡半径增大而增大，底层物理受壁面限制对射流的非对称压力贡献显著。

图表定性明确气泡-壁面相互作用是射流速度的关键控制因素之一。[page::4]

---

图6 最大射流速度随超声频率变化（第5页）

• 曲线显示不同气泡初始半径下，随着超声频率 $f$ 增大（10kHz 到 60kHz），最大射流速度总体下降。

- 该趋势符合超声周期缩短，气泡能量积累不足导致溃灭动力学减弱的物理解释。

图6确认了超声频率调整对气泡溃灭效应调控的可行性和方向性。[page::5]

---

图7 最大射流速度随不同初始半径变化（第5页）

• 曲线描绘在固定近壁距离和不同频率下，最大射流速度与气泡初始半径的正相关关系。

- 半径越大能量潜势越大，动能释放越猛烈，故射流速度随之提升。

图7揭示微气泡尺寸作为动力学驱动参数的重要性，具备实验设计指导意义。[page::5]

---

四、估值分析

本报告属于物理数值模拟领域，无涉及传统金融估值分析，故无相关估值方法。

---

五、风险因素评估

报告中未对风险因素进行专门评述，作为数值模拟科研文献，主要关注模型准确性和计算复杂度：

• 数值误差与网格敏感性问题通过网格无关性验证进行了有效管控，剩余风险较低。

- 模型假设局限：二维轴对称模型对真实三维气泡行为的简化，气泡内气相为可压缩流体、液相不可压缩，VOF模型假定单一压缩相，仅允许一相压缩。

• 超声参数范围有限，实际临床环境参数复杂多变，实验验证及三维模拟是后续工作必要方向。

---

六、批判性视角与细微差别

• 本研究完备地展示了数值模拟方法及结果，逻辑清晰。但模型为二维轴对称，真实三维效应可能存在偏差，特别是复杂形态溃灭和射流的非对称性。

- 参数范围集中在一定超声频率与气泡初始半径中，缺少对更宽或更细参数空间的探索或敏感性分析。

• 二次溃灭的机制解析较深入，但未具体讨论二次溃灭对组织或细胞的实际生物学影响，这一桥接仍需实验和理论补充。

- 算法迭代和残差收敛严格，但缺少针对不同流体模型、表面张力和可压缩性探讨的稳定性讨论。

• 结论中所述“射流速度与超声频率成正比”与图6趋势（频率增大时速度下降）矛盾，文中理应是“反比”表述，注意该处表述不一致，应核实原文页码确认准确性。[page::0,5]

---

七、结论性综合

本报告基于VOF模型的二维轴对称有限体积数值模拟系统地研究了超声波驱动下近壁微气泡的溃灭过程。模拟结果详细揭示：

• 超声激发的近壁微气泡溃灭为复杂的非对称过程，包括气泡远壁端受压凹陷、射流穿透、气泡一次分裂及后续多个碎小气泡的二次溃灭；

- 射流速度峰值高达数百米每秒，边界条件与气泡参数显著影响射流最大速度，反映气泡动力学在微观尺度的复杂反馈机制；

• 最大射流速度与近壁距离无量纲参数 $\gamma$ 呈负相关（壁面限制越强，速度越高），与超声频率呈负相关（频率越高，射流速度越低），与初始气泡半径呈正相关（气泡越大，动能越高）；

- 二次溃灭的最大射流速度高于一次溃灭，提示溃灭全程动力学机理不能忽视二次效应；

• 本研究数值模型与经典理论解高度吻合，网格无关性充分验证了计算的准确性和可靠性；

- 图3与图4中直观展示了复杂的流场结构和物理参数演化，对理解溃灭机制及其医疗应用中的物理本质至关重要。

总体来说，该报告为超声辅助微气泡的医学应用，特别是靶向给药和基因转染提供了坚实的理论和数值基础及操作指导价值。报告严谨详实且图文并茂，展示了超声微气泡物理动力学的关键规律及其优化潜力，具有一定的创新性和应用前景。[page::0-5]

---

参考引证溯源

• 文章各结论及图表内容出自报告原文第0至5页，尤其基于摘要、计算方法、结果讨论与结论章节内容作出分析。[page::0,1,2,3,4,5]
• 图像与数据引用均详见报告附图（图1至图7）及相应页码说明。[page::1,2,3,4,5]

---

总结

该报告系统阐述了超声波作用下近壁微气泡溃灭的过程和关键影响因素，通过严谨的数值模拟验证了物理机制的合理性，确认了二次溃灭在整体能量释放中的主导地位。研究结果为未来基于超声的精准医疗技术优化提供了理论依据和数值支持，是微气泡动力学领域的高质量学术贡献。

报告