Comparison of different numerical methods used to handle the open boundary of a regional ocean circulation model of the Bay of Biscay

Sylvain Cailleau (1) , Veronika Fedorenko (2) , Bernard Barnier (1) , Eric Blayo (2) , Laurent Debreu (2)

1 LEGI - Laboratoire des Écoulements Géophysiques et Industriels [Grenoble]
2 MOISE - Modelling, Observations, Identification for Environmental Sciences

Several methods for specifying boundary conditions at the limits of a regional model are compared. The methods investigated are those using clamped and radiation boundary conditions, one-way and two-way nesting in a model for a more extensive area, and "full" coupling based on domain decomposition techniques. These methods are compared in the realistic framework of interactions between a 1/15° model of the Bay of Biscay and a 1/3° model of the North Atlantic, over a 3-year simulation (1996-1998). The clamped and radiation boundary conditions systematically lead to energy accumulation and problematic recirculations along the boundary, and can disturb the internal dynamics of the regional domain. The one-way or two-way-nesting and the full-coupling methods result in far more satisfactory behaviour. For long periods of integration, the two-way mode improves both the fine and coarse-grid solutions. The full coupling method provides the most regular solution at the boundary, and also opens interesting new perspectives since it should enable the coupling of models with different physics. However it requires much more computation time.

Mots clés: Regional modelling; Open boundary conditions; Nesting; Coupling; Bay of Biscay

Domaines: Analyse numérique [math.NA]; Océan, Atmosphère

Sylvain Cailleau, Veronika Fedorenko, Bernard Barnier, Eric Blayo, Laurent Debreu. Comparison of different numerical methods used to handle the open boundary of a regional ocean circulation model of the Bay of Biscay. Ocean Modelling, 2008, 25 (1-2), pp.1-16. ⟨10.1016/j.ocemod.2008.05.009⟩. ⟨hal-00387446⟩