Numerical Simulations of the Inviscid Primitive Equations in a Limited Domain

Antoine Rousseau (1) , Roger M. Temam (2, 3) , Joe Tribbia (4)

1 IDOPT - System identification and optimization in physics and environment
2 LM-Orsay - Laboratoire de Mathématiques d'Orsay
3 ISC - Institute for Scientific Computing and Applied Mathematics
4 CGD - Climate and Global Dynamics Division [Boulder]

This work is dedicated to the numerical computations of the primitive equations (PEs) of the ocean without viscosity with nonlocal (mode by mode) boundary conditions. We consider the 2D nonlinear PEs, and firstly compute the solutions in a "large" rectangular domain D with periodic boundary conditions in the horizontal direction. Then we consider a subdomain D', in which we compute a second numerical solution with transparent boundary conditions. Two objectives are achieved. On the one hand the absence of blow-up in these computations indicates that the PEs without viscosity are well-posed when supplemented with the boundary conditions. On the other hand they show a very good coincidence on the subdomain D' of the two solutions, thus showing also the computational relevance of these new boundary conditions. We end this study with some numerical simulations of the linearized primitive equations, which correspond to the theoretical results established by the authors, and evidence the transparent properties of the boundary conditions.

Domaines: Equations aux dérivées partielles [math.AP]; Analyse numérique [math.NA]; Dynamique des Fluides [physics.flu-dyn]; Géophysique [physics.geo-ph]; Géophysique [physics.geo-ph]; Bibliothèque électronique [cs.DL]

Antoine Rousseau, Roger M. Temam, Joe Tribbia. Numerical Simulations of the Inviscid Primitive Equations in a Limited Domain. Mathematical Methods for Hydrodynamics, Université Lille 1, Jun 2005, Lille, France. pp.163-181, ⟨10.1007/978-3-7643-7742-7_10⟩. ⟨inria-00172561⟩