高速排水型船舶通过采用更先进的船体设计、使用轻量化材料和先进的推进系统,可以减少船舶的阻力,提高航速和燃油效率,是一种非常具有潜力的新型船舶类型。
高速排水型船舶的尾浪对于环境和自身运行都会产生不同程度的影响,一方面,高速船舶航行产生的尾浪可能会冲刷堤岸、危害水工建筑。另一个方面,高速排水型船舶在产生尾浪时,由于波浪的作用力,船舶的航行稳定性可能会下降,这可能会对船舶的操纵和控制带来一定的挑战,特别是在恶劣海况下。因此,有些航行区域会限制高速船舶的行驶速度。
本文基于流体动力学理论[1],对高速排水型船舶的水动力特性进行分析,并通过建立流体有限元模型,结合广义边界元计算理论,对高速排水型船舶的尾浪特性进行研究,对于指导高速排水型船的优化设计有一定作用。
1 高速排水型船舶的控制方程针对高速排水型船舶的运动特性,本文建立随动坐标系如图1所示:
![]() |
图 1 高速排水型船舶的随动坐标系 Fig. 1 The servo coordinate system of the high speed displacement ship |
图中,原点o在船舶质心,oz垂直向上,船舶前进方向为oy轴正向。
结合随动坐标系,建立高速排水型船舶的运动模型如下:
m0∂2∂t2Xg=→F,→T→Ig→T−1→w+→w×→T→Ig−1=→Mg。 |
式中:
高速排水型船舶的摩擦系数
μf=0.075(logRn−3)2,Ff=12ρV02S0μf。 |
其中:
结合海水流场的特性,建立高速排水型船舶的流体连续性方程如下:
ρ∗dVdt+δ(ρl)δy+δ(ρm)δz+δ(ρn)δx=0。 |
式中:
建立高速排水型船舶的流体动量方程如下:
∂(ρux)∂t+∂(ρuxuy)∂x=−∂p∂x+∂∂x[μ(∂ux∂y+∂uy∂x)]+Fi。 |
式中:
建立高速排水型船舶的流体湍流模型如下:
{∂k∂t+ˉu∂k∂x=∂∂y[(v+σkvt)∂k∂xj]+Pk−β∗kω,∂ω∂t+ˉu=∂∂x[(v+σkvt)∂ω∂xj]+2(1−Fi)σ。 |
式中:
高速排水型船舶在流场中的水动力性能与流体的基本特性密切相关,如图2所示,流体根据内部分子的运动状态分为层流、湍流和页面边界。
![]() |
图 2 流体的特性示意图 Fig. 2 Diagram of fluid characteristics |
图中,页面边界又可分为静力学边界和动力学边界,本文在分析高速排水型船舶的流体特性时,采用动力学边界条件[3]。
建立边界面方程为:
F(x,y,z,t)=η(x,y,t)−z=0, |
质点的法向速度建模为:
vn=vxFx+vyFy+vzFz√F2x+F2y+F2z, |
进而可得边界面上的质点法向速度建模为:
ϕn=ϕxFx+ϕyFy+ϕzFz√F2x+F2y+F2z, |
由于
∂η∂t+∂ϕ∂x∂η∂x+∂ϕ∂y∂η∂y−∂ϕ∂z=0, |
根据等压条件和拉格朗日积分原理[4],可得:
∂ϕ∂t+12∇ϕ⋅∇ϕ+gη=0,η(x,y,t)=z。 |
与传统边界元法相比,广义边界元法能够处理更加复杂的流体流动现象。广义边界元法基于边界积分方程,通过将流体领域划分为边界和内部2个区域,将流体力学问题转化为求解边界上的积分方程。在边界上,通过离散化边界上的节点,并利用边界元法的基本原理,将积分方程转化为代数方程组。在内部区域,通过使用有限元法或有限差分法等数值方法,求解流体的速度和压力分布。
广义边界元法的优点在于能够处理复杂的流体流动问题,如自由表面流动、多相流动、湍流等。由于只需要在边界上离散化节点,相对于有限元法等方法,广义边界元法的计算量较小,适用于大规模问题的求解。此外,广义边界元法还可以方便地处理边界条件的变化和流体领域的变形等问题。
流体广义边界元的原理模型如图3所示。
![]() |
图 3 流体广义边界元的原理模型 Fig. 3 Principle model of fluid generalized boundary element |
流体广义边界元的边界包括粘性边界、流固边界和流气边界3种,建立流体广义边界元模型为:
12P(rP)=∫(G(rP,rQ)∂2P(rP)∂n−1)ds, |
流体边界的法向导数与速度有下式关系:
∂P(rQ)∂nQ=−jρ0ω(rQ), |
流体边界元模型可简化为:
−12P(rP)=∫(G(rP,rQ)jρ0ω(rQ)−1)ds。 |
式中:
流体广义边界元法的计算流程如图4所示。
![]() |
图 4 流体广义边界元法的计算流程 Fig. 4 The calculation flow of fluid generalized boundary element method |
1)网格划分。将流体领域划分为边界和内部2个区域。在边界上离散化节点,并根据问题的要求选择适当的节点分布方式。
2)边界积分方程。根据流体力学问题的边界条件和物理方程,建立边界积分方程,通常使用格林公式将流体力学问题转化为边界上的积分方程。
3)离散化。将边界上的积分方程离散化,将积分方程中的未知量和积分变量用离散的节点和权重表示。通常使用基函数展开边界上的未知量,并将积分方程转化为代数方程组。
4)内部区域求解。在内部区域使用其他数值方法,如有限元法或有限差分法,求解流体的速度和压力分布。通常将流体的速度和压力表示为基函数的线性组合,并根据流体的动力学方程和边界条件建立代数方程组。
5)耦合求解。将边界上的代数方程组和内部区域的代数方程组耦合起来,通过迭代方法求解整个流体力学问题。
6)后处理。根据求解得到的速度和压力分布,可以计算流体的其他物理量,如流量、涡量等,并进行结果的分析和评估。
2.2 基于广义边界元法的高速排水型船舶尾浪数值计算本文基于流体力学有限元分析软件Fluent进行高速排水型船舶的尾浪数值仿真,关键步骤如下:
1)广义边界元模型
建立船舶的广义边界元模型如下:
12P(t)=∫(G(x,t)∂2P(t)∂t−1)dt。 |
式中:
压力与船舶角速度关系满足:
∂P(t)∂t=−jρ0ω(t), |
可得简化后的广义边界元模型:
f(t)=12∫(G(x,t)jρ0ω(t)−1)dt。 |
式中,
2)有限元建模
为了提高高速排水型船舶的尾浪仿真精度,在进行有限元建模时,重点对以下2个方面进行网格的特殊处理:
1)网格细化。对船舶尾浪区域的网格进行细化处理,以更准确地捕捉尾浪的细节。尾浪区域通常是船舶尾部周围的区域,这是尾浪形成和传播的主要区域。通过细化网格,可以更好地模拟尾浪的形状和传播过程。尾浪区域的网格细化精度为10~30 μm逐渐过渡。
2)自适应性网格。Fluent软件[5]支持网格的自适应性调整,在尾浪区域,可以使用网格适应性技术,根据流场的变化自动调整网格的密度。网格适应性可以根据流场的梯度、涡度等参数来判断网格的密度是否足够,如果不够,则自动细化网格;如果足够,则可以保持网格的密度不变或进行适当的稀疏。通过网格适应性,可以更好地捕捉尾浪形成和传播过程中的细节和变化。
建立高速排水型船舶尾浪仿真的有限元模型如图5所示。
![]() |
图 5 高速排水型船舶尾浪仿真的有限元模型 Fig. 5 A finite element model for the simulation of stern waves of high speed drainage ships |
3)入口边界参数施加
入口边界参数决定了有限元仿真是否发散和收敛,入口边界条件参数主包括海水速度,压力等,仅考虑沿水平ox方向的平面波,模型为:
ζ(x,t)=acos(kx−ω⋅t)。 |
式中:
ω2=kgtanh, |
式中,h为水深。此时,平面波的速度为:
Vp=ω/k=√g/k√tanh(kh)。 |
4)求解
给定流场的来流速度为15 m/s,在Fluent求解器中高速排水型船舶的尾浪特性仿真曲线如图6所示。
![]() |
图 6 高速排水型船舶的尾浪特性仿真曲线 Fig. 6 Simulation curve of tailwave characteristics of high speed vessel |
图中可见,随着与船舶距离的增加,尾浪的波高从0.7 m逐渐增加至3 m左右。
3 结 语高速排水型船舶的尾浪对环境和自身动力学性能均有影响,本文针对高速排水型船舶的尾浪特性仿真进行研究,通过建立高速排水型船舶的流动动力学模型,结合广义边界元法,在有限元仿真软件Fluent中进行有限元建模、边界模型建立、求解等,完成了高速排水型船舶的尾浪数值模拟与仿真。
[1] |
王怀兵, 徐新, 陈姚节. 基于粒子系统的船行波三维仿真[J]. 计算机技术与发展, 2021, 31(8): 134-138. WANG Huai-bing, XU Xin, CHEN Yao-jie. Three-dimensional simulation of ship traveling wave based on particle system[J]. Computer Technology and Development, 2021, 31(8): 134-138. |
[2] |
郭亚飞, 姚怡超, 陆颖等. 基于单视点视频的舰艇尾浪三维重建[J]. 系统仿真学报, 2015, 27(10): 2475-2482. GUO Ya-fei, YAO Yi-chao, LU Ying, et al. 3D reconstruction of ship wake based on single-viewpoint video[J]. Journal of System Simulation, 2015, 27(10): 2475-2482. |
[3] |
殷宏, 叶伟, 张宏军等. 破障艇作业可视化仿真与效能评估[J]. 系统仿真学报, 2009, 21(4): 1066-1070. YIN Hong, YE Wei, ZHANG Hong-jun, et al. Visual simulation and efficiency evaluation of wreck-breaking boat[J]. Journal of System Simulation, 2009, 21(4): 1066-1070. |
[4] |
任俊生, 杨盐生, 刘秀文. 高速水翼船操纵模拟器中运动数学模型的研究[J]. 系统仿真学报, 2005(2): 316-318+336. REN Jun-sheng, YANG Yan-sheng, LIU Xiu-wen. Research on Mathematical Model of Motion in High-speed Hydrofoil Maneuvering Simulator[J]. Journal of System Simulation, 2005(2): 316-318+336. |
[5] |
卢本卓, 陈慰祖, 王存新. 边界元法与广义Langevin动力学相结合的模拟方法[J]. 科学通报, 2000(14): 1482-1486. LU Ben-zhuo, CHEN Wei-zu, WANG Cun-xin. Simulation method combining boundary element method and generalized Langevin dynamics[J]. Chinese Science Bulletin, 2000(14): 1482-1486. DOI:10.3321/j.issn:0023-074X.2000.14.004 |