[论文解读] Three-Dimensional Simulations of Deep-Water Breaking Waves

Kyle A. Brucker, Thomas T. O’Shea+2 more|arXiv (Cornell University)|Oct 7, 2014

Ocean Waves and Remote Sensing参考文献 7被引用 16

一句话总结

本论文提出了一种基于体积分数法的三维数值模拟框架，采用一种新颖的法向应力强迫技术，以生成完全非线性的深水推进波。主要发现包括：识别出不同的破碎波类型（从弱到极强的冲浪型破碎），量化了空气卷入与能量耗散（在强事件中超过50%的峰值能量被耗散），并证明Favre平均相比雷诺平均能更准确地跟踪自由表面。

ABSTRACT

The formulation of a canonical deep-water breaking wave problem is introduced, and the results of a set of three-dimensional numerical simulations for deep-water breaking waves are presented. In this paper fully nonlinear progressive waves are generated by applying a normal stress to the free surface. Precise control of the forcing allows for a systematic study of four types of deep-water breaking waves, characterized herein as weak plunging, plunging, strong plunging, and very strong plunging.

研究动机与目标

• 开发一个用于在三维空间中研究深水破碎波的基准模拟框架。
• 实现对各种破碎波类型（从弱冲浪型到极强冲浪型）的系统性模拟，并精确控制激励条件。
• 量化破碎波中空气卷入与能量耗散机制。
• 比较雷诺平均与Favre平均在自由表面流统计分析中的性能。
• 评估在下一代硬件上使用NFA代码对包含50亿至100亿未知数的大规模模拟的可行性。

提出的方法

• 在移动计算框架中，通过在自由表面施加法向应力强迫技术，生成完全非线性的推进波。
• 采用体积分数法（VOF）追踪气水界面，其中体积分数α表示每个单元中流体的占据比例。
• 求解不可压Navier-Stokes方程，采用包含压力和表面张力效应的动量方程保守形式。
• 数值方法采用有限体积法，基于非结构化非交错网格，并结合压力-速度耦合算法。
• 采用一种改进的MPI通信策略，使用一维子通信器替代全局集体操作，从而实现对数千个处理器的高效扩展。
• 对速度场和体积分数场应用Favre平均与雷诺平均，以比较其在捕捉自由表面动态中的统计性能。

实验结果

研究问题

• RQ1 法向应力强迫技术如何实现在深水环境中对完全非线性推进波的可控生成？
• RQ2 不同破碎波类型（从弱冲浪型到极强冲浪型）所对应的三维流场结构与涡量分布有何特征？
• RQ3 在破碎过程中能量耗散的程度如何？不同破碎强度下的耗散程度有何差异？
• RQ4 在不同破碎类型中，自由表面下方的空气脱气速率有何比较？
• RQ5 在模拟中，哪种统计平均方法——雷诺平均还是Favre平均——能更准确地反映真实自由表面位置？

主要发现

• 在强和极强冲浪型破碎事件中，超过50%的峰值能量被耗散，表明强烈破碎过程伴随显著的能量损失。
• 在主要冲浪型破碎事件中，空气卷入后不久即形成强烈的顺流向涡量，如三维涡量等值面所示。
• 在初始破碎事件之后，大量空气仍长期滞留在自由表面下方，尤其在强和极强冲浪型破碎情况下更为显著。
• 所有情况下空气脱气速率均为线性且相同，表明其物理机制与破碎强度无关。
• 对速度和体积分数场进行Favre平均，相比雷诺平均，能更准确地匹配真实自由表面位置，尤其在捕捉界面动态方面表现更优。
• 在下一代硬件上，使用NFA代码在最优条件下，对包含50亿至100亿未知数的模拟已具备可行性，计算速度可达每小时约100次迭代。