HAL : hal-00665826, version 2

DOI : 10.1016/j.ocemod.2012.04.007

Fiche détaillée

Récupérer au format BibTeX - EndNote - TEI - RefWorks -

Ocean Modelling 52-53 (2012) 9-35

Versions disponibles :

v1 (02-02-2012)

v2 (26-04-2012)

On the Stability and Accuracy of the Harmonic and Biharmonic Isoneutral Mixing Operators in Ocean Models

Florian Lemarié (  Auteur correspondant ) 1, Laurent Debreu 2, Alexander Shchepetkin 1, James C. Mcwilliams 3

(08/2012)

Ocean models usually rely on a tracer mixing operator which diffuses along isoneutral directions. This requirement is imposed by the highly adiabatic nature of the oceanic interior, and a numerical simulation needs to respect these small levels of dianeutral mixing to maintain physically realistic results. For non-isopycnic models this is however non-trivial due to the non-alignment of the vertical coordinate isosurfaces with local isoneutral directions, rotated mixing operators must therefore be used. This paper considers the numerical solution of initial boundary value problems for the harmonic (Laplacian) and biharmonic rotated diffusion operators. We provide stability criteria associated with the conventional space-time discretizations of the isoneutral Laplacian operator currently in use in general circulation models. Furthermore, we propose and study possible alternatives to those schemes. A new way to handle the temporal discretization of the rotated biharmonic operator is also introduced. This scheme requires only the resolution of a simple one-dimensional tridiagonal system in the vertical direction to provide the same stability limit of the non-rotated operator. The performance of the various schemes in terms of stability and accuracy is illustrated by idealized numerical experiments of the diffusion of a passive tracer along isoneutral directions.

1 :	Department of Atmospheric Sciences [Los Angeles]
	Institute of Geophysics and Planetary Physics, Los Angeles – University of California, Los Angeles
2 :	MOISE (INRIA Grenoble Rhône-Alpes / LJK Laboratoire Jean Kuntzmann)
	CNRS : UMR5224 – INRIA – Laboratoire Jean Kuntzmann – Université Joseph Fourier - Grenoble I – Institut Polytechnique de Grenoble - Grenoble Institute of Technology
3 :	Institute of Geophysics and Planetary Physics [Los Angeles] (IGPP)
	University of California at Los Angeles

Domaine

:

Mathématiques/Equations aux dérivées partielles Mathématiques/Analyse numérique Planète et Univers/Océan, Atmosphère

Mots Clés : Isoneutral Mixing – Diffusion Equation – Time Schemes – Von-Neumann stability analysis – Initial-Boundary Value Problems

hal-00665826, version 2

http://hal.archives-ouvertes.fr/hal-00665826

oai:hal.archives-ouvertes.fr:hal-00665826

Contributeur : Florian Lemarié <>

Soumis le : Jeudi 26 Avril 2012, 16:43:43

Dernière modification le : Lundi 25 Février 2013, 19:56:28