Gravity currents entering a two-layer fluid

1980 ◽  
Vol 100 (4) ◽  
pp. 739-767 ◽  
Author(s):  
Judith Y. Holyer ◽  
Herbert E. Huppert

This paper presents a study of steady gravity currents entering a two-layer system, with the current travelling either along the boundary to form a boundary current, or between the two different layers to form an intrusion. It is shown that, at the front of an intrusion, the streamlines meet at angles of 120° at a stagnation point. For an energy-conserving current the volume inflow rate to the current, the velocity of propagation and the downstream depths are determined. In contrast to the pioneering study of Benjamin (1968), it is found that the depth of the current is not always uniquely determined and it is necessary to use some principle additional to the conservation relationships to determine which solution occurs. An appropriate principle is obtained by considering dissipative currents. In general, if the volume inflow rate to a current is prescribed, the current loses energy in order to maintain a momentum balance. We thus suggest the criterion that the energy dissipation is a maximum for a fixed volume inflow rate. It is postulated that the energy which is lost will go to form a stationary wave train behind the current. A nonlinear calculation is carried out to determine the amplitude and wavelength of these waves for intrusions. Such waves have been observed on intrusions in laboratory experiments and the results of the calculation are found to agree well with the experiments. Similar waves have not been observed on boundary currents because the resulting waves have too much energy and break.

2014 ◽  
Vol 740 ◽  
pp. 97-113 ◽  
Author(s):  
Joseph J. Kuehl ◽  
V. A. Sheremet

AbstractThe problem of oceanic gap-traversing boundary currents, such as the Kuroshio current crossing the Luzon Strait or the Gulf Stream traversing the mouth of the Gulf of Mexico, is considered. Systems such as these are known to admit two dominant states: leaping across the gap or penetrating into the gap forming a loop current. Which state the system will assume and when transitions between states will occur are open problems. Sheremet (J. Phys. Oceanogr., vol. 31, 2001, pp. 1247–1259) proposed, based on idealized barotropic numerical results, that variation in the current’s inertia is responsible for these transitions and that the system admits multiple states. Generalized versions of these results have been confirmed by barotropic rotating-table experiments (Sheremet & Kuehl, J. Phys. Oceanogr., vol. 37, 2007, 1488–1495; Kuehl & Sheremet,J. Mar. Res., vol. 67, 2009, pp. 25–42). However, the typical structure of oceanic boundary currents, such as the Gulf Stream or Kuroshio, consists of an upper-layer intensified flow riding atop a weakly circulating lower layer. To more accurately address this oceanic situation, the present work extends the above findings by considering two-layer rotating table experiments. The flow is driven by pumping water through sponges and vertical seals, creating a Sverdrup interior circulation in the upper layer which impinges on a ridge where a boundary current is formed. The $\beta $ effect is incorporated in both layers by a sloping rigid lid as well as a sloping bottom and the flow is visualized with the particle image velocimetry method. The experimental set-up is found to produce boundary currents consistent with theory. The existence of multiple states and hysteresis, characterized by a cusp topology of solutions, is found to be robust to stratification and various properties of the two-layer system are explored.


2008 ◽  
Vol 38 (10) ◽  
pp. 2294-2307 ◽  
Author(s):  
Hristina G. Hristova ◽  
Joseph Pedlosky ◽  
Michael A. Spall

Abstract A linear stability analysis of a meridional boundary current on the beta plane is presented. The boundary current is idealized as a constant-speed meridional jet adjacent to a semi-infinite motionless far field. The far-field region can be situated either on the eastern or the western side of the jet, representing a western or an eastern boundary current, respectively. It is found that when unstable, the meridional boundary current generates temporally growing propagating waves that transport energy away from the locally unstable region toward the neutral far field. This is the so-called radiating instability and is found in both barotropic and two-layer baroclinic configurations. A second but important conclusion concerns the differences in the stability properties of eastern and western boundary currents. An eastern boundary current supports a greater number of radiating modes over a wider range of meridional wavenumbers. It generates waves with amplitude envelopes that decay slowly with distance from the current. The radiating waves tend to have an asymmetrical horizontal structure—they are much longer in the zonal direction than in the meridional, a consequence of which is that unstable eastern boundary currents, unlike western boundary currents, have the potential to act as a source of zonal jets for the interior of the ocean.


2008 ◽  
Vol 616 ◽  
pp. 327-356 ◽  
Author(s):  
BRIAN L. WHITE ◽  
KARL R. HELFRICH

A steady theory is presented for gravity currents propagating with constant speed into a stratified fluid with a general density profile. Solution curves for front speed versus height have an energy-conserving upper bound (the conjugate state) and a lower bound marked by the onset of upstream influence. The conjugate state is the largest-amplitude nonlinear internal wave supported by the ambient stratification, and in the limit of weak stratification approaches Benjamin's energy-conserving gravity current solution. When the front speed becomes critical with respect to linear long waves generated above the current, steady solutions cannot be calculated, implying upstream influence. For non-uniform stratification, the critical long-wave speed exceeds the ambient long-wave speed, and the critical-Froude-number condition appropriate for uniform stratification must be generalized. The theoretical results demonstrate a clear connection between internal waves and gravity currents. The steady theory is also compared with non-hydrostatic numerical solutions of the full lock release initial-value problem. Some solutions resemble classic gravity currents with no upstream disturbance, but others show long internal waves propagating ahead of the gravity current. Wave generation generally occurs when the stratification and current speed are such that the steady gravity current theory fails. Thus the steady theory is consistent with the occurrence of either wave-generating or steady gravity solutions to the dam-break problem. When the available potential energy of the dam is large enough, the numerical simulations approach the energy-conserving conjugate state. Existing laboratory experiments for intrusions and gravity currents produced by full-depth lock exchange flows over a range of stratification profiles show excellent agreement with the conjugate state solutions.


2006 ◽  
Vol 19 (8) ◽  
pp. 1531-1544 ◽  
Author(s):  
Naoki Sato ◽  
Masaaki Takahashi

Abstract Statistical features of quasi-stationary planetary waves were examined on the subtropical jet in the midsummer Northern Hemisphere by using objectively analyzed data and satellite data. As a result, a quasi-stationary wave train that is highly correlated with the midsummer climate over Japan was identified. A clear phase dependency of the appearance of waves was also confirmed. An analysis of temporal evolution and wave activity flux revealed that the eastward propagation of the wave packet starts in the Middle East, passes over East Asia, and reaches North America. The anomaly pattern is strengthened through kinetic energy conversion near the entrance of the Asian jet over the Middle East. The interaction between the anomaly pattern and the basic field contributes to the appearance of the anomalous wavelike pattern. Although the wave train is correlated with the anomaly of convective activity over the western North Pacific and the Indian Ocean, it is implied that internal dynamics are important in determining the statistical features of the appearance of anomalous quasi-stationary waves on the subtropical jet.


1968 ◽  
Vol 31 (2) ◽  
pp. 209-248 ◽  
Author(s):  
T. Brooke Benjamin

This paper presents a broad investigation into the properties of steady gravity currents, in so far as they can be represented by perfect-fluid theory and simple extensions of it (like the classical theory of hydraulic jumps) that give a rudimentary account of dissipation. As usually understood, a gravity current consists of a wedge of heavy fluid (e.g. salt water, cold air) intruding into an expanse of lighter fluid (fresh water, warm air); but it is pointed out in § 1 that, if the effects of viscosity and mixing of the fluids at the interface are ignored, the hydrodynamical problem is formally the same as that for an empty cavity advancing along the upper boundary of a liquid. Being simplest in detail, the latter problem is treated as a prototype for the class of physical problems under study: most of the analysis is related to it specifically, but the results thus obtained are immediately applicable to gravity currents by scaling the gravitational constant according to a simple rule.In § 2 the possible states of steady flow in the present category between fixed horizontal boundaries are examined on the assumption that the interface becomes horizontal far downstream. A certain range of flows appears to be possible when energy is dissipated; but in the absence of dissipation only one flow is possible, in which the asymptotic level of the interface is midway between the plane boundaries. The corresponding flow in a tube of circular cross-section is found in § 3, and the theory is shown to be in excellent agreement with the results of recent experiments by Zukoski. A discussion of the effects of surface tension is included in § 3. The two-dimensional energy-conserving flow is investigated further in § 4, and finally a close approximation to the shape of the interface is obtained. In § 5 the discussion turns to the question whether flows characterized by periodic wavetrains are realizable, and it appears that none is possible without a large loss of energy occurring. In § 6 the case of infinite total depth is considered, relating to deeply submerged gravity currents. It is shown that the flow must always feature a breaking ‘head wave’, and various properties of the resulting wake are demonstrated. Reasonable agreement is established with experimental results obtained by Keulegan and others.


A simplified analysis is made of the stability of long waves in a sand bed under deep, slow, and steady (or slowly varying) water flow. Allowing for vertical variation in density and shear, the linearized hydrodynamical equations yield a symmetrical flow whose only action is to impart a slow phase velocity to existing sand-waves without altering their amplitude. The only mechanism found under which sand-waves could grow under the assumed conditions is that of a stationary wave train in the lee of a permanent obstacle. The lee-waves require density gradients greater than a certain minimum, independently of any stability due to shear. Application of this model to tidal flow in the Summer thermocline over the Continental Shelf west of Brittany yields a spectrum of wave-building increments which agrees in general wavelength and shape with that of sand-waves measured on La Chapelle Bank (47° 41' N, 7° 13' W). Changes in amplitude of a few sand-grain diameters per year would be expected. Thermal stratification would be insufficient for the same mechanism to generate waves in the North Sea, but the possibility of density gradients due to suspended sediment is suggested as a likely factor of importance.


2014 ◽  
Vol 142 (10) ◽  
pp. 3860-3880 ◽  
Author(s):  
Thomas Dubos ◽  
Marine Tort

Abstract The curl form of equations of inviscid atmospheric motion in general non-Eulerian coordinates is obtained. Narrowing down to a general vertical coordinate, a quasi-Hamiltonian form is then obtained in a Lagrangian, isentropic, mass-based or z-based vertical coordinate. In non-Lagrangian vertical coordinates, the conservation of energy by the vertical transport terms results from the invariance of energy under the vertical relabeling of fluid parcels. A complete or partial separation between the horizontal and vertical dynamics is achieved, except in the Eulerian case. The horizontal–vertical separation is especially helpful for (quasi-)hydrostatic systems characterized by vanishing vertical momentum. Indeed for such systems vertical momentum balance reduces to a simple statement: total energy is stationary with respect to adiabatic vertical displacements of fluid parcels. From this point of view the purpose of (quasi-)hydrostatic balance is to determine the vertical positions of fluid parcels, for which no evolution equation is readily available. This physically appealing formulation significantly extends previous work. The general formalism is exemplified for the fully compressible Euler equations in a Lagrangian vertical coordinate and a Cartesian (x, z) slice geometry, and the deep-atmosphere quasi-hydrostatic equations in latitude–longitude horizontal coordinates. The latter case, in particular, illuminates how the apparent intricacy of the time-dependent metric terms and of the additional forces can be absorbed into a proper choice of prognostic variables. In both cases it is shown how the quasi-Hamiltonian form leads straightforwardly to the conservation of energy using only integration by parts. Relationships with previous work and implications for stability analysis and the derivation of approximate sets of equations and energy-conserving numerical schemes are discussed.


2008 ◽  
Vol 65 (2) ◽  
pp. 536-556 ◽  
Author(s):  
George H. Bryan ◽  
Richard Rotunno

Abstract This study presents analytic results for steady gravity currents in a channel using the deep anelastic equations. Results are cast in terms of a nondimensional parameter H/H0 that relates the channel depth H to a scale depth H0 (the depth at which density goes to zero in an isentropic atmosphere). The classic results based on the incompressible equations correspond to H/H0 = 0. For cold gravity currents (at the bottom of a channel), assuming energy-conserving flow, the nondimensional current depth h/H is much smaller, and nondimensional propagation speed C/(gH)1/2 is slightly smaller as H/H0 increases. For flows with energy dissipation, C/(gH)1/2 decreases as H/H0 increases, even for fixed h/H. The authors conclude that as H/H0 increases the normalized hydrostatic pressure rise in the cold pool increases near the bottom of the channel, whereas drag decreases near the top of the channel; these changes require gravity currents to propagate slower for steady flow to be maintained. From these results, the authors find that steady cold pools have a likely maximum depth of 4 km in the atmosphere (in the absence of shear). For warm gravity currents (at the top of a channel), h/H is slightly larger and C/(gH)1/2 is much larger as H/H0 increases. The authors also conduct two-dimensional numerical simulations of “lock-exchange flow” to provide an independent evaluation of the analytic results. For cold gravity currents the simulations support the analytic results. However, for warm gravity currents the simulations show unsteady behavior that cannot be captured by the analytic theory and which appears to have no analog in incompressible flow.


The linearized theory of unsteady wind-driven currents in a horizontally stratified ocean is applied to the northern part of the Indian Ocean. This is argued to be a suitable area for detailed application and evaluation of the theory because (i) the theory has certain advantages near the equator (for example, influence of detailed bottom topography is reduced, thermoclines are somewhat less variable in character, and speeds of baroclinic propagation are enhanced relative to current speeds), and (ii) the wind-stress pattern undergoes a well marked change with onset of the Southwest Monsoon, a change to which the pattern of currents shows a more or less identifiable, and rather quick, response which may be compared with theoretical predictions. Response is predicted to be found principally in two modes as far as vertical distribution of current is concerned; to a somewhat lesser extent in the barotropic mode with uniform distribution, and to a somewhat greater extent in the first baroclinic mode with current distribution as in figure 7, concentrated predominantly in the uppermost 200 m (see Appendix for detailed analysis of the modes appropriate to the equatorial Indian Ocean). Of particular interest is the strong Somali Current, that flows northward along the Somali coast only during the northern hemisphere summer (after monsoon onset) but during that time is comparable in volume flow (about 5 x 107 m3/s) to other western boundary currents such as the Gulf Stream. Detailed discussion of the application of linearized theory to equatorial oceans with western boundaries leads the author to conclude, both in the barotropic (§ 2) and baroclinic (§ 4) cases, that c wave packets5 of current pattern reaching such a boundary deposit the c flux5 they carry (velocity normal to the boundary integrated along it) in a boundary current which rather rapidly takes a rather concentrated form. Linear theory with horizontal transport neglected indicates that such flux requires of the order of 10 days to become concentrated in a current of 100 km width, but that thereafter it continues to become still thinner; however, with horizontal transport included, a steady-state finite thickness of current is reached. In reality, nonlinear effects would play an important additional part in limiting steady-state current thickness to the observed 100 km or thereabouts, but the time scale required to bring the thickness down to this value is probably given reasonably well by linear theory. Calculations for a zonal distribution of winds, which rather rapidly make a reversal of direction and increase of strength somewhat north of the Equator characteristic of the onset of the Southwest Monsoon, predict westward propagation of both barotropic and baroclinic wave energy at comparable speeds of the order of 1 m/s; the marked contrast here with other oceans (in the comparability of speeds) is given particularly detailed study. Calculations indicate that the barotropic signal is considerably distorted (figure 3) by the fact that low-wavenumber components reach the western boundary first. Baroclinic propagation takes the form of special planetary-wave modes concentrated near the equator (§3), of which perhaps four, delivering flux patterns depicted in figure 5, and possessing wave velocities of 0.9, 0.55, 0.4 and 0.3 m/s towards the west, are specially relevant to generation of the Somali Current. Peak surface flows in that current are predicted to be influenced about three times as much by this baroclinic propagation as by the barotropic. Theory indicates 1 month (of which two-thirds is needed for propagation of current patterns and one-third for their concentration in a boundary current) as characteristic time scale for formation of the Somali Current (see figure 6 in particular for the calculated baroclinic component) in contradistinction to the ‘decades’ predicted by the same type of theory in mid-latitude oceans (Veronis & Stommel 1956). Observations do, indeed, make clear that the time scale is not significantly more than 1 month, although the possibility that it might be still less cannot yet be decided on the basis of observational evidence. The flow is calculated as reaching 40 % of a typical maximum value (observed in August) already within 1 month of monsoon onset (May), even though no effect of wind stress acting within 500 km of the coast has been taken into account. The linearized theory predicts the current as reaching as far north as 6° N or 7°N, but nonlinear terms are generally found in computational studies (Bryan 1963; Veronis 1966) to bring about some ‘ inertial overshoot ’ in concentrated boundary currents, which may explain why the current does not in fact separate until about 9°N.


2005 ◽  
Vol 51 (172) ◽  
pp. 125-138 ◽  
Author(s):  
Perry Bartelt ◽  
Othmar Buser ◽  
Martin Kern

AbstractWe derive work dissipation functionals for granular snow avalanches flowing in simple shear. Our intent is to apply constructive theorems of non-equilibrium thermodynamics to the snow avalanche problem. Snow chute experiments show that a bi-layer system consisting of a non-yielded flow plug overriding a sheared fluidized layer can be used to model avalanche flow. We show that for this type of constitutive behaviour the dissipation functionals are minimum at steady state with respect to variations in internal velocity; however, the functionals must be constrained by subsidiary mass- continuity integrals before the equivalence of momentum balance and minimal work dissipation can be established. Constitutive models that do not satisfy this equivalence are henceforth excluded from our consideration. Fluctuations in plug and slip velocity depend on the roughness of the flow surface and viscosity of the granular system. We speculate that this property explains the transition from flowing avalanches to powder avalanches. Because the temperature can safely be assumed constant, we demonstrate within the context of non-equilibrium thermodynamics that granular snow avalanches are irreversible, dissipative systems, minimizing – in space – entropy production. Furthermore, entropy production is linear both near and far from steady-state non-equilibrium because of the mass-continuity constraint. Finally, we derive thermodynamic forces and conjugate fluxes as well as expressing the corresponding phenomenological Onsager coefficients in terms of the constitutive parameters.


Sign in / Sign up

Export Citation Format

Share Document