基于RANS方程与大涡模拟法的畸形波数值造波对比研究.pdfVIP

基于RANS方程和大涡模拟法的 畸形波数值造波对比研究 柏君励1，马 宁1’，顾解忡1，一，彭 勃1 (1．上海交通大学船舶海洋与建筑工程学院；2．上海交通大学海洋T程国家重点实验室) 摘要 极端海洋环境下出现的畸形波具有异常大的波高和波陡，具有很强的非线性特征，因此水流黏性流场 流动信息对畸形波的预报十分重要。相对于势流理论，采用CFD的方法能更准确地模拟出畸形波特征。 在CFD方法中，雷诺平均法求解RANS方程，可以获得平均的湍流场；而大涡模拟法对大尺度涡直接进 行数值模拟，可以获得大涡的脉动信息，比RANS方程的模拟更精细。为研究RANS方程和大涡模拟法 对畸形波的模拟效果，本文利用Fluent软件中的这两种模块对畸形波进行模拟。文中采用动网格技术模拟 造波板造波，采用VOF法追踪自由液面。研究采用过渡水波手法实现深水聚焦波，首先对过渡水波波陡 的变化对畸形波的影响进行数值模拟研究，在此基础上选择合适的波陡研究流对畸形波传播及破碎过程的 影响，并对比分析RANS方程模型和大涡模拟法在破碎畸形波的模拟方面的异同。 关键词：畸形波；大涡模拟；RANS方程 1前言 畸形波通常被定义为波高大于两倍的有义波高，且波峰尖锐而波谷平坦的波浪。在各种恶劣海 洋环境中，畸形波以其异常大的波高、巨大的能量以及瞬问出现和消失的特性而成为海洋结构物的 重要威胁之一。近几十年来已有许多学者采用势流理论对畸形波的生成机理和物理特性等进行了深 入的研究【l，2】，但是采用势流理论研究畸形波时，只能模拟弱非线性的波浪，而对于强非线性波浪 问题，例如波浪的翻卷、破碎等，则无法给出合理的结果。实海域中，畸形波的波高极大，往往伴 随着波浪翻卷、破碎、上浪抨击等现象，而在这些复杂过程中，流体的黏性因素起重要作用，因此 采用CFD的方法能更准确地模拟出畸形波特征，并为波浪与结构物相互作用研究奠定基础。 ANSYS CFX软件模拟了聚焦波的产生和变化过程，计算了聚焦波与水平圆柱和垂直圆柱相互作用 问题，给出了圆柱体的受力以及自由面的变化情况，并与试验值进行了对比验证。Zhao等[4]采用 包含七模型的RANS方程模拟了规则的短波和长波，并在此基础上进一步地对畸形波进行了模拟研究。 和水质点速度并与试验结果进行了对比验证。给出了破碎过程中的紊动强度分布、涡场分布以及速 度场分布情况，分析了破碎过程中的能量耗散过程。Cui等[6]采用包含忌模型的RANS方程模拟了具 有变化地形的水槽(包括底部为斜坡以及底部有半圆柱形凸起)中畸形波的传播过程，分析了地形 变化对畸形波的影响，以及高频能量和低频能量的变化情况。Ojieh等【7】采用Fluent软件模拟了一 个带斜坡的深水水槽中的畸形波，给出了质点的速度和加速度沿水深方向的分布图。还模拟了畸形 波与水平圆柱相互作用的问题，给出了惯性力、拖曳力和倾覆力矩沿水深方向的对比图。 387 目前针对畸形波的CFD模拟研究中，基本上都是以RANS方程作为理论基础的，而采用大涡模 拟方法研究畸形波的鲜有报道。RANS方程是经过雷诺平均的N．S方程，它只能提供湍流的平均信 息，且雷诺平均模型不具有普适性。而大涡模拟是通过滤波处理，将小于某个尺度的漩涡从流场中 过滤掉，仅直接计算大涡，然后通过求解附加的方程得到小涡的解。大涡模拟方法在湍流模拟方面， 方程模型和大涡模拟法两种方式对畸形波进行数值模拟研究，探讨畸形波的破碎过程以及流对畸形 波的影响等，并对比分析两种方法的特点。 2控制方程 上述两种模型理论介绍如下【8】。 2．1 雷诺平均N．S方程模型介绍 iop+兰(肛，)=o (1) c2， 昙(腭)+专(以¨=，一嘉+苦心堕axj+簧一吾色筹]]+毒(一p巧) 上式中U，7和U，’表示脉动速度分量，方程中其余各量(p、P、U，、U，)为雷诺平均量。式中 最后一项一PUt’Uj7为雷诺应力项，该项采用