非静压模型NHWAVE学习（9）——波浪破碎RANS紊流模型算例学习（Wave Breaking on a Planar Beach）

本blog介绍了NHWAVE模型自带算例 Regular Depth-Limited Wave Breaking on a Planar Beach 的编译运行；并分享了该算例配置文件、模拟结果的学习记录。
关于该算例的详细描述见 Morteza Derakhti 等人的论文（Derakhti et. al, 2016-a¹, 2016-b²）。

此外，本人在NHWAVE原版源代码基础上，更新了垂向网格分层百分比可指定、被动示踪剂输移、温度场输入输出等模块，并修正了植被水流模型中的一些BUG。有兴趣的朋友请上我的GitHub下载代码。

欢迎各位朋友对此提出指导意见；大家可以一起交流，共同进步！

算例简介

该算例模拟的是孤立波靠近斜面时发生的变形与崩破（Spilling Breaking）。模型中运用了 k-ε 紊流模型，可用于检验NHWAVE中紊流模型的模拟性能。相应的物理实验详见 Ting and Kirby (1994) ³。
模拟所需参数文件和初始水深文件在代码包 /examples/tingkirby 目录下。

地形与网格

本算例的计算域是一个垂向二维的水槽，如下图所示。在初始时刻，孤立波从左侧边界进入水域（Wavemaker）；计算域右侧是一个斜坡面，其立面为一直角梯形，坡度比为1：35。水槽总长22.0m，斜坡坡脚据左边界 5.6m。初始时刻水深 h₀和坐标系如下图中所示。

在模型中，我们设置880个水平均匀网格（Δx = 0.025 m）和20个均匀的垂向σ层。以上网格的部分我们在input.txt中设定。

! --------------------DIMENSION---------------------------------
! cell numbers
Mglob = 880
Nglob = 1
Kglob = 20

! ------------------------GRID----------------------------------
! grid sizes
DX = 0.025
DY = 0.025
! ---------------------VERTICAL GRID OPTION--------------------
! IVGRD = 1: uniform; 2: exponential
IVGRD = 1
GRD_R = 1.1

运行参数配置

首先需要确定边界条件。在本例中，左侧为波浪入射边界，底面边界设置阻力系数；右侧及地边界均为滑移边界（free-slip）。边界条件设定见 input.txt 的 BOUNDARY_TYPE 部分。

! -------------------BOUNDARY_TYPE--------------------------------
! bc_type=1: free-slip
!         2: no-slip
!         3: influx
!         4: outflux (specified eta)
!         5: bottom or wall friction
!         6: radiation bc
BC_X0 = 3
BC_Xn = 1
BC_Y0 = 1
BC_Yn = 1
BC_Z0 = 5
BC_Zn = 1

其中的波浪参数在 input.txt 的 WAVEMAKER 部分中设定，底部边界的摩阻力参数在 input.txt 的 BOTTOM ROUGHNESS 部分中设定。

计算区域左侧进入的孤立波参数为：

1. 波高 H = 0.125 m；
2. 周期 T = 2.0 s；

上述参数的设定如下（对于孤立波，其控制作用的参数只有AMP和DEP）：

! ---------------------WAVEMAKER------------------------------
! wavemaker
! AMP - wave height; PER - wave period; DEP - incident water depth
! THETA - incident wave angle
! LEF_SOL - left boundary solitary wave, need AMP,DEP
! LEF_LIN - left boundary linear wave, need AMP,PER,DEP
! LEF_CON - left boundary cnoidal wave, need AMP,PER,DEP
! LEF_STK - left boundary stokes wave, need AMP,PER,DEP
! LEF_TID - left boundary tide wave, has to specify in subroutine
! LEF_JON - left boundary for JONSWAP spectrum
! RIG_LIN - right boundary linear wave, need AMP,PER,DEP,THETA
! INI_ETA - initial surface elevation specified in subroutine initial
! INT_LIN - internal wavemaker for linear wave
! INT_CON - internal wavemaker for cnoidal wave
! INT_SOL - internal wavemaker for solitary wave
! INT_JON - internal wavemaker for JONSWAP spectrum
! INT_SPC - internal wavemaker for 2D spectrum (need spc2d.txt)
! INT_IRR - internal wavemaker for irregular wave (need irr2d.txt)
! FLUX_LR - impose flux at both left and right boundaries
! FOCUSED - left boundary focusing wave packet (isolated whitecap)
! WAV_CUR - left boundary coexisting waves and currents
WAVEMAKER = LEF_CON
AMP = 0.125
PER = 2.0
DEP = 0.4
THETA = 0.0
CUR = 0.0
sd_return = 0.0

底摩阻采用拖曳力公式计算（详细公式见手册），拖曳力系数C_d设定为0.006。

! ----------------------BOTTOM ROUGHNESS-------------------
! Ibot = 1: given the drag coefficient Cd0
! Ibot = 2: given the bottom roughness height Zob
Ibot = 1
Cd0 = 0.006
Zob = 0.0001
Dfric_Min = 0.0

注：上面的 Ibot = 1 表示本模型的摩阻力系数通过给定拖曳力系数C_d来设定；而 Ibot = 2 表示本模型的摩阻力系数通过给定底面粗糙高度Z_ob来设定。

此外，本算例启用了紊流模型；相关的参数设置见 input.txt 的 VISCOSITY 部分。本算例采用垂向二维 RNG k-ε 模型，分子粘度为 1.0×10^-6，Schmidt数为1.0。

VISCOUS_FLOW = T
IVTURB = 10
IHTURB = 10
PRODTYPE = 3
VISCOSITY = 1.e-6
Schmidt = 1.0
Chs = 0.001
Cvs = 0.001
RNG = T

以上各参数的含义详见模型手册。

模型的编译运行

首先，将代码包 /examples/tingkirby目录下的 input.txt 和 depth.txt 复制到 nhwave 所在的运行文件夹内。之后设定并行相关的参数，编译运行程序。

模拟结果

物理实验与NHWAVE模拟的对比结果详见 Derakhti 等人的论文^1&2。

---

1. Derakhti M , Kirby J T , F Shi, et al. Wave breaking in the surf zone and deep-water in a non-hydrostatic RANS model. Part 1: Organized wave motions[J]. Ocean Modelling, 2016, 107:125-138. ↩︎ ↩︎

2. Derakhti M , Kirby J T , Shi F , et al. Wave breaking in the surf zone and deep water in a non-hydrostatic RANS model. Part 2: Turbulence and mean circulation[J]. Ocean Modelling, 2016:139-150. ↩︎ ↩︎

3. Ting F , Kirby J T . Observation of undertow and turbulence in a laboratory surf zone[J]. Coastal Engineering, 1994, 24(1-2):51-80. ↩︎