An ALE approach to reduce spurious numerical mixing through variational minimizers: application to internal waves

TL;DR

This paper introduces an ALE-based variational method to minimize spurious numerical mixing in ocean models, specifically applied to internal waves, improving physical accuracy in challenging nonlinear wave simulations.

Contribution

The paper presents a novel variational approach to define vertical grid motion in ALE ocean models, reducing numerical mixing and enhancing simulation fidelity for internal waves.

Findings

Effective reduction of spurious mixing in internal wave simulations.

Method maintains physical relevance in nonlinear wave breaking and overturning.

Applicable to any ocean model with ALE vertical coordinates.

Abstract

Spurious numerical mixing is a frequent phenomenon in ocean models. In the present paper, we present an efficient and robust methodology that defines the vertical grid motion so that this mixing is reduced. This motion is defined through the solution of an optimization problem that -- using the ideas of the calculus of variations -- results in an elliptic equation. This framework is generally applicable to any ocean model that uses an ALE vertical coordinate and can be tuned to fit the modeler's specific needs based on the guidelines presented herein. The method is applied to the nonhydrostatic solver presented by the authors in [Alexandris-Galanopoulos et al., 2024] and its applicability in fully nonlinear internal waves is investigated for the demanding test cases of wave breaking and overturning. These numerical benchmarks show the ability of the method to reduce spurious mixing, while attaining the physical relevancy of the results.