力学学报, 2019, 51(6): 1712-1719 DOI: 10.6052/0459-1879-19-321

海洋工程专题

波浪破碎的一种混合湍流模拟模式1)

詹杰民,2), 李熠华

中山大学应用力学与工程系,广州 510275

A HYBRID TURBULENCE MODEL FOR WAVE BREAKING SIMULATION1)

Zhan Jiemin,2), Li Yihua

Department of Applied Mechanics and Engineering, Sun Yat-sen University, Guangzhou 510275,China

通讯作者: 2) 詹杰民,教授,主要研究方向:流体力学. E-mail:stszjm@mail.sysu.edu.cn

收稿日期: 2019-10-23   接受日期: 2019-11-5   网络出版日期: 2019-11-18

基金资助: 1) 国家重点项目.  6140206040301
广东省科技计划项目资助.  2016A050502022

Received: 2019-10-23   Accepted: 2019-11-5   Online: 2019-11-18

作者简介 About authors

摘要

近岸波浪的变形与破碎,一方面影响水体标识码运输,另一方面对消波护岸具有指导意义.本文提出一种三维混合湍流模拟模式方法,将求解区域分为造波区、波浪传播区和消波区. 造波区采用层流模式,通过基于 Fluent 的二次开发的 UDF 方法在边界进行速度造波.这种方法在给定入口速度的条件标识码据已知波高进行 精准插值,从而控制水的体积分数. 在波浪传播区域,采用大涡标识码行模拟研究,在消波区,采用 RANS 模型并利用多孔介质消波法进行消波. 模标识码 VOF 方法捕捉波浪破碎过程中的自由面变化.本文对波高为 5.5cm 的规则波 (M1)、波高为 13.5cm 的 规则波 (M3)、有效标识码 7.75cm 的 TMA 谱单向不规则波 (U1) 和有效波高为 19cm 的 TMA 谱多向不规则波 (B5) 展开标识码研究,并与前人的相关标识码果作比较,各条件下模拟结果与实验结果吻合.模拟结果说明本文提出的模型能够准确模拟出波浪传播过程中的折射和标识码象,并且能够捕捉的波浪破碎过程中的自由面变化,为三维波浪的传播与破碎的数值模拟提供一种模拟方法.

关键词: 波浪破碎 ; 湍流 ; 数值模拟 ; 混合模式 ; 大涡模拟

Abstract

On the one hand, the deformation and breaking of near-shore waves affect the transportation of water and sediment, and on the other hand, it is of guiding significance for wave elimination and revetment. In this paper, a 3-dimensional hybrid turbulence simulation model is proposed. Laminar flow mode is adopted in the wave-making region, and velocity wave-making is carried out at the boundary through the User Defined Function(UDF) developed based on fluent. This method can control the volume fraction of water by precise interpolation according to the known wave height under the condition of given inlet velocity. In the wave propagation region, large eddy simulation (LES) is used for simulation research, and in the wave elimination region, RANS model with the porous media wave absorber is used for wave elimination. The model uses the VOF method to capture the free surface changes in the process of wave breaking. This paper carried out simulation studies on the regular wave (M1) with a wave height of 5.5cm, the regular wave (M3) with a wave height of 13.5cm, the unidirectional irregular wave (U1) with a TMA spectrum with an effective wave height of 7.75cm and the multidirectional irregular wave (B5) with an effective wave height of 19cm. Simulation results show that the proposed model can accurately simulate the wave propagation in the process of refraction and diffraction phenomenon, and be able to catch the waves of the free surface in the process of change, broken for 3D wave propagation and broken numerical simulation provides a method of simulation.

Keywords: wave breaking ; turbulence ; numerical simulation ; hybrid model ; large eddy simulation

PDF (5979KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

詹杰民, 李熠华. 波浪破碎的一种混合湍流模拟模式1). 力学学报[J], 2019, 51(6): 1712-1719 DOI:10.6052/0459-1879-19-321

Zhan Jiemin, Li Yihua. A HYBRID TURBULENCE MODEL FOR WAVE BREAKING SIMULATION1). Chinese Journal of Theoretical and Applied Mechanics[J], 2019, 51(6): 1712-1719 DOI:10.6052/0459-1879-19-321

引 言

波浪破碎是重要的水动力学现象之一,在海洋环境与近岸工程存在不可忽视的影响.许多学者对波浪破碎的影响开展了深入的研标识码孙宝楠等[1]将波浪破碎引起的油水混合模型加入海洋模型,研究了溢油事件后的海面残油量变化,分析了波浪破碎对残油量的影响; Kim 等[2]研究了陆架地形和波浪破碎对 N 型海啸的影响. Ravindar 等[3]通过实验研究不同类型的破碎波对竖直墙标识码击压力的差异,结果表明不同类型的破碎波造成的冲击压力在大小与位置上均有显著差异.可见波浪破碎的研究对环境和工程都有重大意义.

国内外学者使用了多种不同的方法研究波浪破碎问题.这些方法大致可分类为数值模拟和水槽试验两类.在数值模拟方标识码 20 世纪 60 年代,Peregrine[4]对非均匀水深条件下的波浪问题提出了数学模型,许多研究波浪问题的手段都在标识码以无黏的势流理论为基础.随着理论与技术的发展,近年来学者运用新的手段对波浪破碎进行了研究.郑茜和张陈浩[5]、 张陈浩和郑茜[6]在 OpenFOAM 软件中将RANS模型加入动量阻尼方程,对斜坡波浪破碎实验进行模拟,数值模拟结果与相应的 PIV 标识码果互相吻合,捕捉到了波浪破碎过程中发生的涡旋扩散现象.许璐璐和郄禄文[7]通过 SPH 方法分析了波浪破碎后的波速标识码对基床的波压等水动力性能的影响,研究表明随着海底坡度增大,波浪破碎形态有从崩破波转变为卷破波的趋势.标识码[8] 应用不忽略动压的NS方程,以 S-A 模型为湍流模型,基于有限体积法,建立了基于非静压模型的波浪破碎模式,分别模拟了斜标识码台地形附近的波浪破碎情况,模拟结果表明该模式精度与计算效率良好. Liu 等[9]则采用特定的浅水波波动色散方标识码 Camassa-Holm 系统进行了波浪破碎分析. Luo 等[10]通过引入涡黏性的非线性 Boussinesq 方程,应用非结构有限单元标识码拟了多向波的传播与破碎.在水槽试验方面,诸裕良等[11]在复合坡度的珊瑚礁上开展实验标识码波浪破碎的类型、波高、位置等特性开 展了深入的研究. 柳淑学等[12]对特定坡度的珊瑚礁地形波浪破碎过程进行了研究并给出了该地形下波浪破碎后波高的计算公式. 对波浪破碎的研究较多以二维单向波为研究对象,柳淑学等[13]通过水槽实验表明不同波浪入射方向对波浪破碎的波高存在影标识码明三维的波浪破碎问题仍有较大研究 空间.

本文拟在三维角度上对波浪经过椭圆浅滩的破碎问题进行数值模拟研究. Vincent 和 Briggs[14]对该地形就多向波的传播与破标识码进行了实验.该实验研究了水底椭圆浅滩附近的折射与绕射问题,同时提供了破碎与不破碎情况的实验结果. 不少数值模拟的研究通标识码浅滩的案例来验证数值模型的效果.詹杰民等[15]基于 Chebyshev 广义有限谱建立了求解 Boussinesq 方程的高精度模式,良标识码拟了该实验非破碎的情况. Gouin 等[16]也针对该实验结果用高阶谱方法进行了研究.有关混合湍流模式的研究, Zhan 等[17]尝试在 Boussinesq 方程的水平方向上加入大涡模拟模式,该混合模式成功捕捉到了波面的能量耗散情况,但是标识码面的详细破碎情形模拟则仍有不足. Zou 和 Zhang[18]改善了波 浪-标识码型的消波边界条件并在椭圆浅滩案例上标识码验证,研究结果表明改善的消波条件能使模拟结果更好. Kang 和 Guo[19]等用一种交错方向隐格式的非静力模型模拟了椭标识码的算例,结果表明该模型对波浪的折射与集中模拟良好,但是会低估绕射的效果. Zhan 等[20]在对椭圆浅滩算例的模拟中标识码了不同湍流模式并对比了模拟效果. Griffiths 和 Porter[21]对该地形引起的波浪聚焦效应展开了重点研究. Choi 和 Yoon[22]通过 fluent 运用 RANS 模型模拟了相关波浪经过椭圆浅滩时的变化,并研究了波浪传播过程中的速度场和湍动能. Yin 等[23]提出一种混合有限体积和有限差分方法,并以该地形上的波浪传播情况对方法进行了验证. Smit 和 Janssen[24]通 过辐射传播方程模型模拟了波浪在椭圆浅滩上的传播,讨论了该方法对波浪研究的适用性. 胡克林和丁平兴[25]应用SWAN模型研标识码浪在椭圆浅滩地形上的衍射效应.陈小文等[26]在椭圆浅滩地形上运行了 WABED 模型和 SWAN 模型,并评估了两种模型对衍射的模拟效果. Hajime[27]对该地形上多向传播的波浪进行了模拟. Janssen 和 Herbers[28]则研究了非线性波在椭圆浅滩地形上的聚焦效应. 金红和邹志利[29]建立具有精确色散性的非线性波浪模型,并在对椭圆浅滩地形上的波浪模拟中验证模拟效果.

本文针对三维波浪的传播与破碎问题建立一种混合湍流模型.将模拟结果与 Vincent 和 Briggs[14]的浅滩实验结果进行对标识码验混合湍流模拟模式的模拟效果,为三维波浪传播与破碎的数值模拟提供一种新的方法.

1 数学模型与方法

本文采用基于 fluent 进行二次开发的一种包含造波、模拟和消波模块的混合数值模拟方法,将求解区域分为造波区、传播标识码波区 3 个区域,通过 VOF 模型捕捉波浪的自由面.其中,造波区采用层流模式并通过 UDF 模块进行边界造波.波浪传播区为数值标识码核心区域,采用大涡模拟模型,消波区采用 RANS 湍流模型,在相应区域中增加多孔介质,以达到消波的效果.通过对不同算例标识码模拟结果与实验结果进行对比,以确定模型结果的准确度,同时在破碎算例结果中研究数值模型的波浪破碎效果.

1.1 基本方程

本文研究波浪经过椭圆浅滩的波浪变化和破碎问题,以 $xy$ 平面为水平面,$x$ 轴正方向为波浪行进方向,$z$ 轴负方向为重力方向. 控制方程为笛卡尔坐标系下以速度和压力为变量的不可压缩黏性流体的 Navier-Stokes 方程和连续性方程

$\left. \begin{array}{l@{\quad}l} \rho \Big ( \dfrac{\partial u}{\partial x}+\dfrac{\partial v}{\partial y} \Big ) = S_{\rm m} \\ \rho \Big ( \dfrac{\partial u}{\partial t}+u\dfrac{\partial u}{\partial x}+u\dfrac{\partial u}{\partial y} \Big )=- \dfrac{\partial p}{\partial x}+\mu \Big ( \dfrac{\partial^2u}{\partial x^2}+\dfrac{\partial^2 u}{\partial y^2} \Big ) +S_x \\ \rho \Big ( \dfrac{\partial v}{\partial t}+u\dfrac{\partial v}{\partial x}+u\dfrac{\partial v}{\partial y} \Big )=- \dfrac{\partial p}{\partial y}-\rho g +\mu \Big ( \dfrac{\partial^2v}{\partial x^2}+\dfrac{\partial^2 v}{\partial y^2}\Big ) +S_y \end{array} \right\}$

其中,$\rho$ 为流体密度,$u$ 和 $v$ 分别为 $x$ 和 $y$ 方向的速度分量, $p$ 为压强,$\mu$ 为动力黏性系数,$g$ 为重力加速度, $S_{\rm m}$, $S_x$ 和 $S_y$ 分别为附加质量源项、$x$ 和 $y$ 方向的动量源项.本文将求解区域分为造波区、传播区与消波区,分别对应的层流模式, 大涡模拟模型和 RANS 型 Relizeable $k$-$\varepsilon $ 湍流模式,消波方法采用多孔介质消波法,具体详见文献[20].本文不具体展开相应表达形式.

自由面的追踪采用 VOF 方法.体积分数函数 $\alpha_q$ 定义为单元内第 $q$ 相流体所占有体积与该单元的体积之比. 即若 $\alpha_q=1$,则表示单元内全部为第 $q$相流体;若 $\alpha_q=0$,则说明单元内没有第 $q$ 相流体;若 $0<\alpha_q<1$,则该单元称为交界面单元. 对于水波问题,由于只有气液两相,故 $q=1$ 或 2. $\alpha_{2}$ 满足方程

$\left. \begin{array}{l@{\quad}l}\dfrac{\partial \alpha_q}{\partial t}+\dfrac{\partial(u\alpha_q )}{\partial x}+ \dfrac{\partial(v\alpha_q )}{\partial y}+\dfrac{\partial(w\alpha_q )}{\partial z}=0\\ \sum^2_{q=1} \alpha_q=1, \ q=1,2 \end{array} \right\}$

1.2 造波方法

无论是线性规则波,非线性规则波还是不规则波,本文基于 Fluent 造波方法的实现方式是相同的,就基于自定义函数 (UDF) 模式,采用速度造波方法,进行入口边界造波.在技术实现的层面,根标识码理论给定入口面水体部分速度分布,为了保证计算的协调性,入口面的空气部分的速度采用水气交界面水的速度,这样就保标识码体速度的连续性.由于波浪理论能够完整给出自由面的瞬时状态,本文采用 UDF,在水气交界面采用线性插值,准确给出了入标识码面单元的体积分数,从而准确解决的入口处波高变化对波浪场的影响.

(1) 对于线性规则波,波高和速度分量由线性理论给出.

(2) 对于非线性规则波,按非线性波理论 (如 Stokes 波和椭余波理论) 给出入射波的波高和速度分量.

(3) 对于不规则波,水位和速度用规则波叠加求得,例如,对于水位,有

$\varsigma _I (x,y,t) = \sum_{m = 1}^{M_f } \sum_{n = 1}^{N_\theta } a_{mj} \cos \Big [k_{mj} (x\cos \theta _j + y\sin \theta _j ) -\\ 2\pi f_{mj} t + \varepsilon _{mj} \Big]$

其中

$f_{mj} = \hat {f}_m - \dfrac{1}{2}\Delta f + (j - 1 + RAN_{mj} )\Delta f / N_\theta$
$a_{mf} = \sqrt {2\hat {S}(\hat {f}_m ,\theta )\Delta f\Delta \theta }$
$\Delta f = (f_H - f_L ) / M_f , \ \ \Delta \theta = (\theta _{\max} - \theta _{\min } ) / N_\theta$
$\hat {S}(f,\theta ) = S(f)G(f,\theta )$

式中,$S$ 为频率谱,$G$ 为方向分布函数,其余变量含义参见文献 [14,30].

2 模拟案例与计算结果

数值模型参照 Vincent 和 Briggs[14]的浅滩实验建立.模型长 23m,宽 23.76m,高 0.6097m,初始网格 $\Delta x=0.05$m, $\Delta y=0.2$m, $\Delta z=0.02$m,在靠近边界与空气-水交界面的地方存标识码应加密网格. 时间步长为 0.02s,对于有效波高为 19cm的多向不规则波算例 (B5),时间步长取 0.005s. 椭圆浅滩中心位于 $x=6.1$m, $y=11.88$m 处,椭圆浅滩曲面方程为

$\Big ( \dfrac{X’}{3.05} \Big )^2+\Big (\dfrac{Y’}{3.96} \Big )^2 =1$

式中,$X’=x-6.1$, $Y’=y-11.88$. 数值水槽平坦处底部 $z$ 坐标为 $-0.4572$,浅滩范围内水槽底部 $z$ 坐标 $Z_{\rm b}$ 由方程 (9) 确定

$Z_{\rm b}=-0.4572+0.762\times \Big [1-\Big ( \dfrac{X’}{3.81} \Big )^2-\Big ( \dfrac{ Y’}{4.95} \Big )^2 \Big ]^{0.5} $

在与实验对应截面位置设置报告监测点,模型中浅滩位置与测量端面位置如图1 ,将相应截面所测得的相标识码与实验结果对比,以验证混合模式模型结果的可靠性.

图1

图1   椭圆浅滩的位置和测量断面

Fig.1   Locations of eillptical shoal and measurement transects


模型中造波区长 2.8m 位于模型前端,消波区长 3m 位于模型末端,其余区域为模拟波浪传播和破碎的湍流区.造波区标识码区、消波区的分布如图2 所示.

图2

图2   模型建立与分区

Fig.2   Model set up and area divided


运用混合模式分别对 Vincent 和 Briggs[14]的浅滩实验中相应编号的实验进行了数值模拟,本文模拟的实验条标识码为规则波初始算例 M1,规则波破碎算例 M3,不规则波初始算例 U1 以及不规则波破碎案例 B5,各算例编号与对应条件见表1.表中, $\gamma $ 与 $\sigma_{\rm m}$ 分别为实验中不规则波能量谱与方向谱中对应的参数,具体的能量谱与方向谱可参见文献[14] 与文献[30].

表1   各算例条件

Table 1  Conditions for case series

新窗口打开| 下载CSV


2.1 规则波计算结果与分析

Vincent 和 Briggs[14]的实验中,在 $x=3.05$m, $y=y_{\rm m}+7.62$m 处设置了波高计用以记录实验中的参考波高,$y_{\rm m}$ 为 $y$ 方向的中间点,图1中各个标号截面均匀分布 9 个监测点,将每个监测点记录得到的平均或标识码高 $H$ 与参考波高 $H_0$ 的比值 $H/H_0$ 作为相对波高.在以往对椭圆浅滩的模拟研究中,大多对算例 M1 或 M2 进行验证. 将 M1 的模拟结果与实验结果标识码与邹志利[29]的提出的非线性波浪模型在文献[31] 中的模拟结果以及 Vincent 和 Briggs[14]的实验结果标识码比,以验证混合湍流模型的模拟改进效果.

算例 M1 截面 4 实验与数值模拟结果对比如图3 所示.

图3

图3   算例M1截面4相对波高

Fig.3   Relative wave height at plane 4 (Case M1)


图3 可以看出,相较于文献[29] 的非线性波浪模型,混合湍流模型更加贴近实验结果.从相对波高分布看出,波标识码过地形后到达截面 4 标识码中间的相对波高能抬高到 2 倍以上,而在正中间往两侧偏 1.535m 处出现相对波高的低谷,

低谷处相对波高大约为 1/2. 混合标识码规则波不破碎案例的模拟结果与实验结果吻合良好,较为准确地计算出波浪在经过地形后相对波高的变化情况.为研究分析标识码过椭圆地形时的流场变化情况,将算例 M1 运行至 16s 时位于 $x \in [3.15, 12.2]$, $y \in [8.83, 14.93]$ 波浪自由面附近的流线图如图4 所示.

图4

图4   指定自由面流线俯视图(算例 M1)

Fig.4   Vertical view of streamlines from specified surface (Case M1)


图4 左边界为 $x=3.05$ (浅滩起始位置),右边界为 $x=12.2$ (截面 4 位置),波浪中水的质点在 $xz$ 平面内做环形运动标识码面俯视图中的流线呈由左到右的平行线条,由流线变化可以看出,波浪在行进过程中逐渐向中间聚拢,最终在截面 4 标识码图3 相对波高的结果.算例 M1 的结果表明,混合模式能够良好地模拟波浪经过浅滩地形时的折射与绕射情况.

算例M3实验与数值模拟结果对比如图5所示.

图5

图5   算例 M3 截面 4 相对波高

Fig.5   Relative wave height at plane 4 (Case M3)


图5 实验数据结果表明破碎案例中截面 4 正中间处相对波高相较不破碎案例大大减小,相对波高值大约为 1,在正中间标识码两侧偏 0.7625m 处出现相对波高的一个低谷,随后相对波高沿外逐渐增加,在距中心点 2.2875m 处达到相对标识码极大值,大约为 1.3.对比结果表明数值模型对规则波破碎的案例的计算结果也能较为准确地贴近实验结果,但在截面4两 侧模拟结果较实验结果偏小,数值模型能够捕捉波浪发生破碎时的自由面形态,如图6 所示.

图6

图6   波浪破碎处自由面速度云图(a) 和 $y=10$ 处相云图(b)

Fig.6   Contours of velocity on free surface at breaking site (a) and contours of phases at $y=10$ (b)


图6(a) 为在波浪破碎处 ($x \in [7.15, 7.6]$, $y \in [9, 14.76 ]$) 速度云图的俯视图与斜视图, 图6(b) 则为 $x \in [3.05, 9.15]$, $y=10$ 处的气液相云图,图6(a) 能从三维角度观测到模式捕捉到的自由面所展示的波浪破碎效果,图6(b) 则在二维角度下表现标识码经过浅滩前后发生的变化,表明混合模式能够计算并捕捉到波浪越过浅滩后破碎时扩散、卷折的形态.

为分析波浪破碎过程的自由面流场变化情况,将与图4 相同位置的自由面附近的流线图绘制如图7 所示.

图7

图7   指定自由面流线俯视图(算例M3)

Fig.7   Vertical view of streamlines from specified surface (Case M3)


图7 流线图的左半部分,即波浪在跨越浅滩时,相较于不破碎的算例,流线出现了分团的现象,在图6 的中间为波浪破碎的位标识码度在这里变大,流线变长且集中向后延伸,在图7 的方框标记处,流线向上下两侧有明显扩散,且有些许流线沿 $z$ 方向螺旋标识码伸,与不破碎情况不同,在方框处能看到投影在 $xy$ 平面上呈环形的流线.这些流线特征能表现出波浪破碎过程中的涡旋扩散现标识码一定程度上解释了算例M3截面4中间相对波高不像算例M1般集中抬高的原因.

2.2 不规则波计算结果

单向不规则波算例 U1 和多向不规则波破碎算例 B5 的数值模拟结果与实验结果对比如图8 所示.

图8

图8   算例 U1(a) 和算例 B5(b) 截面 4 相对波高

Fig.8   Relative wave height at plane 4(Case U1 and Case B5)


对于图8(a) 的算例 U1,相对波高在截面 4 正中间位置有极大值约 1.5,沿两侧逐渐减小到 1 左右.数值模拟结果分布与实验结果相似,但在中间位置处的结果相较于实验结果偏小. 图8(b) 的算例 B5. 实验结果相对波高分布呈现中间低两边高的特征,中间极小值大约为 0.55,两侧逐渐向 1 逼近,由于不规则波存在一定偶然性,导致数值模拟结果在截面 4 中其中一个点出现较大误差,其余点的结果与实验结果吻合. 不规则波波浪破碎自由面如图9.

图9

图9   $t=36$s 的波面传播规律 (Case B5)

Fig.9   Perspective of computed wave surface at $t=36$s (Case B5)


图9 可以看出,波的破碎特征主要呈现陡峭的前倾状态,在浅滩的破碎区,波的聚焦特征已由于波浪破碎而而消标识码比规则波破碎的情况,不规则波破碎的形态显得更为复杂,更具随机化、多样化.另外,3 个区域的模拟对接效果良好. 图8图9 表明混合模式对不规则波的模拟整体可行.

3 结 论

本文提出一种研究波浪破碎的混合湍流模式方法,该方法将流体域分为层流造波区、大涡模拟区和消波区. 运用该方法分别模拟了规则与不规则波下破碎与不破碎的情况,并与相应的文献结果进行了对比分析,得出结论如下:

(1) 对于规则波不破碎的案例 M1,模拟结果与实验结果吻合,相应自由面出的流线图表明混合模式能够准确模拟出标识码过浅滩地形时发生的折射与绕射.

(2) 对于规则波破碎的案例 M3,模拟结果与实验结果大致吻合,混合模式能够捕捉到波浪破碎处的自由面形态,流线图表明标识码碎后存在投影于 $xy$ 平面上呈环形的流线,一定程度上能够解释波浪破碎后相对波高的分布原因.

(3) 对于不规则波的案例 U1 和 B5,模拟结果与实验结果整体相似,表明模式对于不规则波的模拟也有的可行性.

本文提出的混合湍流模式能准确模拟波浪的折射、反射、破碎等现象,为波浪破碎的研究提供一种方法选择. 本文提出的方法是能够较准确模拟波浪与空气运动的三维模拟模式,只要在入口处和部分顶部加上风边界条件,就可以用于标识码对于波浪破碎的影响,方法具有应用前景.

参考文献

孙宝楠, 连展, 杨永增 .

波浪破碎过程对油滴垂直混合影响的数值模拟及实际验证

船海工程, 2018,47(2):15-19

[本文引用: 1]

( Sun Baonan, Lian Zhan, Yang Yongzeng .

Numerical simulation and verification of the impact of wave breaking to the vertical mixing of oil drops

Ship Ocean Engineering, 2018,47(2):15-19 (in Chinese))

[本文引用: 1]

Kim DH, Son SY .

Role of shelf geometry and wave breaking in single N-type tsunami runup under geophysical-scale

Ocean Modelling, 2019,138:13-22

DOI      URL     [本文引用: 1]

Sriram R, Stagonas S .

Characterization of breaking wave impact on vertical wall with recurve

ISH Journal of Hydraulic Engineering, 2019,25(2):153-161

DOI      URL     [本文引用: 1]

Peregrine DH .

Long waves on a beach

Journal of Fluid Mechanics, 1967,27(4):815-827

DOI      URL     PMID      [本文引用: 1]

Family-centered prevention programs for couples with children are being increasingly disseminated, with the hope that improving couples' romantic relationships will lead to other benefits for families. To date, however, it is unclear whether these interventions do in fact yield these benefits. The current study addressed this gap by examining whether postintervention improvements in couples' relationship functioning following family-centered prevention predicted longer-term change in coparenting, and whether postintervention improvements in coparenting predicted longer-term change in relationship functioning. We used four waves of data collected over 2 years from 346 rural African American couples with an early adolescent child who participated in a randomized controlled trial of the Protecting Strong African American Families (ProSAAF) program, an intervention designed to promote strong couple, coparenting, and parent-child relationships in two-parent African American families. Results indicated that ProSAAF had significant short-term positive effects on both romantic relationship functioning and coparenting and that these effects did not differ in magnitude. Over time, however, only romantic relationship functioning postintervention was positively associated with long-term changes in coparenting; coparenting postintervention was not associated with long-term changes in relationship functioning and this association was significantly weaker than the other pathway. These findings support a key premise underlying relationship enhancement programs for parents, indicating that improving couples' romantic relationship functioning can have longer-term benefits for the coparenting relationship as well. Further research examining long-term parent, child, and family outcomes following family-centered prevention for couples and the mechanisms of change underlying these outcomes is needed.

郑茜, 张陈浩 .

近岸带珊瑚礁台波浪破碎变形流场结构的试验和数值模拟

中国港湾建设, 2019,39(2):36-40

[本文引用: 1]

( Zheng Qian, Zhang Chenhao .

Experimental and numerical simulation on flow field structure of wave breaking deformation at coral reefstation in the nearshore

China Harbour Engineering, 2019,39(2):36-40 (in Chinese))

[本文引用: 1]

张陈浩, 郑茜 .

基于OpenFOAM的浅滩波浪破碎区气液两相流耦合模拟

水运工程, 2018(11):38-42

[本文引用: 1]

( Zhang Chenhao, Zheng Qian .

Simulation of gas-liquid two-phase flow coupling in wave surf zone in shallow water based on OpenFOAM

Port & Waterway Engineering, 2018(11):38-42 (in Chinese))

[本文引用: 1]

许璐璐, 郄禄文 .

基于SPH法波浪破碎数值模拟分析

山西建筑, 2018,44(12):218-220

[本文引用: 1]

( Xu Lulu, Qie Luwen .

Numerical simulation analysis of broken waves based on SPH technique

Shanxi Architecture, 2018,44(12):218-220 (in Chinese))

[本文引用: 1]

张景新 .

基于非静压模型的波浪破碎模拟

水科学进展, 2017,28(3):438-444

[本文引用: 1]

( Zhang Jingxin .

Wave breaking simulation by non-hydrostatic numerical model

Advances in Water Science, 2017,28(3):438-444 (in Chinese))

[本文引用: 1]

Liu JJ, Patrizia P, Zhang QH .

Wave breaking analysis for the periodic rotation-two-component Camassa-Holm system

Nonlinear Analysis, 2019,187:214-228

DOI      URL     [本文引用: 1]

Luo L, Liu SX, Li JX , et al.

Numerical simulation of oblique and multidirectional wave propagation and breaking on steep slope based on FEM model of Boussinesq equations

Applied Mathematical Modelling, 2019,71:632-655

DOI      URL     [本文引用: 1]

诸裕良, 宗刘俊, 赵红军 .

复合坡度珊瑚礁地形上波浪破碎的试验研究

水科学进展, 2018,29(2):717-727

[本文引用: 1]

( Chu Yuliang, Zong Liujun, Zhao Hongjun , et al.

Experiment of wave breaking on coral reef topography with compound slope

Advances in Water Science, 2018,29(2):717-727 (in Chinese))

[本文引用: 1]

柳淑学, 刘宁, 李金宣 .

波浪在珊瑚礁地形上破碎特性试验研究

海洋工程, 2015,33(2):42-49

[本文引用: 1]

( Liu Shuxue, Liu Ning, Li Jinxuan , et al.

Study of wave breaking characteristic on coral reef topography

The Ocean Engineering, 2015,33(2):42-49 (in Chinese))

[本文引用: 1]

柳淑学, 胡书义, 李金宣 .

斜向和多向波浪的破碎指标研究

水动力学研究与进展(A辑), 2017,32(4):423-432

[本文引用: 1]

( Liu Shuxue, Hu Shuyi, Li Jinxuan , et al.

Research of diagonal and multidirectional wave breaking index

Chinese Journal of Hydrodynamics, 2017,32(4):423-432 (in Chinese))

[本文引用: 1]

Vincent CL, Briggs MJ .

Refraction - diffraction of irregular waves over a mound

Journal of Waterway Port Coastal and Ocean Engineering-ASCE, 1989,115(2):269-284

DOI      URL     [本文引用: 8]

詹杰民, 林东, 李毓湘 .

线性与非线性波的Chebyshev广义有限谱模拟

物理学报, 2007(4):3649-3654

[本文引用: 1]

( Zhan Jiemin, Lin Dong, Li Yok-Sheng .

Chebyshev generalized finite spectral method for linear and nonlinear waves

Acta Physica Sinica, 2007(4):3649-3654 (in Chinese))

[本文引用: 1]

Gouin M, Ducrozet G, Ferrant P .

Propagation of 3D nonlinear waves over an elliptical mound with a high-order spectral method

European Journal of Mechanics / B Fluids, 2017,63(Complete):9-24

[本文引用: 1]

Zhan JM, Li YS, Wai OWH .

Numerical modeling of multi-directional irregular waves incorporating 2-D numerical wave absorber and subgrid turbulence

Ocean Engineering, 2003,30(1):23-46

DOI      URL     [本文引用: 1]

Zou GL, Zhang QH .

Improvement of absorbing boundary conditions for non-hydrostatic wave-flow model SWASH

Applied Mechanics and Materials, 2013, 353-356:2676-2682

DOI      URL     [本文引用: 1]

Kang L, Guo X .

Depth-integrated, non-hydrostatic model using a new alternating direction implicit scheme

Journal of Hydraulic Research, 2013,51(4):368-379

DOI      URL     [本文引用: 1]

The hydrostatic pressure assumption has been widely applied in simulating rivers, lakes and reservoirs, but it has been found inappropriate in various cases where the vertical acceleration is significant. This paper presents a depth-integrated, non-hydrostatic model using a new alternating direction implicit scheme. Using the proposed scheme, each step is split into two half steps. In the first half step, the dynamic pressure and the x-direction velocity in the continuity equation and the momentum equations in the x-direction and z-direction are expressed implicitly, and the others explicitly; in the second half step, the dynamic pressure and the y-direction velocity in the continuity equation and the momentum equations in the y-direction and z-direction are discretized implicitly, and the others explicitly. The Thomas algorithm is applied to solve the tri-diagonal linear system at each half step. The model is developed and validated with several analytical solutions and laboratory experiments. The results show that the model can provide comparable results at very low computational cost.

Zhan JM, Dong Z, Jiang W , et al.

Numerical simulation of wave transformation and runup incorporating porous media wave absorber and turbulence models

Ocean Engineering, 2010,7(14-15):1261-1272

[本文引用: 2]

Griffiths LS, Porter R .

Focusing of surface waves by variable bathymetry

Applied Ocean Research, 2012,34(none):150-163

DOI      URL     [本文引用: 1]

Choi J, Yoon SB .

A three-dimensional numerical simulation of wave and current evolving over a submerged shoal

Coastal Engineering Journal, 2011,53(2):63-85

DOI      URL     [本文引用: 1]

Yin J, Sun JW, Wang XG , et al.

A hybrid finite-volume and finite difference scheme for depth- integrated non-hydrostatic model

China Ocean Engineering, 2017,31(1):261-271

DOI      URL     [本文引用: 1]

Smit PB, Janssen TT .

The evolution of inhomogeneous wave statistics through a variable medium

Journal of Physical Oceanography, 2013,43(8):1741-1758

DOI      URL     [本文引用: 1]

Hu KL, Ding PX .

Numerical study of wave diffraction effect introduced in the SWAN model

China Ocean Engineering, 2007(3):128-139

[本文引用: 1]

Chen XW, Zheng JH, Zhang C , et al.

Evaluation of diffraction predictability in two phase averaged wave models

China Ocean Engineering, 2010,24(2):235-244

[本文引用: 1]

Hajime M .

Multi-directional random wave transformation model based on energy balance equation

Coastal Engineering Journal, 2001,43(4):317-337

DOI      URL     [本文引用: 1]

Janssen TT, Herbers THC .

Nonlinear wave statistics in a focal zone

Journal of Physical Oceanography, 2009,39(8):1948-1964

DOI      URL     [本文引用: 1]

金红, 邹志利 .

具有精确色散性的非线性波浪数学模型

力学学报, 2010,42(1):23-34

URL     [本文引用: 3]

以完全非线性的自由表面边界条件为基础,以波面升高\eta和自由表面速度势\phi _\eta为待求变量,建立了新的波浪方程.方程在色散性上是完全精确的,非线性近似至三阶.与缓坡方程相比较,两者都具有精确的色散性,但该方程属于非线性模型,可模拟波浪的非线性效应,且适用于不规则波.方程的特点是属于微分-积分方程,对如何处理方程中积分项进行了讨论,并数值模拟了不同周期的线性波和二阶Stokes波,也模拟了波群的非线性演化,以对模型进行验证.

( Jin Hong, Zou Zhili .

Non-linear wave model with accurate dispersion

Chinese Journal of Theoretical and Applied Mechanics, 2010,42(1):23-34 (in Chinese))

URL     [本文引用: 3]

以完全非线性的自由表面边界条件为基础,以波面升高\eta和自由表面速度势\phi _\eta为待求变量,建立了新的波浪方程.方程在色散性上是完全精确的,非线性近似至三阶.与缓坡方程相比较,两者都具有精确的色散性,但该方程属于非线性模型,可模拟波浪的非线性效应,且适用于不规则波.方程的特点是属于微分-积分方程,对如何处理方程中积分项进行了讨论,并数值模拟了不同周期的线性波和二阶Stokes波,也模拟了波群的非线性演化,以对模型进行验证.

Bouws E, Günther H, Rosenthal W , et al.

Similarity of the wind wave spectrum in finite depth water: 1. Spectral form

Journal of Geophysical Research: Oceans, 1985,90:975-986

[本文引用: 2]

金红 .

具有精确色散性的线性和非线性波浪模型

[博士论文]. 大连:大连理工大学, 2008

[本文引用: 1]

( Jin Hong .

Linear and non-linear wave models with accurate dispersion

[PhD Thesis]. Dalian: Dalian University of Technology, 2008 (in Chinese))

[本文引用: 1]

/