Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa 239-8686, Japan May 9, 2014 Abstract This paper presents a theory escribing the energy buget of a flui uner the Boussinesq approximation: the theory is evelope in a manner consistent with the conservation law of mass. It shows that no potential energy is available uner the Boussinesq approximation, an also reveals that the work one by the buoyancy force ue to changes in temperature correspons to the conversion between kinetic an internal energy. This energy conversion, however, makes only an ignorable contribution to the istribution of temperature uner the approximation. The Boussinesq approximation is, in physical oceanography, extene so that the motion of seawater can be stuie. This paper consiers this extene approximation as well. Uner the extene approximation, the work one by the buoyancy force ue to changes in salinity correspons to the conversion between kinetic an potential energy. It also turns out that the conservation law of mass oes not allow the conition u = 0 on the flui velocity u to be impose uner the extene approximation; the conition to be impose instea is presente. 1 Introuction The Boussinesq approximation is frequently use to stuy the motion of a flui with a non-uniform temperature istribution. Uner this approximation, Winters et al. (1995) presente a theoretical framework for analyzing the energy buget of a ensity-stratifie flow. Since then, in physical oceanography, the energetics of the overturning circulation of the oceans has been stuie within this theoretical framework (Huang 1998; Hughes, Hogg & Griffiths 2009). This framework, however, is physically inamissible as explaine below. Winters et al. assume in their paper that the ensity ρ appearing in their efinition of the potential energy of a flui fulfille the following equation: ρ/ t+u ρ = j, (1) where u is the flui velocity satisfying the conition u = 0, an j enotes the iffusion flux ensity of mass allowe in their paper to exist. This equation can be rewritten, using u = 0, in the form ρ/ t+ (ρu+j ) = 0. (2) 1
It follows from this expression that ρu+j is the momentum of unit volume of flui (see Lanau & Lifshitz 1987, 49). Hence, so long as j 0, we arrive at the irrational conclusion that u is not the momentum of unit mass of flui. The origin of this irrational conclusion is as follows. To erive the equation of motion for a flui uner the Boussinesq approximation, the ensity ρ of the flui is customarily assume to be given by (see Lanau & Lifshitz 1987, 56) ρ = ρ 0 ρ 0 βt, (3) where ρ 0 is a constant reference ensity, β is the thermal expansion coefficient, an T is the temperature of the flui measure from a reference value. This assumption, combine with the conition u = 0, gives rise to the iffusion flux of mass; we are consequently le to the above irrational conclusion. The assumption (3), however, shoul be interprete as a mere expeient to erive the equation of motion for the following reason. Suppose that the flui is containe in a omain with a fixe volume V. The mass M of the flui is then given by M = ρ 0 V ρ 0 βvt, (4) where T enotes the average temperature of the flui. From this expression, it follows that the mass M varies when the flui is heate or coole. Eviently, this conclusion violates the conservation law of mass: the law requires, instea of (3), that ρ = ρ 0. It is now evient that the efinition by Winters et al. of the potential energy of a flui is physically unacceptable. Since their theoretical framework contains no mathematical error, it may be use for analyzing numerical simulations (see e.g. Gayen et al. 2013); it cannot be use, nevertheless, for analyzing real flows. Thus the aim of this paper is to present a physically reasonable theory that can escribe the energy buget of a flui uner the Boussinesq approximation. It is also extene so that the motion of seawater can be ealt with. 2 Energetics uner the Boussinesq approximation Let us consier the motion uner gravity of a viscous flui containe in a fixe omain. We set up in the omain a system of rectangular coorinates (x 1,x 2,x 3 ) with the x 3 -axis taken vertically upwars. The unit vectors in the positive x 1 -, x 2 -, an x 3 -irections are respectively enote by e 1, e 2, an e 3. In the following, Latin inices are use to represent the numbers 1, 2, an 3; the summation convention is also implie. As explaine in 1, in view of the conservation law of mass, the ensity ρ of the flui shoul be regare uner the Boussinesq approximation as constant: ρ = ρ 0. (5) This is consistent with the usual assumption uner the approximation that the velocity u of the flui is solenoial: u = 0. (6) 2
Nevertheless, the thermal expansion coefficient β of the flui is assume not to vanish: β = v 1 ( v/ T) p 0. (7) Here v = 1/ρ enotes the specific volume of the flui; T an p are respectively the temperature an the pressure of the flui. We first wish to fin an equation for the rate of change of the internal energy of the flui. This can be erive from the general equation of heat transfer (see Lanau & Lifshitz 1987, 49): ρt = τ u i ij q, (8) in which D/ enotes the material erivative, s is the specific entropy of the flui, τ ij are the components of the viscous stress tensor, u i are those of u, an q is the heat flux ensity; the first term on the right represents the heating ue to viscous issipation. Uner the Boussinesq approximation, the temperature an the pressure are assume to vary only slightly in the flui (see Lanau & Lifshitz 1987, 56). Accoringly, T can be written in the form T = T 0 +T, (9) where T 0 is a constant reference temperature, an T is the small variation from T 0. We also write, enoting by p the small perturbation pressure, p = p 0 +p. (10) Here p 0, which satisfies the hyrostatic equation, is efine by where g is the acceleration ue to gravity. Now let us express / in terms of T an p: ( ) ( ) s = DT s T + p p 0 = ρ 0 gx 3 +constant, (11) p T Dp. (12) Here ( s/ T) p an ( s/ p) T are given by (see e.g. Batchelor 1967, 1.5) ( ) s = C ( ) ( ) p s v T T, = = vβ, (13) p T p in which C p is the specific heat at constant pressure. Since T an p vary only slightly, C p an β may be taken to be constant; we may also put T = T 0 in (13). Thus we have, since v = v 0 = 1/ρ 0, T = C p T 0 DT v 0β p ( ) Dp0 + Dp, (14) where (9) an (10) have been use. It is also possible to express / in terms of the specific energy e an the specific volume v: approximating T an p by T 0 an p 0, we get = 1 T 0 De + p 0 T 0 Dv. (15) 3
We note, however, that Dv/ = v u = 0. Hence it follows that Using (5), (9), (14), an (16), we fin ρt = ρ 0T 0 +ρ 0T = ρ 0 De +ρ 0T = 1 De T 0. (16) { Cp T 0 DT v 0β ( )} Dp0 + Dp. (17) Since Dp 0 / = ρ 0 gu 3, we obtain, to the first orer of prime variables, The substitution of this result into (8) yiels ρt = ρ De 0 +ρ 0βT gu 3. (18) De ρ 0 = τ u i ij q ρ 0 βt gu 3. (19) Integrating (19) over, whose bounary is enote by, we have u i ρ 0 ev = τ ij V q ns ρ 0 βt gu 3 V. (20) Here n is the unit outwar normal on, an we have assume that u n = 0 on. This is the esire equation for the rate of change of the internal energy. Let us next consier the kinetic energy of the flui. Uner the Boussinesq approximation, the equation of motion is given by (see Lanau & Lifshitz 1987, 56) Du ρ 0 = p + τ ij e i +ρ 0 βt ge 3. (21) From this equation, we can obtain the following equation for the rate of change of the kinetic energy: 1 2 ρ 0 u 2 u i V = u i τ ij n j S τ ij V + ρ 0 βt gu 3 V, (22) where n j are the components of n. On the other han, the potential energy of the flui is invariable: ρ 0 gx 3 V = 0. (23) Hence, uner the Boussinesq approximation, no potential energy is available. Aing (20), (22), an (23), we get ( ρ 1 0 2 u 2 +gx 3 +e ) V = u i τ ij n j S q ns. (24) This expression states that the total energy of the flui changes owing to the work one by the viscous force acting on an owing also to the heat transfer across. We have thus obtaine a physically reasonable result. 4
We now focus our attention on the equation (22). Of the three terms on its right-han sie, the last term represents the work one by the buoyancy force: this force is represente by the last term of (21). The last term of (22) appears, with an opposite sign, in (20) too. Hence the work one by the buoyancy force is seen to correspon to the conversion between kinetic an internal energy. Finally, let us return to the equation (19). Rewriting this equation with the help of (14) an (16), we obtain DT { ( } ρ 0 C p = τ u i Dp0 ij q + βt 0 ) ρ + Dp 0 βt gu 3. (25) The sum of the terms in the braces, however, is neglecte uner the Boussinesq approximation together with the term representing the heating ue to viscous issipation (see Lanau & Lifshitz 1987, 56). Thus we have DT = 1 ρ 0 C p q. (26) This equation, combine with e.g. Fourier s law, etermines the istribution of temperature uner the Boussinesq approximation. However, as is apparent from the erivation, this temperature equation takes no account of the contribution of the conversion between kinetic an internal energy iscusse above. 3 Energetics uner the oceanographic Boussinesq approximation In physical oceanography, the Boussinesq approximation is extene to stuy the motion ofseawater(see e.g. McWilliams 2006, 2.2). Ouraim in this section is to investigate the energy buget of a flui uner this extene approximation which is calle in this paper the oceanographic Boussinesq approximation. We consier the same situation as that in 2, but the flui to be stuie is now regare as seawater. The salinity c, the mass of issolve salt per unit mass of flui, is assume to change only slightly in the flui. Thus we write c = c 0 +c, (27) where c 0 is some constant reference salinity, an c is the small eviation from c 0. The equation for c is as follows (see Lanau & Lifshitz 1987, 58): ρ Dc = i, (28) where i enotes the iffusion flux ensity of salt; it shoul be note that the flux ensity of salt an that of water are given by ρcu +i an ρ(1 c)u i, respectively. We also nee to introuce the chemical potential µ of the flui (see Lanau & Lifshitz 1987, 58). It is regare as a function of T, p, an c: µ = µ(t,p,c). (29) When c is uniform throughout the flui, the ensity ρ shoul be equate to a constant as in 2. This constant, however, may be a function of c. Hence, in view of the smallness of c, we put ρ = ρ 0 +ρ 0 β s c, (30) 5
in which the haline contraction coefficient β s = ρ 1 ( ρ/ c) T,p is regare as a constant. Then, from the equation of continuity ρ 1Dρ + u = 0, (31) it follows that in general u 0. Using (30) an (31), we can obtain ρgx 3 V = (ρ 0 gu 3 +ρ 0 β s c gu 3 )V. (32) Here we have assume, for simplicity, that u n = 0 on. This is the potential energy equation for the present problem. The internal energy equation can now be foun from the following general equation of heat transfer (see Lanau & Lifshitz 1987, 58): ρt = τ u i ij q +µ i. (33) When s is regare as a function of T, p, an c, we get, in place of (14), = C p DT ( ) ( ) T 0 vβ Dp0 + Dp µ Dc T. (34) Here the thermoynamic relation ( s/ c) T,p = ( µ/ T) p,c has been use (see Lanau & Lifshitz 1987, 59); note also that v 0 in (14) has been replace by v in view of (30). Since ( s/ c) e,v = µ/t, we also have, in place of (16), = 1 De T 0 + p 0 Dv T 0 µ 0 Dc T 0, (35) where µ 0 = µ(t 0,p 0,c 0 ). Thus, neglecting terms above the first orer in prime variables, we can obtain the following equation corresponing to (19): ρ De = p 0 u+τ ij u i q ρ 0 βt gu 3, (36) where ρdv/ = u an (28) have been use. We note in aition that p 0 u = (p 0 u)+u p 0 = (p 0 u) ρ 0 gu 3. (37) The internal energy equation is therefore obtaine as follows: u i ρev = τ ij V q ns (ρ 0 gu 3 +ρ 0 βt gu 3 )V. (38) Next, the equation of motion is now given by (see McWilliams 2006, 2.2) ( ) Du ρ 0 = ρ0 τij p + e i +ρ 0 βt ge 3 ρ 0 β s c ge 3, (39) ρ in which the last term represents the buoyancy force ue to changes in salinity, an the factor ρ 0 /ρ is usually ientifie with unity. Neglecting terms above the first orer in prime variables, we can fin from (39) the kinetic energy equation 1 2 ρ u 2 V = u i τ ij n j S τ ij u i V + p,c (ρ 0 βt gu 3 ρ 0 β s c gu 3 )V. (40) 6
In consequence, aing (32), (38), an (40), we get ρ ( 1 2 u 2 +gx 3 +e ) V = u i τ ij n j S q ns. (41) This total energy equation has exactly the same form as (24); we have obtaine a physically reasonable result again. The secon integran in the last term of (40) represents the work one by the buoyancy force ue to changes in salinity. Comparing (40) with (32), we observe that this work, in contrast to the work one by the buoyancy force ue to changes in temperature, correspons to the conversion between kinetic an potential energy. On the other han, we recognize that the first integran on the right-han sie of (32) correspons to the conversion between potential an internal energy; it appears, with an opposite sign, in the last term of (38). 4 Summary an iscussion The energy buget of a flui has been investigate, in a manner consistent with the conservation law of mass, uner the Boussinesq approximation. It has been shown that no potential energy is available uner the approximation. It has also turne out that the work one by the buoyancy force ue to changes in temperature correspons to the conversion between kinetic an internal energy. This energy conversion, however, is altogether neglecte in the etermination of the istribution of temperature uner the approximation. In contrast, uner the oceanographic Boussinesq approximation, the work one by the buoyancy force ue to changes in salinity represents the conversion between kinetic an potential energy. In aition, energy conversion also arises between potential an internal energy uner the approximation. In the iscussion of the oceanographic Boussinesq approximation, we foun that in general u 0. If we ha assume that u = 0 together with (30), then the iffusion flux of mass woul have arisen. The occurrence of this flux is physically unacceptable, as state in 1. Thus it follows that we cannot take the conition u = 0 uner the oceanographic Boussinesq approximation. Let us examine what conition shoul be taken then in place of u = 0. The substitution of (27) an (30) into (28) leas to ρ 0 Dc = i, (42) where terms above the first orer in prime variables have been neglecte. We can similarly fin from (31) the following equation: From (42) an (43), we get β s Dc = u. (43) (u β s i/ρ 0 ) = 0. (44) We have thus obtaine the conition to be impose in place of u = 0 uner the oceanographic Boussinesq approximation; this conition guarantees that u is the momentum of unit mass of flui. 7
References [1] Batchelor, G. K. 1967 An Introuction to Flui Dynamics. Cambrige University Press. [2] Gayen, B., Griffiths, R. W., Hughes, G. O. & Saenz, J. A. 2013 Energetics of horizontal convection. J. Flui Mech. 716, R10. [3] Huang, R. X. 1998 Mixing an available potential energy in a Boussinesq ocean. J. Phys. Oceanogr. 28, 669 678. [4] Hughes, G. O., Hogg, A. McC. & Griffiths, R. W. 2009 Available potential energy an irreversible mixing in the meriional overturning circulation. J. Phys. Oceanogr. 39, 3130 3146. [5] Lanau, L. D. & Lifshitz, E. M. 1987 Flui Mechanics, 2n e. Pergamon. [6] McWilliams, J. C. 2006 Funamentals of Geophysical Flui Dynamics. Cambrige University Press. [7] Winters, K. B., Lombar, P. N., Riley, J. J. & D Asaro, E. A. 1995 Available potential energy an mixing in ensity-stratifie fluis. J. Flui Mech. 289, 115 128. 8