船舶力学

工况	吃水/m	速度/（m∙s⁻¹）	气泡直径/μm
1	0.25	1.055	200
2	0.25	1.055	500
3	0.25	1.055	800
4	0.25	0.946	500
5	0.25	1.164	500
6	0.138	1.055	500
7	0.29	1.055	500

工况	吃水/m	速度/（m∙s⁻¹）	气泡直径/μm
1	0.25	1.055	200
2	0.25	1.055	500
3	0.25	1.055	800
4	0.25	0.946	500
5	0.25	1.164	500
6	0.138	1.055	500
7	0.29	1.055	500

航速/（m∙s⁻¹）	船舶计算摩擦阻力/N	ITTC公式计算摩擦阻力/N	误差（%）
0.946	7.49	7.57	1.1
1.055	9.22	9.19	0.3
1.164	11.04	11.05	0.1

航速/（m∙s⁻¹）	船舶计算摩擦阻力/N	ITTC公式计算摩擦阻力/N	误差（%）
0.946	7.49	7.57	1.1
1.055	9.22	9.19	0.3
1.164	11.04	11.05	0.1

低速船舶微气泡减阻数值研究

PDF下载

赵晓杰 ¹ , 宗智 ²^,³ , 王加夏 ¹ , 洪智超 ¹^,⁴ , 胡俊明 ¹

船舶力学 | 流体力学 2024,28(3): 368-378

收起

船舶力学 | 流体力学 2024, 28(3): 368-378

低速船舶微气泡减阻数值研究

全屏

赵晓杰¹, 宗智²^,³, 王加夏¹, 洪智超¹^,⁴, 胡俊明¹

作者信息

^1.江苏科技大学　船舶与海洋工程学院，江苏　镇江　212000

^2.大连理工大学　船舶工程学院，辽宁　大连　116024

^3.工业装备结构分析国家重点实验室，辽宁　大连　116024

^4.江苏省船舶与海洋工程装备技术创新中心，江苏　南通　226100

赵晓杰（1992-），男，讲师，E-mail：xjzhao@just.edu.cn。

Numerical study on microbubble drag reduction of low-speed ships

Xiao-jie ZHAO¹, Zhi ZONG²^,³, Jia-xia WANG¹, Zhi-chao HONG¹^,⁴, Jun-ming HU¹

Affiliations

^1.School of Naval Architecture and Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang 212000, China

^2.School of Shipbuilding Engineering, Dalian University of Technology, Dalian 116024, China

^3.State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian 116024, China

^4.Jiangsu Marine Technology Innovation Center, Nantong 226100, China

出版时间: 2024-03-20 doi: 10.3969/j.issn.1007-7294.2024.03.005

文章导航

摘要

收起

为了研究船舶微气泡减阻规律，本文基于OpenFOAM中两相欧拉数值模型，对低速散货船进行微气泡减阻数值研究。对气液两相分别建立控制方程，考虑五种相间作用力及气泡聚合和破碎，采用考虑气泡影响的改进k-ε湍流模型，忽略自由面影响，采用叠模模型研究喷气量、气泡直径、航速及吃水等因素对船舶微气泡减阻的影响，分析气体体积分数、湍流粘度和气泡直径分布等。结果表明：微气泡可以同时减少船舶摩擦阻力、粘压阻力和总阻力；喷气量直接影响减阻率，喷气量越大，减阻率越高；较小气泡的平均气体体积分数较大且气体分布更均匀，同时湍流运动粘度较小，可以更有效减阻；气泡沿着流向会聚并，气泡越小聚并越剧烈；较高航速和小吃水更有利于减阻。

关键词

两相欧拉模型 / 微气泡减阻 / 船舶 / 摩擦阻力 / 数值模拟

Abstract

收起

In order to investigate the microbubble drag reduction (MBDR) of ships, numerical studies of MBDR on a low-speed bulk carrier model were conducted based on the two-phase Euler model in OpenFOAM. The governing equations were established for the gas and liquid phases, respectively, considering five kinds of interfacial forces and bubble coalescence and breakup. The modified k-ε turbulence model considering the effects of bubbles was also used, and the superimposed model was adopted to ignore the influence of free surface. The effects of air flow rate, bubble size, ship speed and draft on MBDR of the ship were investigated while the distributions of air volume fraction, turbulent viscosity and bubble size around the ship were analyzed. The numerical results show that: micro bubbles can simultaneously reduce the frictional drag, viscous pressure drag and total drag of the ship; air flow rate can directly influence the drag reduction and more air flow rate can lead to more drag reduction; smaller micro bubbles can lead to higher average volume fraction, more uniform gas distribution and smaller turbulent viscosity, resulting in drag reduction more effectively;bubble coalescence will occur along the direction of the flow and the coalescence effects are more intense for smaller bubbles; higher ship speed and lower draft are more conducive to drag reduction.

Key words

two-phase Euler model / micro bubble drag reduction / ship / frictional resistance / numerical simulation

引用本文

赵晓杰, 宗智, 王加夏, 洪智超, 胡俊明. 低速船舶微气泡减阻数值研究. 船舶力学, 2024 , 28 (3) : 368 -378 . DOI: 10.3969/j.issn.1007-7294.2024.03.005

Xiao-jie ZHAO, Zhi ZONG, Jia-xia WANG, Zhi-chao HONG, Jun-ming HU. Numerical study on microbubble drag reduction of low-speed ships[J]. Journal of Ship Mechanics, 2024 , 28 (3) : 368 -378 . DOI: 10.3969/j.issn.1007-7294.2024.03.005

正文

收起

0　引言

收起

航行中的船舶会受到阻力。为了实现减阻，研究者们提出了许多减阻技术^[1]，如为了减少兴波阻力，采用球鼻艏或优化船体型线；为了降低摩擦阻力，在船舶表面设置超疏水涂层，或在船底喷入高分子聚合物、微气泡等。其中，微气泡减阻技术受到人们的广泛关注，因为它的减阻效率较高并且无污染。对于低速船舶，摩擦阻力占主要部分，微气泡减阻技术更具应用前景。

几十年来，国内外学者对微气泡减阻进行了许多试验研究，研究对象主要有平板、回转体、船体等。McCormick等^[2]首先通过试验发现了微气泡减阻现象，然后研究学者又相继对平板和回转体的微气泡减阻进行了试验研究^[3]。对于船舶微气泡减阻，减阻情况更加复杂，因为减阻效果会受船型、喷气形式、复杂工况等多个因素影响。Sayyaadi等^[4]对一双体船进行了试验研究，发现喷气量对减阻效果有显著影响，并存在一个临界喷气量；王家楣等^[5]采用多孔硅板生成微气泡，对不同吃水的船模微气泡减阻进行了试验研究，发现船舶在浅吃水下减阻更好。

前人也对微气泡减阻做了许多数值研究。Kunz等^[6]利用两相欧拉模型研究了平板的微气泡减阻，通过和试验结果对比验证了模型的准确性；傅慧萍^[7]利用不同多相流模型对平板微气泡减阻进行了数值模拟并对影响因素进行了分析；对于船舶微气泡减阻，郑晓伟等^[8]利用相间滑移算法研究了二维船的微气泡减阻，所得结果与试验定性符合；吴浩与欧勇鹏^[9]将船底简化为平板并利用Mixture模型研究喷气减阻，较好地模拟出了气液混合流。可以看出，目前船舶微气泡减阻模拟都对数值模型进行了一定简化，几何模型简单，并且没有充分考虑气液间及气泡间相互作用。

因此，为了进一步研究船舶微气泡减阻规律，本文基于OpenFOAM中两相欧拉模型对低速散货船的微气泡减阻展开数值研究，考虑五种相间作用力以及气泡的聚合和破碎，并考虑气泡对液体湍流的影响。区别于二维平板，船舶具有不同阻力成分，本文研究气泡对船舶摩擦阻力、粘压阻力和总阻力的减阻率影响，同时研究喷气量、气泡直径、航速及吃水等因素对减阻率的影响，分析船舶附近气体体积分数、湍流粘度和气泡直径的分布等，为微气泡减阻技术在船舶上的应用提供理论指导。

1　数值模型

收起

1.1　两相欧拉模型

在OpenFOAM的两相欧拉模型中，对气液两相（本文指水和空气）分别进行建模，建立的质量和动量守恒方程^[10]分别为

(1)

(2)

式中：i为相指数，这里i =l或g分别表示液体和气体；α_i表示i相体积分数；U_i表示i相流体速度；ρ_i表示i相流体密度；p为压力；g为重力加速度；

为i相的相间作用力，用于考虑相间动量传递；ν_i表示i相运动粘度系数；ν_t，i表示i相湍流运动粘度系数。

1.2　相间作用力

为了充分考虑气液间动量传递作用，本文考虑五种相间作用力，公式如下：

(3)

式中：F_d为相间阻力，F_l为升力，F_vm为虚质量力，F_td为湍流离散力，F_wl为壁面润滑力。

相间阻力是气泡在液体中运动时受到的阻力，与相对速度方向相反，计算公式为

(4)

式中：α_g为气体体积分数，d_g为气泡直径，ρ_l为液体密度，U_l和U_g分别表示液体速度和气体速度，C_d为相间阻力系数。本文采用Ishii和Zuber^[11]提出的阻力模型，该模型考虑三种气泡形态，对于球形气泡，阻力系数计算公式为

(5)

对于椭球型气泡，阻力系数计算公式为

(6)

式中，Re是气泡雷诺数，Eo表示无量纲Eötvös数，f（α_g）表达式为

(7)

对于帽子形气泡，阻力系数计算公式为

(8)

不同气泡形状区域划分基于以下公式：

(9)

升力是由于液体速度梯度引起的剪切力，它是除阻力外最重要的作用力，可以影响气体垂向分布，计算公式为

(10)

式中：C_l是升力系数。采用Tomiyama升力模型^[12]，可以考虑不同直径的气泡，计算公式为

(11)

式中，函数f（Eo）计算表达式为f（Eo） = 0.00105Eo³ - 0.0159Eo² - 0.0204Eo + 0.474。

虚拟质量力是指气泡相对于流体变速时带动附近流体而产生的作用力，表达式为

(12)

式中，C_vm为虚拟质量力系数，这里假设气泡是球形，根据势流理论取C_vm =0.5。

湍流耗散力为流体湍流脉动作用对气泡产生的作用力，采用Bertodano模型^[13]，计算公式为

(13)

壁面润滑力是为了防止气泡接触壁面而产生的使气泡远离壁面的作用力，计算公式^[14]为

(14)

式中，壁面润滑力系数C_wl的计算公式为

(15)

式中：y_w为气泡离壁面的垂直距离；C_w1和C_w2是校准系数，根据文献^[14]推荐，取值分别为C_w1 =-0.01和C_w2 =0.05，表示壁面润滑力在距离壁面5倍直径距离时变为0。

1.3　气泡群体平衡模型

为了考虑气泡直径变化，采用基于界面面积密度方程的气泡群体平衡模型，控制方程^[15]为

(16)

式中：a表示单位体积气泡界面面积密度；ψ表示气泡形状因子，忽略气泡变形，ψ取值为1/36π；S_j表示气泡聚并和破碎的源项，可以分成三部分：

(17)

其中，S_RC和S_WE分别表示随机碰撞和尾迹夹带作用，它们会引起气泡聚并，表达式分别为

(18)

(19)

式中：n_g为单位体积气泡数密度，u_t与u_r分别为湍流速度和气液相对速度。根据Ishii等^[15]实验研究，对于随机碰撞，常数取值C_RC =0.04和C =3。对于尾迹夹带，常数取值C_WE =0.04和α_max =0.75。S_TI为湍流作用，它导致气泡发生破碎，表达式为

(20)

式中，C_TI为湍流作用参数，We_cr为临界韦伯数，取值分别为C_TI =0.06和We_cr =6.0。

1.4　湍流模型

对于湍流模型，本文采用Lahey^[16]提出的改进k-ε模型。为了考虑气泡影响，在控制方程中加入源项，并增加额外湍流粘度来计算气泡导致的湍流增加。湍流模型的控制方程为

(21)

(22)

式中，湍流粘度系数ν_t的计算公式为

(23)

源项

与

的计算公式为

(24)

式中，常数C_p取值为0.25，其它参数为C_μ =0.09，C_ε，1 =1.44，C_ε，2 =1.92，σ_k =1.0，σ_ε =1.3。

2　计算设置

收起

2.1　几何模型

本文研究的是95 000 DWT散货船，缩尺比为1：50，实船垂线间长为231 m，型宽为38 m，型深为20.7 m，设计吃水为12.5 m，排水体积为91 893.0 m³，湿表面积为12 543.0 m²，设计航速为14.5 kn，船体型线见图1。在船模17站处设置长方形喷气板，长400 mm，宽90 mm，如图2所示。

2.2　计算域和网格划分

计算区域如图3所示，不考虑自由面影响，采用叠模模型，将自由面当作对称滑移壁面。为了减少计算量，采用一半船体进行计算。速度入口距离船体1倍船长，出口边界距离船体3倍船长，底面和侧面距离船体1倍船长。

网格划分利用切割体网格，网格分布如图4所示。最大网格尺寸为231 mm，在船体附近进行网格加密，最小网格尺寸约为7.2 mm，船体表面设置三层边界层，第一层网格厚度满足y⁺ >30。

2.3　边界条件和计算方法

对于边界条件，船体对液体速度采用无滑移壁面条件，对于气泡速度采用滑移边界。入口为速度入口，出口为压力出口。对于喷气面，本文假设采用微米多孔板产生微气泡，为了匹配实际物理条件，假定孔隙率为0.3，与Kunz等^[6]数值研究中条件一致。对于其他边界，设置为对称边界条件。

控制方程采用有限体积法进行离散。对于空间离散，梯度项和拉普拉斯项采用高斯线性格式（Gauss linear），对流项采用限制高斯线性格式（Gauss limited linear），在相输运控制方程中，采用Gauss vanLeer离散格式。空间离散为二阶精度。对于时间离散，采用一阶精度的隐式欧拉格式。

控制方程求解采用隐式速度压力耦合算法（PISO），在OpenFOAM^[17]中利用PIMPLE循环进行设置，外迭代次数设为1，压力矫正次数为2。时间步采用自适应时间步，最大库朗数设为0.9，从而保证计算收敛。

2.4　计算工况设置

为了研究航速、喷气量、气泡直径和吃水等因素对船舶微气泡减阻的影响，设置了不同计算工况。考虑三种航速，分别为0.946 m/s、1.055 m/s和1.164 m/s，其中1.055 m/s为设计航速。考虑三种吃水，分别为0.138 m、0.25 m和0.29 m，其中0.25 m为模型船的设计吃水。考虑三种气泡直径，分别为200 μm、500 μm和800 μm。计算工况设置总结见表1。

采用无量纲喷气量C_v作为参数，表达式为C_v =Q_g/（Q_g+Q_l），式中Q_g是气体喷气流量，Q_l是边界层内液体流量。Q_l计算公式为Q_l=U（δ-δ^∗）b，式中U为船舶航速，b为船宽，δ为边界层厚度，δ^∗为边界层排挤厚度，表达式为δ^∗ =0.125δ，根据边界层理论^[18]，δ计算公式为δ =0.37x/Re^0.2。无量纲喷气量C_v范围为0~0.05，间隔为0.01。

本文研究三种阻力成分的减阻率，计算公式为

，式中R_i表示未喷气时阻力，

表示喷气后阻力，下标i取f、p和t时分别表示摩擦阻力、粘压阻力和总阻力。

3　计算结果与分析

收起

3.1　计算模型验证

为了验证本文模型，采用两种算例：一是计算不喷气下船舶摩擦阻力，并和ITTC推荐的摩擦阻力计算公式进行对比；二是计算喷气情况下平板微气泡减阻，并和前人试验数据进行对比。其中ITTC推荐的摩擦阻力系数公式为

(25)

表2表示船舶在三种速度下摩擦阻力和通过ITTC推荐公式计算得到的摩擦阻力结果对比，可以看出计算误差在2%之内，说明本文的数值模型可以准确计算船舶阻力。

平板微气泡减阻试验验证采用Merkle等^[3]的实验，通过多孔介质板产生微气泡实现减阻，采用实验工况，考虑流速为9.6 m/s和14.2 m/s，本文气泡群体平衡模型需要先给定平均气泡直径，根据实验观测^[3]以及Mohanarangam等^[19]在相同工况计算所得气泡直径分布，两种流速下平均气泡直径分别取600 μm和400 μm。图5和图6分别表示流速为9.6 m/s和14.2 m/s时的数值结果，横坐标为喷气量，纵坐标为平板摩擦阻力减阻率，可以看出数值结果和试验值吻合良好，说明本文微气泡减阻数值模型是可行的和准确的。

3.2　不同气泡直径下船舶微气泡减阻率

气泡直径是影响微气泡减阻率的重要因素。图7（a）~（c）表示在不同气泡直径下船舶摩擦阻力、粘压阻力和总阻力的微气泡减阻率，对应计算工况1~3。从图中可以看出喷气情况下，摩擦阻力、粘压阻力和总阻力都实现了减阻，并且随着喷气量的增大，减阻率也都随之增大，说明喷气量是影响船舶微气泡减阻率的决定因素。另外，摩擦阻力和总阻力的减阻率增大速度一开始很快，但随后会变得缓慢，直到达到最大值后几乎不再改变，说明对于摩擦阻力和总阻力，存在一个最佳临界喷气量，这和Sayyaadi等^[4]试验研究所得的结论相一致。对于粘压阻力，减阻率增加存在波动，这是因为微气泡对流场存在扰动、局部压力出现不稳定波动造成的，但总体上减阻率单调增加。

对于不同气泡直径，从图7（a）~（c）可以看出，较小的气泡可以获得较大的摩擦阻力、粘压阻力和总阻力减阻率，说明小气泡可以更有效地减阻，这和Song等^[20]三维回转体的微气泡减阻试验结论相一致。由于船舶是低速散货船，摩擦阻力占比较大（约80%），而粘压阻力占比较小（约20%），尽管粘压阻力减阻率存在不稳定波动，但是摩擦阻力减阻率和总阻力减阻率的变化趋势还是彼此接近。

3.3　不同气泡直径下船舶附近气体分布

图8（a）~（c）给出了无量纲喷气量为0.05时不同气泡直径下船舶附近气体体积分数分布。可以看出在最小气泡（d=200 μm）时，船底气体分布均匀，计算得到表面气体体积分数均值为0.327，此时减阻率也较大；在最大气泡（d=800 μm）时，靠近喷气板处气体覆盖较多，随着远离喷气板气体逐渐减少，导致船底气体分布不均匀，计算得到表面气体体积分数均值为0.308，此时减阻率也较小，说明近壁面气体体积分数分布是影响船舶微气泡减阻的关键因素。

不同气泡直径导致船体附近气体体积分数分布不同，原因是不同气泡直径的气液相互作用力不同。由公式（4）可以看出，在相同气体体积分数时，气泡直径越小，气液间相间阻力越大，使得小气泡从喷气板喷出后难以远离壁面，更容易沿着液体流动方向延伸至船尾，从而造成气体分布均匀，实现较高减阻率；而对于大气泡，气液相间阻力较小，使得大气泡难以沿着液体流动方向向船尾延伸，从而导致气体在喷气口附近聚集，船尾附近气体体积分数较低，导致减阻率也较低。

3.4　不同气泡直径下湍流运动粘度分布

为了进一步分析不同气泡直径下船舶微气泡减阻率不同的原因，图9（a）~（d）给出了不喷气时及喷气时（C_v=0.05）不同气泡直径的船舶附近湍流运动粘度分布。可以看出，喷气后船底液体湍流运动粘度都会变大，说明气泡对液体湍流具有增强作用，这和Mohanarangam等研究得到的结论相一致。根据公式（23），可以得到在相同气体体积分数下，气泡直径越大湍流运动粘度也越大的结论。图9也表明，800 μm、500 μm和200 μm直径气泡产生的湍流运动粘度依次减小，因此这也是小气泡可以更有效减阻的原因之一。

3.5　船舶附近气泡直径变化

图10（a）~（c）为无量纲喷气量为0.05时不同初始气泡直径下船舶附近气泡直径分布。从图中可以看出，在船底沿着流动方向气泡直径都有所增大，这和Merkle等^[3]平板微气泡减阻实验的气泡变化趋势相同，原因是气泡沿流向会发生随机碰撞，再加上相邻气泡尾迹夹带的作用，气泡会发生聚并，导致气泡直径变大。在船舶尾部附近，由于存在较大的几何形状突变，船尾湍流作用剧烈，就会容易导致气泡破碎，气泡直径又会变小。

另外，从图中可以看出不同初始气泡直径下的气泡直径变化情况不同。对于较小初始直径气泡如200 μm直径气泡，在喷气板不远处就开始发生聚并，最大气泡直径可达约260 μm，气泡直径增加约30%，大气泡在船底分布也多；而对于较大初始直径气泡如800 μm直径气泡，在距离喷气板较远处才开始发生聚并，最大气泡直径约为910 μm，气泡直径增加约13%，增大的气泡仅分布在船尾附近。可见较小气泡的聚并作用更强，原因在于不同气泡直径导致的气泡聚并和破碎作用不同，根据公式（18）~（19），气泡随机碰撞和尾迹夹带作用与气泡数密度有关，相同喷气量下小气泡的气泡数密度更大，所以气泡聚并作用更强。又根据3.4节可得，小气泡的湍流增强作用较小，所以湍流作用导致的气泡破碎作用也较弱，较小气泡的聚并作用更加剧烈。

3.6　航速与吃水对船舶微气泡减阻的影响

图11为不同船舶航速下不同阻力成分的微气泡减阻率，对应工况1、4和5。可以看出，不同航速下摩擦阻力和总阻力减阻率都随喷气量的增加而增加，到最后增加速度会变慢，减阻率都会接近一个恒定值；粘压阻力减阻率也随喷气量的增加而增加，但也存在波动，这是因为气泡导致流场局部压力不稳定造成的。此外，航速较高时船舶减阻率较高，说明此时船舶附近气体体积分数较大。这是因为不同航速下气液相互作用力不同，根据公式（4）可知流速较高时，气泡和水的相对速度较大，此时气液相间阻力较大，气泡喷出后远离壁面趋势较小，气泡在较强的气液相互作用下更容易延伸至船尾，造成船底气体体积分数较大且气体分布均匀，因此船底获得更大气体覆盖面积，从而获得较大减阻率。

图12为不同船舶吃水下不同阻力成分的微气泡减阻率，对应工况1、6和7。可以看出，对于摩擦阻力、粘压阻力和总阻力，小喷气量时减阻率差别较小，随着喷气量增加，减阻率差别逐渐变大。另外，小吃水时减阻率较大，这和王家楣等^[5]试验的结论相一致，这是因为船舶减阻主要原因是船底摩擦阻力的减少，小吃水时船底摩擦阻力在整个船体总阻力中占比较大，更容易获得较大的减阻率。

4　结论

收起

本文基于OpenFOAM中两相欧拉模型对三维低速散货船的微气泡减阻进行了数值研究，考虑五种相间作用力、气泡聚合和破碎以及气泡对液体湍流的影响，研究了气泡对不同阻力成分减阻率的影响，研究了喷气量、气泡直径、航速及吃水等因素的影响，分析了船舶附近气体体积分数、湍流粘度和气泡直径的分布，得出的主要结论如下：

（1）气泡可以同时减少低速船舶摩擦阻力、粘压阻力和总阻力，并且随着喷气量增大，减阻率也都随之增大，说明喷气量是影响船舶微气泡减阻率的决定因素。

（2）船舶附近气体体积分数分布是影响船舶微气泡减阻的关键因素。气泡较小时船舶附近的平均气体体积分数较大，且气体分布更加均匀，减阻率较高。另外，较小气泡产生的湍流运动粘度也较小，这也是小气泡可以更有效减阻的原因之一，因此在工程应用中应尽量保证气泡直径足够小。

（3）气泡会沿着液体流动方向发生聚并，导致气泡直径增大。较小气泡的聚并作用更剧烈，原因是相同喷气量下它的气泡数密度更大，随机碰撞和尾迹夹带作用更强，另外较小气泡的湍流增强作用较小，导致气泡破碎作用也较弱。

（4）较高航速和小吃水更有利于船舶微气泡减阻。较高航速时，气液相间阻力较大，船底气体体积分数较大且气体分布均匀，从而获得较大减阻率。小吃水时，船底摩擦阻力在整个船体总阻力中占比较大，更容易获得较大的减阻。

基金

收起

江苏省自然科学基金资助项目(BK20220652)
国家自然科学基金重大项目(52192692)
面上项目(52171294)
国家自然科学基金青年基金项目(52101359; 52301321)
江苏省自然科学基金资助项目(BK20230669)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

Perlin

, Ceccio

. Mitigation of hydrodynamic resistance: Methods to reduce hydrodynamic drag[M]. World Scientific, 2015.

[2]

Mccormick

M E

, Bhattacharyya

. Drag reduction of a submersible hull by electrolysis[J]. Naval Engineers Journal, 1973, 85(2): 11-16.

[3]

Merkle

C L

, Deutsch

. Microbubble drag reduction in liquid turbulent boundary layers[J]. Applied Mechanics Reviews, 1992, 45(3): 103-127.

[4]

Sayyaadi

, Nematollahi

. Determination of optimum injection flow rate to achieve maximum micro bubble drag reduction in ships: An experimental approach[J]. ScientiaIranica, 2013, 20(3): 535-541.

[5]

王家楣，郑晓伟，姜曼松. 船舶吃水对微气泡减阻影响的水池试验研究[J]. 船舶工程，2004，26(6): 9-12.

Wang

Jiamei

, Zheng

Xiaowei

, Jiang

Mansong

. Experimental study of the effects of ship draft on microbubble drag reduction[J]. Ship Engineering, 2004, 26(6): 9-12.

[6]

Kunz

R F

, Gibeling

H J

, Maxey

M R

, et al. Validation of two-fluid Eulerian CFD modeling for microbubble drag reduction across a wide range of Reynolds numbers[J]. Journal of Fluids Engineering, 2007, 129(1): 66-79.

[7]

傅慧萍. 平板微气泡减阻数值模拟及影响因素分析[J]. 哈尔滨工程大学学报，2015，36(1306): 1297-1301.

Huiping

. Numerical simulation of microbubble drag reduction in a plate and factors influencing its practicality process[J]. Journal of Harbin Engineering University, 2015, 36(1306): 1297-1301.

[8]

郑晓伟，王家楣，曹春燕. 二维船舶微气泡减阻数值模拟[J]. 船舶工程，2005，27(6): 15-18.

Zheng

Xiaowei

, Wang

Jiamei

, Cao

Chunyan

. Numerical simulation of friction reduction by micro-bubble for a two-dimension vessel[J]. Ship Engineering, 2005, 27(6): 15-18.

[9]

吴浩，欧勇鹏. 平板喷气粘性流场数值计算方法研究[J]. 舰船科学技术，2016，38(15): 47-51.

Hao

, Ou

Yongpeng

. Numerical study of method of flat plate viscous flow field with bubble[J]. Ship Science and Technology, 2016, 38(15): 47-51.

[10]

Rusche

. Computational fluid dynamics of dispersed two-phase flows at high phase fractions[D]. Imperial College London (University of London), 2003.

[11]

Ishii

, Zuber

. Drag coefficient and relative velocity in bubbly, droplet or particulate flows[J]. AIChE Journal, 1979, 25(5): 843-855.

[12]

Tomiyama

. Struggle with computational bubble dynamics[J]. Multiphase Science & Technology, 1998, 10(4): 369-405.

[13]

Bertodano

M A L D

. Two fluid model for two-phase turbulent jets[J]. Nuclear Engineering and Design, 1998, 179(1): 65-74.

[14]

CD-adapco. User guide: Star-ccm+. Version 13.02[M]. 2018.

[15]

Ishii

, Kim

, Kelly

. Development of interfacial area transport equation[J]. Nuclear Engineering and Technology, 2005, 37(6): 13-29.

[16]

Lahey

R T

. The simulation of multidimensional multiphase flows[J]. Nuclear Engineering and Design, 2005, 235(10-12): 1043-1060.

[17]

Jasak

, Jemcov

, Zeljko

. OpenFOAM: A C++ library for complex physics simulation[C]//International Workshop on Coupled Methods in Numerical Dynamics, Dubrovnik, Croaia, 2007.

[18]

Marusic

, Monty

J P

, Hultmark

, et al. On the logarithmic region in wall turbulence[J]. Journal of Fluid Mechanics, 2013, 716: R3.

[19]

Mohanarangam

, Cheung

S C P

, Tu

J Y

, et al. Numerical simulation of micro-bubble drag reduction using population balance model[J]. Ocean Engineering, 2009, 36(11): 863-872.

[20]

Song

, Wang

, Wei

, et al. Experimental study of microbubble drag reduction on an axisymmetric body[J]. Modern Physics Letters B, 2018(32): 1850035.

2024年第28卷第3期

PDF下载

108

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2024.03.005

接收时间：2023-09-23
首发时间：2026-03-21
出版时间：2024-03-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2023-09-23

基金

江苏省自然科学基金资助项目(BK20220652)

国家自然科学基金重大项目(52192692)

面上项目(52171294)

国家自然科学基金青年基金项目(52101359; 52301321)

江苏省自然科学基金资助项目(BK20230669)

作者信息

^1.江苏科技大学　船舶与海洋工程学院，江苏　镇江　212000

^2.大连理工大学　船舶工程学院，辽宁　大连　116024

^3.工业装备结构分析国家重点实验室，辽宁　大连　116024

^4.江苏省船舶与海洋工程装备技术创新中心，江苏　南通　226100

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2024.03.005

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT