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 ...
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
- Page 89 and 90: 4.3. A NON-HYDROSTATIC FLAT-BOTTOM
- Page 91 and 92: 4.3. A NON-HYDROSTATIC FLAT-BOTTOM
- Page 93 and 94: 4.3. A NON-HYDROSTATIC FLAT-BOTTOM
- Page 95 and 96: 4.3. A NON-HYDROSTATIC FLAT-BOTTOM
- Page 97 and 98: 4.4. TILTED CONVECTIVE PLUMES IN NU
- Page 99 and 100: 4.4. TILTED CONVECTIVE PLUMES IN NU
- Page 101 and 102: 4.4. TILTED CONVECTIVE PLUMES IN NU
- Page 103 and 104: 4.4. TILTED CONVECTIVE PLUMES IN NU
- Page 105 and 106: 4.4. TILTED CONVECTIVE PLUMES IN NU
- Page 107 and 108: 4.4. TILTED CONVECTIVE PLUMES IN NU
- Page 109 and 110: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 111 and 112: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 113 and 114: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 115 and 116: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 117 and 118: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 119 and 120: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 121 and 122: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 123 and 124: 4.5. MEAN CIRCULATION AND STRUCTURE
- Page 125 and 126: 4.6. ESTIMATION OF FRICTION PARAMET
- Page 127 and 128: 4.6. ESTIMATION OF FRICTION PARAMET
- Page 129 and 130: 4.6. ESTIMATION OF FRICTION PARAMET
- Page 131 and 132: 4.6. ESTIMATION OF FRICTION PARAMET
- Page 133 and 134: 4.6. ESTIMATION OF FRICTION PARAMET
- Page 135 and 136: 4.6. ESTIMATION OF FRICTION PARAMET
- Page 137 and 138: 4.7. ON THE BASIC STRUCTURE OF OCEA
- Page 139: 4.7. ON THE BASIC STRUCTURE OF OCEA
- Page 143 and 144: 4.7. ON THE BASIC STRUCTURE OF OCEA
- Page 145 and 146: 4.7. ON THE BASIC STRUCTURE OF OCEA
- Page 147 and 148: 4.7. ON THE BASIC STRUCTURE OF OCEA
- Page 149 and 150: 4.7. ON THE BASIC STRUCTURE OF OCEA
- Page 151 and 152: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 153 and 154: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 155 and 156: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 157 and 158: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 159 and 160: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 161 and 162: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 163 and 164: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 165 and 166: 4.8. ESTIMATION OF FRICTION LAWS AN
- Page 167 and 168: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 169 and 170: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 171 and 172: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 173 and 174: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 175 and 176: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 177 and 178: 4.9. ON THE NUMERICAL RESOLUTION OF
- Page 179 and 180: Quatrième partie tel-00545911, ver
- Page 181 and 182: 175 tel-00545911, version 1 - 13 De
- Page 183 and 184: 177 Contents 1 Preface 5 tel-005459
- Page 185 and 186: 179 Chapter 1 Preface tel-00545911,
- Page 187 and 188: 181 Chapter 2 Observing the Ocean t
- Page 189 and 190: 183 Chapter 3 Physical properties o
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