EGU22-10204
https://doi.org/10.5194/egusphere-egu22-10204
EGU General Assembly 2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

Finding better numerical solutions for circulation along piecewise-constant coastlines in ocean models

Antoine-Alexis Nasser1,2, Gurvan Madec1,2, Laurent Debreu1, and Casimir De Lavergne2
Antoine-Alexis Nasser et al.
  • 1INRIA Grenoble Rhônes-Alpes
  • 2LOCEAN IPSL

Using a 1/4° shallow water model, Adcroft & Marshall (1998) show that the solution obtained on a rotated mesh depends on the angle between the numerical grid and the physical coastline. They conclude that piecewise constant coastline exert a spurious form stress on the model boundary current.

Here, we reproduce these results and test the combination of free-slip and no-slip boundary conditions with the stress tensors’ vorticity-divergence and symmetric forms. We show that, for all angles of rotation, the equilibrium solutions actually remain similar, provided that the system has numerically converged (i.e. higher resolution with physics unchanged). It appears that having at least 4 grid-points by internal radius of deformation is necessary to accurately represent flows along coastlines. Surprisingly, the symmetric tensor combined with the free-slip condition leads to dynamics that are insensitive to resolution and akin to no-slip. We argue that the way the free-slip condition is implemented with the symmetric tensor is not suitable. We propose an alternative implementation and discuss its advantages and limitations.

How to cite: Nasser, A.-A., Madec, G., Debreu, L., and De Lavergne, C.: Finding better numerical solutions for circulation along piecewise-constant coastlines in ocean models, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10204, https://doi.org/10.5194/egusphere-egu22-10204, 2022.

Comments on the display material

to access the discussion