[论文解读] A penalization method for calculating the flow beneath travelling water waves of large amplitude

Adrian Constantin, Konstantinos Kalimeris+1 more|arXiv (Cornell University)|Aug 8, 2014

Ocean Waves and Remote Sensing参考文献 13被引用 3

一句话总结

本文提出一种惩罚方法，用于数值模拟具有恒定涡度的大振幅行进水波，将控制方程重新表述为约束优化问题。该方法准确捕捉了波面下方的流动特性——包括速度、压力和流线形态，涵盖正涡度或负涡度的旋转流以及无旋情况。

ABSTRACT

A penalization method for a suitable reformulation of the governing equations as a constrained optimization problem provides accurate numerical simulations for large-amplitude travelling water waves in irrotational flows and in flows with constant vorticity.

研究动机与目标

• 开发一种稳健的数值方法，用于模拟具有恒定涡度的流动中的大振幅行进水波，涵盖无旋和旋转情况。
• 解决在陡峭、非线性水波中准确计算波面下方流动动力学的挑战，传统方法在此类情况下可能失效。
• 确保数值解满足底部无通量和自由表面运动学条件的物理约束。
• 提供一种稳定且收敛的迭代算法，能够从更广泛的潜在解中筛选出真实的波解。

提出的方法

• 采用流函数公式将稳态水波问题重新表述为约束优化问题，以减少未知数数量。
• 引入一种惩罚技术，通过变分公式强制执行控制偏微分方程和边界条件。
• 使用有限差分法对流函数 $ h(q,p) $ 进行离散化，以表示流体域并迭代求解该系统。
• 实施一种迭代算法，交替更新高度函数 $ h^{(k)} $ 并最小化偏微分方程与边界条件的残差。
• 基于残差范数和接近停滞点（即 $ 1/h_p o 0 $）的条件设定停止准则，并在 $ h_p^{(k+1)} o 0 $ 时进行自适应重初始化。
• 采用从具有恒定涡度的大振幅非层流波的解析公式导出的初始近似，其中无旋情况为特例（$ u = 0 $）。

实验结果

研究问题

• RQ1 尽管存在非线性且缺乏小振幅假设，如何能准确地数值模拟具有恒定涡度的大振幅行进水波？
• RQ2 在无停滞点的旋转流条件下，惩罚方法能否有效强制执行复杂的边界条件和偏微分方程约束？
• RQ3 在无旋和旋转流中，速度分量和压力分布在大振幅波下方的行为如何？
• RQ4 涡度的存在——无论是正还是负——如何影响自由表面的形状以及波下方的流线结构？

主要发现

• 惩罚方法成功计算出具有恒定涡度的无旋和旋转流中大振幅行进水波的精确数值解。
• 在无旋情况（$ u = 0 $）下，该方法重现了已知的波剖面和与现有理论及数值研究一致的速度分布模式。
• 在正涡度情况（$ u = 2.95 $）下，波剖面变得不对称，前坡更陡峭，波峰更宽广，且近底部的水平速度 $ c - u $ 减小。
• 在负涡度情况（$ u = -1 $）下，波剖面表现出更强的不对称性，前坡更陡峭，压力分布显示出在波峰下方更显著的最大值。
• 竖直速度 $ v $ 在自由表面附近最强，并向底部减小，所有情况下最大值均出现在波峰附近。
• 压力偏离大气压的偏差从自由表面到底部单调增加，在波峰下方略近处达到最大值，与静水假设一致。