# (p.450) Appendix I A brief review of classical fluid dynamics

# (p.450) Appendix I A brief review of classical fluid dynamics

In this appendix, we summarize some of the basic notions of classical fluid dynamics, with particular emphasis on those concepts that are important for TDCDFT. This review is by no means exhaustive; more information and practical examples can be found in the numerous excellent books on this subject, such as the volume by Landau and Lifshitz (1987).

# I.1 Basics and ideal fluids

Fluid dynamics is concerned with the forces on and within a continuous fluid body and the resulting motion. The state of the system is conventionally described by three variables: the mass density *ϱ*(**r**, *t*), the velocity distribution **u**(**r**, *t*), and a third ther-modynamic variable such as the pressure distribution *p*(**r**, *t*). The relation between density and pressure is established via an appropriate thermodynamical equation of state, in which other quantities of interest such as the temperature also appear.^{1}

The most basic relation in fluid dynamics is the *continuity equation*,

*ϱ*in a given spatial region can result only from a net mass current density

*ϱ*

**u**entering or leaving that region.

Next, we wish to establish an equation of motion for the fluid under the influence of internal and external forces. According to Newton's second law, this equation should have the form (mass density) × (acceleration) = (force density).

Let us talk first about the proper treatment of acceleration. Figure I.1 shows the path of a fluid element moving along a streamline. As it travels from position 1 to position 2, the velocity vector of this fluid element changes in accordance with

The last term on the right-hand side is simply the change in the velocity vector at a particular point in space. But we want the rate of change of the velocity vector riding (p.451)

**u**·∇)

**u**+∂

**u**/∂

*t*. The general equation of motion in fluid dynamics can therefore be written in the following way:

**f**

_{total}is the total force density, which describes the sum of all internal and external forces acting on a fluid element.

We distinguish between *ideal* and *viscous* fluids. There is no energy dissipation in an ideal fluid, which means that there are no processes of internal friction or heat exchange between different fluid elements. Under such circumstances, the equation of motion (I.3) becomes

*p*is the pressure force on the fluid element and

**f**

_{ext}is caused by external force fields, usually a gravitational field. Equation (I.4) is known as the

*Euler equation*of fluid dynamics.

If one combines the Euler equation (I.4) and the continuity equation (I.1), it only takes a few elementary steps to show that the rate of change of the momentum density is given by

*μ*denotes one of the Cartesian components

*x*,

*y*, or

*z*, and

_{μν}of this tensor gives the amount of the

*μ*th component of the momentum flowing through a unit area element perpendicular to the

*r*

_{ν}-axis per unit time.

Equation (I.5) can be written in integral form, which gives the rate of change of the total momentum in a volume *V* bounded by a surface *S*:

(p.452) where we have used Green's theorem to transform the volume integral over the diver-gence of the momentum flux density tensor into a surface integral.

Now, assume that the system is finite and let the boundary surface go to infinity: the surface integral in eqn (I.7) then vanishes (since the density and pressure are zero outside the system), and only the term involving *f* _{ext,μ} survives. Thus, as expected, only external forces can cause a change in the total momentum of a finite system.

# I.2 Viscous fluids and dissipation

At the end of the previous section, we learned an important lesson: if we describe all internal forces or interactions within a system as the divergence of a momentum flux density tensor, it is automatically guaranteed that Newton's third law is satisfied. This will provide us with a useful guideline for how to describe viscous effects in a fluid.

Viscous fluids are subject to energy dissipation caused by internal frictional forces between fluid elements in motion. We would like to extend the Euler equation (I.4) or, equivalently, eqn (I.5) to include these internal forces and thus arrive at a more general description of the dynamics of real fluids. To do this, let us add a new term to the momentum flux tensor:

Frictional forces arise between neighboring fluid elements that have different veloc-ities, causing an irreversible tendency for neighboring fluid elements to reach the same final velocity, at which point the frictional forces disappear. Therefore, ${\sigma}_{\mu v}^{\text{VISC}}$ can only depend on the *derivatives* of the velocity, *∇* _{μ} *υ* _{ν}, and not on the fluid velocity itself.

Two generic situations are illustrated in Fig. I.2. The left-hand side shows a fluid flowing in the *x*-direction, with a velocity that is increasing along the *y*-direction, which causes shear stresses between neighboring fluid elements. For a flat cell of cross-sectional area ∆*S*, we find the shear stress—measured as a force per unit area—as

*η*is the coefficient of shear viscosity. This leads to the following shear contribution to the viscous stress tensor:

The same coefficient *η* is used for all *μ*, *ν*, since the fluid is isotropic. Furthermore, eqn (I.10) must contain the symmetric combination *∇* _{ν} *υ* _{μ} + *∇* _{μ} *υ* _{ν} due to the requirement that there can be no friction if the system is in a state of uniform rotation.^{2}

*ζ*′ is often referred to as the second viscosity coefficient. Taking the contributions together, we have

This, however, is not the standard way in which the viscous stress tensor is written. Instead, we redefine the second viscosity coefficient as $\zeta \text{'}=\zeta -\frac{2}{3}\eta $ so that

This expression has the property that the trace of the first term (associated with the shear viscosity) vanishes.

We thus obtain the following generalization of the Euler equation (I.4):

In general, the viscosity coefficients *η* and *ζ* can be complicated functions of space owing to their dependence on pressure and temperature. However, it is often a good approximation to assume that they are constants, and eqn (I.14) then simplifies to

Equation (I.15) is known as the *Navier-Stokes* equation and represents a cornerstone of classical fluid dynamics.

(p.454) Lastly, let us discuss the energy dissipation caused by viscoelastic stresses. The total kinetic energy in the fluid is given by

This expression can be analyzed further by making use of the continuity equation [eqn (I.1)] and the equation of motion (I.14). We will not carry out a full discussion here (see Landau and Lifshitz, 1987), but merely state the obvious fact that in the absence of external forces, the sum of the total kinetic energy and the total internal energy is conserved. In general, exchange between kinetic and internal energy can take place, even without any viscous forces, owing to the compressibility of the fluid. The total kinetic energy remains separately conserved only for incompressible ideal fluids without external forces.

What is of interest to us here is the irreversible change in the total kinetic energy caused by the viscosity of the fluid. We obtain

The change in kinetic energy per unit time due to the fluid viscosity corresponds to a dissipative power loss, and we define

Equation (I.18) can then be recast into

Inserting the explicit form of the viscoelastic stress tensor gives

This expression for the dissipated power must always be negative to indicate loss of mechanical energy due to internal friction, which places constraints on the allowable values of the viscosity coefficients. In the special case of an incompressible liquid (characterized by ∇ · **u** = 0), this implies that the shear viscosity coefficient *η* must always be a positive number.

## Notes:

(^{1})
If the medium is conducting, quantities such as the electric current and the magnetic field will be of interest; this falls within the realm of magnetohydrodynamics.

(^{2})
If the system is rotating with a uniform angular velocity Ω, then the velocity is given by **u** = Ω×**r**, and it can be easily verified that the combination ∇_{ν} *υ* _{μ} + ∇_{μ} *υ* _{ν} vanishes.