Main Body

II. Mathematical Tools

In this chapter we introduce a few mathematical tools that we will use in formulating some of the analysis of fluid flow problems for both inviscid and viscous flows. We introduce the use of tensor notation which is widely used in expressing fluid mechanics governing equations. We also show some mathematical manipulations that will help to provide some physical insight into the governing equations.

Tensor Notation

Most students are very familiar with vector notation (or Gibbs notation) for describing (usually) three component vectors in fluid mechanics. Vector notation implies the existence of components of the vector. A scalar, as is known is fully described by a single number, its magnitude. Of course a vector quantity may have more than three components in general, but here we are using vectors to describe components in three dimensional space. Each component may be a function of both space and time, and as such the entire vector is a function of space and time. The number of components in the vector is determined by the “dimensionality” of the vector. For instance a velocity field may only have two orthogonal directions, say (x,y), and as such the flow becomes two dimensional — the assumption is that there is no velocity variation in the third, z, direction. We can think of the three components as a set of three independent variables that the variable in question is a function of.

Tensor notation is an alternative approach and is a very powerful way of expressing any dimensional vector, as well as what are known as higher order tensors — variables that have several sets of independent variables to be considered. We say that a scalar is a zero order tensor and a vector is a first order tensor, such as velocity. An example of a second order tensor is stress. It requires two sets of indices to define its local value. A third order tensor requires a set of three indices to be fully defined. In our three dimensional world then, a first order tensor has three components. A second order tensor has nine components, determined by the possible combinations of the three elements within each set of indices. To make this clear note that for some vector V in three dimensions:

[latex]\textbf{V}= u\textbf{i}+v\textbf{j}+w\textbf{k }={u{}_{i}}[/latex]

Where u,v, and w are the values of the three components each a function of (x,y,z), [latex]\textbf{i},\textbf{j}[/latex] and [latex]\textbf{k}[/latex] are the unit vectors associated with x,y and z coordinates respectively and [latex]u{}_{i}[/latex] is tensor notation for a vector [latex]V[/latex]. The tensor representation uses a subscript which can have three possible values, each representing one of the three components the vector, or [latex]u{}_{1} = u, u{}_{2}[/latex] = v and [latex]u{}_{3} = w[/latex]. Each component is denoted by a particular subscript that is defined based on the numerical value of the index. Here index number “1” is representative of the x coordinate, etc. But this is not limited to Cartesian coordinates. For example subscripts 1, 2 and 3 could represent r, [latex]\theta[/latex] and z in cylindrical coordinates.

Also, we can define a stress tensor as:

[latex]\tau_{ij}[/latex]

and the number of possible combinations of the subscripts i,j are nine. Do not confuse these subscripts with the unit vectors since they can take on one of the possible three directions. For fluid flow problems we will define the stress tensor and use it to arrive at our conservation of momentum equation. A third order tensor could be written as [latex]T{}_{ijk}[/latex] and it would have 27 elements (3x3x3). As stated above, a single subscript denotes a vector and has 3 possible values, one for each coordinate.

It is possible to write all of our vector equations using tensors instead of the vector notation. This has some advantages and is used widely in fluid mechanics. Table 2.1 illustrates a number of examples of the equivalence of vector and tensor notation. A few of these are worth discussing and we will be using them as we write out the various equations that we will be formulating and using. It is suggested that you examine this notation and look at the various operators listed.

Operations among vectors carry with them some special consideration and rules. The gradient operator is [latex]\mathrm{\nabla }[/latex] in vector notation. If the quantity of interest is the gradient of scalar “a”, written as [latex]\mathrm{\nabla }\mathrm{a\ }[/latex], then this operation results in a vector whose components are partial derivatives in each of the three coordinates, [latex]\frac{\partial a}{\partial x}[/latex] is the x component, etc. In tensor operation this is expressed as [latex]\frac{\partial a}{\partial x_i}[/latex] which is a vector whose components correspond to each of the three possible values of the index “[latex]{i}[/latex]”. If the gradient is being taken of a vector, say [latex]{u{}_{j}}[/latex], the the resulting gradient expression is a second order tensor with the two indices, [latex]{i }[/latex] for the partial derivatives and [latex]\textit{j}[/latex] for the vector quantity [latex]{u{}_{j}}[/latex]. This is then written as [latex]\frac{\partial u_j}{\partial x_i}[/latex]. Notice that the gradient operator index comes first since it operates on the vector [latex]{u{}_{j}}[/latex]. There are nine combinations of [latex]{i,j}[/latex] for this second order tensor. We will see this tensor in our discussions of viscous flows and vorticity.

Table 2.1 Listing of a few vector, tensors or scalars in vector and tensor notation
Examples of Vector VS Tensor Notation
Quality Vector Notation Tensor Notation
Scalar a a
Vector [latex]\bar{a}[/latex] [latex]a_i[/latex]
2nd order tensor [latex]\overline{\overline{a}}[/latex]

[latex]a_{ij}[/latex]
dot product (scalar) [latex]\bar{a}\cdot\bar{b}[/latex] [latex]a_{i} b_{i}[/latex]
cross product (vector) [latex]\bar{a}\times\bar{b}[/latex] [latex]\varepsilon_{ijk}a_{j}b_{k}[/latex]
del operator (vector) [latex]\bigtriangledown[/latex] [latex]\frac{\partial}{\partial x_i}[/latex]
gradient (vector; tensor) [latex]\bigtriangledown a;\bigtriangledown\bar{a}[/latex] [latex]\frac{\partial{a}}{\partial {x_{i}}};\frac{\partial{a_{i}}}{\partial{x_{j}}}[/latex]
divergence (scalar) [latex]\bigtriangledown\cdot\bar{a}[/latex] [latex]\frac{\partial{a_{i}}}{\partial{x_{i}}}[/latex]
curl (vector) [latex]\bigtriangledown\times\bar{a}[/latex] [latex]\varepsilon_{ijk}\frac{\partial{a_k}}{\partial{x_j}}[/latex]
Laplace operator (scalar) [latex]\bigtriangledown^2[/latex] [latex]\frac{\partial^2}{\partial{x_i}\partial{x_i}}[/latex]
divergence of a 2nd order tensor [latex]\bigtriangledown\cdot \overline{\overline{a}}[/latex] [latex]\frac{\partial{a_y}}{\partial{x_i}}(vector)[/latex]
dot product: vector & del [latex]\bar{a}\cdot\bigtriangledown[/latex] [latex]a_j\frac{\partial}{\partial{x_j}}(scalar)[/latex]
dot product: vector & gradient [latex]\bar{a}\cdot\bigtriangledown\bar{b}[/latex] [latex]a_j\frac{\partial{b_i}}{\partial{x_j}}(vector)[/latex]
Material Derivative [latex]\frac{\partial}{\partial{t}}+\bar{u}\cdot\bigtriangledown or \frac{D}{Dt}[/latex]
Vorticity: (pseudovector) [latex]\omega_k=\bigtriangledown x\bar{u}=\left (\frac{\partial{u_i}} {\partial{u_j}}-\frac{\partial{x_j}}{\partial{x_i}} \right )=\varepsilon_{ijk}\frac{\partial {u_k}}{\partial{x_j}}[/latex]
Rate of strain tensor: (2nd order tensor) [latex]\frac{1}{2}\left (\frac{\partial{u_i}} {\partial{u_j}}+\frac{\partial{u_j}}{\partial{u_i}} \right )[/latex]

Consider now the divergence operation, in vector notation written as [latex]\mathrm{\nabla }\mathrm{\cdot }\mathrm{V}[/latex]. This operation is the dot, or scalar, product between the vector operator [latex]\mathrm{\nabla }[/latex] and vector V. The result is a scalar. The tensor notation for this operation is [latex]\frac{\partial u_i}{\partial x_i}[/latex] where the indices for the vector [latex]{u{}_{i}}[/latex] are the same as the index for the partial derivative operator. This implies that the same component of each vector, [latex]\mathrm{\nabla }[/latex] and V, should be used simultaneously. The rule in tensor notation then is to add all of these component operations to result in one final scalar. This becomes for the three coordinates, [latex]x{}_{1}, x{}_{2}[/latex] and [latex]x{}_{3}{}_{,\ }[/latex] which in Cartesian coordinates are x, y and z: [latex]\frac{\partial u_i}{\partial x_i}=\frac{\partial u_1}{\partial x_1}+\frac{\partial u_2}{\partial x_2}+\frac{\partial u_3\ }{\partial x_3}\ =\ \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}[/latex]

The result is a scalar formed from the sum of each of the partial derivatives. As is apparent from this equation, when a term has a “repeated” index (such as the index “[latex]{i}[/latex]” in the above equation for both u and x) then the rule is to sum over all values of the index. This applies to a particular term in an equation.As an additional example consider the term [latex]\frac{\partial a_{ij}}{\partial x_i}[/latex]in this case the index “[latex]{i}[/latex]” is repeated and j is left as a “free index” (which could take on any of the possible three coordinate values.) Consequently this expression is a vector, with one free index, whose components are represented by the value of “[latex]{j}[/latex]”.There is also an operator, [latex]{\varepsilon }_{ijk}[/latex]. This is given the name “permutation operator”. Notice in Table 2.1 this is listed as used in the cross product and curl operation (the curl is the cross product of the vector operator [latex]\mathrm{\nabla }[/latex], which in tensor notation is written as [latex]\frac{\partial }{\partial x_j}[/latex], and the vector [latex]{u{}_{k}})[/latex]. The permutation operator definition results in the cross product between two vectors. For this each of the indices can take on three different values, and there are 27 different possible combinations. The rules associated with the different values are as follows:

  • if any of the three indices have the same value then [latex]{\varepsilon }_{ijk}{}_{\ }= 0[/latex]
  • if the order of the indices is cyclic (1,2,3 or 2,3,1 or 3,1,2) then [latex]{\varepsilon }_{ijk} ={\ }1[/latex]
  • if the order of the indices is anti-cyclic (1,3,2 or 2,1,3 or 3,2,1) then [latex]{\varepsilon }_{ijk\ }={ }_{\ }-1[/latex]

The result is three values for the permutation operator, with values being +1, -1, or 0. Consequently, the following occurs for the cross product: [latex]{a{}_{j}}[/latex] x [latex]{b{}_{k}}[/latex][latex]a_j\times b_k={\varepsilon }_{ijk}a_jb_k=\left(a_2b_3-a_3b_2\right)\boldsymbol{i}+\left(a_3b_1-a_1b_3\right)\boldsymbol{j}\boldsymbol{+}\left(a_1b_2-a_2b_1\right)\boldsymbol{k}[/latex]

The reader should perform this operation for the curl of a vector where [latex]a_i[/latex] is replaced by [latex]\frac{\partial }{\partial x_j}[/latex] and [latex]{b{}_{k}}[/latex] is replaced by [latex]{u{}_{k}}[/latex]

There is another often used operator in tensor notation called the Kronecker delta, [latex]{\delta}_{ij}[/latex]. The rule used is that when [latex]{i=j}[/latex] the value of [latex]{\delta}_{ij}[/latex] is one, otherwise it is set equal to zero, that is:

[latex]{\delta }_{ij}=1\ for\ i=j\\[/latex]

[latex]{\delta }_{ij}=0\ for\ i\neq j\\[/latex]

This comes in handy when one wants to include terms where there may be a nonzero value only when [latex]{i=j}[/latex]. We will see this in our expression in the Navier-Stokes equation for the pressure contribution to the total force on a fluid element. For example suppose we want to include a gradient of a scalar, [latex]{S}[/latex], as [latex]\frac{\partial S}{\partial x_j}[/latex] in an equation, note that this is a vector quantity. However, if we are evaluating different components of the gradient we only want to include those components that align with the vector component “[latex]{i}[/latex]”, then we would write this term [latex]\frac{\partial S}{\partial x_j} {{}δ_{ij}}[/latex]. In doing this we see that when we evaluate this term we only include the gradient component in the [latex]{i}[/latex] direction, since when operating on the gradient with [latex]{\delta}_{ij}[/latex] results in a nonzero possible value only when [latex]{i=j}[/latex]. We will use this for the pressure term in the Navier-Stokes equation as is explained later. We will use gradients, divergence and curl operations in the expressions used in the Navier-Stokes equations. The reader should become familiar with these.

Gradient, Divergence and Curl Operators

The gradient is a vector operator that for instance, when operated on a scalar produces a vector. It takes the partial derivative with respect to each component of the coordinate system. It is denoted as [latex]\bigtriangledown[/latex] and can be expressed in Cartesian coordinates as:

[latex]\bigtriangledown ( ) = \frac {\partial ( )} {\partial x}i + \frac {\partial ( )} {\partial y}j + \frac {\partial ( )} {\partial z} k[/latex]

Where a scalar of a vector or tensor could be inserted into the parenthesis. A gradient of a vector then has 9 components, where you take x,y,z derivatives of each of the three vector components. This would form the secord order tensor:

[latex]\bigtriangledown V =\begin{bmatrix} \frac{\partial u}{\partial x}&\frac{\partial u}{\partial y} &\frac{\partial u}{\partial z} \\ \frac{\partial v}{\partial x}&\frac{\partial v}{\partial y} &\frac{\partial v}{\partial z} \\ \frac{\partial w}{\partial x}&\frac{\partial w}{\partial y} &\frac{\partial w}{\partial z} \end{bmatrix}[/latex]

If one takes the dot product of [latex]\bigtriangledown[/latex] with a vector we get the divergence of that vector which is a scaler:

[latex]\bigtriangledown\cdot V = \frac {\partial u} {\partial x} + \frac {\partial v} {\partial y} + \frac {\partial w} {\partial z}[/latex]

Note that the divergence of the second order tensor yields a vector, such as [latex]\bigtriangledown\cdot\bigtriangledown{V}[/latex]. The x component of this vector is:

[latex]\bigtriangledown\cdot{\bigtriangledown u} = \left ( \frac{\partial^2u}{\partial{x^2}}+\frac{\partial^2u}{\partial{y^2}}+\frac{\partial^2u}{\partial{z^2}} \right )i[/latex]

The y and z components should be apparent where the x component of velocity is replaced by the y and z components of velocity in the numerators for the y and z components of the resultant vector.

 

Vector Identities

Vector Identity

Below are shown a few vector identities that are used in fluid mechanics in the manipulation of the governing equations of motion. A,B are vectors, [latex]phi[/latex], is a scalar.

[latex]\bigtriangledown\left ( A \cdot B \right ) = \left (A\cdot\bigtriangledown \right )B+\left ( B\cdot\bigtriangledown \right )A+A\cdot\left ( \bigtriangledown\cdot B \right )+ B \cdot\left ( \bigtriangledown\cdot A \right ) = A\bigtriangledown{B} + B\bigtriangledown{A}[/latex]

[latex]\frac{1}{2}\bigtriangledown\left ( A\cdot{A} \right )=\left ( A\cdot\bigtriangledown \right )A+A\cdot\left ( \bigtriangledown + A \right ) = A\bigtriangledown{A}[/latex]

[latex]\bigtriangledown\cdot\left ( \bigtriangledown\phi \right )=0[/latex]

[latex]\bigtriangledown\cdot\left ( \bigtriangledown\cdot{A} \right )=0[/latex]

[latex]\bigtriangledown^2\phi=\bigtriangledown\cdot\left ( \bigtriangledown\phi \right )[/latex]

[latex]\bigtriangledown\cdot\left ( \bigtriangledown\cdot A \right )=\bigtriangledown\left ( \bigtriangledown\cdot A \right )-\bigtriangledown^2A[/latex]

Divergence Theorem: What is the name of this symbol? [latex]B\cdot{d}A = \int\int\int\left ( \bigtriangledown\cdot B \right )d V[/latex] (where [latex]dA[/latex] is an elemental area around the boundary of volume [latex]V[/latex])

Stokes’ Theorem: [latex]\oint B\cdot d l = \int \int \left ( \bigtriangledown\cdot B \right )\cdot d A[/latex] (where [latex]dl[/latex] is an element of the line integral around the boundary of area [latex]A[/latex]).

There are a number of vector identities and manipulations that are very useful in fluid mechanics. We list a few of these here for future reference. Unless stated otherwise these are valid for any vector quantity, V.

The reader may want to visit the web site: https://en.wikipedia.org/wiki/Vector_calculus_identities

Material Derivative

The material derivative is a shorthand notation expressing the rate change of some property or parameter, which could be a scalar or a tensor, such as temperature, pressure or momentum (mass times velocity). This time derivative is based on a Lagrangian frame of reference of a fixed mass of fluid. So in a sense one is following a fixed mass of fluid and describing how the given property changes. Since the fluid mass can move we say that the fixed mass depends on both space and time. For example if we pick the scalar pressure, [latex]{P}[/latex], and say that our fixed mass has a pressure that is a function of space and time using Cartesian coordinates we have for the variable pressure:

[latex]{P }= {f}(x,y,z,t\textit{)}[/latex]

Then if we want its time rate of change for a fixed mass at a given time we use the notation that [latex]{DP}[/latex] is the change of quantity [latex]{P}[/latex] and [latex]{DP/Dt}[/latex] is its change over time accounting for changes in both space and time (note we apply the chain rule for each of the possible coordinate directions:

[latex]\frac{DP}{Dt}=\frac{\partial P}{\partial t}+\frac{\partial P}{\partial x}\frac{dx}{dt}+\frac{\partial P}{\partial y}\frac{dy}{dt}+\frac{\partial P}{\partial z}\frac{dz}{dt}[/latex]

Noting that [latex]\frac{dx}{dt}=u[/latex] and similarly for each coordinate direction, the equation becomes:

[latex]\frac{DP}{Dt}=\frac{\partial P}{\partial t}+u\frac{\partial P}{\partial x}+v\frac{\partial P}{\partial y}+w\frac{\partial P}{\partial z}[/latex]

(2.1)

Things to pay attention to in this equation are (i) the velocity components [latex]{u,v}[/latex] and [latex]{w}[/latex] are not vector quantities, but components of the velocity vector, V; (ii) the time variation is specified by the first term and is evaluated at a particular point in space (as are all of the other terms) but all other terms may in fact be time dependent; this equation is valid for scalars or tensors, and the tensor sense is determined by the quantity being evaluated (in Equation (2.1) it is a scalar since [latex]{P}[/latex] is a scalar). The material derivative of the velocity vector V, is:

[latex]\frac{DV}{dt}=\frac{\partial V}{\partial t}+u\frac{\partial V}{\partial x}+v\frac{\partial V}{\partial y}+w\frac{\partial V}{\partial z}[/latex]

(2.2)

Since [latex]{V}[/latex] is a vector then [latex]{DV/Dt}[/latex] is also a vector, and each term on the right hand side is a vector and the components of each of these vectors is determined by inserting the component for the vector [latex]{V}[/latex]. In other words the x component is:

[latex]\frac{Du}{Dt}=\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}+w\frac{\partial u}{\partial z}[/latex]

(2.3)

The above is for Cartesian coordinates. This can be written for any selected coordinate system. For example the result in cylindrical coordinates in [latex]{r,\theta}[/latex] is shown below.

Material Derivative in Cylindrical Coordinates

We would like to find an expression for DV/DT in cylindrical coordinates that we can use to help interpret streamline coordinates. We will only examine a two dimensional situation, [latex]r, \theta[/latex] since [latex]z[/latex] is similar to Cartesian coordinates.

Given the vector [latex]V = u_r\hat{r}+u_\theta\hat{\theta}[/latex] we write this as a material derivative:

[latex]\frac{DV}{Dt}=\frac{\partial V}{\partial t}+u_r\frac{\partial V}{\partial r}+\frac{u_\theta}{r}\frac{\partial V}{\partial \theta}[/latex]

When the angular position changes then there can be a change in both the unit vectors [latex]\hat{r}[/latex] and [latex]\hat{θ}[/latex] since their orientation changes. So we can write the following:

[latex]\frac{\partial V}{\partial\theta}=\frac{\partial\left ( u_r\hat{r}+u_\theta\hat{\theta} \right )}{\partial\theta}=\hat{r}\frac{\partial u_r}{\partial\theta}+u_r\frac{\partial\hat{r}}{\partial\theta}+\hat{\theta}\frac{\partial u_\theta}{\partial\theta}+u_\theta\frac{\partial\hat{\theta}}{\partial\theta}[/latex]

Now we must interpret the derivatives of the unit vectors with respect to θ. Refer to the figure below showing a blob of fluid moving [latex]\partial\theta[/latex] that results in a change of the unit vector [latex]\hat{r}[/latex], that we can write as: [latex]d\hat{r}=\hat{\left | r \right |}d\theta\hat{\theta}=d\theta\hat{\theta}[/latex], which we can write as: [latex]\frac{\partial\hat{r}}{\partial\theta}=\hat{\theta}[/latex].

Similarly the change in the unit vector, [latex]d\hat{\theta}=-rd\theta[/latex], as illustrated above on the right as well. This is negative because if the change in the unit vector is counterclockwise the change in the r unit vector is in the negative r direction. The result is [latex]\frac{\partial\hat{\theta}}{\partial\theta}=-r[/latex].

All of this can now be combined and used in the material derivative above resulting in:

[latex]\frac{DV}{Dt}=\frac{\partial V}{\partial t}+u_r\frac{\partial V}{\partial r}+\frac{u_\theta\partial V}{r\partial\theta} = \left ( \frac{\partial u_r}{\partial t}+u_r\frac{\partial u_r}{\partial r}+\frac{\partial_\theta\partial u_r}{r\partial\theta}-\frac{u_\theta^2}{r} \right )\hat{r}+\left ( \frac{\partial u_\theta}{\partial t}+u_r\frac{\partial u_\theta}{\partial r}+\frac{u_\theta\partial u_\theta}{r \partial\theta}+\frac{u_r u_\theta}{r} \right )\hat{\theta}[/latex]

This applies to Cylindrical Coordinates, but as we shall see helps to interpret streamline coordinates when there is curvature to the streamline. The net result is added terms that account for the change of the unit vectors due to the curving nature of the flow.

Acceleration

If one finds the material derivative of the velocity vector of a fixed mass of fluid, namely Eqn. 2.2, the result is the acceleration of that fixed mass of fluid in time and space. The general expression for this vector quantity can be written in tensor notation as:

[latex]\frac{Du_i}{Dt}=\frac{\partial u_i}{\partial t}+u_j\frac{\partial u_i}{\partial x_j}[/latex]

 (2.4)

Note here that the vector, V is replaced with the tensor notation for a vector written as [latex]{u{}_{i}}[/latex]. Also the summation rule has been used in the second term on the right hand side, this implies that there are really three terms being represented by the one term shown. The first term on the right hand side is denoted as the “local acceleration” and physically represents the rate change of velocity at a fixed point in space. The second expression on the right hand side is denoted as the “convective acceleration” and physically represents how the velocity field changes over space. If a flow is steady the local acceleration is by definition zero. However a fluid mass as it moves through space where the velocity changes in space will experience acceleration (or deceleration). An example is steady flow through a nozzle, or diffuser, where the fluid velocity will increase, or decrease, along the flow direction. For these two cases a fixed mass of fluid will accelerate, or decelerate due to moving in space. The reader is reminded that the components of the acceleration depend on the component of velocity in the derivatives.The convective acceleration can be rewritten using a vector identity listed above. Namely, the convective acceleration in vector notation is [latex]\left(V\cdot \mathrm{\nabla }\right)V[/latex] and can be expressed as:[latex]\left(V\cdot \mathrm{\nabla }\right)V=\frac{1}{2}\mathrm{\nabla }\left(\mathrm{V}\mathrm{\cdot }\mathrm{V}\right)\mathrm{-}\mathrm{V\times (}\mathrm{\nabla }\mathrm{\times V)}[/latex]. In tensor notation this becomes:

[latex]u_j\frac{\partial u_i}{\partial x_j}=\frac{1}{2}\frac{\partial \left(u_ju_j\right)}{\partial x_j}-{\varepsilon }_{ijk}\left(u_j{\varepsilon }_{klm}\frac{\partial u_m}{\partial x_l}\right)[/latex]

(2.5)

Consequently, the convective acceleration can be recast as the expression on the right hand side. Although this may not seem like it is of any benefit and just makes things more complicated we will show that it does indeed have advantages. The reader should pay attention to the subscripts in the above expression. Notice that each term is a vector with the free index of “[latex]{i}[/latex]”, all other indices are repeated and therefore summed. In the summation of vectors we basically need to sum components, this implies that they have the same index, in this case “[latex]{i}[/latex]”.

Streamline Coordinates

Streamline Coordinates define the flow direction of the velocity vector. It is defined as the locus of points tangent to the velocity vector at some instant in time. If the flow is steady these lines do not change. For our purposes we will consider steady flow here, but that is not necessary. Fig. 2.1 illustrates the streamline as well as a streamline coordinate system we will define. Our discussion for simplicity will be two dimensional such that flow is within the [latex]{s-n}[/latex] plane. We retain an orthogonal coordinate system with “s” along the streamline and “n” normal to the streamline. The third coordinate would be normal to both [latex]{s}[/latex] and [latex]{n}[/latex], into the page.

If we take the cylindrical coordinate system as an example and consider flow along a curve, with some non-infinite radius of curvature, then we expect additional terms to appear in the acceleration. These terms account for the changing unit vectors as discussed in the material derivative in cylindrical coordinates. We make note of the fact that the velocity in the [latex]{n}[/latex] direction is zero, by definition of the streamline. However, there can still be acceleration in the [latex]{n}[/latex] direction since the unit vectors are changing. Using the cylindrical coordinate expression as a guide and denoting [latex]{r}[/latex] as normal to the streamline direction, [latex]{n}[/latex], and [latex]{θ}[/latex] as along the streamline direction, [latex]{s}[/latex], then the results for acceleration in the [latex]{s}[/latex] and [latex]{n}[/latex] directions are:

[latex]a_s=\frac{\partial u_s}{\partial t}+u_s\frac{\partial u_s}{\partial s}=\frac{\partial u_s}{\partial t}+\frac{\partial {u_s}^2}{2}[/latex]

(2.6)

[latex]a_n=\frac{{u_s}^2}{{R}}[/latex]

Here [latex]{R}[/latex] is the local value of the radius of curvature of the streamline. The value of [latex]{R}[/latex] goes to infinity for a straight streamline and there is no acceleration in the [latex]{n}[/latex] direction. Recall here that [latex]{u{}_{r}}[/latex] (or [latex]{u{}_{n}}[/latex]) is identically zero at the streamline and that the expression for [latex]a_n[/latex]only contains the last term associated with the changing unit vector [latex]{θ}[/latex], or in this case [latex]{s}[/latex]. Also we define in the streamline coordinates the direction [latex]{n}[/latex] to be inward towards the radius of curvature as shown in Fig. 2.1. This is opposite to the convention used in cylindrical coordinates where [latex]{r}[/latex] is radially outward, hence the sign is positive for the value of [latex]a_n[/latex]. That is, if there is a curving streamline the acceleration is inward towards the direction of the center of the radius of curvature.

 

 

Fig. 2.1 Streamline coordinate illustration; [latex]{s}[/latex] is along the velocity vector direction and tangent to the streamline and [latex]{n}[/latex] is normal to [latex]{s}[/latex] while directed towards the center of radius, denoted as [latex]{R}[/latex].

 

 

Velocity Gradient Tensor

The velocity gradient tensor describes how the velocity varies near a specified location within the flow field. It is represented as [latex]\frac{\partial u_j}{\partial x_i}[/latex] which is a second order tensor, and therefor has nine components in three dimensional space. A second order tensor can be decomposed into symmetric and antisymmetric parts given as follows:

[latex]\frac{\partial u_j}{\partial x_i}=\frac{1}{2}\left(\frac{\partial u_j}{\partial x_i}\mathrm{\ }\mathrm{+}\frac{\partial u_i}{\partial x_j}\mathrm{\ }\right)+\frac{1}{2}\left(\frac{\partial u_j}{\partial x_i}\mathrm{\ }\mathrm{-}\frac{\partial u_i}{\partial x_j}\mathrm{\ }\right)[/latex]

(2.7)

Notice here that the first term adds half of the transpose of the original tensor and the second term subtracts the same quantity completing this identity. The first term is symmetric that is by inverting the [latex]{i}[/latex] and [latex]{j}[/latex] subscripts there is no change in the value of the element. The second term is antisymmetric such that inverting the subscripts the value of the element has the opposite sign. The ½ multiplier is included in the definition of the symmetric and antisymmetric parts such that their sum results in the original tensor.

In fluid mechanics applications the velocity gradient tensor is a measure of the velocity changes experienced in the fluid flow as one moves away (infinitesimally) from a given point. That is to say this tensor is a point function, having a value at any given point in the flow. We refer to that as a local function in space. It may be time dependent as well. So we can say that this tensor has a symmetric part and an antisymmetric part. As we will see when we discuss the viscous forces within a fluid the symmetric part is representative of the deformation rate, or strain rate, experienced by a fluid element. The antisymmetric part represents the rotation rate of a fluid element. If one examines the curl operator on [latex]{u{}_{j}}[/latex] given previously as [latex]{\varepsilon }_{kij}\frac{\partial u_j}{\partial x_i}[/latex], where the subscripts indicate that this is the [latex]{k}th[/latex] element, it can be shown by writing out the terms of the curl operation on [latex]{u{}_{j}}[/latex] that this is equivalent to twice the antisymmetric part of [latex]\frac{\partial u_j}{\partial x_i}[/latex] (or the part in parentheses of the second term in Eqn. (2.7)). Therefore the rotation part of the velocity gradient tensor is one half the vorticity of the flow. To be clear, the antisymmetric part of [latex]\frac{\partial u_j}{\partial x_i}[/latex] is a second order tensor, but from observation if i=j then the value of this is identically zero (all zero components along the diagonal of the tensor). The off diagonal terms are six in total but those with inverted indices are the negatives of the noninverted components ([latex]\frac{\partial u_j}{\partial x_i}=-\frac{\partial u_i}{\partial x_j}[/latex]). Consequently there are only three independent values, so we identify this second order tensor as a pseudo-vector, having three components. It is a so called pseudo-vector since the sign is arbitrary. In other words do we set [latex]\frac{\partial u_1}{\partial x_2}[/latex] as a positive or negative quantity? The sign convention is typically selected so that a counter clockwise rotation is a positive value. This means that for [latex]\frac{\partial u_1}{\partial x_2}[/latex] being positive then as [latex]{u{}_{1}}[/latex] increases in the positive [latex]{x{}_{2}}[/latex] direction then the flow will tend to rotate clockwise and is a negative rotation direction. This is illustrated in Fig. 2.2. Note that if one considers a negative value of [latex]\frac{\partial u_1}{\partial x_2}[/latex] coupled with a positive value of [latex]\frac{\partial u_2}{\partial x_1}[/latex] such that

[latex]{\omega_3}=\left (\frac {\partial u_2}{\partial x_1} -\frac {\partial u_1}{\partial x_2} \right)[/latex] > [latex]0[/latex]

(2.8)

then rotation (or vorticity) about the 3 axis (out of the page in Fig. 2.2) is counterclockwise. Obviously, the signs of the two components of the velocity gradient that go into any of the vorticity components can be positive and/or negative. If in the above example for Eqn. 2.8 the magnitudes of [latex]\frac{\partial u_2}{\partial x_1}[/latex] and [latex]\frac{\partial u_1}{\partial x_2}[/latex] are equal and their signs are the same then the vorticity component [latex]{\omega }{_3}[/latex] will be zero. One can think of the net vorticity or rotation rate vector component to be the combined value from both velocity derivative components, which contribute to rotation about an axis. The same interpretation of the combined effects of the two velocity gradient elements can be made for the strain rate, or the symmetric part of the velocity gradient tensor. When the rotation rates of the vertical and horizontal axes are identical (both elements causing, say, counterclockwise rotation at the same rate) then the strain rate will be zero (the two elements in the symmetric part cancel) and the flow is in pure rotation, as a solid body.

 

Fig. 2.2 Illustration of rotation rate as determined by the velocity gradient tensor components; the original square element has its sides experiencing rotation at a rate of [latex]\frac{\partial u_2}{\partial x_1}[/latex] for the horizontal axis and [latex]\frac{\partial u_1}{\partial x_2}[/latex] for the vertical axis.

Basic Equations

The basic equations generally used in inviscid fluid mechanics are shown in Fig. 2.3. The viscous terms we will add later. But this chart shows various forms of the governing momentum equation which include pressure and gravitational forces contributing to the overall acceleration of the fluid. At the top is the basic Cauchy form of the momentum equation, where viscous effects would be contained in the stress term, [latex]τ[/latex]. The inclusion of viscous terms will be developed in a later chapter. At this point the stress only contains the pressure forces, which is compressive and normal. The streamline coordinate formulation is shown on the right where integration is along the streamline. The introduction of the variable B is a consequence of using the last of the vector identities listed above as equation (2.5) for the convective acceleration term, and then combining the terms shown in the bottom box on the left. Here B is essentially the Bernoulli constant for steady, compressible flow. This is discussed further in the development of the Bernoulli equation. At this point the equations are presented for reference later on when the various terms are developed.

Fig 2.3 Chart Illustrating the basic governing equation in inviscid flow.

 

License

Share This Book