Etudes et évaluation de processus océaniques par des hiérarchies ...

Etudes et évaluation de processus océaniques par des hiérarchies ... Etudes et évaluation de processus océaniques par des hiérarchies ...

tel.archives.ouvertes.fr
from tel.archives.ouvertes.fr More from this publisher
27.12.2013 Views

134 CHAPITRE 4. ETUDES DE PROCESSUS OCÉANOGRAPHIQUES 554 Ocean Dynamics (2009) 59:551–563 tel-00545911, version 1 - 13 Dec 2010 I used the maximal height H = 200 m for the definition of the vein Richardson number and a vein Froude number. A vertical Reynolds number is given by Re = 2(Ek Ri tan α) −1 and depends on the Ekman number, the Richardson number and the slope. This Reynolds number is the geostrophic velocity multiplied by the layer thickness and divided by the vertical viscosity. The influence of a second scalar, salinity and a nonlinear equation of state, leading to, e.g. the thermobaric effect (see Emms 1998), is not considered here. 2.2 The mathematical model The mathematical model for the gravity current dynamics are the Navier–Stokes equations in a rotating frame with a buoyant scalar (temperature). I neglect variations in the stream wise (y) direction of all the variables but include the temperature difference and all the three components of the velocity vector, this type of model is usually revered to as 2.5 dimensional. This leads to a four-dimensional state vector depending on two space and the time variable, (T(x, z, t), u(x, z, t), v(x, z, t), w(x, z, t)). The domain is a rectangular box that spans 51.2 km in the x direction and is 492 m deep (z direction). On the bottom, there is a no-slip and on the top a free-slip boundary condition. The horizontal boundary conditions are periodic. The initial condition is a temperature anomaly which has a parabolic shape, which is 200 m high and 20 km large at the bottom, as described in the previous subsection. The magnitude of the temperature anomaly is varied in the experiments. The initial velocities in the gravity current are geostrophically adjusted, the fluid outside the gravity current is initially at rest. The buoyancy force is represented by an acceleration of strength g ′ cos(α) in the z direction and g ′ sin(α) in the negative x direction to represent the inclination of angle α = 1 ◦ . This geometry represents a rectangular box that is tilted by an angle of one degree. Such implementation of a sloping bottom simplifies the numerical implementation and allows for using powerful numerical methods (see next subsection). 2.3 Numerical implementation of mathematical model The numerical model used is HAROMOD (Wirth 2005). HAROMOD is a pseudo-spectral code, based on Fourier series in all the spatial dimensions, that solves the Navier–Stokes equations subject to the Boussinesq approximation, a no-slip boundary condition on the floor and a free-slip boundary condition at the rigid surface. The time stepping is a third-order, low-storage, Runge–Kutta scheme. A major difficulty in the numerical solution is due to the large anisotropy in the dynamics and the domain, which is roughly 100 times larger than deep. There are 896 points in the vertical direction. For a density anomaly larger than 0.75 K, the horizontal resolution had to be increased from 512 to 2,048 points (see Table 1) to avoid a pile up of small-scale energy caused by an insufficient viscous dissipation range, leading to a thermalised dynamics at small scales as explained by Frisch et al. (2008). The horizontal viscosity is ν h = 5 m 2 s −1 and the horizontal diffusivity is κ h = 1 m 2 s −1 . The vertical viscosity is ν v = 10 −3 m 2 s −1 and the vertical diffusivity is κ v = 10 −4 m 2 s −1 . The anisotropy in the turbulent mixing coefficients reflects the strong anisotropy of the numerical grid. I checked that the results presented here show only a slight dependence on ν h , κ h and κ v by doubling these constants in a control run. This is no surprise as the corresponding diffusion and friction times are larger than the integration time of the experiments. There is a strong dependence on ν v , as it determines the thickness of the Ekman layer and the Ekman transport, which governs the dynamics of the gravity current as will be shown in Section 3. The vertical extension of the Ekman layer is a few metres, while the horizontal extension of the gravity current is up to 50 km. In a fully turbulent gravity current, the turbulent structures within the well mixed gravity current will be isotropic and will therefore measure only a few metres in size. To simulate a fully turbulent gravity current, 10 5 grid points would be necessary in the horizontal direction to obtain an isotropic grid. This is far beyond our actual computer resources. 2.4 Experiments performed The density anomaly was varied in the experiments. The integration was stopped when the downslope side Table 1 Physical and numerical parameters varied in the numerical experiments Exp. T (K) v G (10 −2 )ms −1 Ri N x Integration time (h) G00 0.25 8.30 14.2 512 360 G01 0.5 16.6 7.12 512 192 G03 0.75 24.9 4.75 512 132 G12 1.0 33.2 3.56 2,096 96 G14 1.1 36.5 3.24 2,096 86 G15 1.25 41.5 2.85 2,096 76 G17 1.5 49.8 2.37 2,096 66

4.7. ON THE BASIC STRUCTURE OF OCEANIC GRAVITY CURRENTS 135 Ocean Dynamics (2009) 59:551–563 555 tel-00545911, version 1 - 13 Dec 2010 of the friction layer attained the boundary of the domain. All the experiments discussed in this publication are listed in Table 1, which includes information concerning the physics of the experiments: temperature anomaly T, mean geostrophic velocity v G and Richardson number Ri, as well as numerical values: horizontal resolution and time of integration. 3 Results I start with a qualitative description of the dynamics of the gravity current before quantitatively determining the key parameters. 3.1 Qualitative description The gravity current is initially in a geostrophic state as expressed by Eq. 1. The velocity near the boundary is rapidly reduced due to friction. This rapid decrease of the velocity leads to inertial oscillations throughout the gravity current. Near the ocean floor, in the bottom Ekman layer, frictional forces reduce the long-slope velocity so the flow is no longer geostrophically adjusted, as the Coriolis force no longer balances the hydrostatic pressure gradient and the water in the Ekman layer Fig. 1 Structure of the gravity current after 60 h in exp. G00 (top) and G15 (bottom), see Table 1 for details. The inclination of the floor is exaggerated in the figure, for pedagogical reasons only; the real angle is only 1 ◦ . Please note that the vertical extension is given in metres and the horizontal extension in kilometres. The temperature anomaly, with respect to the ambient water, is shown (given in Kelvin). Isolines of the y component of the velocity (positive into the plane) are given every 0.05 ms −1 in black for positive values and red for negative values, zero line is omitted flows down the pressure gradient. At the downslope border of the vein, a friction layer forms, which consists of gravity current water moving downslope (see Fig. 2). This feature is observed in laboratory experiments by Smith (1977) and Lane-Serf and Baines (1998) (and references therein) and measured by Wirth and Sommeria (unpublished manuscript). In Fig. 1, one sees a typical result of our numerical integrations for two experiments with different density and temperature differences. In the interior of the gravity current, the isolines of the y component are vertical due to geostrophy, which says that this velocity component is determined by hydrostatic pressure gradient, that is, the angle of the surface of the gravity current to the horizontal. In the upslope part of the gravity current, the y component of the velocity is negative due to the negative angle of the surface of the gravity current to the horizontal. Figure 1 clearly shows that the gravity current consists of two parts, a “vein” which is the thick part of the gravity current (thicker than about 20 m) and a “friction layer”, the thin layer of dense fluid that extents downstream of the vein. The vein of the gravity current water detrains through the Ekman layer into the friction layer, which is responsible of the major part of the downslope transport of gravity current water. Please

4.7. ON THE BASIC STRUCTURE OF OCEANIC GRAVITY CURRENTS 135<br />

Ocean Dynamics (2009) 59:551–563 555<br />

tel-00545911, version 1 - 13 Dec 2010<br />

of the friction layer attained the boundary of the domain.<br />

All the experiments discussed in this publication<br />

are listed in Table 1, which inclu<strong>de</strong>s information<br />

concerning the physics of the experiments: temperature<br />

anomaly T, mean geostrophic velocity v G and<br />

Richardson number Ri, as well as numerical values:<br />

horizontal resolution and time of integration.<br />

3 Results<br />

I start with a qualitative <strong>de</strong>scription of the dynamics<br />

of the gravity current before quantitatively d<strong>et</strong>ermining<br />

the key <strong>par</strong>am<strong>et</strong>ers.<br />

3.1 Qualitative <strong>de</strong>scription<br />

The gravity current is initially in a geostrophic state as<br />

expressed by Eq. 1. The velocity near the boundary<br />

is rapidly reduced due to friction. This rapid <strong>de</strong>crease<br />

of the velocity leads to inertial oscillations throughout<br />

the gravity current. Near the ocean floor, in the bottom<br />

Ekman layer, frictional forces reduce the long-slope velocity<br />

so the flow is no longer geostrophically adjusted,<br />

as the Coriolis force no longer balances the hydrostatic<br />

pressure gradient and the water in the Ekman layer<br />

Fig. 1 Structure of the<br />

gravity current after 60 h in<br />

exp. G00 (top) and G15<br />

(bottom), see Table 1 for<br />

d<strong>et</strong>ails. The inclination of the<br />

floor is exaggerated in the<br />

figure, for pedagogical<br />

reasons only; the real angle is<br />

only 1 ◦ . Please note that the<br />

vertical extension is given in<br />

m<strong>et</strong>res and the horizontal<br />

extension in kilom<strong>et</strong>res. The<br />

temperature anomaly, with<br />

respect to the ambient water,<br />

is shown (given in Kelvin).<br />

Isolines of the y component<br />

of the velocity (positive into<br />

the plane) are given every<br />

0.05 ms −1 in black for<br />

positive values and red for<br />

negative values, zero line is<br />

omitted<br />

flows down the pressure gradient. At the downslope<br />

bor<strong>de</strong>r of the vein, a friction layer forms, which consists<br />

of gravity current water moving downslope (see Fig. 2).<br />

This feature is observed in laboratory experiments by<br />

Smith (1977) and Lane-Serf and Baines (1998) (and references<br />

therein) and measured by Wirth and Sommeria<br />

(unpublished manuscript).<br />

In Fig. 1, one sees a typical result of our numerical<br />

integrations for two experiments with different <strong>de</strong>nsity<br />

and temperature differences. In the interior of the gravity<br />

current, the isolines of the y component are vertical<br />

due to geostrophy, which says that this velocity component<br />

is d<strong>et</strong>ermined by hydrostatic pressure gradient,<br />

that is, the angle of the surface of the gravity current<br />

to the horizontal. In the upslope <strong>par</strong>t of the gravity<br />

current, the y component of the velocity is negative<br />

due to the negative angle of the surface of the gravity<br />

current to the horizontal.<br />

Figure 1 clearly shows that the gravity current consists<br />

of two <strong>par</strong>ts, a “vein” which is the thick <strong>par</strong>t of<br />

the gravity current (thicker than about 20 m) and a<br />

“friction layer”, the thin layer of <strong>de</strong>nse fluid that extents<br />

downstream of the vein. The vein of the gravity current<br />

water d<strong>et</strong>rains through the Ekman layer into the friction<br />

layer, which is responsible of the major <strong>par</strong>t of the<br />

downslope transport of gravity current water. Please

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!