力学学报, 2021, 53(4): 944-961 DOI: 10.6052/0459-1879-20-357

流体力学

水中高压脉动气泡与浮体流固耦合特性研究1)

胡振宇, 曹卓尔, 李帅,2), 张阿漫

哈尔滨工程大学船舶工程学院, 哈尔滨 150001

FLUID-STRUCTURE INTERACTION BETWEEN A HIGH-PRESSURE PULSATING BUBBLE AND A FLOATING STRUCTURE1)

Hu Zhenyu, Cao Zhuoer, Li Shuai,2), Zhang Aman

College of Shipbuilding, Harbin Engineering University, Harbin 150001, China

通讯作者: 2)李帅, 副教授, 主要研究方向: 气泡动力学, 流固耦合动力学. E-mail:lishuai@hrbeu.edu.cn

收稿日期: 2020-10-16   接受日期: 2021-03-5   网络出版日期: 2021-04-08

基金资助: 1)国家自然科学基金.  51709056
国家自然科学基金.  51979049
中央高校基本业务费.  3072020CFJ0105
黑龙江省博士后科研启动金.  LBH-Q20016

Received: 2020-10-16   Accepted: 2021-03-5   Online: 2021-04-08

作者简介 About authors

摘要

本文针对水中放电气泡与水面浮体流固耦合作用开展实验和数值研究, 采用边界积分法对气泡运动进行数值模拟, 利用辅助函数法提高非线性流固耦合问题的计算精度, 同时运用双节点法保证气-液-固三相交界线的计算稳定性. 实验中, 采用水下放电技术生成气泡, 使用高速摄影捕捉气泡动力学行为与浮体运动响应. 首先对比数值与实验结果, 二者吻合良好, 验证了数值计算模型的有效性和正确性. 然后通过对气泡与浮体的无量纲距离$\gamma_{s} $ (气泡最大半径为特征长度)进行系统研究发现: (1) $\gamma_{s} $从0.2增大至2时, 气泡在坍塌阶段分别形成了颈缩型环状射流($0.2\leqslant \gamma_{s} \leqslant 0.3)$、接触射流($0.4\leqslant \gamma_{s} \leqslant 0.6)$、非接触射流($0.7\leqslant \gamma _{s} \leqslant 1)$、对射流($1.1\leqslant \gamma_{s} \leqslant 1.3)$和反射流($1.4\leqslant \gamma_{s} \leqslant 2)$等5种典型射流模式; (2)正射流速度随$\gamma_{s} $先增大后减小再增大, 并且当$0.7\leqslant \gamma_{s} \leqslant 0.9$时, 速度可达约1000 m/s; 反射流速度随$\gamma_{s} $增大而增大; (3)在本文实验条件下, $\gamma_{s} <1.5$时浮体对气泡的Bjerknes吸引力强于自由液面的Bjerknes排斥力导致气泡在坍塌阶段向浮体迁移; 当$\gamma_{s} \geqslant 1.5$时自由液面对气泡的排斥作用更强, 气泡在坍塌阶段远离自由液面.

关键词: 气泡 ; 流固耦合 ; 浮体 ; 边界积分法 ; 高速摄影

Abstract

This paper experimentally and numerically investigates the fluid-structure interaction between a spark-induced bubble and a floating structure. The boundary integral method is adopted to simulate the bubble dynamic behaviors and the auxiliary function method is used to improve the computational accuracy of the nonlinear fluid-structure interaction. The double-node method is employed to maintain the computational stability of the gas-liquid-solid interaction line. Besides, we use the underwater electric discharge technique to generate bubbles and the high-speed photography to record the bubble dynamics and the structural responses. Firstly, we compare the numerical result with the experimental data and favorable agreement is achieved which validates this numerical model. Through parametric study with respect to the dimensionless distance $\gamma _{s} $ from the initial bubble center to the floating structure (the reference length is the maximum bubble radius), we then find that (1) as $\gamma_{s} $ increases from 0.2 to 2, five types of jetting pattern such as necking together with annular jet ($0.2\leqslant \gamma_{s} \leqslant 0.3)$, contacting jet ($0.4\leqslant \gamma_{s} \leqslant 0.6)$, non-contacting jet ($0.7\leqslant \gamma_{s} \leqslant 1)$, collision of a jet directed towards the floating body and a counter-jet ($1.1\leqslant \gamma_{s} \leqslant 1.3)$ and individual counter-jet ($1.4\leqslant \gamma_{s} \leqslant 2)$ can be formed; (2) it is also found that the velocity of the jet directed towards the structure first increases, then decreases and finally increases again as $\gamma_{s} $ increases; additionally, it may be in the order of $\sim$1000m/s when $\gamma _{s} $ varies from 0.7 to 0.9; as $\gamma_{s} $ increases, the counter-jet velocity increases; (3) under the conditions of the presented experiments, the bubble migrates towards the floating structure when $\gamma_{s} <\mbox{1.5}$ due to the stronger Bjerknes attraction of the floating structure than the Bjerknes repellence of the free surface on the bubble during the collapsing phase. When $\gamma_{s} \geqslant \mbox{1.5}$, however, the free surface has stronger effects on the migratory behavior of the bubble than the floating structure which causes the bubble to migrate away from the free surface at the collapse stage.

Keywords: bubble ; fluid-structure interaction ; floating structure ; boundary integral method ; high-speed photography

PDF (19341KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

胡振宇, 曹卓尔, 李帅, 张阿漫. 水中高压脉动气泡与浮体流固耦合特性研究1). 力学学报[J], 2021, 53(4): 944-961 DOI:10.6052/0459-1879-20-357

Hu Zhenyu, Cao Zhuoer, Li Shuai, Zhang Aman. FLUID-STRUCTURE INTERACTION BETWEEN A HIGH-PRESSURE PULSATING BUBBLE AND A FLOATING STRUCTURE1). Chinese Journal of Theoretical and Applied Mechanics[J], 2021, 53(4): 944-961 DOI:10.6052/0459-1879-20-357

引言

脉动气泡广泛地存在于自然界中, 并且被应用于水下爆炸[1-2]、空化空蚀[3-4]、超声清洗[5]、海底资源勘探[6-7]和医疗[8-9]等领域. 迄今为止, 前人对气泡动力学的研究主要集中于自由场[10]、壁面[11-12]和自由液面[13-14]等简单边界条件下的气泡; 对复杂边界下的气泡尤其是自由液面-浮体-气泡三者非线性耦合作用的研究相对匮乏. 气泡与近水面结构物的非线性耦合作用是海洋工程和军事研究领域的重要科学问题. 浮体在气泡脉动压力的影响下运动并伴随着自由液面大变形, 这种运动响应反过来又对气泡动力学特性产生影响. 目前气泡-浮体耦合作用中所蕴含的力学机理仍未完全揭示, 因此本文针对该问题进行数值研究, 同时开展放电气泡实验验证数值模型, 旨在增强对气泡与近水面结构物非线性流固耦合效应的认识, 为相关领域的研究和应用提供基础性支撑.

Rayleigh[10]最早基于球状假设提出著名的Rayleigh-Plesset(RP)方程, 随后Herring[15]和Keller等[16]在RP方程的基础上进一步考虑流体的黏性和可压缩性的影响, 得到各种不同形式的改进RP方程. 近边界(自由液面、刚/弹性壁面、其他气泡界面)脉动气泡往往无法保持球状, 这种非球型脉动极大地增加了理论研究的难度. 在实验研究方面, 研究者使用水下爆炸实验[17-19]、放电气泡实验[20-21]和激光气泡实验[22-23]来研究气泡与各种边界的耦合作用, 高速摄影技术也被用于捕捉气泡的动力学行为. Lauterborn和Bolle[24]使用激光聚焦产生单个气泡以研究气泡与刚性壁面的耦合作用, 发现在气泡坍塌后期产生指向壁面的高速射流, 其速度可达120 m/s. 高速水射流与结构表面的砰击被认为是空蚀的重要成因[11,25]. 然而由于实验摄影中时空分辨率的限制, 精细地捕捉射流内部形态仍然是一个挑战. 通过数值模拟, 能得到更多的流场信息以分析气泡与流场边界的非线性耦合作用.

Chahine等[26]最早针对放电气泡与柱形浮体的流固耦合作用进行了实验和数值研究, 分别使用刚固浮体与可动浮体进行研究发现浮体的运动响应对气泡坍塌模式、射流形成和气泡脉动周期具有重要影响. Klaseboer等[27-28]对水下爆炸气泡与漂浮船体的耦合作用进行了数值研究, 发现气泡射流方向受到自由液面、船体表面曲率和重力的影响. 在数值计算中使用了镜像法来模拟自由液面, 该方法对于$\gamma_{f}>1.4(\gamma_{f} $为无量纲气泡与自由液面距离, 取气泡最大半径为特征长度)的工况, 精度较好, 但对于近距离的强耦合工况($\gamma_{f} <1.4)$需要考虑全非线性自由液面条件. Zong等[29]使用边界积分法(BIM)联合有限元法建立了气泡与刚性/可变形浮体耦合模型, 考虑全非线性自由表面条件以及结构的运动和弹性变形, 并采用传统的差分法求解结构表面压力(伯努利方程中${{\partial \phi }/{\partial t}}$项), 然而这种方法可能不适用于结构运动很快的强耦合工况[30-31]. Li等[30]首次将辅助函数法应用于气泡与水中悬浮球体的非线性耦合效应研究中, 相比传统差分法, 计算更为准确, 本文将该方法拓展至气泡与浮体的流固耦合特性研究中. 近年来, 有限体积法[32]、有限元法[33]等方法也被成功的应用于气泡动力学的数值研究中, 相比于这些域离散方法, BIM仅需对流场边界进行离散, 这种降维处理使BIM计算效率具有极大的优势.

本文联合BIM与辅助函数法建立脉动气泡与浮体全非线性耦合模型. 首先开展了一系列放电气泡与柱形浮体耦合作用机理实验, 然后对比数值与实验结果, 验证计算模型与方法的有效性. 采用本文的数值模型系统地研究无量纲气泡与浮体距离参数$\gamma_{s} $对气泡动力学特性和结构响应的影响规律. 本文的研究内容安排如下; 第二章给出本文数值模型的基本理论; 第三章首先验证数值模型然后分析讨论实验与数值结果, 主要关注气泡射流模式、射流速度、气泡迁移以及浮体运动响应; 第四章给出本文研究的结论.

1 基本理论和数值模型

1.1 边界元法

本文采用基于势流理论的轴对称边界元法, 建立脉动气泡与浮体非线性流固耦合作用数值模型. 基于流体无黏、无旋、不可压缩假设, 流场(水域)可由拉普拉斯方程描述

$\begin{eqnarray} \nabla^2 \phi = 0 \end{eqnarray} $

上式可转化为积分形式的边界积分方程[34-35]

$\begin{eqnarray} \label{eq2} &&A({p},t)\phi ({p},t)=\iint_S {\frac{\partial \phi ({q},t)}{\partial n}} \frac{1}{\left| {{p}-{q}} \right|}{d}S-\\ && \iint_S {\phi ({q},t)\frac{\partial }{\partial n}} \left( {\frac{1}{\left| {{p}-{q}} \right|}} \right){d}S \end{eqnarray} $

式中, ${p}$为控制点, ${q}$为源点, $A$为${p}$点观察流场(水域)的立体角, $S$为包括气泡表面$S_{b} $、自由液面$S_{f} $、浮体湿表面$S_{w} $以及无穷大球面$S_{\infty } $的流体边界, $n$表示边界法向, 以指向流场外为正. 图1所示为气泡与浮体非线性流固耦合作用数值模型示意图, 引入柱坐标系$O(r,\theta ,z)$, 其坐标原点位于初始气泡中心, $z$轴指向重力的反方向且为对称轴, 定义浮体的入水深度为$d$,初始气泡中心至初始自由液面的距离为$h$; 本文考察圆柱型浮体与气泡的耦合作用, 圆柱浮体半径为$R_{s} $, 气泡最大等效半径为$R_{m} $.

图1

图1   气泡与浮体非线性流固耦合作用数值模型示意图

Fig.1   Sketch of the numerical model for investigating the nonlinear fluid-structure interaction between a bubble and a floating structure


离散气泡、浮体和自由液面边界, 边界积分方程转化为一个线性方程组, 求解此线性方程组后可获得气泡和自由液面节点的法向速度与浮体湿表面节点速度势, 通过有限差分可获得气泡与自由液面的切向速度从而合成总速度[36].

对浮体附近脉动气泡动力学特性的时域模拟需显式更新边界表面节点位置和速度势. 给出气泡表面与自由液面满足的运动学边界条件如下

$\frac{{d}{x}}{{d}{t}}=\nabla \phi $

其中$x$为节点位置矢量.

气泡表面和自由液面满足的动力学边界条件可由非定常伯努利方程给出, 分别为

$\frac{{d}\phi }{{d}t}=\frac{p_{\infty } -p_{b}}{\rho }+\frac{1}{2}\left| {\nabla \phi } \right|^{2}-gz $
$\frac{{d}\phi }{{d}t}=\frac{1}{2}\left| {\nabla \phi } \right|^{2}-g(z-h) $

式中, $p_{\infty } $为$z=0$平面内流场中无穷远处静水压力, $p_{b} $为气泡内压, $g$表示重力加速度, $\rho $为流体密度. 气泡与自由液面满足的动力学边界条件用于更新节点速度势. 由于气泡内部不可冷凝气体在脉动过程中满足理想气体绝热定律[37], 忽略表面张力作用, $p_{{b}} $可由下式计算

$\begin{eqnarray} \label{eq6} p_{b} =p_{{v}} +p_{0} \left( {\frac{V_{0}}{V}} \right)^{\kappa } \end{eqnarray} $

式中, $p_{{v}} $为气体的饱和蒸汽压, 本文取2337 Pa; $p_{0}$为初始气泡内部不可冷凝气体压力, $V_{0} $为初始气泡体积, $V$为瞬时气泡体积; $\kappa $为气泡内部不可冷凝气体的比热比, 本文取$\kappa =1.4$.

速度势$\phi $需满足浮体湿表面不可穿透条件如下

$\nabla \phi \cdot {n}={V}_{s} \cdot {n} $

式中, ${n}$表示边界单位法向量, ${V}_{s}$表示浮体运动速度.

1.2 辅助函数法

通过牛顿第二定律

$\begin{eqnarray} \label{eq8} m{a}={F}_{{h}} +{F}_{{a}} +{G} \end{eqnarray} $

可求得浮体加速度${a}$. 式中$m$为浮体质量, ${G}$为浮体重力, 空气压力${F}_{{a}} $通过在浮体干表面上对压力(大气压$p_{{atm}} )$积分即可获得, 水动力${F}_{{h}} $通过伯努利方程获得[29]

$\begin{eqnarray} \label{eq9} F_{{h}} =-\rho \iint_{S_{{w}}} {\left( {\frac{\partial \phi }{\partial t}+gz+\frac{1}{2}\left| {\nabla \phi } \right|^{2}-\frac{p_{\infty }}{\rho }} \right)}n{d}S \end{eqnarray} $

传统方法通过对相邻时间步的节点速度势采取向前差分近似求解${{\partial \phi }/{\partial t}}$项, 但是时间步过小或浮体运动剧烈时, 精度不足. 理论上这种方法得到的是${{{d}\phi }/{{d}t}}$, 可通过物质导数与局部导数的关系得到修正的差分法即${{\partial \phi }/{\partial t}}={{{d}\phi }/{{d}t}}-{V}_{s} \cdot \nabla \phi $. 本文采用的辅助函数法可精确求解结构响应, 后面将对比修正后的传统差分法、未修正的差分法和辅助函数法的计算精度.

本文中浮体仅存在$z$向运动, 作用在浮体上的外力沿$z$轴, 故可将式(8)改写为

$\begin{eqnarray} \label{eq10} && (m+\rho \iint_{S_{w}} {\chi n_{z} {d}S} )a_{z} +mg=\\&& -\rho \iint_{S_{w}} \left(\xi -V_{s} \cdot \nabla \phi +\frac{1}{2}\left| {\nabla \phi } \right|^{2}+gz- \frac{p_{\infty }}{\rho }\right)n_{z}{d}S-\\&& p_{{atm}} \iint_{S_{d}} {n_{z} {d}S} \end{eqnarray} $

上式左侧第一项代表惯性力项(其中积分项代表附加惯性力), 左侧第二项代表重力, $\chi $与$\xi $为两个辅助函数, $S_{{d}} $为浮体与大气接触的表面.

在气泡表面与自由液面, $\chi $满足边界条件[30]

$\begin{eqnarray} \label{eq11} \chi =0 \end{eqnarray} $

在结构湿表面, $\chi $满足边界条件如下[38]

$\begin{eqnarray} \label{eq12} \frac{\partial \chi }{\partial n}=n_{z} \end{eqnarray} $

$\xi $在气泡表面与自由液面满足的边界条件可分别表示为[30]

$\xi ={V}_{s} \cdot \nabla \phi +\frac{p_{\infty } -p_{{v}}}{\rho }-\frac{p_{0}}{\rho }\left( {\frac{V_{0}}{V}} \right)^{\kappa }-\frac{\left| {\nabla \phi } \right|^{2}}{2}-gz $
$\xi ={V}_{s} \cdot \nabla \phi -\frac{\left| {\nabla \phi } \right|^{2}}{2}-g(z-h) $

$\xi $在浮体湿表面满足边界条件[31]

$\begin{eqnarray} \label{eq15} \frac{\partial \xi }{\partial n}=0 \end{eqnarray} $

辅助函数$\chi $与$\xi $满足拉普拉斯方程, 并已知其在气泡和自由液面满足的狄利克莱边界条件与在浮体湿表面上满足的诺伊曼边界条件, 采用类似于求解边界积分方程的方式可以得到在浮体湿表面上的值, 代入式(10)中可求解出浮体加速度.

1.3 环状气泡模型

射流砰击时, 原来的单连通气泡转变为双连通气泡, 即变为环状气泡. 此时气泡射流砰击处的速度势不连续. 为保证计算的稳定性, 需要解决砰击处节点的速度势突跳问题. 本文采用Wang等[35]提出的涡环模型, 在气泡内部布置一个涡环, 其环量为砰击处的速度势突跳. 将气泡表面速度势分解为[34]

$\begin{eqnarray} \label{eq16} \phi =\phi_{{vor}} +\phi_{{res}} \end{eqnarray} $

其中, $\phi_{{vor}} $为涡环诱导速度势, $\phi_{{res}}$为残余速度势. 残余速度势$\phi_{{res}} $在环状气泡表面连续, 采用边界积分法进行计算, 涡环诱导速度势采用一种半解析法求解. 详细的计算过程可参考文献[39], 本文不再赘述, 这里仅对涡环模型的速度势分解思想进行介绍.

1.4 三相接触模型

对于水面漂浮结构存在气-液-固三相接触线, 如图2所示, 接触线需要特殊处理以同时满足自由液面狄利克莱边界条件与浮体表面诺伊曼边界条件[40]. 本文采用双节点法追踪交界线, 这种方法同样适用于气泡与结构接触问题. 将接触点$k$分解为$k_{1}$和$k_{2}$两点, 分别属于自由液面和柱体; 使用下标f$_{k_{1}} $和s$_{k_{2}} $分别表示属于自由液面、浮体的交界点. 下标b, f和s分别表示除了交界点外气泡、自由液面和浮体上的节点. 下标$n$表示速度势的法向导数.

图2

图2   三相接触模型示意图

Fig.2   Sketch of the solid-liquid-gas contact model


离散边界积分方程(2)后可得到如下线性方程组

$\begin{eqnarray} \label{eq17} &&\left[ {{\begin{array}{*{20}c} {H_{bb}} & {H_{bf}} & {H_{{bf}_{k_{1}}}} & {H_{{bs}_{k_{2}}}} & {H_{bs}} \\ {H_{fb}} & {H_{{ff}}} & {H_{{ff}_{k_{1}}}} & {H_{{fs}_{k_{2}}}} & {H_{fs}} \\ {H_{{f}_{k_{1}} {b}}} & {H_{{f}_{k_{1}} {f}}} & {H_{{f}_{k_{1}}{f}_{k_{1}}}} & {H_{{f}_{k_{1}} s_{k_{2}}}} & {H_{{f}_{k_{1}} {s}}} \\ {H_{{s}_{k_{2}} {b}}} & {H_{{s}_{k_{2}} {f}}} & {H_{{s}_{k_{2}} f_{k_{1}}}} & {H_{{s}_{k_{2}} {s}_{k_{2}}}} & {H_{{s}_{k_{2}} {s}}} \\ {H_{sb}} & {H_{{sf}}} & {H_{{sf}_{k_{1}}}} & {H_{{ss}_{k_{2}}}} & {H_{{ss}}} \\ \end{array}}} \right]\left[ {{\begin{array}{*{20}c} {{\phi }_{b}} \\ {\phi_{f}} \\ {\phi_{{f}_{k_{1}}}} \\ {\phi_{{s}_{k_{2}}}} \\ {\phi_{s}} \\ \end{array}}} \right]=\\&& \left[ {{\begin{array}{*{20}c} {{G}_{bb}} & {{G}_{bf}} & {{G}_{{bf}_{k_{1}}}} & {{G}_{{bs}_{k_{2}}}} & {{G}_{bs}} \\ {{G}_{fb}} & {{G}_{{ff}}} & {{G}_{{ff}_{k_{1}}}} & {{G}_{{fs}_{k_{2}}}} & {{G}_{fs}} \\ {{G}_{{f}_{k_{1}} {b}}} & {{G}_{{f}_{k_{1}} {f}}} & {{G}_{{f}_{k_{1}} f_{k_{1}}}} & {{G}_{{f}_{k_{1}} s_{k_{2}}}} & {{G}_{{f}_{k_{1}} {s}}} \\ {{G}_{{s}_{k_{2}} {b}}} & {{G}_{{s}_{k_{2}} {f}}} & {{G}_{{s}_{k_{2}} {f}_{k_{1}}}} & {{G}_{{s}_{k_{2}} {s}_{k_{2}}}} & {{G}_{{s}_{k_{2}} {s}}} \\ {{G}_{sb}} & {{G}_{{sf}}} & {{G}_{{sf}_{k_{1}}}} & {{G}_{{ss}_{k_{2}}}} & {{G}_{{ss}}} \\ \end{array}}} \right]\left[ {{\begin{array}{*{20}c} {{\phi }_{{b}{n}}} \\ {{\phi }_{{f}n}} \\ {{\phi }_{{f}_{k_{1}}n}} \\ {{\phi }_{{s}_{k_{2}}n}} \\ {{\phi }_{{s}n}} \\ \end{array}}} \right] \end{eqnarray} $

由于$k_{1}$和$k_{2}$两点在空间位置上完全相同, 影响系数矩阵$H$和$G$第3和第4行完全相同, 整理可得

$\begin{eqnarray} \label{eq18} &&\left[ {{\begin{array}{*{20}c} {{H}_{bb}} & {H_{bf}} & {H_{{bf}_{k_{1}}} +H_{{bs}_{k_{2}}}} & {-G_{{bs}_{k_{2}}}} & {-G_{bs}} \\ {H_{fb}} & {H_{{ff}}} & {H_{{ff}_{k_{1}}} +H_{{fs}_{k_{2}}}} & {-G_{{fs}_{k_{2}}}} & {-G_{fs}} \\ {H_{{f}_{k_{1}} {b}}} & {H_{{f}_{k_{1}} {f}}} & {H_{{f}_{k_{1}} f_{k_{1}}} +H_{{f}_{k_{1}} s_{k_{2}}}} & {-G_{{f}_{k_{1}} s_{k_{2}}}} & {-G_{{f}_{k_{1}} {s}}} \\ {H_{sb}} & {H_{{sf}}} & {H_{{sf}_{k_{1}}} +H_{{ss}_{k_{2}}}} & {-G_{{ss}_{k_{2}}}} & {-G_{{ss}}} \\ \end{array}}} \right]\left[ {{\begin{array}{*{20}c} {{\phi }_{b}} \\ {{\phi }_{f}} \\ {{\phi }_{{f}_{k_{1}}}} \\ {{\phi }_{{s}_{k_{2}} n}} \\ {{\phi }_{{s}n}} \\ \end{array}}} \right]=\\&& \left[ {{\begin{array}{*{20}c} {G_{bb}} & {G_{bf}} & {G_{{bf}_{k_{1}}}} & {-H_{bs}} \\ {G_{fb}} & {G_{{ff}}} & {G_{{ff}_{k1}}} & {-H_{fs}} \\ {G_{{f}_{k_{1}} {b}}} & {G_{{f}_{k_{1}} {f}}} & {G_{{f}_{k_{1}} {f}_{k_{1}}}} & {-H_{{f}_{k_{1}} {s}}} \\ {G_{sb}} & {G_{{sf}}} & {G_{{sf}_{k_{1}}}} & {-H_{{ss}}} \\ \end{array}}} \right]\left[ {{\begin{array}{*{20}c} {{\phi }_{{b}n}} \\ {{\phi }_{{f}n}} \\ {{\phi }_{{f}_{k_{1}} n}} \\ {{\phi }_{s}} \\ \end{array}}} \right] \end{eqnarray} $

求解线性方程组后可得到$k_{1}$点的法向速度.

当气泡产生位置距浮体表面很近时, 由于气泡膨胀或者气泡朝壁面迁移, 导致气泡附着到浮体表面. 若气泡节点与浮体节点距离小于单元尺度, 传统BIM产生严重的数值不稳定性; 如图3(a)所示, 部分气泡节点穿过结构表面, 产生了计算错误. 因此本文将气泡与浮体间的薄水膜去除, 连接气泡与浮体节点, 进而保证数值计算的稳定性, 如图3(b)所示. Wang等[41]的研究表明, 这层薄水膜压力趋近气泡内压, 将其去除未对气泡的动力学特性产生可见的影响. 节点去除的范围和接触判定时刻由气泡最小单元长度$\Delta s$控制, 当气泡节点与浮体节点最小距离小于$\Delta s$时, 即认为气泡即将附着在浮体表面, 并将距离小于$2.5\Delta s$[38]的气泡与浮体节点删除并连接气泡与浮体, 采用双节点法处理产生的气-液-固三相交界线. 气泡附着在浮体后, 气泡内压直接作用在浮体表面上, 此时浮体运动方程如下

$\begin{eqnarray} \label{eq19} && \left(m+\rho \iint_{S_{w}} {\chi n_{z} {d}S} \right)a_{z} +mg=\\&& -\rho \iint_{S_{w}} {\left(\xi -V_{s} \cdot \nabla \phi +\frac{1}{2}\left| {\nabla \phi } \right|^{2}+gz-\frac{p_{\infty }}{\rho }\right)n_{z} {d}S} +\\&& p_{b} \iint_{S_{g}} {n_{z} {d}S}- p_{atm} \iint_{S_{d}} {n_{z} {d}S} \end{eqnarray} $

式中, $S_{g} $为浮体与气泡接触表面.

图3

图3   传统BIM计算结果(a)与采用三相接触模型计算结果(b)对比

Fig.3   Comparison between the traditional BIM (a) and the gas-liquid-solid contact model (b)


1.5 无量纲化

本文将所有物理量进行无量纲化, 并使用无量纲量进行计算与讨论. 取$R_{m}$, $\Delta p=p_{\infty } -p_{{v}} $和$\rho $分别作为长度、压力和密度的特征量, 其他物理量如时间、速度、速度势和加速度的无量纲特征量分别为$R_{m} \left( {{\rho /{\Delta p}}} \right)^{1/2}$, $\left( {{{\Delta p}/\rho }} \right)^{1/2}$, $R_{m} \left( {{{\Delta p}/\rho }} \right)^{1/2}$和$\Delta {p/{(\rho R_{m} )}}$[34-35,37]. 本文后续出现的物理量如无特殊说明均为无量纲形式.

定义浮力参数$\delta $

$\begin{eqnarray} \label{eq20} \delta =\left( {{{\rho gR_{m}}/{\Delta p}}} \right)^{1/2} \end{eqnarray} $

强度参数$\varepsilon $

$\begin{eqnarray} \label{eq21} \varepsilon ={{p_{0}}/{\Delta p}} \end{eqnarray} $

初始气泡中心位置与自由液面距离参数$\gamma_{f} $

$\begin{eqnarray} \label{eq22} \gamma_{f} ={h/{R_{m}}} \end{eqnarray} $

初始气泡中心位置与浮体距离参数$\gamma_{s} $

$\begin{eqnarray} \label{eq23} \gamma_{s} ={{(h-d)}/{R_{m}}} \end{eqnarray} $

尺度比$\zeta $

$\begin{eqnarray} \label{eq24} \zeta =R_{s} /R_{m} \end{eqnarray} $

由于放电气泡实验中很难得到气泡内压, 并且Turangan等[42]的研究发现强度参数$\varepsilon $在100 $\sim$ 500范围内变化, 气泡射流速度变化很小, 气泡的动力学特性变化也很小. 经过计算发现, 本文选取$\varepsilon =300$, 实验与数值结果吻合良好. 计算中气泡初始半径与强度参数相关[43], 本文选取$R_{0} =0.113$. 计算中所用时间步长由下式计算得到[37]

$\begin{eqnarray} \label{eq25} \Delta t=\dfrac{C}{\max \left| {1+0.5\left| {\nabla \phi } \right|^{2}-\varepsilon \left( {\dfrac{V_{0}}{V}} \right)^{\kappa }-\delta ^{2}z} \right|} \end{eqnarray} $

其中$C$为单步最大的速度势改变量.

2 结果分析与讨论

2.1 收敛性分析

为了验证本文数值模型的正确性和有效性, 首先进行收敛性分析, 包括网格与时间步的收敛性分析. 分别使用气泡节点数量$N$等于80, 120, 160, 200和240的模型进行计算. 图4(a)展示了不同网格密度下射流砰击瞬时气泡形态, 节点数为80, 120和160的结果中气泡形态具有一定的差异, 当节点数$N\geqslant 200$时气泡轮廓基本重合, 证明当气泡网格节点取200时能获得收敛的结果.

图4

图4   模型的收敛性分析(a)计算网格的收敛性(b)计算时间步的收敛性

Fig.4   Convergence study with respect to (a) the mesh and (b) the time step


本文采用四阶Runge-Kutta法对脉动气泡与浮体相互耦合的时域运动过程进行模拟, 显式算法的收敛性依赖于时间步长. 取$C$分别为0.04, 0.02, 0.01和0.005分别进行计算. 图4(b)描述了浮体加速度随时间的变化, 4种方案取得了一致的结果, 进一步证明了本文提出的数值模型的有效性. 考虑程序计算的精度以及所需的计算成本, 本文后续计算选取$N=200$, $C=0.01$.

2.2 典型工况对比分析

本节将分别对3个典型距离参数$\gamma_{s} $下气泡与浮体非线性耦合作用进行研究和讨论, 展示实验中的物理现象, 利用数值模拟分析现象背后的力学机理. 本文实验方法详见已发表文献[31, 44].

2.2.1 弱耦合工况

第一个实验中初始气泡深度$h=74.7$ mm, 初始浮体入水深度$d=45.0$ mm, 浮体半径$R_{s} =16.0$ mm, 气泡最大等效半径$R_{m} =16.8$ mm. 根据实验设定计算初始参数为: $\varepsilon =300$, $R_{0} =0.113$, $\kappa =1.4$, $\gamma_{f} =4.45$, $\gamma_{s} =1.77$, $\zeta =0.95$, $\delta =0.04$. 图5所示为本工况数值与实验结果对比图, 在此对本文的数值与实验对比图进行说明; 图中左侧为高速摄影捕捉的气泡图像, 右侧为对应时刻气泡运动模拟结果. 每一张图片左上方显示了实验中气泡脉动时间, 右侧数字为数值计算结果对应的气泡脉动时间, 图5(a)左下方给出了比例尺, 图5(b)中箭头指出了实验图像中自由液面位置, 为了更好地分析浮体与气泡的非线性耦合效应, 数值结果中还给出无量纲压力与速度场.

图5

图5   弱耦合工况数值与实验结果对比

Fig.5   Comparison between the numerical results and the experimental images for the weakly coupling case


图5(a)中, 气泡在内部高压的推动下猛烈膨胀. 图5(b)所示为气泡膨胀达到其最大体积时刻, 此时气泡近似球状, 自由液面和浮体未对气泡脉动产生明显的影响. 由于过度膨胀气泡周围压力比静水压力小, 在接下来的脉动过程中在流体惯性的作用下气泡开始坍塌. 浮体对气泡界面收缩具有延时作用, 在气泡坍塌前期, 气泡顶部略微凸起, 而底部具有较大坍塌速度因此表面较为平坦, 如图5(c)所示. 这种非球状坍塌导致气泡顶部在坍塌中期获得更大的曲率, 如图5(d)所示. 在气泡坍塌后期, 气泡周围流体压力迅速增大, 但气泡顶部压力梯度更大, 使顶部界面向下(逆压力梯度方向)快速收缩, 如图5(e)所示. 第一个周期结束, 并未形成射流, 而在回弹阶段气泡顶部产生远离浮体的射流并穿透另一侧气泡表面, 如图5(f)所示. 气泡顶部射流现象与Lauterborn[45]的预测一致, 即气泡具有高曲率的部分更容易产生射流. 此工况中气泡产生背离浮体的射流(称为"反射流"), 说明浮体对气泡的吸引力较小, 无法诱导指向壁面的射流, 故将本工况定义为弱耦合工况. 总体上, 数值计算得到的气泡运动特性与实验吻合良好.

2.2.2 中度耦合工况

第二个实验中气泡产生位置距离浮体更近, 耦合作用更强, 气泡初始深度$h=60.0$ mm, 初始浮体入水深度$d=42.6$ mm, 浮体半径$R_{s} =16.0$ mm, 气泡最大等效半径$R_{m} =16.4$ mm. 依据实验设定计算初始参数如下; $\varepsilon =300$, $R_{0} =0.113$, $\kappa =1.4$, $\gamma_{f} =3.66$, $\gamma _{s} =1.06$, $\zeta =0.98$, $\delta =0.04$.

图6给出了本工况中几个典型时刻数值与实验结果对比. 铜丝燃烧生成一个高压气核, 在内部高压气体作用下迅速膨胀, 如图6(a)所示. 图6(b)表示气泡达到最大体积时刻, 气泡上方流体运动受到浮体的阻碍, 气泡上表面被压平, 此时气泡上方流体继续向上运动推动浮体远离气泡. 在气泡坍塌前期, 远离浮体的气泡底部近似球形坍塌, 如图6(c)所示. 柱体阻碍了流体从顶部流向气泡中心, 而对侧方的流体的运动影响较小, 因此形成了"梨"形气泡. 图6(d) $\sim\!$图6(f)描述了气泡坍塌中后期的运动, 气泡两侧以及底部收缩的更快, 顶部运动速度很小, "梨"形气泡逐渐转变为"梭"形气泡, 气泡在垂向被拉长. 在坍塌后期气泡底部压力梯度变大, 产生高压区, 使底部气泡表面快速收缩形成射流. 图6(f)的实验图像中气泡底部留下铜丝燃烧产物, 表明形成了一束向上的射流. 此工况中浮体对气泡的动力学行为影响更大, 并能诱导气泡产生指向壁面的射流(称为"正射流"), 因此将此工况命名为中度耦合工况, 从实验中也能观察到浮体在膨胀阶段被气泡排斥而在收缩阶段被吸引的现象.

图6

图6   中度耦合工况数值与实验结果对比

Fig.6   Comparison between the numerical results and the experimental images for the intermediately coupling case


2.2.3 强耦合工况

为研究浮体与脉动气泡强烈的非线性耦合作用, 第三个实验进一步减小了浮体与气泡间的距离. 此时气泡在膨胀阶段附着到浮体底面, 具有明显的流固耦合效应, 固将此工况命名为"强耦合工况". 实验中参数如下; 初始气泡深度$h=42.6$ mm, 浮体入水深度$d=36.0$ mm, 浮体半径$R_{s} =14.0$ mm, 气泡最大等效半径$R_{m} =13.1$ mm. 依据实验设定计算初始参数如下; $\varepsilon =300$, $R_{0} =0.113$, $\kappa =1.4$, $\gamma_{f} =3.25$, $\gamma _{s} =0.5$, $\zeta =1.07$, $\delta =0.04$.

图7给出了强耦合工况中几个典型时刻实验与数值计算结果对比. 如图7(a)所示, 在气泡膨胀初期, 由于柱形浮体的阻碍作用, 气泡顶部变平. 在膨胀阶段气泡上表面逐渐靠近浮体底面, 需启动三相接触模型以保证数值稳定性. 图7(b)所示为气泡与浮体表面接触时刻, 气泡即将附着到浮体表面. 图7(c)表示气泡最大体积时刻, 此时气泡顶部仍有一定的膨胀速度, 而底部与侧面在收缩, 进一步说明浮体对气泡运动具有延时的作用. 侧面流体在气泡坍塌时能自由的流向气泡中心, 顶部流体受到阻碍, 如图7(d)所示. 在接下来的坍塌过程中, 气泡底部产生高压区, 驱动射流直接砰击在浮体上. 实验中还观察到射流前进阶段气泡撕裂为顶部大气泡和底部微小气泡, 计算中并未产生撕裂现象, 如图7(e)和图7(f)所示. 图中最后两帧数值与实验结果具有较大偏差, 这可能是由于气泡脉动过程中浮体底面附着的空化气泡(见图7(c))与放电气泡产生耦合作用所致. 并且三相接触线附近气泡表面产生翻卷和破碎的现象, 这可能与固体材料的润湿性能相关, 而本文计算中未考虑接触角的影响. 总体上气泡形态与周期与实验结果吻合良好. 值得注意的是, 柱形浮体在此强耦合工况中产生了显著的运动.

图7

图7   强耦合工况数值与实验结果对比

Fig.7   Comparison between the numerical results and the experimental images for the strongly coupling case


图8给出了强耦合工况中浮体位移和浮体速度时历曲线, 分别对比了本文辅助函数法模型、传统修正后的差分模型以及实验结果. 实验结果中误差棒的中心代表4次测量的平均值, 误差棒的范围代表标准差; 菱形标记表示气泡最大体积时刻. 由图8(a)可知, 浮体在气泡膨胀阶段具有向上的运动, 由于惯性作用气泡达到最大体积后浮体仍继续向上运动. 在坍塌阶段浮体向下加速运动再逐渐减速. 可以看出, 相比于传统方法本文采用的辅助函数方法与实验结果更加一致. 在气泡膨胀早期, 数值结果中浮体运动较快, 这是由于数值计算中忽略了铜丝燃烧放热的过程; 而在后期随着时间步减小, 传统差分法结果逐渐偏离实验结果. 从图8(b)可以看出在气泡膨胀前期, 两种方法求得的结构速度基本一致, 而在后期传统差分法精度降低, 在气泡附着于浮体表面后(对气泡拓扑结构进行"数值切割"可能进一步增大了不稳定性), 速度时历曲线出现剧烈震荡(图中虚线表示气泡与浮体接触时刻), 而本文方法能得到稳定且准确的结果. 故对于研究近距离的脉动气泡与浮体非线性耦合作用, 本文提出的数值模型更为有效.

图8

图8   柱形浮体(a)位移与(b)速度时历曲线(菱形标记代表最大体积时刻, 虚线表示气泡与浮体接触时刻)

Fig.8   Comparison of (a) the displacement and (b) the velocity histories of the cylindrical floating structure (the diamonds denote the moment of reaching the maximum bubble volume, the dashed line denotes the moment of the bubble-structure attachment)


2.3 无量纲参数讨论

2.3.1 $\gamma_{s} $对气泡射流模式的影响

影响气泡与浮体非线性耦合效应的参数众多, 主要有无量纲气泡与浮体距离参数$\gamma_{s} $、无量纲气泡与自由面距离参数$\gamma_{f} $与尺度比$\zeta $. 众多的影响参数导致了这一问题的复杂性. 限于篇幅, 本文仅研究$\gamma_{s} $的影响, 其余参数与实验参数一致, 本文提出的数值模型也可研究其余参数的影响. 依据实验设定后续数值研究的初始参数; $\varepsilon =300$, $R_{0} =0.113$, $\gamma =1.4$, $d=2.6$, $\zeta =1$, $\delta =0.04$, $\gamma _{s} =0.2\sim 2$.

图9展示了不同距离参数$\gamma_{s} $下, 射流砰击时刻的气泡形态, 明显气泡射流模式分为5种类型, 分别为颈缩型环状射流($0.2\leqslant \gamma_{s} \leqslant 0.3$)、接触射流($0.4\leqslant \gamma_{s} \leqslant 0.6$)、非接触射流($0.7\leqslant \gamma_{s} \leqslant 1$)、对射流($1.1\leqslant \gamma_{s} \leqslant 1.3$)和反射流($1.4\leqslant \gamma_{s} \leqslant 2$), 图中虚线代表浮体.

图9

图9   不同无量纲气泡与浮体距离$\gamma_{s}$下气泡射流模式

Fig.9   Jetting patterns for different dimensionless distances $\gamma_{s}$


图10给出了$\gamma_{s}=0.2$ (颈缩型环状射流)工况中典型时刻的压力和速度场. 气泡收缩时牵引浮体侧面流体向下运动, 同时水平方向气泡能自由收缩, 如图10(a)所示. 水平方向流体与垂向流体在柱体底部最大半径处汇聚, 在此处形成高压区(见图10(b)和图10(c)). 在压力梯度作用下流体向斜下方运动, 形成颈缩现象. 类似于Cui等[44]的分析, 浮体的向下运动可能进一步改变流体运动, 使其偏离形成环状向外的射流导致气泡撕裂为一个大单连通气泡与小环状气泡.

图10

图10   $\gamma _{s} =0.2$工况中典型时刻压力与速度场(分别对应最大体积时刻$t=0.910$, 气泡坍塌时刻$t=1.132$, 环向射流砰击时刻$t=1.403$)

Fig.10   Pressure and velocity fields at some typical moments for $\gamma _{s} =0.2$ (corresponding times are $t=0.910$ at the moment of reaching the maximum bubble volume, $t=1.132$ during the bubble collapse stage and $t=1.403$ at the moment of the annular jet impact respectively)


图9第3 $\sim$ 5帧中气泡附着在浮体表面, 此时产生的射流称为接触射流. 图11给出了$\gamma _{s} =0.6$(接触射流)工况中典型时刻压力和速度场. 在气泡坍塌后期, 底部附近形成了局部高压区域, 驱动气泡界面快速收缩形成高速细射流直接砰击在浮体上, 如图11(a)和图11(b)所示. 射流砰击到浮体底部后, 在浮体底面产生很高的压力. 流体受到壁面的作用变为径向流动, 最终导致气泡脱离浮体表面, 如图11(c)所示. 在气泡脱离瞬时, 浮体底面中心附近压力比射流砰击时刻更高, 这是射流砰击产生的局部高压与气泡收缩产生的脉动压力叠加导致的结果.

图11

图11   $\gamma _{s} =0.6$工况中典型时刻压力与速度场(分别对应气泡坍塌时刻$t=1.792$,射流砰击时刻$t=1.810$和气泡与浮体脱离时刻$t=1.838$)

Fig.11   Pressure and velocity fields at some typical moments for $\gamma _{s} =0.6$ (corresponding times are $t=1.792$ during the bubble collapse stage, $t=1.810$ at the moment of the jet impact and $t=1.838$ at the moment of the bubble-floating structure detachment respectively)


图9第6 $\sim$ 9帧气泡射流模式称为非接触射流. 气泡底部形成的高速射流首先砰击在气泡上表面, 穿过水膜后砰击在浮体底面. 图12给出了$\gamma_{s} =0.9$(非接触射流)工况中典型时刻压力和速度场. 在坍塌阶段气泡下方压力梯度最大, 导致正射流的形成. 射流砰击在气泡上表面后, 在气泡与浮体间形成局部高压区域, 如图12(a)和图12(b)所示. 环状气泡继续收缩同时向浮体运动, 在顶部形成突出物最终与浮体底面接触. 值得注意的是射流砰击后浮体底面中心附近压力比上文接触射流工况更大. 环状气泡继续坍塌, 壁面附近压力超过140 MPa, 如图12(c)所示.

图12

图12   $\gamma _{s} =0.9$工况中典型时刻压力与速度场(分别对应气泡坍塌时刻$t=1.806$, 射流砰击时刻$t=1.807$和气泡与浮体接触时刻$t=1.809$)

Fig.12   Pressure and velocity fields at some typical moments for $\gamma _{s} =0.9$ (corresponding times are $t=1.806$ during the bubble collapse stage, $t=1.807$ at the moment of the jet impact and $t=1.809$ at the moment of the toroidal bubble-floating structure attachment respectively)


图9第10 $\sim$ 12帧展示了反射流与正射流相撞的现象(称为"对射流"). 如图13(a)所示, 在气泡坍塌末期, 气泡顶部首先产生一个宽度较大的凹陷. 气泡底部形成很高的曲率, 在其下方产生一个高压区域. 高压区使气泡底部迅速收缩形成一束指向壁面的正射流. 最终正射流与宽而慢的反射流相撞(图13(b)和图13(c)). 射流相互撞击后在撞击区域附近产生了很高的压力. 反射流内部未与正射流砰击的部分继续向下运动, 导致砰击区域环状凹陷的形成, 如图13(d)所示. 这种由射流相互撞击产生的环状凹陷无法发展成射流, 在回弹阶段逐渐消失, 如图13(e)和图13(f)所示. 在坍塌和回弹阶段, 气泡和高压区逐渐向浮体迁移.

图13

图13   $\gamma _{s} =1.3$工况中典型时刻压力与速度场(分别对应气泡坍塌时刻$t=1.806$, $t=1.807$, 射流砰击时刻$t=1.809$, 气泡坍塌时刻 $t=1.810$ 和气泡回弹时刻$t=1.820$, $t=1.853$)

Fig.13   Pressure and velocity fields at some typical moments for $\gamma _{s} =1.3$ (corresponding times are $t=1.806$, $t=1.807$ during the bubble collapse stage; $t=1.809$ at the moment of the jet impact; $t=1.810$ during the bubble collapse stage; $t=1.820$, $t=1.853$ during the bubble rebounding stage respectively)


图9第13帧中气泡产生了反向射流与径向射流, 径向射流砰击将气泡撕裂为两个单连通的气泡. 图9第14 $\sim$ 15帧中气泡只产生径向凹陷, 无法穿透气泡. 图14给出了$\gamma_{s} =1.5$(反射流)工况中典型时刻压力和速度场. 在气泡坍塌后期, 气泡顶部压力比周围压力更高导致气泡形成向下的凹陷. 然而第一个周期内, 向下的凹陷无法形成射流穿透气泡表面, 如图14(a)所示. 随后气泡两侧收缩较快形成径向的气泡凹陷. 在气泡回弹阶段, 上方凹陷发展为一束向下的反射流. 由于气泡膨胀射流变得越来越细而径向凹陷开始消失, 如图14(b)所示. 在反射流尖端形成了非常奇特的压力分布, 即在射流尖端具有较高的压力, 而在两侧却形成较低的环形压力区域. 反射流与下方气泡界面砰击后, 在气泡底部形成了局部高压, 而其余方向压力由于气泡过度膨胀降至较低水平. 此时气泡呈现为 "蝴蝶"形, 如图14(c)所示.

图14

图14   $\gamma _{s} =1.5$工况中典型时刻压力与速度场(分别对应气泡坍塌时刻$t=1.803$, 气泡回弹时刻$t=1.833$, $t=1.860$)

Fig.14   Pressure and velocity fields at some typical moments for $\gamma _{s} =1.5$ (corresponding times are $t=1.803$ during the bubble collapse stage; $t=1.833$ and $t=1.860$ during the bubble rebounding stage respectively)


图9第16 $\sim$ 19帧中气泡仅产生反射流. 图15给出了$\gamma_{s} =1.9$工况中典型时刻的压力与速度场. 与图14工况相同, 在第一个脉动周期内射流无法穿透气泡表面, 而在回弹阶段射流变细并加速向下运动, 如图15(a)和图15(b)所示. 射流砰击产生后气泡下方形成局部高压, 并且形成了向下的突出物, 如图15(c)和图15(d)所示.

图15

图15   $\gamma _{s} =1.9$工况中典型时刻压力与速度场(分别对应气泡最小体积时刻$t=1.816$, 气泡回弹时刻$t=1.825$, $t=1.837$, $t=1.891$)

Fig.15   Pressure and velocity fields at some typical moments for $\gamma _{s} =1.9$ (corresponding times are $t=1.816$ at the moment of reaching the minimum bubble volume; $t=1.825$, $t=1.837$, and $t=1.891$ during the bubble rebounding phase respectively)


当$\gamma_{s} $很小时($\gamma_{s} \leqslant 0.3)$, 浮体底面两侧产生的局部高压使气泡产生颈缩, 最后导致气泡撕裂. $\gamma _{s} $增大至($0.4\leqslant \gamma_{s} \leqslant 1.3)$, 气泡与浮体具有较强的流固耦合效应. 由于Bjerknes力[46-47]作用, 在坍塌阶段气泡产生指向浮体的正射流, 而当气泡与浮体较远时($1.1\leqslant \gamma _{s} \leqslant 2.0)$, 在坍塌阶段气泡顶部受到浮体的延时作用具有较高曲率并且在气泡上方压力梯度更大, 使气泡上表面产生向下的反射流. 自由液面对气泡具有排斥作用, 进一步促进了反射流的形成[48-49].

2.3.2 $\gamma_{s} $对射流速度的影响

由于本文实验装置和条件的限制, 很难精确测量射流速度, 为验证本文模型能预测射流速度, 首先对比本文模型计算结果与Cui等[44]放电气泡与不同厚度浮冰耦合作用的实验. 图16所示为最大射流速度$V_{jet} $(正射流)与浮冰厚度的关系($\gamma_{s} =0.93)$. 由图可知, 当浮冰厚度大于28.1 mm时, 计算所得最大射流速度基本在实验测量范围内, 而随着浮冰厚度减小或者射流速度增大, 两者偏差增大, 最大误差为16.8${\%}$(定义为实验与计算偏差/实验值). 实验中由于气泡壁面作为一个"凹透镜"折射光线导致测量的射流速度可能偏小[43,50]. 并且当射流速度很大时, 需要极高的拍摄帧率才能精确测量最大射流速度值, Cui等[44]的实验中拍摄帧率可能不足. 这两种因素导致了数值计算结果相比实验结果偏大. 浮冰厚度增大, 射流速度减小, 并且数值与实验结果吻合程度更高, 一定程度上验证了本文计算模型能较准确的预测射流速度. 本文后续进一步研究$\gamma_{s} $对脉动气泡与浮体非线性耦合效应的影响.

图16

图16   最大射流速度数值与Cui等[44]实验结果对比

Fig.16   Comparison of the maximum jet velocity between the present model and the experiment[44]


图17给出了3.3.1节中工况$(\gamma_{s} =0.4\sim2)$最大射流速度$V_{jet} $和最大反射流速度$V_{counter} $与$\gamma_{s}$的关系. 正射流速度远高于反射流且随$\gamma_{s} $先增大后减小再增大, 在$\gamma_{s} =0.8$时达到最大值136.8 (1368 m/s); 反射流速度随$\gamma_{s} $增大而增大, 最大值为27.2(272 m/s). 图12给出的$\gamma_{s} =0.9$工况中, 高速细射流砰击在另一侧气泡表面, 导致浮体底面附近具有很高的压力(超过140 MPa). 而当$0.7\leqslant \gamma_{s} \leqslant 0.9$时, 射流速度均大于100, 可能在浮体底部产生巨大的射流砰击载荷, 造成结构局部损伤. 与Lechner等[51]对超近壁面附近脉动气泡动力学特性的数值研究结果相似, 在气泡坍塌阶段与浮体底面平行的流动比垂直于底面的流动更快, 导致气泡在远端形成很高的曲率, 从而形成了速度约1000 m/s的高速细射流(Lechner等[51]发现的射流速度超过2000 m/s). 本文数值模型中最大射流速度远大于Philipp和Lauterborn[50]通过激光气泡与固定刚性板耦合作用实验测量的射流速度(150 m/s), 可见浮体与脉动气泡的非线性流固耦合效应能显著影响气泡的动力学特性, 在实际应用中需要考虑这种耦合作用的影响. 但实际中如此高速的细射流可能存在不稳定性[51], 并且流体可压缩性具有重要影响, 这不在本文的讨论范围内.

图17

图17   射流速度与无量纲气泡与浮体距离$\gamma_{s} $的关系

Fig.17   Maximum jet velocity versus the dimensionless bubble-floating structure distance $\gamma_{s} $


2.3.3 $\gamma_{s} $对气泡迁移和结构运动的影响

图18展示了不同$\gamma_{s}$下气泡形心的位移时历曲线(计算至射流穿透前). 当$\gamma_{s} <1.5$时, 在气泡膨胀阶段, 由于浮体对气泡上方流体的运动具有很强的阻碍, 导致气泡下表面相比上表面膨胀得更快, 类似于气泡被浮体推开, 形心远离结构. 而在坍塌阶段气泡受到浮体Bjerknes吸引[46]作用强于自由液面的Bjerknes排斥[35,48]作用, 气泡又被浮体吸回. 在气泡坍塌阶段近浮体与远离浮体侧不平衡压力使气泡整体受到一个朝向壁面的力, 从而导致了气泡朝浮体迁移, 如图6(d) $\sim\!$图6(f)所示. 由图18可知随着$\gamma_{s} $的增大, 这种相互排斥、吸引作用减小. 当$\gamma_{s} \geqslant 1.5$时(此时$\gamma _{f} \geqslant 4.1)$, 气泡在膨胀阶段基本保持其形心位置不变, 而在坍塌和回弹阶段气泡逐渐远离浮体. 这可能是因为自由液面对气泡的排斥力比浮体的吸引力更强.

图18

图18   不同$\gamma_{s} $下气泡形心的位移时历曲线

Fig.18   Time histories of the displacement of the bubble centroid for different $\gamma_{s} $


图19给出了浮体最大运动速度$V_{s}^{\max}$与$\gamma_{s} $的关系. 浮体最大速度随$\gamma_{s} $的增大而减小, 显然浮体距离气泡越远, 流固耦合效应越弱, 结构运动速度更小. 图中还给出了本文计算模型与修正差分法、未修正差分法的对比. 对于$\gamma _{s} \geqslant 0.4$的工况, 3种方法吻合较好. 相比未修正差分法, 修正差分法与本文方法吻合更好. 从局部放大图可知, 未修正差分法结果随$\gamma _{s} $减小逐渐偏离本文模型与修正差分法结果. 由图8可知本文模型计算结果更趋近于实验结果, 故假定本文辅助函数法计算结果为真实值, $\gamma_{s}=0.4$时修正差分法误差为0.27${\%}$, 未修正法误差为1.3${\%}$. 对于$\gamma _{s} <0.4$的超近工况, 气泡与浮体接触之后两种差分法获得的加速度时历曲线均出现剧烈震荡(未展示), 导致计算结果不准确, 而辅助函数法在本文参数研究范围内能获得光滑的加速度时历曲线, 证明本文模型计算更为准确.

图19

图19   浮体最大速度与无量纲气泡与浮体距离$\gamma _{s} $的关系

Fig.19   Maximum velocity of the floating body versus the dimensionless bubble-floating structure distance $\gamma _{s} $


3 结论

本文采用边界积分法联合辅助函数法, 建立了高压脉动气泡与浮体非线性流固耦合数值模型, 同时开展了柱形浮体与放电气泡流固耦合作用实验, 对比数值与实验结果, 二者吻合良好. 采用该数值模型, 本文对气泡与浮体无量纲距离参数$\gamma _{s} $进行了系统研究, 得到的主要结论如下;

(1) 将本文流固耦合模型、传统流固耦合模型(差分法计算速度势偏导数)与实验结果进行对比分析, 发现传统方法在处理气泡与浮体接触工况时, 精度明显不足, 而本文的辅助函数法大大提高了计算稳定性与精度.

(2) 在本文所讨论的参数范围内, 当$0.2\leqslant \gamma_{s} \leqslant 2$, 气泡射流模式可分为5类: ① $0.2\leqslant \gamma_{s} \leqslant 0.3$时, 气泡坍塌时产生颈缩, 并且产生环状射流将气泡撕裂为大单连通气泡与小环状气泡(颈缩型环状射流); ② $0.4\leqslant \gamma_{s} \leqslant 0.6$时, 气泡在膨胀阶段与浮体底面接触, 在坍塌阶段射流直接砰击在浮体底面(接触射流); ③ $0.7\leqslant \gamma_{s} \leqslant 1$时, 在坍塌阶段射流先穿透另一侧气泡表面, 然后穿过一层薄水膜后砰击在浮体底面(非接触射流); ④ $1.1\leqslant \gamma_{s} \leqslant 1.3$时产生射流相互撞击的现象, 由于浮体对气泡上表面的坍塌的延时作用, 顶部具有较高曲率, 形成向下的宽而慢的射流, 而在坍塌后期, 气泡底部形成向上的细而快的射流, 两者相互撞击(对射流); ⑤ $1.4\leqslant \gamma_{s} \leqslant 2$时产生反射流同时伴随径向凹陷(反射流).

(3) 正射流速度随着$\gamma_{s} $增大先增大后减小再增大, 而反射流速度随$\gamma_{s} $增大而增大. 正射流速度能达到约1000 m/s的量级, 说明脉动气泡与浮体的流固耦合效应在实际应用中影响重大, 其形成机理为气泡底部高曲率与浮体的Bjerknes吸引力联合作用. 该现象需要未来更精密的实验来验证.

(4) 总体而言, 在气泡膨胀阶段, 浮体排斥气泡, 气泡形心向下运动; 在坍塌阶段, 当$0.2\leqslant \gamma_{s} \leqslant 1.4$时, 浮体对气泡的Bjerknes吸引力强于自由液面Bjerknes排斥力, 导致气泡向浮体迁移; 当$1.5\leqslant \gamma_{s} \leqslant 2$, 浮体对气泡的Bjerknes吸引力比自由液面Bjerknes排斥力弱, 导致气泡远离自由液面.

参考文献

Li SM, Zhang AM, Wang QX, et al.

The jet characteristics of bubbles near mixed boundaries

Physics of Fluids, 2019,31:107105

[本文引用: 1]

Rungsiyaphornrat S, Klaseboer E, Khoo BC, et al.

The merging of two gaseous bubbles with an application to underwater explosions

Computers & Fluids, 2003,32:1049-1074

[本文引用: 1]

黄彪, 王国玉, 王复峰 .

非定常空化流场结构的实验研究

实验力学, 2011,26(4):417-424

[本文引用: 1]

(Huang Biao, Wang Guoyu, Wang Fufeng, et al.

Experimental investigation on flow field structure of unsteady cavitation flow

Journal of Experimental Mechanics, 2011,26(4):417-424 (in Chinese))

[本文引用: 1]

季斌, 程怀玉, 黄彪 .

空化水动力学非定常特性研究进展及展望

力学进展, 2019,49:201906

[本文引用: 1]

(Ji Bin, Cheng Huaiyu, Huang Biao, et al.

Research progresses and prospects of unsteady hydrodynamics characteristics for cavitation

Advances in Mechanics, 2019,49:201906 (in Chinese))

[本文引用: 1]

Song WD, Hong MH, Lukyanchuk B, et al.

Laser-induced cavitation bubbles for cleaning of solid surfaces

Journal of Applied Physics, 2004,95:2952-2956

[本文引用: 1]

Graaf KLd, Brandner PA, Penesis I.

The pressure field generated by a seismic airgun

Experimental Thermal and Fluid Science, 2014,55:239-249

[本文引用: 1]

Li S, van der Meer D, Zhang AM, et al.

Modelling large scale airgun-bubble dynamics with highly non-spherical features

International Journal of Multiphase Flow, 2020,122:103143

[本文引用: 1]

Ohl CD, Arora M, lkink R, et al.

Sonoporation from jetting cavitation bubbles

Biophysical Journal, 2006,91:4285-4295

[本文引用: 1]

Curtiss GA, Leppinen DM, Wang QX, et al.

Ultrasonic cavitation near a tissue layer

Journal of Fluid Mechanics, 2013,730:245-272

[本文引用: 1]

Lord Rayleigh.

VIII. On the pressure developed in a liquid during the collapse of a spherical cavity

The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1917,34:94-98

[本文引用: 2]

Benjamin TB, Ellis AT.

The collapse of cavitation bubbles and the pressures thereby produced against solid boundaries

. Philosophical Transactions for the Royal Society of London Series A, Mathematical and Physical Sciences, 1966,260:221-240

[本文引用: 2]

孙姣, 周维, 蔡润泽 .

垂直壁面附近上升单气泡的弹跳动力学研究

力学学报, 2020,52(1):1-11

[本文引用: 1]

(Sun Jiao, Zhou Wei, Cai Runze, et al.

The bounce dynamics of a rising single bubble near a vertical wall

Chinese Journal of Theoretical and Applied Mechanics, 2020,52(1):1-11 (in Chinese))

[本文引用: 1]

Pearson A, Cox E, Blake JR, et al.

Bubble interactions near a free surface

Engineering Analysis with Boundary Elements, 2004,28(4):295-313

[本文引用: 1]

郭文璐, 李泓辰, 王静竹 .

单空泡与自由液面相互作用规律研究进展

力学学报, 2019,51(6):1682-1698

[本文引用: 1]

(Guo Wenlu, Li Hongchen, Wang Jingzhu, et al.

Reaserch progress on interaction between a single cavitation and free surface

Chinese Journal of Theoretical and Applied Mechanics, 2019,51(6):1682-1698 (in Chinese))

[本文引用: 1]

Herring C.

Theory of the Pulsations of the Gas Bubble Produced by an Underwater Explosion

Columbia Univ., Division of National Defense Research, 1941

[本文引用: 1]

Keller JB, Kolodner II.

Damping of underwater explosion bubble oscillations

Journal of Applied Physics, 1956,27(10):1152-1161

[本文引用: 1]

Li S, Zhang AM, Han R, et al.

Experimental and numerical study of two underwater explosion bubbles: Coalescence, fragmentation and shock wave emission

Ocean Engineering, 2019,190:106414

[本文引用: 1]

张阿漫, 王诗平, 彭玉祥 .

水下爆炸与舰船毁伤研究进展

中国舰船研究, 2019,14(3):1-13

(Zhang Aman, Wang Shiping, Peng Yuxiang, et al.

Research progress in underwater explosion and its damage to ship structures

Chinese Journal of Ship Research, 2019,14(3):1-13 (in Chinese))

Brett JM, Krelle A.

A study of bubble collapse pressure pulse waves from small scale underwater explosions near the water surface

Journal of Sound and Vibration, 2018,435:91-103

[本文引用: 1]

Zhang S, Wang SP, Zhang AM.

Experimental study on the interaction between bubble and free surface using a high-voltage spark generator

Physics of Fluids, 2016,28(3):032109

[本文引用: 1]

张阿漫, 王诗平, 白兆宏 .

不同环境下气泡脉动特性实验研究

力学学报, 2011,43(1):71-83

[本文引用: 1]

(Zhang Aman, Wang Shiping, Bai Zhaohong, et al.

Experimental study on bubble pulse features under different circumstances

Chinese Journal of Theoretical and Applied Mechanics, 2011,43(1):71-83 (in Chinese))

[本文引用: 1]

Lindau O, Lauterborn W.

Cinematographic observation of the collapse and rebound of a laser-produced cavitation bubble near a wall

Journal of Fluid Mechanics, 2003,479:327-348

[本文引用: 1]

Brujan EA, Noda T, Ishigami A, et al.

Dynamics of laser-induced cavitation bubbles near two perpendicular rigid walls

Journal of Fluid Mechanics, 2018,841:28-49

[本文引用: 1]

Lauterborn W, Bolle H.

Experimental investigations of cavitation-bubble collapse in the neighbourhood of a solid boundary

Journal of Fluid Mechanics, 1975,72:391-399

[本文引用: 1]

Brujan EA, Takahira H, Ogasawara T.

Planar jets in collapsing cavitation bubbles

Experimental Thermal and Fluid Science, 2019,101:48-61

[本文引用: 1]

Chahine GL, Kalumuck KM, Hsiao CT.

Simulation of surface piercing body coupled response to underwater bubble dynamics utilizing 3DYNAFS, a three-dimensional BEM code

Computational Mechanics, 2003,32:319-326

[本文引用: 1]

Klaseboer E, Khoo BC, Hung KC.

Dynamics of an oscillating bubble near a floating structure

Journal of Fluids and Structures, 2005,21(4):395-412

[本文引用: 1]

Klaseboer E, Hung KC, Wang C, et al.

Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient/rigid structure

Journal of Fluid Mechanics, 2005,537:387-413

[本文引用: 1]

Zong Z, Wang J, Zhou L, et al.

Fully nonlinear 3D interaction of bubble dynamics and a submerged or floating structure

Applied Ocean Research, 2015,53:236-249

[本文引用: 2]

Li S, Zhang AM, Han R, et al.

3D full coupling model for strong interaction between a pulsating bubble and a movable sphere

Journal of Computational Physics, 2019,392:713-731

[本文引用: 4]

Li S, Khoo BC, Zhang AM, et al.

Bubble-sphere interaction beneath a free surface

Ocean Engineering, 2018,169:469-483

[本文引用: 3]

Zeng Q, Gonzalez-Avila SR, Ohl C-D.

Splitting and jetting of cavitation bubbles in thin gaps

Journal of Fluid Mechanics, 2020,896:A28

[本文引用: 1]

Liu YL, Zhang AM, Tian ZL, et al.

Investigation of free-field underwater explosion with Eulerian finite element method

Ocean Engineering, 2018,166:182-190

[本文引用: 1]

Wang QX, Yeo KS, Khoo BC, et al.

Nonlinear interaction between gas bubble and free surface

Computers & Fluids, 1996,25(7):607-628

[本文引用: 3]

Wang QX, Yeo KS, Khoo BC, et al.

Strong interaction between a buoyancy bubble and a free surface

Theoretical and Computational Fluid Dynamics, 1996,8:73-88

[本文引用: 4]

李帅, 张阿漫, 韩蕊.

气泡多周期运动时引起的流场压力与速度

力学学报, 2014,46(4):533-543

[本文引用: 1]

(Li Shuai, Zhang Aman, Han Rui.

Numerical analysis on the velocity and pressure fields induced by multi-oscillations of an underwater explosion bubble

Chinese Journal of Theoretical and Applied Mechanics, 2014,46(4):533-543 (in Chinese))

[本文引用: 1]

Best JP, Kucera A.

A numerical investigation of non-spherical rebounding bubbles

Journal of Fluid Mechanics, 1992,245:137-154

[本文引用: 3]

Li S, Han R, Zhang AM.

Nonlinear interaction between a gas bubble and a suspended sphere

Journal of Fluids and Structures, 2016,65:333-354

[本文引用: 2]

Zhang AM, Liu YL.

Improved three-dimensional bubble dynamics model based on boundary element method

Journal of Computational Physics, 2015,294:208-223

[本文引用: 1]

Ni BY, Zhang AM, Wu GX.

Numerical and experimental study of bubble impact on a solid wall

Journal of Fluids engineering, 2015,137(3):031206

[本文引用: 1]

Wang QX, Liu WK, Zhang AM, et al.

Bubble dynamics in a compressible liquid in contact with a rigid boundary

Interface Focus, 2015,5(5):20150048

URL     PMID      [本文引用: 1]

Turangan CK, Ong GP, Klaseboer E, et al.

Experimental and numerical study of transient bubble-elastic membrane interaction

Journal of Applied Physics, 2006,100(5):054910

[本文引用: 1]

李帅, 张阿漫, 韩蕊.

水中高压脉动气泡水射流形成机理及载荷特性研究

力学学报, 2019,51(6):1666-1681

[本文引用: 2]

(Li Shuai, Zhang Aman, Han Rui.

The mechanism of jetting behaviors of an oscillating bubble

Chinese Journal of Theoretical and Applied Mechanics, 2019,51(6):1666-1681 (in Chinese))

[本文引用: 2]

Cui P, Zhang AM, Wang SP, et al.

Ice breaking by a collapsing bubble

Journal of Fluid Mechanics, 2018,841:287-309

[本文引用: 6]

Lauterborn W.

Cavitation bubble dynamics - new tools for an intricate problem

Applied Scientific Research, 1982, (38):165-178

[本文引用: 1]

Blake JR, Taib BB, Doherty G.

Transient cavities near boundaries. Part 1. Rigid boundary

Journal of Fluid Mechanics, 1986,170:479-497

[本文引用: 2]

Blake JR.

The Kelvin impulse: application to cavitation bubble dynamics

The Journal of the Australian Mathematical Society Series B Applied Mathematics, 1998,30(2):127-146

[本文引用: 1]

Blake JR, Gibson DC.

Growth and collapse of a vapour cavity near a free surface

Journal of Fluid Mechanics, 1981,111:123-140

[本文引用: 2]

Robinson PB, Blake JR, Kodama T, et al.

Interaction of cavitation bubbles with a free surface

Journal of Applied Physics, 2001,89(12):8225-8237

[本文引用: 1]

Philipp A, Lauterborn W.

Cavitation erosion by single laser-produced bubbles

Journal of Fluid Mechanics, 1998,361:75-116

[本文引用: 2]

Lechner C, Lauterborn W, Koch M, et al.

Fast, thin jets from bubbles expanding and collapsing in extreme vicinity to a solid boundary: A numerical study

Physical Review Fluids, 2019,4:021601

[本文引用: 3]

/