Abstract:The design and implementation of boundary conditions for the robust generation and simulation of periodic finite-amplitude internal waves is examined in a quasi two-layer continuous stratification using a spectral-element-method-based incompressible flow solver. The commonly used Eulerian approach develops spurious, and potentially catastrophic small-scale numerical features near the wave-generating boundary in a non-linear stratification when the parameter A/(δc) is sufficiently larger than unity; A and δ are measures of the maximum wave-induced vertical velocity and pycnocline thickness, respectively, and c is the linear wave propagation speed. To this end, an Euler–Lagrange approach is developed and implemented to generate robust high-amplitude periodic deep-water internal waves. Central to this approach is to take into account the wave-induced (isopycnal) displacement of the pycnocline in both the vertical and (effectively) upstream directions. With amplitudes not restricted by the limits of linear theory, the Euler–Lagrange-generated waves maintain their structural integrity as they propagate away from the source. The advantages of the high-accuracy numerical method, whose minimal numerical dissipation cannot damp the above near-source spurious numerical features of the purely Eulerian case, can still be preserved and leveraged further along the wave propagation path through the robust reproduction of the non-linear adjustments of the waveform. The near- and far-source robustness of the optimized Euler–Lagrange approach is demonstrated for finite-amplitude waves in a sharp quasi two-layer continuous stratification representative of seasonally stratified lakes. The findings of this study provide an enabling framework for two-dimensional simulations of internal swash zones driven by well-developed non-linear internal waves and, ultimately, the accompanying turbulence-resolving three-dimensional simulations.

Numerical solutions of the parabolic wave equation: An ordinary‐differential‐equation approach

Three-dimensional boundary fitted parabolic-equation model of underwater sound propagation

Dynamically orthogonal narrow-angle parabolic equations for stochastic underwater sound propagation. Part I: Theory and schemes

Direct numerical simulation of acoustic wave propagation in ocean waveguides using a parallel finite volume solver

A Stable Chebyshev-Padé Rational Approximation of Parabolic Equation Models for Underwater Sound Field Computation

A Robust Numerical Method Based on a Deep-Learning Operator for Solving the 2D Acoustic Wave Equation

A 3D Parabolic Equation Model with Higher-Order Operator Approximation

A robust numerical method based on deep learning operator for solving the 2D acoustic wave equation

A generalization of the split-step Padé method to the case of coupled acoustic modes equation in a 3D waveguide

A numerical method for wave-structure interactions in the Boussinesq regime

High Accuracy Wave Simulation – Revised Derivation, Numerical Analysis and Testing of a Nearly Analytic Integration Discrete Method for Solving Acoustic Wave Equation

Dynamically orthogonal narrow-angle parabolic equations for stochastic underwater sound propagation. Part II: Applications

A numerical method for computing time-periodic solutions in dissipative wave systems

Numerical solution of highly oscillatory ordinary differential equations

Numerical analysis of the diffusive-viscous wave equation

A simple algorithm for numerical solution of nonlinear parabolic partial differential equations

A Nearly Analytic Discrete Method for Solving the Acoustic-Wave Equations in the Frequency Domain

A New Ocean Acoustic Model in Computational Oceanographic Physics

A Finite Difference Method for Effective Treatment of Mild-Slope Wave Equation Subject to Non-Reflecting Boundary Conditions

A robust numerical method for the generation and propagation of periodic finite-amplitude internal waves in natural waters using high-accuracy simulations

Economic numerical method of solving coefficient inverse problem for 3D wave equation