Table of Contents
ISRN Thermodynamics
VolumeΒ 2012, Article IDΒ 609701, 15 pages
Research Article

Thermodynamics/Dynamics Coupling in Weakly Compressible Turbulent Stratified Fluids

Department of Meteorology, University of Reading, Earley Gate, P.O. Box 243, Reading RG6 6BB, UK

Received 6 December 2011; Accepted 22 January 2012

Academic Editors: E.Β Curotto and N.Β Marmier

Copyright Β© 2012 RΓ©mi Tailleux. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


In traditional and geophysical fluid dynamics, it is common to describe stratified turbulent fluid flows with low Mach number and small relative density variations by means of the incompressible Boussinesq approximation. Although such an approximation is often interpreted as decoupling the thermodynamics from the dynamics, this paper reviews recent results and derive new ones that show that the reality is actually more subtle and complex when diabatic effects and a nonlinear equation of state are retained. Such an analysis reveals indeed: (1) that the compressible work of expansion/contraction remains of comparable importance as the mechanical energy conversions in contrast to what is usually assumed; (2) in a Boussinesq fluid, compressible effects occur in the guise of changes in gravitational potential energy due to density changes. This makes it possible to construct a fully consistent description of the thermodynamics of incompressible fluids for an arbitrary nonlinear equation of state; (3) rigorous methods based on using the available potential energy and potential enthalpy budgets can be used to quantify the work of expansion/contraction 𝐡 in steady and transient flows, which reveals that 𝐡 is predominantly controlled by molecular diffusive effects, and act as a significant sink of kinetic energy.

1. Introduction

A large class of fluid flows of interest in traditional and geophysical fluid dynamics are characterised by fluid velocities much smaller than the speed of sound (low Mach number), strong density stratification yet with small relative density variations, and high Reynolds number (i.e., they are turbulent). It has become common practice to regard such flows as incompressible or nearly incompressible, and to describe them by means of particular approximations to the fully compressible Navier-Stokes equations known as the Boussinesq and anelastic approximations. Traditionally, such approximations are derived in the idealised context of purely adiabatic motions, hence excluding irreversible processes such as molecular viscous and diffusive processes. The equation of state is also usually linearised, so that density anomalies become proportional to temperature anomalies. In that case, it is easily shown that the resulting approximations filter out acoustic waves and thus fully decouple the dynamics from the thermodynamics, in such a way that the mechanical energy (the sum of kinetic energy (KE) and gravitational potential energy (GPE)) and the internal energy (IE) obey independent conservation laws.

The coupling between dynamics and thermodynamics is reinstated, however, when the Boussinesq and anelastic approximations are appended with a representation of irreversible processes due to molecular viscous and diffusive processes, and/or when the nonlinearities of the equation of state are no longer neglected. Such an approach is used for instance in the field of numerical ocean modelling. The way the coupling between the dynamics and thermodynamics manifests itself in that case is through the appearance of nonconservative terms in the local balance equation for the mechanical energy, which indicates the need for conversions with internal energy if total energy is to be conserved. Whereas the approximate forms of KE and GPE are close to their nonapproximated forms, the approximate form of IE has remained so far very mysterious, given that only the conversion terms with IE are explicitly represented in the Boussinesq and anelastic approximations, not IE itself. Since the approximate form of IE is in general unknown and left implicit, it is in general not straightforward to determine whether the approximation obtained by adding viscous and diffusive terms has a well-defined energy budget. Moreover, understanding the role of IE in such approximations may give rise to misinterpretations as in the case discussed by Tailleux [1].

The lack of explicit knowledge about the form assumed by the internal energy in the Boussinesq and anelastic approximations makes it difficult to understand the precise role played by internal energy (and hence the dynamics/thermodynamics coupling) in incompressible or nearly incompressible turbulent stratified fluid flows. Until now, this difficulty has been systematically avoided by assuming the conversions to and from internal energy, and hence the dynamics/thermodynamics coupling, to be dynamically unimportant at leading order. Recently, however, Tailleux [1] challenged this assumption by showing from first principles that it was inconsistent with the second law, by suggesting that the coupling between dynamics and thermodynamics should increase both with the strength of the stratification and the degree of turbulence in the fluid.

The main objective of this paper will be to provide insights into the dynamics/thermodynamics coupling in weakly compressible turbulent stratified fluid flows. In particular, it will serve as follows.(1)It will provide a general overview of the fundamental physical and technical issues involved in understanding the coupling between dynamics and thermodynamics in turbulent stratified fluids. It will explain the importance of isolating the part of the total potential energy (gravitational + potential) that is available for reversible conversions into kinetic energy, as well as the part of the internal energy available for irreversible conversions into mechanical energy. This is based on the results by Tailleux [1, 2].(2)It will show how to construct explicitly the full range of known thermodynamic potentials, including entropy, internal energy, enthalpy, that have previously remained implicit in the Boussinesq and anelastic approximations. It will show how it is possible to establish a formal correspondence between the work of expansion/contraction in compressible fluids, and the changes in gravitational potential energy in the Boussinesq and anelastic approximations.(3)It will show how the above results can help clarify and provide new insights into the energetics of turbulent mixing in turbulent stratified fluid flows, as discussed in Tailleux [1].(4)It will discuss the general issue of how to estimate the overall work of expansion/contraction in steady-state buoyancy driven circulations, which has received much attention over the past decade in relation with understanding the relative importance of the surface buoyancy fluxes in driving and stirring the ocean circulation.

2. Energetics and Dynamics/Thermodynamics Coupling in Compressible Fluids

2.1. Compressible Navier-Stokes Equations for a Binary Fluid

Issues pertaining to the coupling between the dynamics and thermodynamics in turbulent-stratified fluids are fundamentally rooted in the description of the energetics of a fully compressible fluid. For most fluids of interest, this can be achieved in the context of the compressible Navier-Stokes equations. To fix ideas, these equations are given here for a binary fluid with a nonlinear equation of state such as seawater, namely,πœŒπ·π―π·π‘‘+βˆ‡π‘ƒ=βˆ’πœŒβˆ‡Ξ¦+βˆ‡β‹…π’,(1)π·πœŒπ·π‘‘+πœŒβˆ‡β‹…π―=0,(2)πœŒπ·πœ‚π·π‘‘=πœŒΜ‡π‘žπ‘‡=βˆ’βˆ‡β‹…ξ€·πœŒπ…πœ‚ξ€Έ+πœŒΜ‡πœ‚irr,(3)πœŒπ·π‘†π·π‘‘=βˆ’βˆ‡β‹…ξ€·πœŒπ…π‘†ξ€Έ,(4)𝑃=𝑃(πœ‚,𝑆,𝜐)=βˆ’πœ•π‘’πœ•πœ,𝑇=𝑇(πœ‚,𝑆,𝜐)=πœ•π‘’πœ•πœ‚,(5) where 𝐯=(𝑒,𝑣,𝑀) is the three-dimensional velocity field, 𝜌 is the density, 𝑃 is the pressure, Ξ¦=𝑔0𝑍 is the geopotential expressed in terms of a constant gravitational acceleration 𝑔0 and geopotential height 𝑍=𝑍(𝑧) and 𝑧 the geometrical height, 𝐒 is the deviatoric stress tensor, πœ‚ is the specific entropy, 𝑆 is salinity, 𝑒=𝑒(πœ‚,𝑆,𝜐) is the specific internal energy regarded as a function of entropy, salinity, and specific volume 𝜐=1/𝜌, 𝑇 is the temperature, πœ… is the molecular diffusivity, 𝑐𝑝 is the specific heat at constant pressure, πœ€πΎ is the viscous dissipation rate, π…πœ‚ and 𝐅𝑆 are molecular diffusive fluxes of entropy and salt, while Μ‡πœ‚irr>0 is the irreversible production of entropy imposed by the second law of thermodynamics. In practice, these equations need to be supplemented with appropriate boundary conditions, for which specific examples will be discussed in the text.

2.2. Energetics and Dynamics/Thermodynamics Coupling

In order to discuss the nature of the coupling between dynamics and thermodynamics for a compressible fluid, we first form a local evolution equation for the mechanical energy πΈπ‘š=𝐯2/2+Ξ¦, that is, the sum of kinetic energy and gravitational potential energy, by multiplying the momentum equation (1); after some rearrangement, one thus obtainsπœŒπ·π·π‘‘ξ‚Έπ―22+Ξ¦ξ‚Ή+βˆ‡β‹…ξ€·π‘ƒπ―βˆ’πœŒπ…π‘˜π‘’ξ€Έ=πœŒπ‘ƒπ·πœπ·π‘‘βˆ’πœŒπœ€πΎ,(6) which can also be rewritten in conservative form as follows:πœ•ξ€·πœŒπΈπ‘šξ€Έπœ•π‘‘+βˆ‡β‹…ξ‚ΈπœŒξ‚΅πΈπ‘š+π‘ƒπœŒξ‚Άπ―βˆ’πœŒπ…π‘˜π‘’ξ‚Ή=πœŒπ‘ƒπ·πœπ·π‘‘βˆ’πœŒπœ€πΎ,(7) where the work against the stress tensor π―β‹…βˆ‡β‹…π’=βˆ‡β‹…(πœŒπ…π‘˜π‘’)βˆ’πœŒπœ€πΎ has been decomposed into the divergence of the viscous flux of kinetic energy minus the viscous dissipation rate πœŒπœ€πΎ. Physically, the presence of the right-hand side in both (6) and (7) indicates that mechanical energy πΈπ‘š is not a conservative quantity. The principle of total energy conservation suggests that the r.h.s. must therefore represent a conversion with internal energy.

In order to make this apparent, we now turn to the derivation of a local evolution equation for the internal energy. From the well-known expression of its total differential form: d𝑒=𝑇dπœ‚+πœ‡dπ‘†βˆ’π‘ƒd𝜐, combined with the evolution equations for entropy and salt (3) and (4), it follows thatπœŒπ·π‘’π·π‘‘=πœŒξ‚Έπ‘‡π·πœ‚π·π‘‘+πœ‡π·π‘†π·π‘‘βˆ’π‘ƒπ·πœπ·π‘‘ξ‚Ή=π‘‡ξ€Ίβˆ’βˆ‡β‹…ξ€·πœŒπ…πœ‚ξ€Έ+πœŒΜ‡πœ‚irrξ€»βˆ’πœ‡β‹…βˆ‡β‹…ξ€·πœŒπ…π‘†ξ€Έβˆ’πœŒπ‘ƒπ·πœπ·π‘‘=βˆ’βˆ‡β‹…ξ€ΊπœŒξ€·π‘‡π…πœ‚+πœ‡π…π‘†ξ€Έξ€»+πœŒξ€Ίπ…πœ‚β‹…βˆ‡π‘‡+π…π‘†β‹…βˆ‡πœ‡+π‘‡Μ‡πœ‚irrξ€»βˆ’πœŒπ‘ƒπ·πœπ·π‘‘.(8) By summing the mechanical energy equation with the internal energy, we arrive at the following result:πœ•(𝜌𝐸)πœ•π‘‘+βˆ‡β‹…ξ‚ΈπœŒξ‚΅πΈ+π‘ƒπœŒξ‚Άπ―βˆ’πœŒπ…π‘˜π‘’+πœŒπ…π‘žξ‚Ή=+πœŒξ€Ίπ…πœ‚β‹…βˆ‡π‘‡+π…π‘†β‹…βˆ‡πœ‡+π‘‡Μ‡πœ‚irrβˆ’πœ€πΎξ€»ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ŒΜ‡β€ŒπΈirr,(9) where 𝐸=πΈπ‘š+𝑒 is the total energy, while we also defined π…π‘ž=π‘‡π…πœ‚+πœ‡π…π‘† as a generalised heat flux. This in turn provides the expression for the entropy flux π…πœ‚=(π…π‘žβˆ’πœ‡π…π‘†)/𝑇. Now, the principle of total energy conservation implies that the total energy within a control volume can only change through energy and mass fluxes through this control volume. This in turn requires that the term Μ‡β€ŒπΈirr vanishes identically, which imposes the following constraint on the irreversible entropy production term:Μ‡πœ‚irr=πœ€πΎβˆ’π…πœ‚β‹…βˆ‡π‘‡βˆ’π…π‘†β‹…βˆ‡πœ‡π‘‡.(10) This in turn implies constraints on the admissible forms of the entropy flux π…πœ‚ and salt flux 𝐅𝑆, which need to be such that Μ‡πœ‚irrβ‰₯0 in all circumstances. We also see that the irreversible production of entropy must vanish at thermodynamic equilibrium, which for a binary fluid is achieved when 𝑇 and πœ‡ become uniform. This is clearly the case here. It follows from these considerations that the internal energy equation can therefore be rewritten asπœŒπ·π‘’π·π‘‘+βˆ‡β‹…ξ€·πœŒπ…π‘žξ€Έ=βˆ’πœŒπ‘ƒπ·πœπ·π‘‘+πœŒπœ€πΎ.(11) Note that this expression is consistent with the classical form of the first law of thermodynamics 𝐷𝑒/𝐷𝑑=Μ‡π‘ž+Μ‡β€Œπ‘€, with Μ‡π‘ž=πœ€πΎβˆ’πœŒβˆ’1βˆ‡β‹…(πœŒπ…π‘ž) as the local rate of heat transfer, and Μ‡β€Œπ‘€=βˆ’π‘ƒπ·πœ/𝐷𝑑 as the local rate of work transfer. According to (11), the coupling between the dynamics and thermodynamics can be regarded as the sum of two effects. The first effect is due to the viscous dissipation of kinetic energy into the β€œheat” part of internal energy, while the second is related to the conversion between the elastic part of internal energy and kinetic energy. It is the latter term that is usually assumed to be of secondary importance in weakly compressible turbulent stratified fluids at low Mach number, and therefore the one under focus in this paper.

2.3. Energetics of Purely Buoyancy-Driven Circulations

As further motivation and justification for investigating the coupling between the dynamics and thermodynamics in turbulent-stratified fluids, let us simply point out that the work of expansion/contraction plays a central role in the energetics of purely buoyancy-driven circulations, as in horizontal convection, for example, Hughes and Griffiths [3], which has received much attention in the oceanography community over the past ten years. Thus, for a fluid heated and cooled at the top, assumed to be a surface of constant geopotential, the mechanical energy balance is as follows:ξ€œπ‘‰π‘ƒπ·πœπ·π‘‘dπ‘šξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ΅=ξ€œπ‘‰πœŒπœ€πΎd𝑉,(12) where dπ‘š=𝜌d𝑉 is the mass of an elementary fluid element, and expresses the fact that in a steady state, the work of expansion/contraction must balance the volume-integrated viscous dissipation. Such a balance was recently discussed in details by Tailleux [2] in the context of the β€œocean heat engine controversy.” The importance of understanding what such a balance tells us about 𝐡 and the viscous dissipation stems from the central role played by πœ€πΎ in the theory of turbulent flows, since it is the primary quantity entering Kolmogorov theory of turbulence for instance. The riddle here comes from the apparent conflict between the traditional belief that compressibility effects are small in weakly compressible fluids such as water and seawater, which on the basis of (12) would tend to suggest that flows in such fluids cannot really be large or be associated with significant amount of viscous dissipation, and the observation that even in water, buoyancy-driven flows can be sometimes vigorous, as in Rayleigh-BΓ©nard convection, for instance, for which (12) is also supposed to hold. The primary purpose of this paper is to convince the reader that compressible effects in turbulent stratified fluids are significantly more important than usually assumed, and that such large compressible effects are actually described both by the incompressible Boussinesq and Anelastic approximations, in contrast to widely held beliefs to the contrary.

3. Dynamics/Thermodynamics Coupling in Boussinesq/Anelastic Fluids

3.1. Thermodynamically Consistent Boussinesq/Anelastic Approximations

Two main classes of β€œsound-proof” approximations are commonly used for the study of weakly compressible fluid flows at low Mach numbers, such as those typical of oceanic and atmospheric flows. These are the Boussinesq approximation, for example, Boussinesq [4], and the anelastic approximation, for example, Ogura [5]. As shown recently by Pauluis [6], improving on earlier ideas by Ingersoll [7], it is possible to derive a thermodynamically and energetically consistent set of anelastic equations for a binary fluid with an arbitrary nonlinear equation of state that also includes the Boussinesq approximation as a particular case. For the present purposes, these can be written as follows:𝐷𝐯𝐷𝑑+βˆ‡β‹…ξ‚΅π›Ώπ‘ƒπœŒ0ξ‚Ά=βˆ’ξ‚΅πœŒβˆ’πœŒ0𝜌0ξ‚Άβˆ‡Ξ¦+1𝜌0βˆ‡β‹…π’,(13)βˆ‡β‹…ξ€·πœŒ0𝐯=0,(14)𝜌0π·πœ‚π·π‘‘=πœŒΜ‡π‘žπ‘‡=βˆ’βˆ‡β‹…ξ€·πœŒ0π…πœ‚ξ€Έ+Μ‡πœ‚irr,(15)𝜌0𝐷𝑆𝐷𝑑=βˆ’βˆ‡β‹…ξ€·πœŒ0𝐅𝑆,(16) where 𝐯=(𝑒,𝑣,𝑀) is the three-dimensional velocity field, 𝛿𝑃=π‘ƒβˆ’π‘ƒ0 is the difference between the total pressure 𝑃 and the reference pressure 𝑃0, 𝜌 and 𝜌0 are the density and reference density respectively, Ξ¦=𝑔0𝑍 is the geopotential, which can be written as the product of a constant gravitational acceleration times the geopotential height, a function of depth alone; 𝐒 is the stress tensor, πœ€πΎ is the viscous dissipation rate, πœ‚ is the specific entropy, and 𝑆 is the salinity. Pauluis [6] regards the second equation as one for the total water content. The present results are applicable to an arbitrary equation of state for a binary fluid in principle.

3.2. Proof of Energetic and Thermodynamic Consistency

Prior to the studies by Ingersoll [7], Pauluis [6], Young [8], Nycander [9], and Tailleux [10], it was generally thought to be impossible to derive energetically and thermodynamically consistent Boussinesq and anelastic approximations for a binary fluid with an arbitrary nonlinear equation of state. It is therefore of interest to verify such consistency, as it is closely related to understanding the nature of the dynamics/thermodynamics coupling in the BA system. To that end, we first derive the kinetic energy equation similarly as for the compressible case addressed above, which yields𝜌0𝐷𝐷𝑑𝐯22+βˆ‡β‹…ξ€Ίπ›Ώπ‘ƒπ―βˆ’πœŒ0π…π‘˜π‘’ξ€»=𝜌0π‘π·π‘π·π‘‘βˆ’πœŒ0πœ€πΎ,(17) where 𝑏=βˆ’π‘”0(πœŒβˆ’πœŒ0)/𝜌0 is the buoyancy, and where π―β‹…βˆ‡β‹…π’=βˆ‡β‹…(𝜌0π…π‘˜π‘’)βˆ’πœŒ0πœ€πΎ was split into a divergent term, and the classical viscous dissipation term, as previously.

The previous approaches to establishing the energetic consistency of the BA system by Ingersoll [7], Pauluis [6], Young [8], and Nycander [9] all similarly relied on linking the term 𝜌0𝑏𝐷𝑍/𝐷𝑑 to the following function:β„Žβ€‘(πœ‚,𝑆,𝑍)=β„Ž0(πœ‚,𝑆)βˆ’ξ€œπ‘π‘0π‘ξ‚€πœ‚,𝑆,π‘ξ…žξ‚d𝑍,(18) where β„Ž0(πœ‚,𝑆)=β„Žβ€‘(πœ‚,𝑆,𝑍0) is a priori an arbitrary function of the adiabatically conserved variables πœ‚ and 𝑆. The evolution equation for such a function is, therefore,𝜌0π·β„Žβ€‘π·π‘‘=𝜌0πΆπœ‚π·πœ‚π·π‘‘+𝜌0πΆπ‘†π·π‘†π·π‘‘βˆ’πœŒ0𝑏𝐷𝑍𝐷𝑑,(19) whereπΆπœ‚=πœ•β„Ž0πœ•πœ‚βˆ’ξ€œπ‘π‘0πœ•π‘πœ•πœ‚ξ‚€πœ‚,𝑆,π‘ξ…žξ‚dπ‘ξ…ž,𝐢𝑆=πœ•β„Ž0πœ•π‘†βˆ’ξ€œπ‘π‘0πœ•π‘πœ•Sξ‚€πœ‚,𝑆,π‘ξ…žξ‚dπ‘ξ…ž.(20) Now, combining the kinetic energy equation with (19) leads to:𝜌0𝐷𝐷𝑑v22+β„Žβ€‘ξ‚Ή+βˆ‡β‹…ξ€Ίπ›Ώπ‘ƒπ―βˆ’πœŒ0π…π‘˜π‘’ξ€»=𝜌0ξ‚ΈπΆπœ‚π·πœ‚π·π‘‘+πΆπ‘†π·π‘†π·π‘‘βˆ’πœ€πΎξ‚Ή.(21) This is an important result, which clearly shows that the energy quantity 𝜌0(𝐯2/2+β„Žβ€‘) is conserved in absence of viscous (πœ€πΎ=0) and diabatic effects (π·πœ‚/𝐷𝑑=𝐷𝑆/𝐷𝑑=0) and hence can be regarded as the relevant energy of the BA system in such conditions. As argued by Tailleux [10], however, it is important to realize that the principle of energy conservation must hold both for adiabatic/inviscid conditions, as well as for diabatic/viscous ones. In the latter case, the r.h.s. of (21) will not vanish in general, which calls for the existence of an additional energy quantity β„Žβˆ— satisfying an evolution equation of the type:𝜌0π·β„Žβˆ—π·π‘‘+βˆ‡β‹…ξ€·πœŒ0π…β„Žξ€Έ=βˆ’πœŒ0ξ‚ΈπΆπœ‚π·πœ‚π·π‘‘+πΆπ‘†π·π‘†π·π‘‘βˆ’πœ€πΎξ‚Ή,(22) where the r.h.s. of (22) exactly cancels out the r.h.s. of (21), making the quantity 𝐯2+β„Žβ€‘+β„Žβˆ— the relevant total energy for the BA system, rather than the quantity 𝐯2/2+β„Žβˆ—. The problem of the latter approach, however, is that it focuses on the coupling between the adiabatic energy 𝐯2+β„Žβ€‘ and the diabatic energy β„Žβˆ—, which because of the function β„Ž0(πœ‚,𝑆) does not appear to be uniquely defined. As a result, this makes it difficult to relate this way of discussing the energetics of the BA system to the classical approach.

Tailleux [10] showed that the energetic and thermodynamic consistency of the BA system can in fact be approached in a way closer to the classical approach by rewriting the kinetic energy equation as follows:𝜌0𝐷𝐷𝑑𝐯22βˆ’π‘ξ€·π‘βˆ’π‘0ξ€Έξ‚Ά+βˆ‡β‹…ξ€Ίπ›Ώπ‘ƒπ―βˆ’πœŒ0π…π‘˜π‘’ξ€»=βˆ’πœŒ0ξ€·π‘βˆ’π‘0ξ€Έπ·π‘π·π‘‘βˆ’πœŒ0πœ€πΎ(23) by using the result that 𝑏d𝑍=d[𝑏(π‘βˆ’π‘0)]βˆ’(π‘βˆ’π‘0)d𝑏, which in conservative form becomesπœ•ξ€·πœŒ0πΈπ‘šξ€Έπœ•π‘‘+βˆ‡β‹…ξ‚ΈπœŒ0ξ‚΅πΈπ‘š+π›Ώπ‘ƒπœŒ0ξ‚Άπ―βˆ’πœŒ0π…π‘˜π‘’ξ‚Ή=βˆ’πœŒ0ξ€·π‘βˆ’π‘0ξ€Έπ·π‘π·π‘‘βˆ’πœŒ0πœ€πΎ(24) by defining the total mechanical energy by 𝜌0πΈπ‘š=𝜌0𝐯2/2βˆ’πœŒ0𝑏(π‘βˆ’π‘0)=𝜌0𝐯2/2+𝑔0(πœŒβˆ’πœŒ0)(π‘βˆ’π‘0), which is consistent with the traditional definition. Now, it is easily seen that (23) and (24) are directly comparable with the compressible mechanical energy equations (6) and (7) provided that it is meaningful to establish the following correspondence:πœŒπ‘ƒπ·πœπ·π‘‘βŸ·βˆ’πœŒ0ξ€·π‘βˆ’π‘0𝐷𝑏𝐷𝑑.(25) Tailleux [10] shows that such a correspondence is legitimate. In the Boussinesq case (𝜌0=constant), this can be shown by expanding the pressure and specific volume as follows: 𝑃=π‘ƒπ‘Žβˆ’πœŒ0𝑔0𝑍+β‹― and 𝜐=1/𝜌=1/𝜌0βˆ’(πœŒβˆ’πœŒ0)/𝜌20+β‹―, so thatπœŒπ‘ƒπ·πœπ·π‘‘β‰ˆβˆ’πœŒ0ξ€·π‘ƒπ‘Žβˆ’πœŒ0𝑔0π‘ξ€Έπ·π·π‘‘βŽ›βŽœβŽπœŒβˆ’πœŒ0𝜌20⎞⎟⎠=βˆ’πœŒ0ξ€·π‘βˆ’π‘0𝐷𝑏𝐷𝑑(26) in agreement with (25), provided that one defines 𝑍0 by 𝑍0=π‘ƒπ‘Ž/(𝜌0𝑔0). Although the anelastic case is somewhat more involved, it turns out that such a correspondence can also be established in the anelastic case, as discussed by Tailleux [10].

3.3. Thermodynamics of Boussinesq and Anelastic Fluids

The possibility to regard the term π›Ώπ‘Šπ‘π‘Ž=(π‘βˆ’π‘0)d𝑏 as the counterpart of the compressible work π›Ώπ‘Š=βˆ’π‘ƒd𝜐 is an important result, because it naturally points to the method by which to construct the thermodynamics of the Boussinesq/Anelastic system. Indeed, let us recall that a basic tenet of classical thermodynamics is that the specific internal energy 𝑒=𝑒(πœ‚,𝑆,𝜐), whose total differential is given byd𝑒=𝑇dπœ‚+πœ‡dπ‘†βˆ’π‘ƒd𝜐,(27) where πœ‡ is the relative chemical potential of seawater in seawater, is a function of state independent of the thermodynamic path followed. As a result, it follows that the generalised heat transfer differential 𝛿𝑄=𝑇dπœ‚+πœ‡d𝑆 is not independent of the work transfer π›Ώπ‘Š=βˆ’π‘ƒd𝜐. In classical thermodynamics, this interdependence is imposed by the so-called Maxwell relationships, for example, Callen [11], which simply expresses the result that the cross-derivatives with respect to two different variables must be equal for twice continuously differentiable functions. In the present case, the assumption that 𝑒 is a function of the thermodynamic state only implies that𝑇=πœ•π‘’πœ•πœ‚,πœ‡=πœ•π‘’πœ•π‘†,𝑃=βˆ’πœ•π‘’πœ•πœ,(28) which in turn implies the following three Maxwell relationships:πœ•π‘‡πœ•πœ=πœ•2π‘’πœ•πœ‚πœ•πœ=βˆ’πœ•π‘ƒπœ•πœ‚,πœ•π‘‡πœ•π‘†=πœ•2π‘’πœ•πœ‚πœ•π‘†=πœ•πœ‡πœ•πœ‚,πœ•πœ‡πœ•πœ=πœ•2π‘’πœ•π‘†πœ•πœ=βˆ’πœ•π‘ƒπœ•π‘†.(29) The above remarks suggest that the simplest way to ensure that the BA system is energetically and thermodynamically consistent is to ensure that the approximation to the generalised heat transfer, denoted here by π›Ώπ‘„π‘π‘Ž=π‘‡π‘π‘Ždπœ‚π‘π‘Ž+πœ‡π‘π‘Ždπ‘†π‘π‘Ž, is similarly linked via relevant Maxwell relationships to the approximation to the work transfer identified previously, namely,π›Ώπ‘Šπ‘π‘Ž=ξ€·π‘βˆ’π‘0ξ€Έd𝑏,(30) where the subscript ba was used to denote the approximation to the generalised β€œheat variables and functions” 𝑇, πœ‚, πœ‡, and 𝑆. If so, this would in turn allows one to regard the following expression:dπ‘’π‘π‘Ž=π‘‡π‘π‘Ždπœ‚+πœ‡π‘π‘Žd𝑆+ξ€·π‘βˆ’π‘0ξ€Έd𝑏,(31) as the natural counterpart of the fundamental relation of thermodynamics (27), and hence as the total differential of the relevant approximation to the internal energy of the β€œBA fluid,” for which the natural variables are entropy πœ‚, salinity 𝑆, and buoyancy 𝑏. In practice, however, it is often more convenient to work with pressure 𝑃 as a dependent variable rather than specific volume, as well as with temperature 𝑇 rather than entropy πœ‚, which motivates the introduction of additional thermodynamic potentials constructed from (27) by means of the Legendre transform, for example, Alberty [12]. The most common thermodynamic potentials that are also the most relevant for the present work are the specific enthalpy β„Ž=𝑒+π‘πœ and the Gibbs free energy 𝑔=𝑒+π‘πœβˆ’π‘‡πœ‚, whose natural dependent variables are (πœ‚,𝑆,𝑃) and (𝑇,𝑆,𝑃), respectively. From (31), it is easy to convince oneself that the corresponding approximations to β„Ž and 𝑔 are given by β„Žπ‘π‘Ž=π‘’π‘π‘Žβˆ’π‘(π‘βˆ’π‘0) and π‘”π‘π‘Ž=π‘’π‘π‘Žβˆ’π‘(π‘βˆ’π‘0)βˆ’π‘‡π‘π‘Žπœ‚π‘π‘Ž, with the following total differentials:dβ„Žπ‘π‘Ž=π‘‘ξ€Ίπ‘’π‘π‘Žβˆ’π‘ξ€·π‘βˆ’π‘0ξ€Έξ€»=π‘‡π‘π‘Ždπœ‚π‘π‘Ž+πœ‡π‘π‘Ždπ‘†π‘π‘Žβˆ’π‘d𝑍(32)dπ‘”π‘π‘Ž=π‘‘ξ€Ίπ‘’π‘π‘Žβˆ’π‘ξ€·π‘βˆ’π‘0ξ€Έβˆ’π‘‡π‘π‘Žπœ‚π‘π‘Žξ€»=βˆ’πœ‚π‘π‘Ždπ‘‡π‘π‘Ž+πœ‡π‘π‘Ždπ‘†π‘π‘Žβˆ’π‘d𝑍.(33) Until now, the above considerations have remained rather formal, and while they indicate that it is in principle possible to construct the full range of known thermodynamic potentials for a BA fluid, they have not addressed the issue of how such potentials might be constructed in practice.

In the previous approaches by Ingersoll [7], Pauluis [6], Young [8], and Nycander [9], the authors have generally assumed the buoyancy 𝑏=𝑏(πœ‚,𝑆,𝑍) to be given as a function of the dependent variables considered, typically entropy, salinity, and geopotential height 𝑍. From (32), this can be integrated with respect to 𝑍, namely,β„Žπ‘π‘Ž(πœ‚,𝑆,𝑍)=β„Žπ‘π‘Žξ€·πœ‚,𝑆,𝑍0ξ€Έβˆ’ξ€œπ‘π‘0π‘ξ‚€πœ‚,𝑆,π‘ξ…žξ‚dπ‘ξ…ž,(34) for some reference geopotential height 𝑍0, traditionally taken at the ocean surface in the oceanic case, but this only provides an expression for the specific enthalpy up to an indeterminate function of entropy and salinity, which cannot be specified without additional thermodynamic information about the fluid. This problem arises simply because the knowledge of density alone is in general never sufficient by itself to specify all possible thermodynamic properties of the fluid, which in general also requires some knowledge of such properties as the heat capacity, speed of sound, for example, see Feistel [13] or Callen [11] for exhaustive discussion of this issue. For this reason, it seems preferable to investigate the possibility of deriving the approximate thermodynamic potentials from the knowledge of the Gibbs function 𝑔=𝑔(𝑇,𝑆,𝑃), for instance, which is currently used to synthesise all thermodynamic properties for seawater as discussed in Feistel [13]. Here, we show that the function,π‘”π‘π‘Ž=π‘”π‘π‘Ž(𝑇,𝑆,𝑍)=𝑔𝑇,𝑆,𝑃0(𝑍)ξ€Έ+𝑔0ξ€·π‘βˆ’π‘0ξ€Έ,(35) satisfies the above differential relationships and hence can be regarded as the relevant Boussinesq/Anelastic approximation to the Gibbs function. Indeed, from the result that d𝑔=βˆ’πœ‚d𝑇+πœ‡d𝑆+𝜐d𝑃, it follows that the total differential of π‘”π‘π‘Ž is given bydπ‘”π‘π‘Ž=βˆ’Μƒπœ‚d𝑇+ξ‚πœ‡d𝑆+𝑔0ξ€·1βˆ’πœŒ0Μƒπœξ€Έd𝑍,(36) where we used d𝑃0=βˆ’πœŒ0𝑔0d𝑍, and where the tilde quantities denote the exact quantities in which the actual pressure is replaced by the reference pressure 𝑃0(𝑍), that is, Μƒπœ‚=πœ‚(𝑇,𝑆,𝑃0(𝑍)), ξ‚πœ‡=πœ‡(𝑇,𝑆,𝑃0(𝑍)), and Μƒπœ=𝜐(𝑇,𝑆,𝑃0(𝑍)). This expression is consistent with the previous one, provided that one choosesπœ‚π‘π‘Ž=Μƒπœ‚,πœ‡π‘π‘Ž=ξ‚πœ‡,𝑏=βˆ’π‘”0ξ€·1βˆ’πœŒ0Μƒπœξ€Έ=βˆ’π‘”0ξ‚΅ΜƒπœŒβˆ’πœŒ0ΜƒπœŒξ‚Ά.(37) This in turn shows that the density entering the definition of the buoyancy is actually related to the β€œtrue” density by 𝜌=𝜌0+𝜌0/ΜƒπœŒ(ΜƒπœŒβˆ’πœŒ0)=2𝜌0βˆ’πœŒ20/ΜƒπœŒ.

Assuming that 𝑔 is known, we can deduce the following thermodynamic potentials as follows:π‘’π‘π‘Ž(𝑆,𝑇,𝑍)=π‘”π‘π‘Ž+π‘ξ€·π‘βˆ’π‘0ξ€Έ+π‘‡π‘π‘Žπœ‚π‘π‘Ž=̃𝑔+𝑔0+π‘ξ€Έξ€·π‘βˆ’π‘0ξ€Έβˆ’π‘‡πœ•Μƒπ‘”πœ•π‘‡,β„Žπ‘π‘Ž(𝑆,𝑇,𝑍)=π‘’π‘π‘Žβˆ’π‘ξ€·π‘βˆ’π‘0ξ€Έ=̃𝑔+𝑔0ξ€·π‘βˆ’π‘0ξ€Έβˆ’π‘‡πœ•Μƒπ‘”πœ•π‘‡,(38) which are, respectively, associated with the total differentials:dπ‘’π‘π‘Ž=𝑇dΜƒπœ‚+ξ‚πœ‡d𝑆+ξ€·π‘βˆ’π‘0ξ€Έd𝑏,dβ„Žπ‘π‘Ž=𝑇dΜƒπœ‚+ξ‚πœ‡dπ‘†βˆ’π‘d𝑍.(39) While the above approach is well suited to the case of seawater whose Gibbs function is now well documented, it is arguably less relevant to the case of a fluid whose thermodynamic properties have been synthesised using a different kind of thermodynamic potential. In theory, such a synthesis can indeed be achieved in many ways, for instance, using internal energy, enthalpy, Helmholtz-free energy, or even the more esoteric Massieu functions as discussed by Callen [11]. While it is beyond the scope of this paper to treat all possible cases, it is useful for illustration to discuss the particular case where the specific enthalpy β„Ž=β„Ž(πœ‚,𝑆,𝑃) is the thermodynamic potential used to synthesise the thermodynamic properties of the fluid considered. In that case, it is left to the reader to check that the relevant approximation β„Žπ‘π‘Ž to the specific enthalpy β„Ž is then given by β„Žπ‘π‘Ž=β„Ž(πœ‚,𝑆,𝑃0(𝑍))+𝑔0(π‘βˆ’π‘0), as can be verified from taking the total differential of the latter expression, namely,dβ„Žπ‘π‘Ž=Μƒβ€Œπ‘‡dπœ‚+ξ‚πœ‡d𝑆+𝑔0ξ€·1βˆ’πœŒ0π‘”Μƒπœξ€Έd𝑍=Μƒβ€Œπ‘‡dπœ‚+ξ‚πœ‡dπ‘†βˆ’π‘d𝑍,(40) which is in agreement with (32). As seen previously, π‘’π‘π‘Ž=β„Žπ‘π‘Ž+𝑏(π‘βˆ’π‘0), which implies that the the internal energy can then be written as π‘’π‘π‘Ž=ξ‚β„Ž+(𝑔0+𝑏)(π‘βˆ’π‘0)=ξ‚β„Ž+𝜌0𝑔0(π‘βˆ’π‘0)Μƒπœ.

The above derivations therefore make it clear that it is possible to use any existing knowledge of the thermodynamic properties of the fluid considered to reconstruct the relevant approximations for a BA fluid, and also that the approximations underlying the construction of the BA system of equations only slightly alter the form of the thermodynamic properties, which remain close to that of a fully compressible fluid.

4. Consequence for Understanding the Energetics of Turbulent Stratified Fluids

4.1. Model Formulation

In this section, we discuss the consequences of the above results for our understanding of the energetics of turbulent mixing in stratified fluids, which were originally pointed out by Tailleux [1]. Until recently, the energetics of turbulent-stratified mixing had been essentially discussed for an incompressible Boussinesq fluid with a linear equation of state governed by the following equations:𝐷𝐯𝐷𝑑+1𝜌0βˆ‡π‘ƒ=βˆ’πœŒπœŒ0𝑔0𝐀+πœˆβˆ‡2𝐯,βˆ‡β‹…π―=0,𝐷𝑇𝐷𝑑=πœ…βˆ‡2𝑇,𝜌=𝜌0ξ€Ί1βˆ’π›Όξ€·π‘‡βˆ’π‘‡0ξ€Έξ€»,(41) where 𝛼 is the thermal expansion coefficient, 𝜈 the molecular viscosity, and πœ… the molecular diffusivity, with all remaining notations as before. The situation considered is that of a freely decaying turbulent flow. For instance, a configuration extensively studied is that of a stably stratified parallel shear flow that becomes unstable. Turbulent behaviour ensues until the flow relaminarises, see Tailleux [1] for references. As a result of turbulence, diffusive mixing is enhanced. Here, mixing refers to the irreversible modification of temperature by molecular diffusion. Winters et al. [14] introduced a conceptual framework able to quantify precisely the amount of mixing taking place during a turbulent mixing event.

4.2. The Classical View of the Energetics of Turbulent Mixing

The classical view on the energetics of turbulent mixing has usually revolved on the idea that during a turbulent mixing event, although most of the initial kinetic energy of the parallel shear flow is lost to viscous dissipation as expected, a significant fraction (often reported to be around 20 percent) appears to be lost to gravitational potential energy. This can be mathematically formulated by considering the evolution equation for the kinetic energy, namely,𝜌0𝐷𝐷𝑑𝐯22+βˆ‡β‹…ξ€·π‘ƒπ―βˆ’πœŒ0π…π‘˜π‘’ξ€Έ=βˆ’πœŒπ‘”π‘€βˆ’πœŒ0πœ€πΎ,(42) integrated over the volume of the fluid, as well as over the duration of a turbulent mixing event, which yieldsΞ”KE=βˆ’ξ€œπ‘‰πœŒπ‘”π‘€dπ‘‰βˆ’ξ€œπ‘‰πœŒ0πœ€πΎd𝑉<0,(43) where the overbar denotes the integration over the relevant time interval, and Ξ”KE is the difference between the final and initial kinetic energy. Looking now at the evolution of the gravitational potential energy,𝐷(πœŒπ‘”π‘§)𝐷𝑑=πœŒπ‘”π‘€+π‘”π‘§π·πœŒπ·π‘‘=πœŒπ‘”π‘€+π‘”π‘§πœ…βˆ‡2𝜌,(44) integrating over the volume, as well as in time over the duration of the mixing event, we obtainΞ”GPE=ξ€œπ‘‰πœŒπ‘”π‘€d𝑉+ξ€œπ‘‰π‘”π‘§πœ…βˆ‡2πœŒξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œ.π‘Šlaminar(45) According to the latter equation, the net change in gravitational potential energy Ξ”GPE is made up of two contributions, one in appearance occurring at the expenses of the kinetic energy of the shear flow, and the other to due to internal energy via the term π‘Šlaminar. In this description, the role of internal energy appears to be very limited and restricted to the weak laminar conversion π‘Šlaminar. On the other hand, a somewhat puzzling feature of the above description is that the term πœŒπ‘”π‘€ is supposed to represent in principle a reversible conversion of KE into GPE. Yet, the net change in GPE is irreversible by construction, as the time integration is achieved over a turbulent mixing event separated by laminar evolution. In other words, while reversible conversions of KE into GPE are physically well understood, irreversible conversion of KE into GPE is significantly harder to rationalise.

4.3. Available Potential Energy and Irreversible Mixing

As first recognised by Lorenz [15], only a small fraction of the total gravitational potential energy is generally available for reversible conversion into kinetic energy. This fraction is usually referred to as the β€œavailable potential energy” (APE) and is defined as the difference in potential energy between the actual state and the reference state of minimum potential energy obtained in an adiabatic rearrangement of the fluid parcels. Physically, it can be shown that the density πœŒπ‘Ÿ(𝑧,𝑑) and pressure π‘ƒπ‘Ÿ(𝑧,𝑑) of the reference state depend on 𝑧 and 𝑑 only. Moreover, they must be in hydrostatic equilibrium at all times, namely,πœ•π‘ƒπ‘Ÿπœ•π‘§=βˆ’πœŒπ‘Ÿπ‘”0.(46) Physically, the reference state can be described in terms of the function π‘§π‘Ÿ=π‘§π‘Ÿ(𝐱,𝑑) to represent the position that a parcel located at (𝐱,𝑑) in the actual state would have in the reference state. As a result, the assumption that the parcel is displaced adiabatically from (𝐱,𝑑) to π‘§π‘Ÿ(𝐱,𝑑) implies that πœŒπ‘Ÿ(π‘§π‘Ÿ,𝑑)=𝜌(𝐱,𝑑), that is, the parcel retains its density in the adiabatic displacement. This makes it possible to derive the following expressions for the background gravitational potential energy and available gravitational potential energy:GPEπ‘Ÿ=ξ€œπ‘‰πœŒπ‘Ÿξ€·π‘§π‘Ÿξ€Έπ‘”0π‘§π‘Ÿd𝑉=ξ€œπ‘‰πœŒπ‘”0π‘§π‘Ÿd𝑉,AGPE=ξ€œπ‘‰πœŒπ‘”0ξ€·π‘§βˆ’π‘§π‘Ÿξ€Έd𝑉.(47) Because the reference state is obtained from the actual state by an adiabatic transformation, it follows that it can only change with time if diabatic effects are present in the system. The work by Winters et al. [14] was the first study to suggest that such a property could be used as a rigorous way to diagnose the amount of irreversible mixing taking place in stably stratified turbulent fluids, and suggested to estimate the degree of turbulent mixing taking place in a fluid from the rate of change of GPEπ‘Ÿ. This idea can be put on a rigorous footing by deriving the following evolution equation for GPEπ‘Ÿ:dGPEπ‘Ÿd𝑑=ξ€œπ‘‰πœŒπ‘”0π·π‘§π‘Ÿπ·π‘‘d𝑉+ξ€œπ‘‰π‘”0π‘§π‘Ÿπ·πœŒπ·π‘‘d𝑉=ξ€œπ‘‰π‘”0π‘§π‘Ÿπœ…βˆ‡2𝜌d𝑉.(48) A key result in Winters et al. [14] was to show that the first term in the r.h.s. of (48) always vanishes, by using an explicit construction of the reference density profile. Since then, the result was generalised to a fully compressible fluid by Tailleux [1], while Pauluis [16] showed that the result follows naturally from the reference state being a constrained potential energy minimum. As a result, (48) becomesπ‘Šπ‘Ÿ,mixing=ξ€œπ‘‰π‘”0π‘§π‘Ÿπœ…βˆ‡2𝜌d𝑉=βˆ’ξ€œπ‘‰π‘”0πœ…πœ•πœŒπ‘Ÿπœ•π‘§π‘Ÿβ€–β€–βˆ‡π‘§π‘Ÿβ€–β€–2d𝑉>0,(49) where the last part was obtained by using an integration by parts, assuming the density flux to vanish at all boundaries, and using the result that 𝜌(𝐱,𝑑)=πœŒπ‘Ÿ(π‘§π‘Ÿ(𝐱,𝑑),𝑑) by construction, so that βˆ‡πœŒ=πœ•πœŒπ‘Ÿ/πœ•π‘§π‘Ÿβˆ‡π‘§π‘Ÿ. The result that π‘Šπ‘Ÿ,mixing>0 follows from the fact that πœ•πœŒπ‘Ÿ/πœ•π‘§π‘Ÿβ‰ 0 everywhere, since the reference state is stably stratified by definition. Equation (49) makes it clear that π‘Šπ‘Ÿ,mixing measures the amount of turbulent mixing taking place in the system, because β€–βˆ‡π‘§π‘Ÿβ€–2 measures the degree to which isopycnal surfaces are distorted by turbulent stirring. As a result, we can therefore writedAGPEd𝑑=dξ€·GPEβˆ’GPEπ‘Ÿξ€Έd𝑑=ξ€œπ‘‰πœŒπ‘”0𝑀d𝑉+π‘Šlaminarβˆ’π‘Šπ‘Ÿ,mixing.(50)

4.4. Energy Budget of a Turbulent Mixing Event

Let us now return to the problem of understanding the energetics of a turbulent mixing event, again defined as an episode of intense turbulent mixing characterised by potentially high values of AGPE preceded and followed by laminar conditions for which AGPE=0. Thus, integrating the above equations over the duration of one such turbulent mixing event yieldsΞ”AGPE=0=ξ€œπ‘‰πœŒπ‘”0𝑀d𝑉+π‘Šlaminarβˆ’π‘Šπ‘Ÿ,mixing,(51)Ξ”GPEπ‘Ÿ=π‘Šπ‘Ÿ,mixing.(52) The first important conclusion that one may infer from (51) is that the buoyancy flux βˆ«π‘‰πœŒπ‘”π‘€d𝑉 can be linked to the time-and volume-integrated irreversible diabatic effects due to molecular diffusion since (51) impliesξ€œπ‘‰πœŒπ‘”0𝑀d𝑉=π‘Šπ‘Ÿ,mixingβˆ’π‘Šlaminar=𝐷(APE),(53) where 𝐷(APE) represents the dissipation of AGPE by means of molecular diffusive processes. The second important conclusion inferred from (52) is that the total amount of turbulent mixing having taking place during the turbulent mixing event can be measured from the net change in GPEπ‘Ÿ. Regarding the latter equation, Winters et al. [14] suggested that its physical meaning is best understood by rewriting it as follows:Ξ”GPEπ‘Ÿ=𝐷(APE)+π‘Šlaminar,(54) which suggests that the fraction of kinetic energy that is consumed by turbulent mixing but not dissipated by viscous friction is first converted reversibly into AGPE, but then irreversibly converted somehow into GPEπ‘Ÿ by molecular diffusion, as illustrated in Panel (b) of Figure 1, with an additional small contribution to the net increase of GPEπ‘Ÿ being due to internal energy.

Figure 1: Evolution of ideas about the energetics of turbulent stratified mixing. (a) Interpretation seemingly following the result of classical energetics arguments based on the three reservoirs KE, GPE, and IE; (b) interpretation proposed by Winters et al. [14] based on the use of Lorenz [15] ideas to separate the gravitational potential energy into an available and unavailable component; (c) new interpretation proposed by Tailleux [1] on the basis of a rigorous analysis of the energetics of compressible Navier-Stokes equations aimed at elucidating the role of internal energy.
4.5. Alternative Interpretation of the Energetics of Turbulent Mixing

It is important to remark that the validity of Winters et al. [14] interpretation relies on the idea that the irreversible conversion of one form of mechanical energy (AGPE) into another form of mechanical energy (GPEπ‘Ÿ) is a physically legitimate one. As discussed by Tailleux [1], however, this idea is difficult to accept from the vantage point of thermodynamics, because the irreversible character of such a conversion necessarily implies that it should be accompanied by a change of entropy, and therefore of internal energy. Indeed, assuming the reverse to be true would lead to a contradiction. This is because for no change in internal energy to occur, an equal and opposite amount of compressible work π›Ώπ‘Š=βˆ’π‘ƒd𝜐 would be needed to exactly balance the heat transfer 𝛿𝑄=𝑇dπœ‚irr implied by the assumed diabatic and irreversible character of the AGPE/GPEπ‘Ÿ conversion. However, this would violate the second law of thermodynamics, because it is well known from the theory of heat engines that only a small fraction of the heat transfer 𝛿𝑄 can be transformed into useful work π›Ώπ‘Š, for example, Callen [11], which implies the impossibility for π›Ώπ‘Š and 𝛿𝑄 to be equal in magnitude.

The above difficulty is easily resolved, however, if one accepts to regard the term (π‘βˆ’π‘0)d𝑏 as a conversion between internal energy and gravitational potential energy in the energetics of the BA system. Indeed, note that the term π‘Šπ‘Ÿ,mixing can be rewritten as follows:π‘Šπ‘Ÿ,mixing=ξ€œπ‘‰π‘”0π‘§π‘Ÿπ·πœŒπ·π‘‘d𝑉=βˆ’ξ€œπ‘‰πœŒ0π‘§π‘Ÿπ·π‘π·π‘‘d𝑉,(55) which is precisely of the form associated with a conversion of internal energy into mechanical energy, with π‘§π‘Ÿ instead of π‘βˆ’π‘0. If so, π‘Šπ‘Ÿ,mixing and 𝐷(APE) should therefore be naturally regarded as conversions between internal energy and gravitational potential energy, suggesting that the internal energy should evolve according tod𝐼𝐸d𝑑=𝐷(KE)+𝐷(APE)βˆ’π‘Šπ‘Ÿ,mixing.(56) Although such an interpretation now appears to be consistent with thermodynamics, the form of (56) suggests that it might still be possible for the APE dissipated into internal energy to be then converted into GPEπ‘Ÿ, in which case the whole argument would arguably look like a minor quibble, since it is hard to understand why it should matter whether APE is directly dissipated into GPEπ‘Ÿ or via prior conversion into internal energy.

The latter issue was addressed by Tailleux [1], who argued that internal energy is best regarded as the sum of distinct types of internal energy, such that the kind of internal energy into which APE is dissipated differs from the kind of internal energy converted into GPEπ‘Ÿ, making it possible for (56) to hold without necessarily implying that the APE dissipated into internal energy is converted into GPEπ‘Ÿ. In particular, Tailleux [1] argue that in the same way that it is possible to decompose the temperature as follows:𝑇(𝐱,𝑑)=𝑇0(𝑑)+π‘‡π‘Ÿ(𝑧,𝑑)+π‘‡ξ…ž(𝐱,𝑑),(57) that is, as the sum of a time-dependent mean temperature, a part accounting for the mean vertical gradient, and a residual, it is similarly possible to decompose the internal energy into three reservoirs:IE=IE0+IEexergy+AIE,(58) where IE0 is the internal energy of the thermodynamic equilibrium obtained from Lorenz [15] state in a process conserving the total enthalpy of the system, IEexergy=IEπ‘Ÿβˆ’IE0 is the difference in internal energy between the internal energy of Lorenz [15] reference state, and AIE=IEβˆ’IEπ‘Ÿ is the available internal energy, that is, the difference in internal energy between the actual state and Lorenz [15] reference state. Tailleux [1] showed how to derive evolution equations for each subcomponents of the internal energy and was able to show that at leading order, neglecting only very small terms, it is possible to decompose (56) into the two following equations for IE0 and IEexergy as follows:dIE0dπ‘‘β‰ˆπ·(KE)+𝐷(APE),(59)dIEexergydπ‘‘β‰ˆβˆ’π‘Šπ‘Ÿ,mixing.(60) In physical terms, Tailleux [1] shows that increasing or decreasing IE0 is equivalent to increasing or decreasing the mean thermodynamic equilibrium temperature 𝑇0, whereas decreasing IEexergy is equivalent to smooth out the vertical temperature gradient of π‘‡π‘Ÿ(𝑧,𝑑). As a result, (59) interpretation is that both the viscous dissipation of KE and the diffusive dissipation of APE increase 𝑇0, whereas (60) interpretation is that turbulent mixing smoothes out the vertical gradient of π‘‡π‘Ÿ(𝑧,𝑑). This is an important result, because both π‘Šπ‘Ÿ,mixing and 𝐷(APE) are widely recognised to represent leading order energy conversions in turbulent-stratified mixing. In the current dominant paradigm based on Winters et al. [14] interpretation, however, these conversions have been rationalised without invoking any significant role for the internal energy, in contrast to our interpretation. Physically, our interpretation also implies that whereas the dissipation of KE and APE should not be associated with any significant overall expansion of the fluid owing to its assumed weakly compressible nature, the release of exergy IEexergy should in contrast be comparatively associated with a significant reduction in the volume of the fluid. By comparison, Winters et al. [14] interpretation call only for a tiny volume expansion due to the viscous dissipation, and for a tiny volume contraction due to the conversion of internal energy into GPEπ‘Ÿ at the laminar rate π‘Šlaminar. Figure 2 illustrates a laboratory experiment that clearly tends to support our interpretation.

Figure 2: Experimental setup of an experiment showed to the author by Peter Rhines during a visit to his GFD Laboratory in Seattle in June 2008, suggesting that compressibility effects increase with the degree of turbulent in stably stratified fluids. An Erlenmeyer flask is filled up with water at room temperature. A vertical stratification is then created by bringing the upper part of the water near boiling point, and a magnetic stirrer is dropped in the fluid. The top of the flask is then sealed up with a cork through which a thin glass tube is inserted in order to magnify the variations of the air-water interface. In laminar conditions (a), the fluid is near resting conditions, and the air-water interface moves downward very slowly as the result of molecular diffusion and/or radiative cooling. Activating the magnetic stirrer generates stratified turbulence (b) that results in the spectacular drop of the air-water interface, most likely as the result of the well-known contraction upon mixing effect that is due to the temperature dependence of the thermal expansion coefficient.

5. Net Compressible Work in Purely Buoyancy-Driven Flows

Having determined that the classical Boussinesq and anelastic approximations formally possess a representation of the coupling between the dynamics and thermodynamics that is comparable to that of the fully compressible Navier-Stokes equations, we now seek to illustrate some aspects of such coupling in the particular case of a purely buoyancy-driven circulation driven by differential heating/cooling imposed at the top boundary of the domain, as in the oceanic case, a configuration generally referred to as β€œhorizontal convection,” see Hughes and Griffiths [3] for a review. The importance of the dynamics/thermodynamics coupling, and of the compressible work, comes from that in such a context the mechanical energy balance reduces to a balance between the compressible work and the total viscous dissipation, namely, ξ€œπ‘‰π‘ƒπ·πœπ·π‘‘dπ‘š=ξ€œπ‘‰πœ€πΎdπ‘š,(61) where dπ‘š=𝜌d𝑉 is the mass of an elementary fluid element. The central question here is to what extent is it possible to quantify the compressible work in (61), and whether the resulting estimates are sensitively affected by using the Boussinesq/Anelastic approximation as compared with using the fully compressible Navier-Stokes equations.

5.1. Antiturbulence Theorem and πœ€-Theorems

The above issue was first addressed for a Boussinesq fluid with a linear equation of state governed by (41), as part of the derivation of the so-called antiturbulence theorem by Paparella and Young [17], whose particular aim was to derive an upper bound for the total viscous dissipation βˆ«π‘‰πœ€πΎdπ‘š. Paparella and Young [17] derived their result for a two-dimensional rectangular ocean forced at its top by a Dirichlet-like surface temperature boundary condition, with insulating boundary conditions everywhere else. For the particular system of equations considered, the local balance equations for the kinetic energy and gravitational potential energy are then, respectively, given by𝜌0𝐷𝐷𝑑𝐯22+βˆ‡β‹…ξ‚΅π‘ƒπ―βˆ’πœˆβˆ‡π―22ξ‚Ά=βˆ’πœŒπ‘”0π‘€βˆ’πœŒ0πœ€πΎ,π·ξ€·πœŒπ‘”0𝑧𝐷𝑑=πœŒπ‘”0𝑀+𝑔0π‘§π·πœŒπ·π‘‘.(62) For a steady state, summing these two equations, and integrating over the whole domain yieldsξ€œπ‘‰g0π‘§π·πœŒπ·π‘‘d𝑉=ξ€œπ‘‰π‘”0π‘§πœ…βˆ‡2𝜌dπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ‘Šlaminar=ξ€œπ‘‰πœŒ0πœ€πΎd𝑉.(63) The diffusive term π‘Šlaminar has been previously encountered. As a result, such a balance can be written:π‘Šlaminar=πœ…π‘”0π΄ξ€ΊβŸ¨πœŒβŸ©bottomβˆ’βŸ¨πœŒβŸ©topξ€»=ξ€œπ‘‰πœŒ0πœ€πΎ.(64) Paparella and Young [17] showed that for the boundary conditions considered, the maximum principle imposes to the quantity ⟨𝜌⟩bottomβˆ’βŸ¨πœŒβŸ©top to be bounded from above by a constant independent of the molecular diffusivity πœ… and viscosity 𝜈. This result implies that βˆ«π‘‰πœ€πΎdπ‘š must converge toward zero in the limit of zero diffusivity by holding the Prandtl number π‘ƒπ‘Ÿ=𝜈/πœ… constant, which is the so-called antiturbulence theorem. For the present purposes, the result is of interest, because it provides an explicit expression of the work of expansion/contraction in terms of quantities that are all observable or measurable. Thus, Wang and Huang [18] used typical oceanic values to estimate that 𝐡=𝑂(15GW) based on (64), which is about two orders of magnitude smaller than the power input due to the wind, generally estimated to be 𝑂(1TW). This result has been very influential, because it has generally been concluded that the observed amount of viscous dissipation in the oceans can only be accounted for by the mechanical stirring due to the wind and tides. It has also been regarded as implying that surface buoyancy fluxes can only drive a weak overturning circulation. The latter inference, however, requires that it is possible to link in some way the volume-integrated viscous dissipation, which only measures the amount of kinetic energy present at the dissipation scale, to the amount of kinetic energy of the large-scale overturning circulation. As far as we know, such a link is yet to be established rigorously.

5.2. Generalisation to a Nonlinear Equation of State

The equation of state of seawater is strongly nonlinear, however, so that it is unclear whether (64) remains valid for the more realistic oceanic case. The following provides a brief overview of various attempts at generalising (64) to a nonlinear equation of state typical of that for seawater.

5.2.1. Unsuitability of the Direct Approach

The most natural approach to estimating the volume-integrated work of expansion/contraction for a compressible ocean with a nonlinear equation of state is perhaps to regard the specific volume as a function of entropy and pressure (we discard salinity temporarily for simplicity), which allows one to expand the rate of expansion/contraction as follows:π·πœπ·π‘‘=πœ•πœπœ•πœ‚π·πœ‚π·π‘‘+πœ•πœπœ•π‘ƒπ·π‘ƒπ·π‘‘=π›Όπ‘‡πœŒπ‘π‘ξ€Ίβˆ‡β‹…ξ€·πœ…πœŒπ‘π‘βˆ‡π‘‡ξ€Έ+πœŒπœ€πΎξ€»πœŒπ‘‡βˆ’1𝜌2𝑐2𝑠𝐷𝑃𝐷𝑑,(65) and therefore the volume-integrated work of expansion/contraction as follows:𝐡=ξ€œπ‘‰π›Όπ‘ƒπœŒπ‘π‘βˆ‡β‹…(πœ…πœŒπ‘π‘βˆ‡π‘‡)dπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ΅1+ξ€œπ‘‰π›Όπ‘ƒπœŒπ‘π‘πœŒπœ€πΎdπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ΅2βˆ’ξ€œπ‘‰π‘ƒπœŒπ‘2𝑠𝐷𝑃𝐷𝑑dπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ΅3,(66) where the expression was obtained by using the results that πœ•πœ/πœ•πœ‚=𝛼𝑇/(πœŒπ‘π‘)=Ξ“ is the so-called adiabatic lapse rate, and πœ•πœ/πœ•π‘ƒ=πœ•(1/𝜌)/πœ•π‘ƒ=βˆ’(1/𝜌2)πœ•πœŒ/πœ•π‘ƒ=βˆ’1/(𝜌2𝑐2𝑠), where 𝑐𝑠 is the speed of sound. We also assumed the evolution equation for entropy to be given byπ·πœ‚π·π‘‘=βˆ‡β‹…ξ€·πœŒπœ…π‘π‘βˆ‡π‘‡ξ€Έ+πœŒπœ€πΎπœŒπ‘‡,(67) which assumes that the only terms contributing to the diabatic heating are the molecular diffusion of heat and the viscous dissipation. Equation (66) regards 𝐡 as being made up of essentially three contributions, where 𝐡1 is the net effect due to molecular diffusive heating/cooling, 𝐡2 is the net effect due to viscous dissipation, and 𝐡3 represents the net effect due to the adiabatic expansion/contraction. The two first terms 𝐡1 and 𝐡2 can be seen to involve the key dimensionless quantity Ξ₯=𝛼𝑃/(πœŒπ‘π‘), which plays the role of a thermodynamic efficiency-like quantity. In the oceans, typical values are 𝛼=10βˆ’4Kβˆ’1, π‘ƒβ‰ˆ4000dbar=4.108Pa (near the ocean bottom), πœŒβ‰ˆ103kgΒ·mβˆ’3, π‘π‘β‰ˆ4.103JΒ·kgβˆ’1Β·Kβˆ’1, leading to Ξ₯β‰ˆΓ—10βˆ’2β‰ͺ1. Since 𝐡 must balance the overall viscous dissipation βˆ«π‘‰πœŒπœ€πΎd𝑉 in a steady-state, it follows immediately that the term 𝐡2 must be negligible compared to the other two terms. Obviously, the term 𝐡1 is the only one that can be compared with Paparella and Young [17] expression, since it is the only one depending on molecular diffusion. By integrating by parts, this term can be rewritten:𝐡1=ξ€œπœ•π‘‰Ξ₯πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…π§dπ‘†βˆ’ξ€œπ‘‰πœ•Ξ₯πœ•π‘‡πœ…πœŒπ‘π‘β€–βˆ‡π‘‡β€–2dπ‘‰βˆ’ξ€œπ‘‰πœ•Ξ₯πœ•π‘ƒπœ…πœŒπ‘π‘βˆ‡π‘‡β‹…βˆ‡π‘ƒd𝑉.(68) Note here that Paparella and Young [17] result can be recovered by assuming 𝜌, 𝑐𝑝, πœ… to be constant, and the pressure to be approximated by the hydrostatic pressure π‘ƒβ‰ˆβˆ’πœŒ0𝑔𝑧, as in that case, πœ•Ξ₯/πœ•π‘‡=0, while πœ•Ξ₯/πœ•π‘ƒ=𝛼/(πœŒπ‘π‘), and the boundary term vanishes owing to the assumption that βˆ«πœ•π‘‰πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…π§d𝑆=0 required by the steady-state assumption. Ultimately, however, such a direct approach appears to be impractical, because it is unclear how to estimate the term 𝐡3.

5.2.2. Approach Based on the Budget of Total Available Energy

A significant breakthrough in understanding how to get around the difficulty posed by the term 𝐡3 above is due to McIntyre [19] and Nycander [9], which both showed that it could be circumvented by focusing on the budget of the total available energy 𝐯2/2+β„Žβ€‘ discussed above. The following briefly summarises the key steps of their approach, which were derived for a Boussinesq ocean with a nonlinear equation of state. Here, we extend their discussion to the more general BA system discussed previously. To that end, let us first recall that in the BA system, the local evolution equation for the kinetic energy takes the form:𝜌0𝐷𝐷𝑑𝐯22+βˆ‡β‹…ξ€Ίπ›Ώπ‘ƒπ―βˆ’πœŒ0π…π‘˜π‘’ξ€»=𝜌0π‘π·π‘π·π‘‘βˆ’πœŒ0πœ€πΎ.(69) In their derivation, both McIntyre [19] and Nycander [9] assumed the buoyancy 𝑏=𝑏(S,Θ,𝑍) to be a given function of salinity 𝑆, conservative temperature Θ, and geopotential height 𝑍. Moreover, Θ and 𝑆 were assumed to satisfy evolution equations of the form:𝜌0π·Ξ˜π·π‘‘=βˆ’βˆ‡β‹…ξ€·πœŒ0π…Ξ˜ξ€Έ+𝜌0Μ‡β€ŒΞ˜irr,𝜌0𝐷𝑆𝐷𝑑=βˆ’βˆ‡β‹…ξ€·πœŒ0𝐅𝑆,(70) with π…Ξ˜ and 𝐅𝑆 being linear functions of the molecular diffusive fluxes of heat and salt, whose forms are traditionally constrained by the second law of thermodynamics, as discussed in de Groot and Mazur [20], for instance. The specific form of these fluxes is not important for what follows, and is therefore not further detailed. The term Μ‡β€ŒΞ˜irr represents the nonconservative production/destruction of the conservative temperature, which McDougall [21] argued is small enough to be regarded as negligible. As a result, this term was neglected by both McIntyre [19] and Nycander [9]. As discussed by Tailleux [22], however, this introduces some inconsistency in the discussion, as the resulting equations are no longer energetically consistent. This is a minor difficulty, however, which could easily be corrected by retaining Μ‡β€ŒΞ˜irr. To proceed, McIntyre [19] and Nycander [9] introduced the following function:β„Žβ€‘(𝑆,Θ,𝑍)=βˆ’ξ€œπ‘π‘0𝑏𝑆,Θ,π‘ξ…žξ‚dπ‘ξ…ž,(71) which was coined β€œdynamic enthalpy” by Young [8]. Differentiating with respect to time yieldsπ·β„Žβ€‘π·π‘‘=βˆ’π‘π·π‘π·π‘‘+𝐢𝑆𝐷𝑆𝐷𝑑+πΆΞ˜π·Ξ˜π·π‘‘,(72) where𝐢𝑆=βˆ’ξ€œπ‘π‘0πœ•π‘πœ•π‘†ξ€·π‘†,Θ,π‘ξ…žξ€Έdπ‘ξ…ž,𝐢Θ=βˆ’ξ€œπ‘π‘0πœ•π‘πœ•Ξ˜ξ€·π‘†,Θ,π‘ξ…žξ€Έdπ‘ξ…ž.(73) As a result, summing (69) and (72) yields𝜌0𝐷𝐷𝑑𝐯22+β„Žβ€‘ξ‚Ά+βˆ‡β‹…ξ€·π›Ώπ‘ƒπ―βˆ’πœŒ0π…π‘˜π‘’ξ€Έ=𝜌0𝐢𝑆𝐷𝑆𝐷𝑑+πΆΞ˜π·Ξ˜π·π‘‘βˆ’πœ€πΎξ‚.(74) The final step consists in integrating the latter equation over the whole ocean domain. In a steady-state, this leads to the following budget:ξ€œπ‘‰πœŒ0πœ€πΎd𝑉=𝐡=βˆ’ξ€œπ‘‰ξ€·πΆπ‘†βˆ‡β‹…ξ€·πœŒ0𝐅𝑆+πΆΞ˜βˆ‡β‹…ξ€·πœŒ0π…Ξ˜ξ€Έξ€Έd𝑉=βˆ’ξ€œπœ•π‘‰πœŒ0(𝐢𝑆𝐅𝑆+CΞ˜π…Ξ˜)⋅𝐧dπ‘†ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œ=0+ξ€œπ‘‰πœŒ0ξ€·π…π‘†β‹…βˆ‡πΆπ‘†+π…Ξ˜β‹…βˆ‡πΆΞ˜ξ€Έd𝑉,(75) where the first term in the second line of (75) relates to the surface fluxes of π…Ξ˜ and 𝐅𝑆, which vanishes because by construction, 𝐢Θ=𝐢𝑆=0 at the ocean surface where 𝑍=𝑍0. In McIntyre [19] and Nycander [9], the focus is on deriving upper bounds for βˆ«π‘‰πœŒ0πœ€πΎd𝑉, but clearly the method also provides a way to discuss the magnitude of 𝐡, since the two quantities must be equal in a steady state. The important point to be noted here is that (75) provides an expression for 𝐡 that only involves the molecular diffusive fluxes, with no term proportional to 𝐷𝑃/𝐷𝑑, similarly as in Paparella and Young [17] study.

5.2.3. Approach Based on the Budget of Potential Enthalpy

In this section, we aim to further clarify the above issues, by showing how the above results can be extended to a fully compressible ocean by considering the budget of potential enthalpy. To simplify the discussion, the effects of salinity are discarded. The potential enthalpy, introduced in the oceanographic context by McDougall [21], is defined as the enthalpy that a fluid parcel would have if brought adiabatically to the reference pressure 𝑃0, usually taken as the standard atmospheric pressure. As a result, the potential enthalpy, denoted by β„Žπœƒ thereafter, can be regarded as being implicitly defined by β„Žπœƒ(πœ‚,𝑃;𝑃0)=β„Ž(πœ‚,𝑃0), with β„Ž being the specific enthalpy. From the definition of its total differential dβ„Ž=𝑇dπœ‚+𝜐d𝑃, it follows thatπ·β„Žπœƒπ·π‘‘=πœƒπ·πœ‚π·π‘‘=πœƒπ‘‡ξ€Ίβˆ‡β‹…ξ€·πœ…πœŒπ‘π‘βˆ‡π‘‡ξ€Έ+πœŒπœ€πΎξ€»,(76) where πœƒ is the potential temperature, that is, the temperature that a parcel would have if brought adiabatically to the ocean surface, and hence defined as the implicit solution of πœ‚(𝑇,𝑃)=πœ‚(πœƒ,P0) or as the solution of πœ‚(𝑇,𝑆,𝑃)=πœ‚(πœƒ,𝑆,𝑃0) if salinity was to be retained. Now, if we take the overall budget in a steady state, we have the result:ξ€œπœ•π‘‰πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…π§dπ‘†ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œ=0βˆ’ξ€œπ‘‰πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…βˆ‡ξ‚€πœƒπ‘‡ξ‚dπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ŒΜ‡β€Œβ„Žirr+ξ€œπ‘‰πœƒπ‘‡πœŒπœ€πΎd𝑉=0,(77) which was obtained by noting that πœƒ=𝑇 at the reference pressure 𝑃=𝑃0. Here, it is useful to introduce the quantity:𝑠=βˆ«π‘‰((π‘‡βˆ’πœƒ)/𝑇)πœŒπœ€πΎdπ‘‰βˆ«π‘‰πœŒπœ€πΎd𝑉.(78) Physically, πœƒ never differs more than of a few degrees Celsius from the in-situ temperature, which implies 𝑠β‰ͺ1. As a result, it follows that (77) implies βˆ«π‘‰πœŒπœ€πΎd𝑉=Μ‡β€Œβ„Žirr/(1βˆ’π‘ )β‰ˆΜ‡β€Œβ„Žirr, which in turn impliesπ΅β‰ˆΜ‡β€Œβ„Žirr=ξ€œπ‘‰πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…βˆ‡ξ‚€πœƒπ‘‡ξ‚d𝑉.(79) To make progress, note that from the definition of potential temperature πœ‚(𝑇,𝑃)=πœ‚(πœƒ,𝑃0), we may regard πœƒ as a function of temperature and pressure. In particular, using the result that dπœ‚=(𝑐𝑝/𝑇)dπ‘‡βˆ’(𝛼/𝜌)d𝑃, it follows that𝑐𝑝𝑇dπ‘‡βˆ’π›ΌπœŒd𝑃=π‘π‘ξ€·πœƒ,𝑃0ξ€Έπœƒdπœƒ,(80) which impliesπœ•πœƒπœ•π‘‡=π‘π‘πœƒπ‘π‘Ÿπ‘π‘‡,πœ•πœƒπœ•π‘ƒ=βˆ’π›ΌπœƒπœŒπ‘π‘Ÿπ‘,(81) where for simplicity, we defined 𝑐𝑝(πœƒ,𝑃0)=π‘π‘Ÿπ‘. Now, using (81) leads toβˆ‡ξ‚€πœƒπ‘‡ξ‚=π‘‡βˆ‡πœƒβˆ’πœƒβˆ‡π‘‡π‘‡2=1π‘‡ξ‚ƒπœ•πœƒπœ•π‘‡βˆ‡π‘‡+πœ•πœƒπœ•π‘ƒβˆ‡π‘ƒξ‚„βˆ’πœƒβˆ‡π‘‡π‘‡2=1π‘‡ξ‚Έπ‘π‘πœƒπ‘π‘Ÿπ‘π‘‡βˆ‡π‘‡βˆ’π›ΌπœƒπœŒπ‘π‘Ÿπ‘βˆ‡π‘ƒξ‚Ήβˆ’πœƒβˆ‡π‘‡π‘‡2=ξ‚΅π‘π‘βˆ’π‘π‘Ÿπ‘π‘π‘Ÿπ‘ξ‚Άπœƒβˆ‡π‘‡π‘‡2βˆ’π›ΌπœƒπœŒπ‘π‘Ÿπ‘π‘‡βˆ‡π‘ƒ,(82) which in turn allows one to express Μ‡β€Œβ„Žirr as follows:Μ‡β€Œβ„Žirr=ξ€œπ‘‰πœƒξ‚΅π‘π‘βˆ’π‘π‘Ÿπ‘π‘π‘Ÿπ‘ξ‚Άπœ…πœŒπ‘π‘β€–βˆ‡π‘‡β€–2𝑇2dπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ»1βˆ’ξ€œπ‘‰π›Όπ‘π‘πœƒπ‘π‘Ÿπ‘π‘‡πœ…βˆ‡π‘ƒβ‹…βˆ‡π‘‡dπ‘‰ξ„Ώξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…ƒξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…€ξ…Œπ»2.(83) Equation (83) is a key result, as it currently represents the most general extension of Paparella and Young [17] result that is valid for a compressible ocean with a nonlinear equation of state. Although salinity was not considered, the result could be easily extended to include it as well. As in Paparella and Young [17], McIntyre [19], and Nycander [9], the expression obtained only involves terms related to molecular diffusive processes. It is composed of mainly two terms 𝐻1 and 𝐻2. The first term involves the quantity 𝑅=πœƒ(π‘π‘βˆ’π‘π‘Ÿπ‘)/π‘π‘Ÿπ‘, which can be shown to be in general negative for most conditions encountered in the oceans. It is zero at the surface, and increase with pressure. By invoking the intermediate value theorem, it follows that there exists a value π‘…βˆ— so that 𝐻1 can be written as follows:𝐻1=π‘…βˆ—ξ€œπ‘‰πœ…πœŒπ‘π‘β€–βˆ‡π‘‡β€–2𝑇2d𝑉<0.(84) Interestingly, under this form, 𝐻1 appears as being proportional to the irreversible entropy production by molecular diffusive heat fluxes. Since 𝐻1 is negative, 𝐻2 must also be negative in order for 𝐡 to positive in order to balance the overall viscous dissipation. Some insight into 𝐻2 can be obtained by rewriting it as follows:𝐻2=ξ€œπ‘‰ξ‚΅π‘π‘πœƒπ‘π‘Ÿπ‘π‘‡ξ‚Άπ›Όπœ…βˆ‡π‘ƒβ‹…βˆ‡π‘‡d𝑉=ξ€œπ‘‰π‘π‘πœƒπ‘π‘Ÿπ‘π‘‡π›Όπœ…βˆ‡β„Žπ‘ƒβ‹…βˆ‡β„Žπ‘‡d𝑉+ξ€œπ‘‰π‘π‘πœƒπ‘π‘Ÿπ‘π‘‡π›Όπœ…πœ•π‘ƒπœ•π‘§πœ•π‘‡πœ•π‘§d𝑉.(85) Interestingly, it is the last term in (85) that allows one to make the connection with Paparella and Young [17] result, as seen by the following manipulation:ξ€œπ‘‰π‘π‘πœƒπ‘π‘Ÿπ‘π‘‡π›Όπœ…πœ•π‘ƒπœ•π‘§πœ•π‘‡πœ•π‘§dπ‘‰β‰ˆβˆ’ξ€œπ‘‰π›ΌπœŒπ‘”0πœ…πœ•π‘‡πœ•π‘§d𝑉=ξ€œπ‘‰πœ…π‘”πœ•πœŒπœ•π‘§d𝑉=βˆ’π‘Šlaminar(86) by approximating π‘π‘πœƒ/(π‘π‘Ÿπ‘π‘‡)β‰ˆ1 and using the hydrostatic approximation πœ•π‘ƒ/πœ•π‘§β‰ˆβˆ’πœŒπ‘”0. It is more difficult, however, to quantity the correlation between the horizontal pressure and temperature gradients βˆ‡β„Žπ‘ƒβ‹…βˆ‡β„Žπ‘‡; presently, even the sign of the later, let alone its magnitude, is challenging to estimate, because the correlation between the two gradients is likely to occur at the molecular diffusive scales. Somehow, understanding how to quantify this term is equivalent to understanding the characteristics of the turbulent pressure and temperature fluctuations at such small diffusive scales. More research is needed to clarify this issue.

5.3. General Remarks on Deriving Expressions for 𝐡

Reflecting in hindsight on the different ideas that have been developed over the past decade to generalise Paparella and Young [17] result to a fully compressible ocean with a nonlinear equation of state, it appears that the most general approach to estimating the viscous dissipation and hence the overall work of expansion/contraction is ultimately rooted in the steady-state budget for some function 𝐹(πœ‚) of the specific entropy. Indeed, if we consider the evolution equation for such a function, it is given by𝜌𝐷𝐹(πœ‚)𝐷𝑑=πœŒπΉξ…ž(πœ‚)π·πœ‚π·π‘‘=πΉξ…ž(πœ‚)βˆ‡β‹…ξ€·πœ…πœŒπ‘π‘βˆ‡π‘‡ξ€Έ+πœŒπœ€πΎπ‘‡.(87) Integrating this equation over the ocean volume and assuming a steady state thus yieldξ€œπ‘‰πΉξ…ž(πœ‚)βˆ‡β‹…ξ€·πœ…πœŒπ‘π‘βˆ‡π‘‡ξ€Έπ‘‡d𝑉+ξ€œπ‘‰πΉξ…ž(πœ‚)πœŒπœ€πΎπ‘‡d𝑉=0.(88) Now, if 𝐹(πœ‚) is chosen as an increasing function of πœ‚ so that πΉξ…ž(πœ‚)>0, then the theorem of intermediate values implies that there must exist some value πœ‚βˆ— and π‘‡βˆ— allowing to rewrite the above formula as follows:ξ€œπ‘‰πœŒπœ€πΎd𝑉=βˆ’π‘‡βˆ—πΉξ…ž(πœ‚βˆ—)ξ€œπ‘‰πΉξ…ž(πœ‚)βˆ‡β‹…ξ€·πœ…πœŒπ‘π‘βˆ‡π‘‡ξ€Έπ‘‡d𝑉=βˆ’π‘‡βˆ—πΉξ…ž(πœ‚βˆ—)βŽ§βŽ¨βŽ©ξ€œπœ•π‘‰πΉξ…ž(πœ‚)πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…π§π‘‡dπ‘†βˆ’ξ€œπ‘‰πœ…πœŒπ‘π‘βˆ‡π‘‡β‹…βˆ‡βŽ›βŽœβŽπΉξ…ž(πœ‚)π‘‡βŽžβŽŸβŽ dπ‘‰βŽ«βŽ¬βŽ­.(89) Clearly, this method appears to systematically provide an expression for the viscous dissipation as the sum of one term involving the surface heat fluxes and one term linked to the irreversible molecular diffusive fluxes. Clearly, the method based on potential enthalpy derived above is based on the particular choice 𝐹(πœ‚)=β„Žπœƒ(πœ‚,𝑃;𝑃0) and has the particularity that the term involving the surface fluxes identically vanish. Two different particular choices corresponding to 𝐹(πœ‚)=πœ‚ and to 𝐹(πœ‚)=β„Ž(πœ‚,π‘ƒπ‘Ÿ(πœ‚)) for some reference pressure π‘ƒπ‘Ÿ(πœ‚) related to Lorenz [15] reference profile were considered by Tailleux [2] as part of a discussion on the nature of the work of expansion/contraction. The current challenge consists in identifying the most sensible choices of 𝐹(πœ‚), and to determine what kind of information can be extracted from making such choices. To that end, however, the main difficulties is in understanding how to quantify the kind of correlation between the pressure and temperature gradients that enter such expressions, which requires understanding the nature of such gradients at the molecular diffusive scales that are likely to control the overall value of such a correlation.

6. Summary and Conclusions

For over a century, the coupling between the dynamics and thermodynamics has been overwhelmingly regarded as being of little or no relevance to the understanding of stratified turbulence in weakly compressible fluids such as water or seawater. As a result, the study of such fluids has been nearly systematically been carried out in the context of dynamical approximations, such as the Boussinesq approximation or anelastic approximation, that decouples the dynamics and thermodynamics at leading order. It is not necessarily fully realized, however, that such a decoupling is achieved only in absence of diabatic effects and for a linear equation of state. Indeed, coupling is unavoidably reintroduced whenever any one of the latter effects is retained. If so, a key issue is what can we say about the nature of such a coupling in that case? How does it compare with that of the fully compressible Navier-Stokes equations?

In order to address these questions, the present paper first showed that the particular anelastic approximation derived by Pauluis [6], which admits the Boussinesq approximation as a particular case, can be endowed with a fully consistent energetics and thermodynamics, even for a binary fluid with an arbitrary nonlinear equation of state. As showed in more details by Tailleux [10], it is not only possible to show that the energetics of such a system possesses a term that can be directly interpreted as the approximation to the work of expansion/contraction, it is also possible to construct explicitly the full range of thermodynamic potentials for such a fluid. This is an important result, because it permits for the first time to clarify the role played by internal energy in the Boussinesq/Anelastic system, which has been previously a large source of confusion, owing to the β€œincompressible” label attached to the Boussinesq approximation, or to the notion that such approximations filter out sound waves, which is often interpreted as the lack of a coupling with internal energy. In particular, the clarification allows one to confirm the previous results by Tailleux [1], who argued that the dissipation of available potential energy, previously interpreted by Winters et al. [14] as an irreversible conversion of APE into background gravitational potential energy, should in fact be regarded as an irreversible conversion of APE into internal energy, in the same way that viscous dissipation is irreversible dissipation of mechanical energy into internal energy. Finally, we also determined how to extend the derivation of so-called πœ€-theorems for horizontal convection to the case of a fully compressible fluid with an arbitrary nonlinear equation of state.

An important feature of the energetics of the BA system is that there is a priori no physical mechanism to prevent the conversion term between internal energy and mechanical energy to become potentially very large. In other words, even though the BA system is constructed under the premises that compressibility effects are small, it is now clear that the BA compressible effects will in general increase with the strength of the stratification and turbulence. In principle, the magnitude of the compressibility effects can be diagnosed a posteriori by estimating the divergent velocity field 𝐯𝑑 from solving the equation βˆ‡β‹…π―π‘‘=βˆ’πœŒβˆ’1𝐷𝜌/𝐷𝑑. For instance, in the case of the Boussinesq model with a diffusive model for the density, this would consist in solving βˆ‡β‹…π―π‘‘=βˆ’πœ…βˆ‡2𝜌/𝜌, which makes it clear that 𝐯𝑑 can become large if the local values of diffusion πœ…βˆ‡2𝜌 become large, which is expected to be the case as the degree of turbulence and strength of the stratification increase. Further work is needed to determine whether it may become possible for 𝐯𝑑 to become sufficiently large as to become locally important in the advection of momentum and tracers. If so, this would also potentially suggest the importance of sound waves. As is well known, acoustic waves can induce mean flows, as in the case of acoustic streaming. Whether such effects may be identified in stratified turbulence opens up an exciting new terrain of investigation.


  1. R. Tailleux, β€œOn the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy,” Journal of Fluid Mechanics, vol. 638, pp. 339–382, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus
  2. R. Tailleux, β€œEntropy versus APE production: on the buoyancy power input in the oceans energy cycle,” Geophysical Research Letters, vol. 37, no. 22, Article ID L22603, 2010. View at Publisher Β· View at Google Scholar
  3. G. O. Hughes and R. W. Griffiths, β€œHorizontal convection,” Annual Review of Fluid Mechanics, vol. 40, pp. 185–208, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus
  4. J. Boussinesq, Théorie Analytique de la Chaleur, vol. 2, Gauthiers-Villars, Paris, France, 1903.
  5. Y. Ogura and N. A. Phillips, β€œScale analysis and shallow convection in the atmopshere,” Journal of the Atmospheric Sciences, vol. 19, pp. 173–179, 1962. View at Google Scholar
  6. O. Pauluis, β€œThermodynamic consistency of the anelastic approximation for a moist atmosphere,” Journal of the Atmospheric Sciences, vol. 65, no. 8, pp. 2719–2729, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus
  7. A. P. Ingersoll, β€œBoussinesq and anelastic approximations revisited: Potential energy release during thermobaric instability,” Journal of Physical Oceanography, vol. 35, no. 8, pp. 1359–1369, 2005. View at Publisher Β· View at Google Scholar Β· View at MathSciNet Β· View at Scopus
  8. W. R. Young, β€œDynamic enthalpy, conservative temperature, and the seawater Boussinesq approximation,” Journal of Physical Oceanography, vol. 40, no. 2, pp. 394–400, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus
  9. J. Nycander, β€œHorizontal convection with a non-linear equation of state: generalization of a theorem of Paparella and Young,” Tellus, vol. 62, no. 2, pp. 134–137, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus
  10. R. Tailleux, β€œThermodynamics/dynamics coupling and thermodynamic consistency of Boussinesq and anelastic binary fluids with an arbitrary nonlinear equation of state,”
  11. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, Wiley, 1985.
  12. R. A. Alberty, β€œUse of legendre transforms in chemical thermodynamics,” Pure and Applied Chemistry, vol. 73, no. 8, pp. 1349–1380, 2001. View at Google Scholar Β· View at Scopus
  13. R. Feistel, β€œA new extended Gibbs thermodynamic potential of seawater,” Progress in Oceanography, vol. 58, no. 1, pp. 43–114, 2003. View at Publisher Β· View at Google Scholar Β· View at Scopus
  14. K. B. Winters, P. N. Lombard, J. J. Riley, and E. A. D'Asaro, β€œAvailable potential energy and mixing in density-stratified fluids,” Journal of Fluid Mechanics, vol. 289, pp. 115–128, 1995. View at Google Scholar Β· View at Scopus
  15. E. N. Lorenz, β€œAvailable potential energy and the maintenance of the general circulation,” Tellus, vol. 7, pp. 157–167, 1955. View at Google Scholar
  16. O. Pauluis, β€œSources and sinks of available potential energy in a moist atmosphere,” Journal of the Atmospheric Sciences, vol. 64, no. 7, pp. 2627–2641, 2007. View at Publisher Β· View at Google Scholar Β· View at Scopus
  17. F. Paparella and W. R. Young, β€œHorizontal convection is non-turbulent,” Journal of Fluid Mechanics, vol. 466, pp. 205–214, 2002. View at Publisher Β· View at Google Scholar Β· View at MathSciNet Β· View at Scopus
  18. W. Wang and R. X. Huang, β€œAn experimental study on thermal circulation driven by horizontal differential heating,” Journal of Fluid Mechanics, vol. 540, pp. 49–73, 2005. View at Publisher Β· View at Google Scholar Β· View at Scopus
  19. M. E. McIntyre, β€œOn spontaneous imbalance and ocean turbulence: generalisations of the Paparella-Young epsilon theorem,” in Proceedings of the International IUTAM/Newton Institution Workshop, D. G. Dritschel, Ed., Turbulence of the Atmosphere and Oceans, Springer, 2010.
  20. S. R. de Groot and P. Mazur, Non-Equilibrium Thermodynamics, North Holland, 1962.
  21. T. J. McDougall, β€œPotential enthalpy: a conservative oceanic variable for evaluating heat content and heat fluxes,” Journal of Physical Oceanography, vol. 33, no. 5, pp. 945–963, 2003. View at Publisher Β· View at Google Scholar Β· View at Scopus
  22. R. Tailleux, β€œIdentifying and quantifying nonconservative energy production/destruction terms in hydrostatic Boussinesq primitive equation models,” Ocean Modelling, vol. 34, no. 3-4, pp. 125–136, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus