基于STARCCM+的船舶尾部流动问题研究公开版.docxVIP

西门子船舶/海工行业CFD专题研讨会上海7.17

吴铁成、邓锐、骆婉珍

中山大学海洋工程与技术学院

2019.07.17

Unrestricted

目录

一、基于CFD-PIV的尾流场速度与旋涡空间分布研究

二、基于CFD-DEM的碎冰航道内船舶碎冰绕流研究

三、某型气垫船绕流流动问题初步研究

一、基于CFD-PIV的尾流场速度与

旋涡空间分布研究

1.所应用的方法

1.1数值方法

某型船尾部流动数值模拟基于STAR-CCM+软件进行URANS模拟：

TurbulenceModel:SSTk-omegaturbulencemodel

TreatmentofFreeSurface:HRIC-VOFofSTAR-CCM+

WaveDampingmodelofSTAR-CCM+byChoiandYoon(2009)

…

网格分布

1.所应用的方法

1.2试验方法-PIV

1.2.1测量设备与2D-3C-SPIV系统

HEU拖曳水池示意图及其主要参数

拖曳水池环境下的拖曳水下2D-3CSPIV测量系统

测试视频

1.所应用的方法

1.2.23D-3C尾流场测量与重构

3D-3C尾流场的多切面CT扫掠式重构过程：

a）2D-3C速度切面分布与位置；b）CT扫掠后形成的

体测量区域；c）原始测量3D-3C体数据集合网格分布；

d）插值后3D-3C体数据集合网格分布

多切面2D-3C速度场和CT扫掠重构后的体空间流场：

a）一半区域的多切面2D-3C速度场；b）CT扫掠重构

后的一半区域体空间流场；c）全区域的多切面2D-3C

速度场；d）CT扫掠重构后的全区域体空间流场

2.船舶尾流场数值与试验结果

2.12D桨盘面区域标称伴流分布（2维切面验证）

2.1.1标称伴流场速度与流线分布

标称伴流场速度与流线分

布图PIVVSCFD（右）

2.1.2标称伴流场涡量分布

标称伴流场涡量分布CFD

（左）VSPIV分布图

2.船舶尾流场数值与试验结果

2.2船艉区域3D-3C尾流场速度特性空间分布（3维空间验证）

2.2.1无量纲轴向速度u/U三维空间分布

无量纲轴向速度u/U三维

空间分布图（PIV）

轴向速度三维空间分布PIVVSCFD

（star-ccm+URANS）

2.船舶尾流场数值与试验结果

2.2.2无量纲垂向速度w/U三维空间分布

垂向速度三维空间分布PIVVSCFD（star-ccm+URANS-右）

2.船舶尾流场数值与试验结果

2.2.3无量纲空间速度长度L/U三维空间分布

无量纲空间速度长度L/U

三维空间分布图（PIV）

空间速度长度三维空间分布PIVVS

CFD（star-ccm+URANS）

2.船舶尾流场数值与试验结果

2.3船艉区域3D-3C尾流场旋涡特性空间分布

2.3.1涡量场Vorticity（x）三维空间分布

涡量场Vorticity（x）的

三维空间分布图（PIV）

涡量场Vorticity（x）的三维空间分PIV测

量结果与CFD数值计算结果对比图

2.船舶尾流场数值与试验结果

2.3.2旋涡空间演化特性

不同切面位置处的涡量Vorticity(x)和轴向速度边

界层分布图舭涡与假毂毂帽涡空间演化分布

二、基于CFD-DEM的碎冰航道船

舶碎冰绕流特性

1.所应用的方法

1.1DEM方法

?DEM方法最早由Cundall和Strack[1]创立，可以看作为包含密集粒子流的拉格

朗日多项流的扩展。用来模拟大量致密堆积的不同尺寸大小、不同形状的颗粒行

为。

?能精确和快速模拟大量固体颗粒的运动和接触

以离散元法（DEM）为基础，采用CFD-DEM耦合的方式对船-冰-水作用进行计算。

1.2DEM粒子的选择

?DEM复合颗粒是由很多球形颗粒组合成的刚性颗

粒，作为一个具有更复杂形状的单一粒子。

?升力、曳力和自旋力作用在复合颗粒的质心，这

些力的计算是基于该质心点的流场求解获得。

?每个复合颗粒有一个球形度属性相关的形状，用

来计算颗粒的曳力；

?复合颗粒形状能基于CAD几何自动创建

Ref：1.Cundall,P.A.,Strack,O.D.L.1979.“Adiscretenumericalmodelforgranularassemblies”,Geotechnique,29,pp.47-65.

1.所应用的方法

1.3碎冰航道CFD-DEM分析流程

?几何?网格?物理条件设定

?求解与后处理

RANS

DEM

VOF

Oversetmesh

…

1.所应用的方法

1.4碎冰航道与计