MA3D1 – Fluid Dynamics and MA4L0 – Advanced Topics in Fluids Lecture notes

Chapter 1 Preliminaries

The theory of fluid dynamics rests on the developments in some greater disciplines such as mathematics and physics, especially mechanics. In this chapter, let us summarize and revise some of these topics.

1.1 Newtonian mechanics

While it is not possible to do justice to Newtonian mechanics here, we present the basic tenets. 11 1 For more detailed treatments, refer to a book on undergraduate mechanics, e.g., Halliday and Resnick, Physics - I, or, if you prefer the economy of Soviet authors, I. E. Irodov, Fundamental Laws of Mechanics. Newtonian mechanics, at its core, describes physical laws governing the motion of point material particles. Consider a point material object of mass m moving through empty space. Let us say that the objects location is given by 𝒙(t) as a function of time t. The velocity of the body is 𝒗(t)=d𝒙/dt. The three laws of motion according to Newtonian mechanics are:

  1. 1.

    First law: In an inertial reference frame, every body continues in its state of rest or uniform motion until and unless acted upon by an external unbalanced force. Essentially, this law defines the notion of an inertial reference frame, i.e. one in which a free body moves uniformly.

  2. 2.

    Second law: When an external unbalanced force, 𝑭, acts on a body it must equal the rate of change of momentum of this body. The momentum (also known as the linear momentum) of this point mass is defined as 𝒑=m𝒗. Therefore, the mathematical statement of the second law is:

    𝑭=d𝒑dt. (1.1)

    Noting that the acceleration of the particle, 𝒂=d𝒗/dt=, if the mass of the body does not change with time, then this law is also written as 𝑭=m𝒂. This law describes the motion of the particle.

  3. 3.

    Third law: Every action has an equal and opposite reaction. In other words, if a body A exerts a force 𝑭 on a body B, then the body B must exert a force of 𝑭 on body A. This law describes the interaction between particles.

As a consequence of Newton’s laws of motion for point particles, it is possible to derive the following “conservation laws” of momentum and angular momentum for extended deformable bodies. A common terminology used in physics is the moment of a quantity, which is obtained by making a cross product of the position vector 𝒙 of a point particle with that quantity. Suppose that an extended material object at time t occupies a region of space given by Ω. Then its total momentum is obtained by breaking it up into infinitesimal pieces of mass dm, determining the momentum of these individual pieces, and summing the momentum. Mathematically, this amounts to integrating as

𝒑=Ω𝑑m𝒗(𝒙,t), (1.2)

where 𝒗 is the velocity of the piece at position 𝒙Ω. The total mass m of the body is simply m=Ω𝑑m. The moment of momentum , also known as the angular momentum 𝒍, is then

𝒍=Ω𝑑m𝒙×𝒗. (1.3)

These conservation laws invoke the “centre of mass” 𝒙c of this object, defined as

𝒙c=Ω𝑑m𝒙Ω𝑑m=1mΩ𝑑m𝒙. (1.4)

Suppose that the external forces 𝑭i, i=1,2,,n act on this object at locations 𝒙i. The net torque 𝝉 on this body is the moment of the external forces,

𝝉=i=1n𝒙i×𝑭i. (1.5)

So long as the object does not exchange material with anything outside it, the laws of conservation state that:

  1. 1.

    Conservation of (linear) momentum: The rate of change of total momentum must equal the total external unbalanced force on the object, i.e.

    d𝒑dt=i=1n𝑭i. (1.6)

    Noting the relation between 𝒑 and 𝒙c, this law is also stated as

    md2𝒙cdt2=i=1n𝑭i, (1.7)

    which means that the mass times the acceleration of the centre of mass must equal the external unbalanced force.

  2. 2.

    Conservation of angular momentum: The rate of change of angular momentum must equal the net torque on the body. Mathematically,

    d𝒍dt=𝝉. (1.8)

While the momentum and angular momentum are not constants of motion, here the term “conservation law” is used in the sense that it is possible to keep an account of these quantities in relation to their respective sources.

1.2 Units and dimensions

Units are used to describe the values of physical quantities. For example, Dr. Shreyas Mandre is 170 cm tall, weigh 10 stones and consume about 100 Watts when resting. The choice of the unit is arbitrary, as reflected in the plethora of units available for the same quantity, e.g. 170 cm 68 inches 1.8 ×1016 light-years. The SI system is popular for the purpose of standardization (and that purpose alone).

The fundamental nature of a quantity attributes it with a unique dimension. Example of dimensions are length L, mass M, time T, force F, etc. Many different units are associated with each dimension, as described above. The basic rules for analysis based on dimensions are as follows. The dimension of product of physical quantities is the product of the dimensions of the individual multiplicands. It is an error to add quantities of different dimensions, i.e. lengths may be added with lengths but not mass. For an equation in physics to be dimensionally consistent, the dimensions on both sides of the equation must be identical.

Dimensions of quantities are not necessarily independent of each other. For example, the dimensions of volume L3 depend on the dimensions of length. This dependence is best expressed in systems of dimensions. In the MLT system, the dimensions M, L and T are considered basic, from which all other dimensions are derived, e.g. F=ML/T2. There is some arbitrariness in constructing systems of dimensions, for example, equally valid is the FLT system where M=FT2/L is a derived dimension.

Exercise: Construct a ρVL system, where ρ is the dimension of density and V of velocity. Express the dimensions of mass and time in terms of ρ, V and L.

Exercise: What are the dimensions of p, μ and ν?

This arbitrariness in the nature of physical quantities underlies some of the most powerful theoretical tools in physics. It was on the basis of such considerations that Sir G. I. Taylor estimated the yield of the atomic bomb, a heavily guarded American secret at the time, from publicly available images. Ignore dimensions at your own peril!

1.3 Cartesian tensors

Vectors and tensors offer a compact representation to tackle the large number of variables whilst accounting for the relationship between them. Vectors are quantities that, in the colloquial sense, have a magnitude and a direction. Mathematically, vectors are elements of 3, where the components commonly represent Cartesian components. Tensors are generalization of vectors which have more than three components in 3-dimensions. An nth rank tensor has 3n components in 3 dimensions. In this way, a vector is a first-rank tensor. Examples of second rank tensors include moment of inertia of a solid body. Here the nine Cartesian components of the moment-of-inertia tensor are

Ixx=Vx2𝑑m, Iyy=Vy2𝑑m, Izz=Vz2𝑑m,
Iyx= Ixy=Vxy𝑑m, Izy= Iyz=Vyz𝑑m, Izx= Ixz=Vxz𝑑m,

where dm represents an infinitesimal mass element spanning the extent of the body V.

Tensors of rank zero are scalars.

Two different notations are available to represent vector and tensor expressions, they are called the vector notation and the index notation. Depending on the circumstance, one notation may be preferred to the other either for clarity or compactness.

We will use boldface symbols to represent vectors and tensors in vector notation. In this notation, it is crucial to declare the rank of the tensor a symbol represents, e.g., let 𝒖 represent first-rank velocity tensor, and let 𝑻 represent the second rank stress tensor. The operations on these symbols are understood as a matter of convention, e.g., 𝒂+𝒃 is the sum of two tensors 𝒂 and 𝒃 of equal rank. Vector notation is not equipped to unambiguously represent dot products and cross products, but it is denoted as 𝒂𝒃 and 𝒂×𝒃, respectively, when no ambiguity exists. Examples will be abundant in the rest of the notes. The advantage of this notation is that it treats tensors in their abstraction without reference to any specific coordinate system.

The index notation explicitly invokes (for this module Cartesian) coordinate system. The Cartesian coordinates are themselves written as x1=x, x2=y and x3=z. Let 𝒆^1, 𝒆^2 and 𝒆^3 be the unit vectors along x1, x2 and x3 coordinates, respectively. Then the components u1, u2 and u3 of the tensor a first rank tensor (i.e., a vector) 𝒖 are written as ui, where the index i ranges from 1 to 3. The components of a second rank tensor are represented using two indices, e.g., the components of the moment-of-inertia tensor are represented as Iij, where

Iij=Vxixj𝑑m,and i and j both independently range from 1 to 3.

It is customary to omit the range of the indices in this notation. In general, nth rank tensors have n indices.

The summation convention: Expressions, such as the dot product between two vectors are written as a sum in index notation, e.g.

𝒂𝒃=a1b1+a2b2+a3b3=i=13aibi. (1.9)

As a matter of convention due to Einstein, it is customary to also omit the summation sign over an index which appears twice in any product. Hence, 𝒂𝒃=aibi, with an implied sum over i from 1 to 3. When a sum on an index appearing twice is not implied, this must be stated explicitly. Similarly, whether a sum is implied on an index that repeats more than two times must be explicitly stated.

Two special tensors: Let us introduce two special and basic tensors. The first is called the Kronecker delta, defined as

δij={0for ij1for i=j. (1.10)

The second is the third-rank alternating tensor, also known as the Levi-Civita, which is defined as

ϵijk={0ifi=j,j=kori=k,1for(i,j,k)=(1,2,3),(2,3,1)or(3,1,2),1for(i,j,k)=(1,3,2),(2,1,3)or(3,2,1). (1.11)

Tensor algebra: Using these rules, the notation for tensor algebra is as follows. Here w is a scalar, 𝒂, 𝒃 and 𝒄 are vectors, 𝑻, 𝑺, 𝑨 and 𝑩 are second rank tensors.

Vector notation Index notation Tensorial rank
Sum of vectors: 𝒄=𝒂+𝒃 ci=ai+bi 1
Dot product of vectors: w=𝒂𝒃 w=aibi 0
Dot product of tensors: 𝑻=𝑨𝑩 Tij=AikBkj 2
Outer product of vectors 𝑨=𝒂𝒃 Aij=aibj 2
Outer product of tensors: 𝑨𝑩 AijBkm 4
Contraction: tr(𝑨) Aii 0
Scalar triple product: 𝒂(𝒃×𝒄) ϵijkaibjck 0
Vector triple product: 𝒂×(𝒃×𝒄) ϵijkϵkmnajbmcn 1

A combination of the above operations is a contraction between two second rank tensors denoted by 𝑨:𝑩 defined as AikBki. Note that 𝑨:𝑩=tr(𝑨𝑩).

A symmetric second-rank tensor is on that satisfies Sij=Sji for all combinations of i and j. An anti-symmetric one is that satisfies Aij=Aji.

Components of vectors transform in a predictable way when represented in a different coordinate system. In other words, if a vector 𝒖 has components ui, i.e. 𝒖=ui𝒆^i (note: sum implied), then the components of 𝒖 change when the basis vectors change to 𝒆^i. Similarly, tensor components also change in a predictable way under coordinate transformations. This can be seen using the representation of tensors in terms of the Cartesian basis vectors as 𝑻=Tij𝒆^i𝒆^j (sum implied). Note that 𝒆^i𝒆^j is the outer product between 𝒆^i and 𝒆^j, and just like 𝒆^i is a unit vector, the term 𝒆^i𝒆^j is one of the nine unit second-rank tensors.

The reason to introduce the Kronecker delta and the alternating tensor are that they are special second and third rank tensors. They are the so-called “unit isotropic” tensors, i.e. the Cartesian components of these tensors do not change under coordinate transformations. (We shall take it for granted here without proof.) By extension, any algebraic expressions constructed using these unit isotropic tensors are also isotropic. In fact, the most general second, third, and fourth rank isotropic tensors are

Iij2 =λδij (1.12a)
Iijk3 =λϵijk (1.12b)
Iijkl4 =λδijδkl+αδikδjl+βδilδjk, (1.12c)

where λ, α and β are arbitrary scalars.

The so-called ϵδ identity is useful in simplifying the inner products of two alternating tensors as

ϵijkϵkmn=δimδjnδinδjm. (1.13)

Here note that the left hand side is an isotropic fourth-rank tensor, which is expressed on the right-hand side as a special case of the most general isotropic fourth-rank tensor.

1.4 Tensor Calculus

Tensorial notation also aids in representing multivariate calculus of fields. Here are the basic elements in terms of scalar field ϕ(𝒙,t), vector field 𝒖(𝒙,t) and tensor fields 𝑻(𝒙,t). A special notation called the comma notation is used the denote derivatives with respect to Cartesian coordinates, ϕ/xi=ϕ,i, i.e. an index following a comma in the subscript of a symbol implies differentiation with respect to the corresponding Cartesian coordinate.

Vector notation Index notation Tensorial rank
Gradient of a scalar: ϕ ϕxi=ϕ,i 1
Divergence of a vector: 𝒖 uixi=ui,i 0
Gradient of a vector: 𝒖 uixj=ui,j 2
Divergence of a tensor: 𝑻 Tijxj=Tij,j 1
Curl of a vector: ×𝒖 ϵijkukxj=ϵijkuk,j 1

Results in integral calculus may also be condensed using the index notation. Consider a volume Ω in space, with its unit normal oriented outwards denoted by 𝒏^. The divergence theorem expressed in vector and index notation is

Ω𝒖𝒏^𝑑A=Ω𝒖𝑑ΩorΩuini𝑑A=Ωui,i𝑑Ω, (1.14)

where Ω represents the bounding surface of Ω and dA is the infinitesimal area element on it. A generalization of this theorem for second-rank tensor is

Ω𝑻𝒏^𝑑A=Ω𝑻𝑑ΩorΩTijnj𝑑A=ΩTij,j𝑑Ω. (1.15)

Stokes theorem applies to closed curves C, the unit tangent to it 𝒕^, the surface it encloses S and the unit normal to the surface 𝒏^. For vectors, it reads

C𝒖𝒕^𝑑s=S(×𝒖)𝒏^𝑑AorCuiti𝑑s=Sϵijkukjni𝑑A, (1.16)

where ds is the infinitesimal arc-length element.