Simulation of acoustic and flexural-gravity waves in ice-covered oceans
Simulation of acoustic and flexural-gravity waves in ice-covered oceans
复制标题
DOI：
10.1016/j.jcp.2018.06.060
复制
发表时间：
2018-11
期刊：
影响因子：
--
通讯作者：
K. Mattsson;E. Dunham;J. Werpers
中科院分区：
文献类型：
--
作者：
K. Mattsson;E. Dunham;J. Werpers
我们介绍了一种可证明稳定的、高精度的有限差分方法，用于模拟部分被薄弹性层覆盖的可压缩、无粘流体中的声波和弯曲重力波。当研究海浪与漂浮的冰架、海冰和漂浮结构物的相互作用时，就会产生这样的海浪。特别强调的是流体-冰耦合的良好界面处理。为了保证数值的稳定性和效率，基于部分求和(SBP)框架的有限差分近似与惩罚技术(同时逼近项，SAT)相结合来施加边界条件和界面条件。得到的SBP-SAT近似是用无条件稳定的有限差分方法进行时间积分的。二维数值模拟证实了预测的效率和稳定性行为。该方法可以用目前的形式来研究海浪和海啸在冰架上的传播，并在耦合到冰和水下面的弹性半空间时，研究与长周期地震面波有关的冰运动。
We introduce a provably stable, high-order-accurate finite difference method for simulation of acoustic and flexural gravity waves in compressible, inviscid fluids partially covered by a thin elastic layer. Such waves arise when studying ocean wave interactions with floating ice shelves, sea ice, and floating structures. Particular emphasis is on a well-posed interface treatment of the fluid-ice coupling. To ensure numerical stability and efficiency, finite difference approximations based on the summation-by-parts (SBP) framework are combined with a penalty technique (simultaneous approximation term, SAT) to impose the boundary and interface conditions. The resulting SBP-SAT approximations are time integrated with an unconditionally stable finite difference method. Numerical simulations in 2D corroborate the predicted efficiency and stability behaviors. The method can be used in its current form to study transmission of ocean waves and tsunamis through ice shelves, and upon coupling to an elastic half-space beneath the ice and water, to study ice motions associated with long-period seismic surface waves.