Fluid Mechanics for Mechanical Engineers/Scalar, Vectors and Tensors

From testwiki
Jump to navigation Jump to search

>back to Chapters of Fluid Mechanics for Mechanical Engineers

Motivation

Each scientific field needs an optimized language. Without it, this field can not develop. For example, calculus could not have developed as easily with Roman numerals:

Template:Center topXXVIII*XXIV=?Template:Center bottom

Without Arabic numbers and the decimal system, even basic operations would not be possible.

Fluid Mechanics deals with the changes of flow and fluid variables in three space dimensions and in time. In order to overcome this complexity, there is also a need for an optimized language in the field of fluid mechanics. Unfortunately, for the same vector operation there are multiple mathematical notations: Template:Center top div a=aorgrad a= a and for a vector a, a,a¯, a¯ or a~ Template:Center bottom Symbols are not unique and somewhat confusing, therefore a unified language accelerates the development of the field.

In this chapter we will be introduced to the fascinating Tensor, a unified mathematical concept which is very important to computational fluid dynamics and other scientific fields. This chapter will give a gentle introduction, but readers are encouraged to read more from the references. Due to this short introduction, Tensors will not be used as a separate chapter.

Tensor

The concept of tensor could be easily understood if we want to imagine a function which takes vector as input and do some operation with it and then give another vector as output. For example, let's assume a vector which is defined in Euclidean reference co-ordinate system. Now for some unknown reason, the co-ordinate system changed via rotation or translation. To find our previous vector defined in the new co ordinate system, we can assume a magical operator or function which will take the initial components of the vector in old co-ordinate system as input and give us back the new vector component in new co-ordinate system as output. In other way, New position vector= f(Old position vector). It is found that if this operator or function is linear, then this function or operator would be a 2nd rank tensor.

A second rank tensor looks like 3×3 matrix. For example,
|153294618| is a 3×3 matrix. If we now multiply this matrix with another matrix with 1×3 matrix we will get as an output, 3×1 matrix where the element of these matrices could be imagined as the vector components of initial and new vector.

So we can say that a 2nd rank Tensor is capable to keep the map or information of how a vector will change between new to old co-ordinate system.

Now, we have small introduction to tensor, we can present the definition of tensor. A tensor 𝑨 is a linear transformation from a vector space 𝒱 to 𝒱.[1] Thus, we can write:

𝑨:𝐮𝒱𝐯𝒱

which means tensor 𝑨 operates on vector 𝐮 to give new vector 𝐯.

This concept is found to be extremely useful where all computational calculations have to be performed in somehow matrix format. A very nice illustration of use of tensor is shown here[2]. Very common examples of Tensor are : Deformation Gradient Tensor, Stress Tensor, Einstein's famous stress-energy tensor. Even dot and cross product of two vectors are also tensor.

Although the concept of linear operation is very useful to understand the application of tensor at first time,it's not all. Let's consider about a cylinder which is stretched along it's length and due to poisson's contraction, we will see contraction in radial direction to keep the volume constant. So we can say that the deformation along different directions is somehow related on the stress applied along different directions.It appears that this relationship is a tensor which is called Stiffness Tensor and this relationship forms like :[Strain Tensor(written in column vector and reduced form)]=[Stiffness Tensor,C][Stress Tensor(same as Strain Tensor)]. The relationship for isotropic material will look like below:


[ε11ε22ε332ε232ε312ε12]=[ε11ε22ε33γ23γ31γ12]=1E[1νν000ν1ν000νν10000002(1+ν)0000002(1+ν)0000002(1+ν)][σ11σ22σ33σ23σ31σ12]

This above matrix looks differently for an anisotropic material[3]. So in other way, we can say this particular tensor can help to express how the Stress Tensor will be influencing Strain Tensor and this 'how' is a signature of type of material which is point of interest. Tensor could be envisioned in more easy way. Let us consider a vector which has components of ax1,ay1,az1. This components are essential to express this vector precisely in a given co-ordinate system.These components together will give value of the vector by ax12+ay12+az12. If this co-ordinate system change also, the vector will be represented with 3 components with different value ax2,ay2,az2 but ax22+ay22+az22 would give the same value as previous which means that something remains as identity of the vector. The components changed keeping something constant. Now lets consider about tensor. Imagine we form an entity where we have 3 vectors as components(meaning actually it has 9 or 3×3 components of vectors now) and they are somehow coupled.This means that when the co-ordinate system will change, the individual components of vector(hence the entity) will be determined keeping 'some other value' constant like ax,ay,az in the case of vector. Of course, this 'other value' has to be calculated by using 9 components of that entity. After transformation, this new 9 components will give the same 'other value'. It would be shown later on actually that vector itself is 1st order tensor.

Scalar: Tensor of 0th order

We may associate some properties of fluid or flow at one point without any preference of direction. Template:Center topP=1.01×105magnitude [Nm2]unit,T=293[K] and ρ=1.01×103[kgm3]Template:Center bottom These are some examples of such fluid variables.

Summation can be only done between variables having the same unit, i.e. Template:Center topa[Pa]+b[Pa]=a+b[Pa]Template:Center bottomTemplate:Center topa[Pa]+b[K]= ?Template:Center bottom

Multiplication is valid for variables having different units. Template:Center topa[Pa]b[K]=ab[PaK]Template:Center bottom

Products of variables having different units, have a new physical meaning. Template:Center topm=ρVTemplate:Center bottomTemplate:Center top[kg]=[kgm3][m3]Template:Center bottom

Power loss (Energy loss per unit time):

File:Scalars Vectors Tensors1.png
Pressure drop measurement in a pipe

Template:Center topΔPV˙=|ΔP||V˙|[Nm2][m3s]Template:Center bottom Template:Center topΔPV˙=|ΔP||V˙|[Nms]Template:Center bottomTemplate:Center topΔPV˙=|ΔP||V˙|[Watt]Template:Center bottom Watt is the unit for power.

Vectors: 1st order tensors

File:Scalars Vectors Tensors2.png
A vector

Variables of flow, associated with a point and a direction, like force and velocity, are vectors.

The direction of a vector should be specified in relation with respect to a given frame of reference. This frame of reference is arbitrary as the units. Let us use the Cartesian coordinate system having 3 mutually orthogonal axes.

U1,U2,U3

are components of this vector and can differ with respect to selected coordinate system but the amplitude not. Components allow us to reconstruct the vector in a particular system of reference. One should distinguish vector as an entity from its components.

File:Scalars Vectors Tensors3.png
Components of a vector

Repeated indices in Sums

Summation of variables, identified with indices can be shortly written with the summation symbol:

Template:Center topa1x1+a2x2+a3x3+...+anxn=i=1naixiTemplate:Center bottom

In short, we will write aixi

In other words, repeated indices means series sum over this index. Thus, Template:Center topU=Ui eiTemplate:Center bottom Template:Center topU=|U|cosαi eiTemplate:Center bottom

Now we can introduce the component representation of a vector U: Template:Center topUi=|U|cosαiTemplate:Center bottom

Let's look some examples: Template:Center topaibi=a1b1+a2b2+a3b3(i,j=1,2,3)Template:Center bottom whereas Template:Center topaibjTemplate:Center bottom indicates no summation over i and j and Template:Center topaijbkTemplate:Center bottom indicates also no summation. Other examples of summation are: Template:Center topaiibj=a11bj+a22bj+a33bjTemplate:Center bottom Template:Center topaijbj=ai1b1+ai2b2+ai3b3Template:Center bottom

Dummy and free index

In the following expression Template:Center topaijbj=aikbk (i,j,k=1 ... n)Template:Center bottom the repeating index j is also a dummy index, i.e. can be replaced by k. The following equality Template:Center topaijbj=akjbjTemplate:Center bottom is valid only when i=k. Then, i is called free index

yi=airxr for i,r=1,2,3

can be expanded as

y1=a1rxr=a11x1+a12x2+a13x3

y2=a2rxr=a21x1+a22x2+a23x3

y3=a3rxr=a31x1+a32x2+a33x3

Einstein summation convention

Any expression involving a repeated index (sub- or superscript) shall automatically considered for its sum over the index range 1,2,3...n.

Remark:

  • All indices have the same range, n, unless stated otherwise.
  • No index appear more than twice in any given expression.

Double sum

Summed over i:

Template:Center topaijxiyj=a1jx1yj+a2jx2yj+a3jx3yjTemplate:Center bottom


Summed over j:

Template:Center topaijxiyj=(a11x1y1+a12x1y2+a13x1y3)+(a21x2y1+a22x2y2+a23x2y3)+(a31x3y1+a32x3y2+a33x3y3)Template:Center bottom

Substitution

Let Q=bijyixj and yi=aijxj.


Then Q=bijaijxjxj. However, we have now more than two repeated indices and it is not correct. In order to have a correct substitution,


1st Find duplicated dummy index:  j.


2nd Change dummy index:  yi=airxr.


3rd Substitute:  Q=bijairxrxj.

Kronecker Delta

Template:Center topδij={1i=j0ijTemplate:Center bottom


Template:Center topδii=δ11+δ22+δ33=3Template:Center bottom


Template:Center topδijxixj=1x1x1+0x1x2+0x1x3+0x2x1+1x2x2+0x2x3+0x3x1+0x3x2+1x3x3Template:Center bottom


Template:Center topδijxixj=x1x1+x2x2+x3x3Template:Center bottom

Template:Center topδijxixj=xixi Template:Center bottom

Template:Center topδij replaces one of the dummy index with the other one.Template:Center bottom

Scalar Multiplication

Remember that,

Template:Center topU=Ui eiTemplate:Center bottom


Template:Center topU=|U|cosαi eiTemplate:Center bottom


Template:Center topαU=αUieichange only in amplitude.Template:Center bottom


Addition of Vectors

Template:Center topW=U+V=Uiei+Viei=(Ui+Vi)eiTemplate:Center bottom


Template:Center topWi=Ui+ViTemplate:Center bottom


  • Commutative:

Template:Center topU+V=V+UTemplate:Center bottom


  • Associative:

Template:Center topU+(V+W)=(U+V)+WTemplate:Center bottom

Scalar Product

Template:Center topUV=UiViTemplate:Center bottom


If m and n are unit vectors in the direction of U and V.

Template:Center topm=U|U|=cosαi eiTemplate:Center bottom


Template:Center topn=V|V|=cosβi eiTemplate:Center bottom


Template:Center topmn=cosαicosβi=cosθTemplate:Center bottom


Template:Center topfor θ=90cosθ=0  i.e.  mn=0UV=|U||V|mn=0UVTemplate:Center bottom


Vector Product

File:Scalars Vectors Tensors5.png
Vector product of two vectors

Consider vector products of unit vectors.

Template:Center tope2×e3=e3×e2=e1Template:Center bottom


Template:Center tope3×e1=e1×e3=e2Template:Center bottom


Template:Center tope1×e2=e2×e1=e3Template:Center bottom


Template:Center top(a1e1+a2e2+a3e3)×(b1e1+b2e2+b3e3)=(a2b3a3b2)e1+(a3b1a1b3)e2+(a1b2a2b1)e3Template:Center bottom


The determinant of the symmetric matrix will be the same: Template:Center top|e1e2e3a1a2a3b1b2b3|Template:Center bottom



Permutation symbol

The permutation symbol is defined as: Template:Center topϵijk={0if any of i,j,k are the same1if i,j,k is an even permutation1if i,j,k is an odd permutationTemplate:Center bottom

Hence, the vector product can be written as:

Template:Center topa×b=ijkaibjek=ijkajbkeiTemplate:Center bottom

Triple Scalar Product

Template:Center topa(b×c)=a[ijkbicjek]=akijkbicj=ϵijkakbicj=ijkaibjckTemplate:Center bottom


Tensor Product (Dyadic product)

Product of two tensors

Template:Center topaibj=cij for i,j=1,2,3Template:Center bottom

where Template:Center topcij=[c11c12c13c21c22c23c31c32c33]Template:Center bottom


Similarly,

Template:Center topcijDklm=Eijklm (2+3=5th order tensor)Template:Center bottom


Thus, generally, the product of a tensor having mth order and another tensor having nth order will result in a tensor having an order of m + n.

Therefore,

Template:Center topδijckl=DijklTemplate:Center bottom


The order is increased by two.


Substitution

Template:Center topδijcjk=DijjkTemplate:Center bottom


Template:Center topDijjk=Di11k+Di22k+Di33k=EikTemplate:Center bottom


Consequently if a tensor has a pair of repeated indices, its order reduces by two. And it can be shown that:

Template:Center topEik=cik i.e.Template:Center bottom


Template:Center topδijcjk=cikTemplate:Center bottom


Index j is replaced by index i. This operation is called substitution. Generally:

Template:Center topδijAkjrs=AkirsTemplate:Center bottom


Contraction

Template:Center topδijAijrs=Bijijrs=DrsTemplate:Center bottom

The total order is reduced by four.

Or equally:

Template:Center topδijAijrs=Aiirs=DrsTemplate:Center bottom


The order of Aijrs is reduced by two.

This operation is called contraction with respect to i and j.


Scalar product of tensors

Remember:

Template:Center topab=aibi=δijaibj=δijcij=cii=aibiTemplate:Center bottom


Contraction of a tensor product of two 1st order tensors results in scalar product of two vectors.

Hence, the scalar product of two tensors are:

Template:Center topδipAijBpqr=AijBiqr=ciqr=ApjBpqr=cjqrTemplate:Center bottom


Remember:

Template:Center topA×B=CTemplate:Center bottom


Template:Center topϵijkAjBk=Dijjkk=CiTemplate:Center bottom


Consider,

Template:Center topA(B×C)=Ai(ϵijkBjCk)=Diijjkk (scalar)Template:Center bottom


In conclusion, using component notation and basic rules the form of the tensor can be deduced.

Important Relations

Template:Center topδii=3Template:Center bottom


Template:Center topδijδjk=δikTemplate:Center bottom


Template:Center topδikϵikm=0Template:Center bottom (repeated index)


Template:Center topϵijkϵlmn=|δilδimδinδjlδjmδjnδklδkmδkn|Template:Center bottom


Template:Center topϵijkϵlmn=δilδjmδkn+δimδjnδkl+δinδjlδkm(δilδjnδkm+δimδjlδkn+δinδjmδkl)Template:Center bottom


Contraction with respect to k n.

Template:Center topδknϵijkϵlmn=ϵijkϵlmk=δilδjmδimδjlTemplate:Center bottom


Twice contraction,

Template:Center topδjmδknϵijkϵlmn=δjm[δilδjmδimδjl]Template:Center bottom


Template:Center topδjmδknϵijkϵlmn=[δilδjjδijδjl]=[3δilδil]=2δilTemplate:Center bottom


Template:Center topδjm[δknϵijkϵlmn]=δjm ϵijk ϵlmk = ϵijkϵljk =2δilTemplate:Center bottom


and

Template:Center topϵijkϵijk=6Template:Center bottom

Tripple Vector Product

File:Scalars Vectors Tensors10.png
Triple vector product

Therefore, it can be a linear sum of two vectors.

Template:Center topa×(b×c)=βb+γc=TTemplate:Center bottom


Use component notation.

Template:Center topTi=ijkaj(klmblcm)=βbi+γciTemplate:Center bottom


Template:Center topTi=ijkklmajblcm=βbi+γciTemplate:Center bottom


Template:Center topTi=(δilδjmδimδjl)ajblcmTemplate:Center bottom


Template:Center topTi=δilδjmajblcmδimδjlajblcmTemplate:Center bottom


Template:Center topTi=ambicmajbjci=bi(amcm)ci(ajbj)=(ac)b(ab)cTemplate:Center bottom


Derivative of a Tensor

Let A scalar derivative with respect to xi:

Template:Center topAxi=BiTemplate:Center bottom


Let A a vector:

Template:Center topAixj=BijTemplate:Center bottom


Generally, the mth derivative of a tensor of the nth order is a tensor of the (n+m)th order.

Special case,

Template:Center topxjxi=1 (if i=j)Template:Center bottom


Template:Center topxjxi=0 (if ij)Template:Center bottom


Template:Center topxjxi=δijxixi=δii=3Template:Center bottom


Let A(xj):

Template:Center topA(xj)xi=Axjxjxi=δijAxj=δijCj=CiTemplate:Center bottom


A(n) n(xj):


Template:Center topAxi=Annxjxjxi=δijAnnxjvektorvektorvektorTemplate:Center bottom

Reference

Template:Reflist