STARCCM+在船舶波浪增阻预报的应用公开版.docxVIP

上海船舶研究设计院

StarCCM+在船舶波浪增阻预报中的应用

张亚晖

2019/7/18

CUSTOMIZEDDESIGN,INNOVATIVETECHNOLOGY

目录CONTENTS

1

背景

2

计算模型

3

算例展示

4

结论

1

背景

CUSTOMIZEDDESIGN,INNOVATIVETECHNOLOGY

背景

绿色运营

?性能指标：静水阻力+波浪中阻力

?油耗预报、航线规划等

关键:波浪增阻的预报

安全运营

?散货船、油船“最小推进功率”

?客船“安全返港”

背景

波浪增阻预报方法

?水池试验

?数值计算

?线性势流方法如Sesam、二维切片方法

?非线性势流方法如ShipFLOW

?CFD方法如StarCCM+、Fluent、OpenFOAM

2

计算模型

CUSTOMIZEDDESIGN,INNOVATIVETECHNOLOGY

计算模型

数值模型

?常规计算设置

湍流模型：Standardk-e运动模拟：DFBI+OVERSET

计算域、边界条件

?背景计算域：船前:1λ(1Lpp)船后：4λ(2Lpp)

水深：1λ(2Lpp)横向：2Lpp

阻尼消波区长度：1.5~2λ

?OVERSET：1.5Lpp*B*2D

计算模型

网格划分波浪是主角，船是配角

?自由面网格

加密区域：静水面上下各一倍波高

一个波高内10~20层网格，一个波长内60~120个网格

网格长宽比：1:4、1:8

?OVSRSET尽量保证与背景网格一一对应

?网格数300~600万

数据处理方法

傅立叶级数方法

短期预报方法

某时间变量写为傅立叶级数的形式：

?

?0

??=+??cos???+??,?=1,2,…

2

?=1

其中

不规则波中平均波浪增阻、运动有义值

?2?

?=2????????

?1

?2

??=2

?1

?(?)2?????

??=??2+??2

其中，Raw为规则波中波浪增阻传递函数

?

?

?1

?=???(

?

)

?

?

?

2

?

??=

??cos?????

0

?

2

??=?

??sin?????

?

0

F为规则波中运动传递函数

S为波浪谱密度，选用Jonswap谱

??

155H944

2

S()1/3exp??

??

?????

TT

4?54?4

??

0101

????2?

??0.191?T1?

exp

????

01

?2??

??

??

3

计算算例

数值波浪场

运营航速下波浪增阻预报

“最小推进功率”所需波浪增阻预报

CUSTOMIZEDDESIGN,INNOVATIVETECHNOLOGY

算例

数值波浪场

?波浪参数：λ=4.0m，H=0.08m

综合考虑计算精度、网格数、计算耗时

dt=1/400Tdz=1/16Hdz/dx=1/8

?计算域：20m*(8+2.5)m*2cell

?波高测点：距入口4m

dz(H)Waveamplitude(m)Err

1/100.0388-2.99%

自由面网格厚度dz

?dt=1/400T，dz/dx=1/8

1/120.0390-2.48%

1/160.0391-2.16%

1/180.0389-2.66%

1/200.0386-3.54%

TimestepWaveamplitude(m)Err

时间步长

1/3000.0383-4.29%

1/4000.0391-2.16%

?dz=1/16H，dz/dx=1/8

1/5000.0393-1.65%

1/6000.0396-1.08%

算例-运营航速

BC1-12kn

Lpp185.0m

计算工况

B32.2m

?/???0.420.560.700.851.01.121.271.421.551.701.832.25

T12.5m

V12kn

H0.02Lpp

?1:52

1.01.2

0.8

1.0

8EXP-

CurrentMethod

V=12kn

0.8

0.6

6

EXP-

CurrentMethod

V=12kn

0.6

0.4

EXP-

CurrentMethod

V=12kn

0.4

4

0.2

2

0.2

0.0

0.0

0

0.51.01.52.02.5

?/Lpp

0.51.01.52.02.5

?/Lpp

0.51.01.52.02.5

?/Lpp

算例