Fluid Mechanics for MAP/Scalar, Vectors and Tensors

From testwiki
Revision as of 21:25, 6 March 2022 by imported>Dave Braunschweig (Remove math tags from headings)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

>back to Chapters

Motivation

Each scientific field needs an optimised language. Without it, this field can not develope. For example, calculus could not develop with Roman numbers.


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


Without arabic numbers and the decimal system, even the 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 tackle with 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 topdiva=a or grad a= a <mrow data-mjx-texclass="ORD">and<mtext> </mtext>for<mtext> </mtext>a<mtext> </mtext>vector</mrow> a¯, a¯ or a~Template:Center bottom

Symbols are not unique and, somewhat, confusing, i.e. unified language is beneficial to speed up the development in this field. In this chapter, we will be particularly introduced to such unified mathematical concept, Tensor which is very important for computational fluid dynamics and also many other field of science. This chapter will give a gentle introduction about fascinating Tensor but readers are encouraged to read more from the references.Due to this short introduction,Tensor will not be used any other 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 A operates on vector u to give new vector v.

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


Those are 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 bottom


 Template:Center topa[Pa]+b[K]= ?Template:Center bottom
 
 

Multiplication is valid for variables having different units.


 Template:Center topa[Pa]*b[K]=a*b[PaK]Template:Center bottom
   
 

Products of variables having different units, have a new physical meaning.


 Template:Center topm=ρ*VTemplate:Center bottom


 Template: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 bottom


 Template: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 fro series:

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

In short, we will write Template:Center topaixiTemplate:Center bottom 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 idex appear more than twice in any given expresion.

Double sum

Sumed over i:

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


Sumed 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. Here we find more than two repeated indices.


1st Find duplicated dumy index:  j.


2nd Change dummy index:  yi=airxr.


3rd Substitude:  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=xixiTemplate:Center bottom


replaces one of the dummy index with the other one.

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 symmetric determinant 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

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 opperation 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