数值模型与求解方法介绍 [page::2][page::3][page::4]

• 开发基于OpenFOAM的compressibleInterPhaseChangeFoam求解器，考虑液-蒸气两相可压缩混合物，采用VOF方法捕捉界面。

- 引入Schnerr–Sauer模型处理相变，使用物态方程闭合方程组，考虑表面张力和热传导。

• 计算域利用轴对称简化，局部网格加密，时间空间数值格式为二阶精度，保证流场及声压场准确。

模型验证与参数选择 [page::5][page::6]

• 通过对比简化Rayleigh-Plesset方程解析解验证气泡半径时间演化，确定径向网格细胞数为150以平衡精度与效率。

- 验证压力峰值与经验公式吻合良好，确保声能计算可靠。

气泡在壁面附近的塌陷形态演化 [page::6][page::7]

• 观察到气泡形态由球形变为心形，射流击穿气泡底部，形成环状空泡及两个震荡冲击波的声能发射过程。

- 射流强度及塌陷速度受无量纲距离γ的影响，壁面约束显著改变气泡动力学行为。

声能随时间及γ的变化规律 [page::8][page::9]

• 声能在气泡塌陷大部分时间处于低水平，气泡达到最小体积后声能迅速上升。

- γ对最大声能影响明显，γ<0.9时声能随γ减小而增加，γ>0.9时声能随γ增大而增加，揭示壁面距离调节空化声能释放的机制。

不同γ值下速度与压力场差异 [page::9]

• γ=0.9时，环状涡核向壁面移动，抑制气泡压缩及声能释放。

- γ=0.3时涡环强度减弱，气泡压缩更剧烈，声能更强。

结论与应用意义 [page::9]

• 本研究为刚性壁附近气泡塌陷音频辐射提供了准确数值模拟工具，揭示了不同泡壁距离对空化声能释放与侵蚀机理的影响，有助于水力机械设备设计中的侵蚀风险早期评估。

对报告《Numerical Study of the Emission of Acoustic Energy of Single Collapsing Vapor Bubble near a Rigid Wall》的全面详尽分析

---

1. 元数据与概览

• 报告标题: Numerical Study of the Emission of Acoustic Energy of Single Collapsing Vapor Bubble near a Rigid Wall

- 作者: Hangbo Duan, Linya Chen, Xiaoyu Liang

• 发布机构: College of Metrology and Measurement Engineering, China Jiliang University

- 发布日期: 2022年2月2日

• 刊物来源: Water, 2022, 第14卷，文章编号455

- 研究主题: 本文聚焦单个蒸气空穴泡在刚性壁面附近塌缩时声能释放的数值模拟研究，重点探讨空穴泡与刚性壁之间的无量纲距离参数$\gamma$对声波辐射能量的影响机理，旨在为防止空穴侵蚀的数值预测与工业应用提供理论基础。

核心论点与目标信息：
本文建立了基于OpenFOAM开源平台的可压缩两相流模型，加入了相变和热力学影响，实现了对单个空穴泡塌缩过程及其声能释放的高精度仿真。在对比实验结果验证可靠性的基础上，分析了不同无量纲距离$\gamma=L/R0$对声能幅度的影响，发现$\gamma$值在0.1~0.9区间内声能幅度随距离减小而升高，在0.9~2.0区间则相反，且声能峰值可占流场总能量的8%~25%。该规律有助于理解空穴侵蚀的物理机理，为工程设备设计提供数值预报支持[page::0,1,2,8,9]。

---

2. 逐节详细分析

2.1 引言

• 章节内容：介绍了空穴侵蚀在水轮机等液压设备中的普遍存在及其危害，阐述了传统实验方法的高成本与局限性。回顾了国内外针对空穴侵蚀机理的两个主流模型：微射流冲击理论和基于能量平衡的冲击波能量理论。作者倾向于能量平衡方法，因其适用于工程中非定常雷诺时均RANS求解器，便于工业设备的侵蚀风险估算[page::0]。
• 论点分析：引言清晰铺垫了研究的背景、问题重要性与已有研究方法的不足。强调数值仿真在早期设计阶段对侵蚀风险预测的成本效益，突显研究的应用价值。
• 关键引用：实验观测确认空穴泡塌缩产生高速微射流及冲击波导致表面损伤[9]，存在不同模型间机理争议[14-22]，综合微射流与冲击波对侵蚀贡献的研究仍在深化[13]。

2.2 空穴泡-刚性壁关系与声能研究现状

• 章节内容：详细分析了无量纲距离$\gamma=L/R0$对空穴泡动力学及冲击压力的影响。综述了一些经典经验规律（如Philipp等发现$\gamma<2$时产生侵蚀且在特定区间侵蚀最强），指出了现有研究对声能演化的数值分析不足，主要体现在忽略相变效应或仅限于球形空穴泡的压缩塌缩模型[20,30-35]。
• 逻辑与假设：作者指出传统方法未充分考虑液体的压缩效应及相变，这些因素对声能释放及能量分布有显著影响，从而为本研究开展压缩两相流模型并引入相变奠定基础[page::1]。

2.3 研究目标及结构

• 内容摘要：阐释本研究主要目的为开发包含相变和压缩效应的数值模型，重点研究塌缩阶段声能释放与$\gamma$参数的关系，不涉及泡生长阶段的墙面影响。文章结构安排包括物理模型描述、数值实现、模型验证及结果分析[page::2]。

---

3. 数值模型部分详解（第2章）

3.1 控制方程与模型假设

• 模型假设：泡的动力学以惯性力和压缩性为主，液体粘性、热力学效应、相变及界面表面张力均被纳入考虑，重力忽略。流体采用VOF方法，液蒸气两相均作为混合均质物质处理[page::2]。
• 控制方程：

- 连续性方程及动量方程详尽列出，含混合密度$\rho$、速度$\mathbf{U}$、压力$p$、粘度$\mu$及界面张力项$\sigma\kappa\mathbf{n}$。
- 流体体积分数$\alphai$定义，密度和粘度均为体积分数加权。
- 质量方程兼顾相变质量交换率$\dot{m}$，将液相蒸发和蒸气凝结纳入动力学描述，实现动态相变过程模拟[page::2-3]。

• 相变模型：引入Schnerr-Sauer模型，具体公式详述了凝结率与蒸发率依赖泡核半径、体积分数及局部压力差。相核密度$n$和Antoine方程用于确定饱和蒸汽压$pv$和温度关系，体现热力学耦合[page::3]。
• 状态方程和能量方程：

- 液相采用线性弹性模型，其中密度与压力正相关，通过声速$cl$连接；
- 蒸气泡压力假定不变，基于Rayleigh-Plesset方程简化假设；
- 能量方程包含热容和导热项，考虑了液、蒸气热物性参数，保留动能影响，较全面描述热交换和动力学耦合[page::3]。

• 能量定义：定义泡的势能$E{pot}=\frac{4}{3}\pi R^3 \Delta p$，声能为积分形式$E{wave}=\int \frac{\Delta p^2}{(\rho cl)^2} dV$，用于后续计算声能释放量[page::4]。

3.2 数值实现

• 方法论：

- 采用有限体积法，VOF方法捕捉两相界面，基于OpenFOAM v1906定制开发“compressibleInterPhaseChangeFoam”求解器；
- 模型假设轴对称，计算域为楔区域，局部网格细化保证关键区域精度，网格生成及细化策略具体说明，最大限度保证计算的准确性与效率[page::4]。

• 数值格式与边界条件：

- 时间导数二阶欧拉隐式格式，空间对流项采用基于van Leer的MUSCL格式，粘性项二阶中心差分；
- 压力-速度耦合用PISO算法；
- 边界条件详述，墙面无滑移，压力与液相体积分数采取零法线梯度，远场压力固定[page::4]。

---

4. 验证与结果分析（第3章）

4.1 网格独立性与模型验证

• 泡半径演化对比：

- 采用雷利-普莱斯特方程(Rayleigh–Plesset)解析解验证，设置$\gamma=50$近似自由场泡塌缩[page::5]。
- 不同网格（45~250个径向单元）结果比较，发现相较解析解，计算结果在塌缩末段有一定时间延迟，主因模型考虑了热传递和相变，真实泡塌缩受阻导致延迟；
- 网格细化至150个单元后延迟趋于稳定，确认选用150网格进行后续模拟[page::5]。

• 压力场验证：

- 利用Supponen等经验表达式$P{max} \propto (\rho c{l}^2 \Delta p)^{1/2} (R0/d)^{\beta} \zeta^{(2\beta/3 -1)}$，其中$\beta=2$适用于近场[page::6]；
- 数值结果中泡心附近不同点峰值压力与理论关系曲线拟合良好，说明模型对压力场的准确捕捉能够有效支持声能分析[page::6]。

4.2 塌缩过程模拟与声能释分析

• 泡塌缩过程呈现（$\gamma=1.1$示例，图4[page::7]）：

- 初期泡形近球形受压开始塌缩，末期因刚性壁影响底部塌缩较缓慢，形成椭圆形态；
- 顶部喷射产生高速微射流，“心形”泡形成，微射流冲击泡下表面产生首次声能释放；
- 泡塌缩至最小体积时，内部压差峰值促使第二次更强的球形声波发射，声波辐射至壁面并反射。

• 声能演化规律（图5）：

- 声能在塌缩绝大多数时间保持低水平，最小体积时剧烈上升[page::8]；
- 不同$\gamma$值下声能峰值及上升时间表现不同；
- 关键点在$\gamma=0.9$附近，$\gamma<0.9$时声能峰值随距离减小增强，$\gamma>0.9$则相反[page::8]。

• $\gamma$与声能峰值的关系（图6）：

- 最大声能随$\gamma$在0.1~0.9区间呈下降趋势，0.9~2.0区间呈上升趋势，曲线中值为临界点[page::9]；
- 归一化气泡塌缩时间及声能上升时间与声能峰值呈反向变化。

• 物理机理解释：

- $\gamma>0.9$时泡距墙面超过临界水平，壁面对气泡动力学干扰减弱，微射流强度减弱，局部径向流动增强，气泡塌缩更剧烈，声波释放更强；
- $\gamma<0.9$时气泡近半球形，环形涡流体累计于泡周边减弱泡压缩性，导致声能释放减弱；
- 但随着距离进一步减小，涡流强度下降，泡压缩加剧，促使声能增强，形成非单调性趋势[page::8-9]。

• 涡流及压力场对比图（图7）：

- $\gamma=0.9$时涡环明显，涡流缓冲效应显著，声能减弱；
- $\gamma=0.3$时涡环减弱，压缩效应增强，造成声能峰值提升。

---

3. 图表深度解读

图1（页4）

• 描述：泡与刚性壁的二维轴对称计算区域示意图，局部网格细化区域从泡的中心向外辐射2R至4R，分辨率高以捕捉界面细节。

- 数据意义：保证关键泡塌缩及喷射过程中速度和压力梯度的准确计算，为声能计算精度奠定基础。

• 文本联系：辅助数值模型实现阶段，确保后续结果有效性[page::4]。

图2（页5）

• 描述：不同网格密度下泡径向半径$R^$随时间演化曲线与Rayleigh-Plesset解析解对比。

- 趋势及意义：
- 初期泡收缩趋势高度一致；
- 塌缩晚期数值结果由于热传递和相变效应显示时延，且细网格减少误差；
- 体现模型物理真实性及计算稳定性。

• 支持文本结论：选用150单元后结果收敛，保证计算均衡[page::5]。

图3（页6）

• 描述：泡塌缩近场不同点最大压力$P{max}$与函数$(\rho cl^2 \Delta p)^{1/2}(R0/d)^\beta \zeta^{(2\beta/3 -1)}$的关系曲线。

- 数据解读：数值点（黑色）与拟合曲线（红色）吻合良好，验证模型压强场捕捉准确。

• 文本联系：压力场正确性直接影响声压和声能计算，确保数值声能可靠[page::6]。

图4（页7）

• 描述：$\gamma=1.1$时泡塌缩的速度场和压力场云图序列展现，从初始均匀到喷射穿透、泡体变形、声波发射及反射全过程。

- 趋势与物理意义：
- 泡形变非对称化导致微射流产生；
- 两段主要声波辐射显著，声波传播和反射机制清晰；

• 关联章节：可视化了泡动力学过程，为声能峰值形成的物理机理提供形象支撑[page::7]。

图5（页8）

• 描述：不同$\gamma$值下声能归一化曲线，纵轴为$E_{wave}^$，横轴为归一化时间$t^*$。

- 趋势说明：
- 初期低声能，塌缩末期声能急剧上升；
- 曲线最大峰值和上升时间随$\gamma$非单调变化，呈现复杂的物理耦合效应。

• 文本支撑：验证声能与距离敏感性及不同泡动态过程的对应性[page::8]。

图6（页9）

• 描述：$\gamma$与最大声能峰值、峰值时间及泡塌缩时间的关系分布图。

- 趋势解读：
- 最大声能（红点）呈“U”形随$\gamma$变化；
- 时间指标（黑三角和蓝倒三角）走势相反；

• 物理关联：对应环流结构与压缩动力学的共同作用规律，揭示泡-壁间相互作用复杂性[page::9]。

图7（页9）

• 描述：对比$\gamma=0.9$和$\gamma=0.3$的气泡周围速度场和压力场云图。

- 差异解读：
- $\gamma=0.9$时涡环明显，压缩减弱；
- $\gamma=0.3$时涡环减弱，压力集中度提高，声能更强。

• 意义：定性说明声能幅度背后复杂的流体动力学和压缩耦合机制[page::9]。

---

4. 估值分析

本研究未涉及传统金融报告中的估值概念，但在数值模拟领域，其“估值”可类比为模型对声能及压力的预测能力。通过引入：

• 基于能量守恒的声能计算公式，结合精细网格和物理本构方程（VOF，Schnerr-Sauer模型等）；

- 压缩流体与相变影响参数偏差分析；

• 以及 与解析解及实验参数经验式相符的压力峰值验证，

形成了声能预测的“估值基础”，保证了最终数据的物理合理性及仿真精度。

---

5. 风险因素评估（模型局限性与假设）

• 相变模型简化假设：泡内压力保持不变，忽略热力学驱动下的复杂内部泡动力学，可能影响极端条件下吸声波产生的预测。

- 轴对称假设的局限：现实泡塌缩具有三维非对称性，轴对称简化可能丢失局部细节喷射行为。

• 无重力假设：适合微小泡的短暂塌缩过程，非适用所有工业条件。

- 声能积分计算基于压力变化平方，对瞬态强冲击声波计算仍有一定数值误差隐患。

• 数值时间延迟分析表明热传递和相变效应可能导致泡塌缩时序与现实存在差异。

尽管存在上述潜在风险，模型通过实验对比验证压力峰值及泡径演化趋势，一定程度缓解了风险，对工程应用仍具较高参考价值[page::5,6,9]。

---

6. 批判性视角与细微差别

• 模型假设较多且理想化，尤其对泡内恒压的假设，可能在强非线性膨胀和塌缩中低估内部动力学复杂性。

- 涡环形成与声能关系的解释趋于定性，缺乏定量机制模型，未来研究可加强流体-声场耦合机理解析。

• 验证案例主要集中在单泡孤立条件，未涵盖多泡群体协同效应，现实工业环境更复杂。

- $\gamma=0.9$临界点的物理机制尚未完全明晰，存在探讨空间。

• 受计算资源限制，时间步长与空间网格虽经过优化，但极端细节可能未充分捕捉。

总体上，报告结构缜密，验证充足，但对复杂现象的解释多基于数值观察，缺少更深层的理论推导。

---

7. 结论性综合

该报告通过基于OpenFOAM的可压缩两相流数值模拟方法，首次系统研究了单个空穴泡在刚性壁面附近塌缩时声能释放的空间距离依赖特性，具体总结如下：

• 模型开发：创新引入相变动力学与热力学耦合，采用VOF方法捕捉界面，确保压力场及声能计算的物理真实性和数值稳定性。
• 验证阶段：针对网格独立、泡径演化与压力峰值进行严格对比，选取合理网格参数保证模型精度，验证数值压力场与经验公式吻合良好。
• 动力学观察：明确泡塌缩过程中的非对称结构形成、微射流生成与两个阶段声波释放过程，形象展示复杂流体-结构相互作用。
• 声能规律揭示：

- 不同无量纲距离$\gamma$对声能释放有显著影响；
- 声能峰值及上升时间呈现非单调关系，临界$\gamma\approx0.9$左右分割两种趋势；
- 物理机制涉及涡环结构对泡压缩性的缓冲及微射流强度的竞争影响。

• 图表深度支持：各图详细说明了模型的实现细节、压力及速度场演变、能量变化过程及$\gamma$ 相关声能差异，直观体现研究成果的科学性和可信度。
• 工程应用价值：研究成果有助于通过数值手段早期评估水轮机叶片等设备的空穴侵蚀风险，提升设计阶段优化的有效性，降低开发成本。

综上，本文为空穴泡塌缩声能释放规律的数值研究提供了坚实的理论与计算基础，对深入理解空穴侵蚀机制和工业应用均具有重要意义[page::0-9]。

---

总体评价

报告内容详实，理论与数值方法结合紧密，验证严谨，图表解析全面，结果呈现明晰，满足科研报告高标准要求。尤其将相变、压缩效应纳入数值建模，填补了该领域关于非球形泡塌缩声能释放的研究空白。报告虽存在部分简化假设的局限，但整体科研价值显著，值得行业研究与应用推广参考。

---

以上为该篇报告的系统性、细致化技术解读与分析，涵盖论文重点段落、理论模型、数值方法、图表数据及结论性论点，力求完整呈现原文全貌，提供严谨学术视角的专业评析。

报告