Next Article in Journal
Modeling the Viscosity of Concentrated Nanoemulsions and Nanosuspensions
Next Article in Special Issue
Fofonoff Negative Modes
Previous Article in Journal
Investigation of Slot-Burner Aerodynamics with Recessed-Type Nozzle Geometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Semicompressible Ocean Thermodynamics and Boussinesq Energy Conservation

1
Department of Earth, Ocean, and Atmospheric Science, Florida State University, Tallahassee, FL 32309, USA
2
School of Mathematics and Statistics, UNSW, Sydney, NSW 2052, Australia
3
Institut für Mathematik, Freie Universität, Berlin 14195, Germany
*
Author to whom correspondence should be addressed.
Submission received: 29 January 2016 / Revised: 21 March 2016 / Accepted: 24 March 2016 / Published: 8 April 2016
(This article belongs to the Collection Geophysical Fluid Dynamics)

Abstract

:
Equations more accurate than the Boussinesq set that still filter out sound were recently introduced. While these equations were shown to have a consistent potential energy, their thermodynamical behavior and associated implications were not fully analyzed. These shortcomings are remedied in the present note that argues both sets are fully consistent from a thermodynamic perspective. It is further argued that both sets remain computationally competitive with the Boussinesq set.

1. Introduction

Dewar et al. [1] recently developed two sets of ocean dynamical equations that include some degree of compressibility and exclude sound waves. They were argued to be simultaneously more accurate than the Boussinesq equations and computationally competitive with them. These were referred to as the Type I and Type II Semicompressible equations.
In subsequent communications with Prof. R. Klein, he correctly pointed out that the thermodynamic consistency of these equations had not yet been demonstrated. Thermodynamic consistency in this context means that total energy is conserved and that an entropy variable exists for the system that obeys the second thermodynamic law. This remark was made in comparison to pseudoincompressible equations derived in [2] for a general equation of state in which a much more thorough thermodynamic analysis had been undertaken. The purpose of this short communication is to analyze similar characteristics of both types of the semicompressible equations.

2. The Semicompressible Equations

As discussed in [1], the Type I equations are
ρ * ( t u + u · u + 2 Ω × u ) = p [ ρ * ρ r ( z ) + p ρ r c s 2 ρ * ] g k + · ρ * ν v i s c u t ρ * + · ( ρ * u ) = 0 ρ * = ρ ( S A , Θ , P * )
d d t Θ = Θ o d d t S A = S A o
where density, ρ * , is evaluated using Conservative Temperature, Θ, Absolute Salinity, S A and the hydrostatic pressure, P * , determined from the reference density profile, ρ r (at most a function of z) according to
z P * = ρ r g
The quantity ν v i s c is viscosity and ν * = 1 / ρ * .
De Szoeke and Samelson [3] have shown that the hydrostatic Boussinesq equations are equivalent to the compressible hydrostatic Navier-Stokes equations written in pressure coordinates. Equation (1a) improves upon the classic hydrostatic balance in the hydrostatic limit through the presence of the sound speed. The density approximation it represents is a more accurate approximation of true density than ρ * , density evaluated at the hydrostatic pressure. Equation (1) is also more flexible in that they can be used in non-hydrostatic settings.
The notation X o in Equation (1) denotes the non-advective contributions to X. Absolute Salinity, being a purely conserved quantity (ignoring chemistry, see [4]), the form is clear
S A o = ν · F S
where F S is the diffusive flux of Absolute Salinity. We purposely leave the form for Θ o unspecified, as part of the present exercise is to determine it. Other notation is standard.
The equations in Equation (1) satisfy the energy equation
t ( ρ * ( K + h I + ) ) + · ( u ( ρ * ( K + h I + ) ) + p u ) = p R + · ρ * ν v i s c K ρ * ϵ
with potential energy played by the quantity
h I + = P o P * ρ ρ r ρ ρ r d P = P o P * b I g ρ r d P
The quantity ϵ = ν v i s c | u | 2 is the viscous dissipation of kinetic energy, b I is the Type I buoyancy, the quantity R
R = Θ ν * d d t Θ S A ν * d d t S A
denotes diffusive terms and Equation (1b) demonstrates the presence of compressibility in these equations.
The Type II equations are
ρ I I ( t u + u · u + 2 Ω × u ) = p ( ρ * ρ I I ) g k p c s I I 2 g k + · ρ I I ν v i s c u t ρ I I + · ( ρ I I u ) = · ( ρ I I u ) = 0 ρ * = ρ ( S A , Θ , P * )
d d t Θ = Θ o d d t S A = S A o
where
ρ I I ( z ) = ρ I I ( 0 ) e o z g c s I I 2 d z
is either a chosen reference density profile yielding the sound speed profile c s I I 2 or a profile constructed from a chosen sound speed profile, c s I I 2 . In either case, density, ρ I I , is a function of depth only and the hydrostatic pressure is computed from it
z P * = ρ I I g
Compressibility appears in these equations in Equation (7b) where the time independence of ρ I I has been recognized. The Type II equations satisfy the energy constraint
t ( ρ I I ( K + h I I + ) ) + · ( u ρ I I ( K + h I I + ) + p u ) = · ρ I I ν v i s c K ρ I I ϵ
where the role of potential energy is now played by
h I I + = P o P * ρ ρ I I ρ I I 2 d P = P o P * b I I g ρ I I d P
The quantity b I I is the Type II buoyancy. The advantages and disadvantages of these equations are discussed in [1]. Extracting Equations (1) and (7) from the full Navier–Stokes equations is akin to substituting the Boussinesq density for the full density in most places, and retaining a reference density profile where needed for potential energy consistency.

3. Semicompressible Thermodynamics

We now consider how to augment (1) and (7) with a consistent thermodynamic behavior. The First Law of Thermodynamics is
d d t h ν d d t P = T d d t η + μ d d t S A = 1 ρ · F Q + ϵ
where T is temperature, h ( η , S A , P ) is specific enthalpy, η is specific entropy, specific volume is denoted by ν = 1 / ρ , P is pressure and μ is the relative chemical potential of salt in seawater. The quantity ϵ is the heat of viscous energy dissipation, and F Q is a generalized heat flux. Enthalpy is defined by
h = e + P ν
where e is specific internal energy. Enthalpy from Equation (12) is seen to depend naturally on the variables η, S A , and P, and in this form plays the role of a thermodynamic potential from which all thermodynamic variables can be obtained by differentiation. For example, derivatives of enthalpy with respect to its natural variables are T, μ, and ν, respectively. The derivation of Equation (12) is discussed extensively in the TEOS-10 manual ([4]; see also [5,6]).
Equation (12) is exact. The semicompressible equations consider modifications to the thermodynamic variables caused mainly by the difference between static and full pressure, and we analyze Equation (12) from that perspective. In what follows, we will adopt entropy, salinity, and pressure as state variables, and will often cast pressure dependency in terms of static pressure. For example,
h = h * ( S A , η , P * ) + h P * ( S A , η , P * ) p + O ( p 2 )
where P * is static pressure and p the dynamic pressure, given by p = P P * . The superscript ( * ) will denote quantities evaluated at the static pressure. Specific volume, density, etc. will be written in a manner similar to Equation (14). Note, according to Equation (12), Equation (14) is equivalent to
h = h * ( S A , η , P * ) + ν * p

3.1. Type I Thermodynamics

We first develop what is a consistent statement of Type I thermodynamics. We begin by considering the form of Equation (12) to O ( p 2 ) ; i.e.,
d d t ( h * + ν * p ) = ( ν * + p P ν ) d d t P * + ν * d d t p + T + d d t η + μ + d d t S A + O ( p 2 )
where
T + = T * + p P T
and similarly for μ + . Evaluating the derivative of the static pressure, Equation (16) becomes
d d t ( h * + ν * p ) = ( ν * + p P ν ) ( w g ρ r ) + ν * d d t p + T + d d t η + μ + d d t S A
While Equation (18) is literally correct to O ( p 2 ) , it will turn out that a less accurate approximation, to O ( p ) , is required on the right hand side of Equation (12)
T + d d t η + μ + d d t S A = ν * F Q + ϵ
The reason for this will become clear later, here we stress only that it assures thermodynamic consistency of the Type I equations.

3.2. Type I Mechanical Energy Equation

The mechanical energy equation derived from Equation (1a) is
ρ * d d t K = u · p w g ( ρ * ρ r + ρ r p ρ * c s 2 ) + · ρ * ν v i s c K ρ * ϵ
where K = ( u · u ) / 2 is the kinetic energy density. Introducing the gravitational potential Φ = g z and using Equation (18)
ρ * ( d d t K + Φ + h * + p ν * ) = u · p + d d t p + · ρ * ν v i s c K ρ * ϵ + ρ * ( T + d d t η + μ + d d t S A )
Using Equation (12)
P h * = ν *
so
h * = h o + P o P * ν d P = h o + P o P * 1 ρ r d P + P o P * b ρ r g d P
where
b = g ( ρ * ρ r ) ρ *
The quantity h o is recognized as the potential enthalpy defined by [6], and the last integral in Equation (23) is the so-called dynamic enthalpy of [7]. The quantity P o is a constant reference surface pressure. The middle integral in Equation (23) reduces upon inspection to Φ .
Thus, Equation (21) becomes
ρ * d d t ( K + h o + h I + ) = · p u + · ρ * ν v i s c K ρ * ϵ + ρ * ( T + d d t η + μ + d d t S A )
Lastly, we exploit Equation (19) to obtain
ρ * d d t ( K + h o + h I + ) = · p u + · ρ * ν v i s c K · F Q
which is in conservative form. This implies the Type I equations have a well-formed energy principle, provided they are augmented by Equations (18) and (19) as their thermodynamic equations. Note that the entropy in Equation (19) can be rewritten as
ρ * d d t η = · F Q μ + F S T + + ( F Q μ + F S ) · ( 1 T + ) F S · ( μ + T + ) + ρ * ϵ T +
which is entropy “non-conservation” in its familiar form. Hence, the Type I equations also have a consistent entropy variable.
At this point, a practical difficulty is that entropy is not a typical ocean modelling variable. Instead, McDougall’s Conservative Temperature
Θ = h o c P o
McDougall ([6]) is preferable and so it is necessary to work out a Conservative Temperature equation to complement (27). From Equation (18), one finds
d d t h * ν * d d t P * = 1 ρ * · F Q p P T d d t η p P μ d d t S A + ϵ
By definition
h o = h ( η , S A , P o )
so
d d t h o = θ d d t η + μ o d d t S A
where θ is potential temperature and μ o = μ ( η , S A , P o ) . Potential enthalpy is thus related to entropy via
d d t η = d d t h o μ o d d t S A θ
From Equation (23)
d d t h o = c p o d d t Θ = 1 ρ * · F Q ( T + θ ) d d t η ( μ + μ o ) d d t S A + ϵ
so using Equation (32)
d d t Θ = ν * · F Q + ϵ ̲ c P o ( 1 + ( T + θ ) θ ̲ ) ( ( T + θ ) μ o μ + ) d d t S A c P o ( 1 + ( T + θ ) θ ) ̲
The underlined terms in Equation (34) are quite small, according to [8]. Neglecting them leads to the statement that Conservative Temperature is very accurately portrayed as a conserved quantity; however, all of those terms must be included in order to ensure full energy conservation.
It is seen, therefore, that the Type I equations yield energy conservation, but the cost is the insertion of a large quantity ( c P o Θ ) into the energy equation. A similar result was labelled in [7] as a “crippling” disadvantage of the Boussinesq-Bernoulli equation. It might thus be more useful to work with a partial equation when considering energy; i.e., one that doesn’t explicitly include Conservative Temperature, but in any case, the full potential enthalpy Equation (34) is required. We have also found the Type I Semicompressible equations have a proper entropy variable and conserve total energy, and so are fully thermodynamically consistent. These equations are very similar to the pseudo-incompressible equation set for general equations of state in [2].

3.3. Type II Equations: Potential, Kinetic, and Thermal Energies

The Type II equations are quite similar to the Type I equations except for the appearance of the reference density profile, rather than the Boussinesq density, in the mass conservation equation and in front of the momentum acceleration. Most of the analysis carries over in a manner similar to the Type I analysis, with the one important distinction that the velocity divergence diagnosed from mass conservation has no diabatic contributions
· u I I = g w c s I I 2
Therefore, the mechanical energy equation for the system analogous to Equation (21) is
ρ I I d d t K = · p u + p w g c s I I 2 ̲ w g ( ρ * ρ I I + p c s I I 2 ̲ ) + · ρ I I ν v i s c K ρ I I ϵ
with the underlined terms cancelling. Taking the material derivative of Equation (11) and substituting in Equation (36) yields
ρ I I d d t K = · p u + ρ I I P o P * S A b I I g ρ r d P d d t S A + ρ I I P o P * η b I I g ρ r d P d d t η ρ I I d d t h I I + + · ρ I I ν v i s c K ρ I I ϵ
To put Equation (37) in conservative form requires the potential enthalpy equation
d d t h * = d d t h o + ν * d d t P * + P o P * ν S A d P d d t S A + P o P * ν η d P d d t η
which is equivalent to Equation (33). This leads to
ρ I I d d t ( K + h I I + + c P o Θ ) = · p u · F Q + · ρ I I ν v i s c K
The Type II Conservative Temperature equation becomes
d d t Θ = ν I I · F Q + ϵ ( ( T + θ ) μ o μ + ) d d t S A c P o T + θ
The most significant difference relative to the Type I set is in Equation (35), where the diabatic contributions to expansion are absent. Like for Type I, it is still necessary to retain all terms for Conservative Temperature, rather than the approximate form recommended in [6]. We conclude that the Type II equations also have a consistent mechanical and thermal energy structure. Our Type II equations are very similar to the thermodynamically consistent anelastic set derived in [9].

3.4. A Comparison to the Boussinesq Set

To ensure full energetic consistency for the semicompressible equations, it has been necessary to carry along some additional non-conservative effects in the calculations. We now analyze the full energetics behavior of the Boussinesq set to gauge any additional computational burden of the semicompressible set relative to them. If we adopt Conservative Temperature, salinity, and pressure as the thermodynamic coordinates and follow the procedure outlined in [7], we arrive at the Boussinesq version of the Bernoulli equation
d d t ( K + p + h B ) = t p + P o P * S A b g ρ o d P d d t S A + P o P * Θ b g ρ o d P d d t Θ ϵ
(see Equation (48) in [7]) where b = g ( ρ * ρ o ) / ρ * is the Boussinesq buoyancy and ρ o is the Boussinesq reference density. Dynamic enthalpy as defined in [7] is
h B = 1 g ρ o P o P * b d P
Returning to a form of the potential enthalpy consistent with the Boussinesq approximation
h * ( S A , Θ , P * ) = h o + 1 ρ o P o P * ( 1 + b g ) d P
where h o is potential enthalphy, and taking its material derivative yields
d d t h * ν * d d t P * = c P o d d t Θ + P o P * S A b g ρ o d P d d t S A + P o P * Θ b g ρ o d P d d t Θ = 1 ρ o · F Q + ϵ
Using Equation (44) to eliminate the pressure integrals in Equation (41) turns it into the conservative Boussinesq energy equation
d d t ( K + p + h B + c P o Θ ) = t p + · ρ I I ν v i s c K 1 ρ o · F Q
The Boussinesq Conservative Temperature equation is given by
d d t Θ = ν o · F Q + ϵ ( ( T * θ ) μ o μ + ) d d t S A c P o T * θ
The point of this section is that to obtain a conservative energy equation out of the Boussinesq set, it is necessary to carry along the same additional diffusive terms as for the semi-compressible equations. A similar result has recently been derived in [10,11]. This is not done in any Boussinesq ocean circulation model of which we are aware.

4. Summary

The Type I and II semicompressible equations discussed in [1], which are analogues of the pseudo-incompressible and anelastic equations known in meteorology, have undergone an examination of their thermodynamic structure. To summarize, the Types I and II equations are fully consistent from an energetic and thermodynamic perspective. To obtain this, it is necessary to include a more complete equation for Conservative Temperature as a part of the equation sets than is normally done. Interestingly, it turns out the Boussinesq equations must also carry analogous quantities to be in purely conservative energetic form. We thus see the semicompressible equations remain computationally competitive with the Boussinesq equations while admitting some degree of compressibility. Other results here include the associated definitions of the dynamic enthalpy variables belonging to the sets.
From a simple computer throughput perspective, this analysis recommends Type II relative to Type I because they are slightly less computationally intensive. However, the Type II equations represent the oceanic density structure by a single vertical profile. Any such profile cannot perform well globally, so the “cost” of the simplicity is a less accurate density representation. This promises to be the most troubling when computing sea level changes via the continuity equation. The additional accuracy provided by the Boussinesq density in the Type I set might well compensate for the (slight) additional complexity.
Any further exploration with these equation sets should bear in mind the above points.

Acknowledgments

This research was make possible by grants from the National Science Foundation Grants OCE-1537304 and OCE-100743 and the Gulf of Mexico Research Initiative through consortium funding to the University of Miami for the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE). William Young is thanked for many interesting conversations that added significantly to the manuscript.

Author Contributions

William K. Dewar and Rupert Klein conceived the research. Trevor McDougall and Joseph Schoonover contributed substantially to the thermodynamical analysis. William K. Dewar wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dewar, W.; Schoonover, J.; McDougall, T.; Young, W. Semicompressible ocean dynamics. J. Phys. Oceanogr. 2015, 45, 149–156. [Google Scholar] [CrossRef]
  2. Klein, R.; Pauluis, O. Thermodynamic consistency of a psuedoincompressible approximation for general equations of state. J. Atmos. Sci. 2012, 69, 961–968. [Google Scholar] [CrossRef]
  3. De Szoeke, R.; Samelson, R. The duality between the Boussinesq and non-Boussinesq hydrostatic equations of motion. J. Phys. Oceanogr. 2002, 32, 2194–2203. [Google Scholar] [CrossRef]
  4. Intergovernmental Oceanographic Commission (IOC); Scientific Committee on Oceanic Research (SCOR); International Association for the Physical Sciences of the Oceans (IAPSO). The International Thermodynamic Equation of Seawater—2010: Calculation and Use of Thermodynamic Properties; 2010 Manuals and Guides No. 56, UNESCO (English); Intergovernmental Oceanographic Commission: West Perth, Australia, 2010; p. 196. [Google Scholar]
  5. Landau, L.M.; Lifshitz, E.M. Fluid Mechanics, 1st ed.; Pergamon Press: Oxford, UK, 1959. [Google Scholar]
  6. McDougall, T. Potential enthalpy: A conservative oceanic variable for evaluating heat content and heat fluxes. J. Phys. Oceanogr. 2003, 33, 945–963. [Google Scholar] [CrossRef]
  7. Young, W. Dynamic enthalpy, conservative temperature and the seawater Boussinesq approximation. J. Phys. Oceanogr. 2010, 40, 394–400. [Google Scholar] [CrossRef]
  8. Graham, F.S.; McDougall, T.J. Quantifying the Nonconservative Production of Conservative Temperature, Potential Temperature, and Entropy. J. Phys. Oceanogr. 2013, 43, 838–862. [Google Scholar] [CrossRef]
  9. Pauluis, O. Thermodynamic consistency of the anelastic approximatio for a moist atmosphere. J. Atmos. Sci. 2008, 65, 2719–2729. [Google Scholar] [CrossRef]
  10. Tailleux, R. On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy. J. Fluid Mech. 2009, 638, 339–382. [Google Scholar] [CrossRef]
  11. Tailleux, R. Irreversible compressible work and available potential energy dissipation in turbulent stratified fluids. Phys. Scr. 2013. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Dewar, W.K.; Schoonover, J.; McDougall, T.; Klein, R. Semicompressible Ocean Thermodynamics and Boussinesq Energy Conservation. Fluids 2016, 1, 9. https://doi.org/10.3390/fluids1020009

AMA Style

Dewar WK, Schoonover J, McDougall T, Klein R. Semicompressible Ocean Thermodynamics and Boussinesq Energy Conservation. Fluids. 2016; 1(2):9. https://doi.org/10.3390/fluids1020009

Chicago/Turabian Style

Dewar, William K., Joseph Schoonover, Trevor McDougall, and Rupert Klein. 2016. "Semicompressible Ocean Thermodynamics and Boussinesq Energy Conservation" Fluids 1, no. 2: 9. https://doi.org/10.3390/fluids1020009

APA Style

Dewar, W. K., Schoonover, J., McDougall, T., & Klein, R. (2016). Semicompressible Ocean Thermodynamics and Boussinesq Energy Conservation. Fluids, 1(2), 9. https://doi.org/10.3390/fluids1020009

Article Metrics

Back to TopTop