1 CL305 Vectors Tensors
1 CL305 Vectors Tensors
𝜹𝒊 × 𝜹𝒋 = 𝜺𝒊𝒋𝒌 𝜹𝒌
𝒌=𝟏
• These two relations enable us to develop analytic expressions for all the common
dot and cross operations
• For developing expressions for vector and tensor operations, we have to break all
vectors up into components and then apply above two relations
• Expansion of a vector in terms of its components:
• Any vector v can be specified by giving values of its projections v1, v2, v3 on
coordinate axes 1,2,3 as shown below
• Vector can be constructed by adding vectorially the components multiplied by their
corresponding unit vectors
• 𝐯 = 𝜹𝟏 𝒗𝟏 + 𝜹𝟐 𝒗𝟐 + 𝜹𝟑 𝒗𝟑 = σ𝟑𝒊=𝟏 𝜹𝒊 𝒗𝒊
• 𝒗𝒊 are called as “components of vector v” and are scalars
• 𝜹𝒊 𝒗𝒊 are vectors, which when added together vectorially give v
= 𝛅𝑖 𝛅𝑙 𝜎𝑖𝑗 𝜏𝑗𝑙
𝑖 𝑙 𝑗
• Not commutative: 𝛻𝑠 ≠ 𝑠𝛻
• Not associative: 𝛻𝑟 𝑠 ≠ 𝛻 𝑟𝑠
• Distributive: 𝛻 𝑟 + 𝑠 = 𝛻𝑟 + 𝛻𝑠
• Divergence of a vector field:
• If vector v is a function of space variables x1, x2, x3, then scalar
product can be formed with operator 𝛻 as
𝜕 𝜕
• 𝛻·𝐯 = σ𝑖 𝛅𝑖
𝜕𝑥𝑖
· σ𝑗 𝛅𝑗 𝑣𝑗 = σ𝑖 σ𝑗 𝛅𝑖 · 𝛅𝑗 𝑣
𝜕𝑥𝑖 𝑗
=
𝜕 𝜕𝑣𝑖
σ𝑖 σ𝑗 δ𝑖𝑗 𝑣 = σ𝑖
𝜕𝑥𝑖 𝑗 𝜕𝑥𝑖
• It is also abbreviated as div v
• Not commutative: 𝛻·𝐯 ≠ 𝐯·𝛻
• Not associative: 𝛻 · 𝑠𝐯 ≠ 𝛻𝑠 · 𝐯
• Distributive: 𝛻· 𝐯+𝐰 = 𝛻·𝐯 + 𝛻·𝐰
• Curl of a vector field:
• Cross product may also be formed between operator 𝛻 and vector v, which is function of
space variables
𝜕
• 𝛻×𝐯 = σ𝑗 𝛅𝑗
𝜕𝑥𝑗
× σ𝑘 𝛅𝑘 𝑣𝑘
𝜕 𝜕
• = σ𝑗 σ𝑘 𝛅𝑗 × 𝛅𝑘 𝑣
𝜕𝑥𝑗 𝑘
= σ𝑖 σ𝑗 σ𝑘 𝜀𝑖𝑗𝑘 𝛅𝑖 𝜕𝑥 𝑣𝑘
𝑗
𝛅1 𝛅2 𝛅3
𝜕 𝜕 𝜕 𝜕𝑣3 𝜕𝑣 𝜕𝑣1 𝜕𝑣 𝜕𝑣2 𝜕𝑣
• = 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
= 𝛅1 𝜕𝑥2
− 𝜕𝑥2 + 𝛅2 𝜕𝑥3
− 𝜕𝑥3 + 𝛅3 𝜕𝑥1
− 𝜕𝑥1
3 1 2
𝑣1 𝑣2 𝑣3
• Thus vector constructed is called curl of v
• Other notations for 𝛻 × 𝐯 are curl v and rot v
• Like divergence, the curl operation is also distributive but not commutative or associative
• Gradient of a vector field:
• In addition to scalar product 𝛻 · 𝐯 and vector product 𝛻 × 𝐯 , one can also form
dyadic product 𝛻𝐯 as:
𝜕 𝜕
𝛻v = 𝛅𝑖 𝛅𝑗 𝑣𝑗 = 𝛅𝑖 𝛅𝑗 𝑣
𝜕𝑥𝑖 𝜕𝑥𝑖 𝑗
𝑖 𝑗 𝑖 𝑗
• This is gradient of vector v and is sometimes written grad v
𝜕
• It is a second-order tensor whose ij-component is 𝜕𝑥 𝑣𝑗
𝑖
𝜕
• Its transpose is: 𝛻v † = σ𝑖 σ𝑗 𝛅𝑖 𝛅𝑗 𝑣
𝜕𝑥𝑗 𝑖
𝜕
• Whose ij-component is 𝜕𝑥 𝑣𝑖
𝑗
• In addition: 𝛻v ≠ 𝐯𝛁 and 𝛻v †
≠ v𝛻
• Divergence of a tensor field:
• If the tensor τ is a function of space variables x1, x2, x3, then a vector product may be
formed with operator 𝛻 as:
𝜕 𝜕
𝛻·𝛕 = 𝛅𝑖 · 𝛅𝑗 𝛅𝑘 τ𝑗𝑘 = 𝛅𝑖 · 𝛅𝑗 𝛅𝑘 τ𝑗𝑘
𝜕𝑥𝑖 𝜕𝑥𝑖
𝑖 𝑗 𝑘 𝑖 𝑗 𝑘
𝜕 𝜕
= δ𝑖𝑗 𝛅𝑘 τ = 𝛅𝑘 τ
𝜕𝑥𝑖 𝑗𝑘 𝜕𝑥𝑖 𝑖𝑘
𝑖 𝑗 𝑘 𝑘 𝑖
• This is called divergence of tensor τ and is written as div τ
𝜕
• kth component of 𝛻 · 𝛕 is σ𝑖 𝜕𝑥 τ𝑖𝑘
𝑖
𝜕
• If τ is the product svw, then 𝛻 · 𝑠𝐯𝐰 = σ𝑘 𝛅𝑘 σ𝑖 𝜕𝑥 𝑠𝑣𝑖 𝑤𝑘
𝑖
• Laplacian of a scalar field:
• It is divergence of a gradient of scalar function s:
𝜕 𝜕𝑠
• 𝛻 · 𝛻𝑠 = σ𝑖 𝛅𝑖
𝜕𝑥𝑖
∙ σ𝑗 𝛅𝑗
𝜕𝑥𝑗
𝜕 𝜕𝑠 𝜕2
• = σ𝑖 σ𝑗 δ𝑖𝑗
𝜕𝑥𝑖 𝜕𝑥𝑗
= σ𝑖 2
𝜕𝑥𝑖
𝑠
• 𝛻 · 𝐯𝐰 = 𝐯 · 𝛻𝐰 + 𝐰 𝛻 · 𝐯
• 𝑠𝛅: 𝛻𝐯 = 𝑠(𝛻 · 𝐯)
• 𝛻 · 𝑠𝛅 = 𝛁𝑠
• 𝛻 · 𝑠𝛕 = 𝛻𝑠 · 𝛕 + 𝑠 𝛻 · 𝛕
• 𝛻 𝐯 · 𝐰 = (𝛻𝐯) · 𝐰 + (𝛻𝐰) · 𝐯
• Example: prove that for symmetric τ 𝛕: 𝛻𝐯 = 𝛻 · 𝛕 · 𝐯 − 𝐯 · 𝛻 · 𝛕
• Solution:
• Write RHS in terms of components:
𝜕 𝜕
• 𝛻 · 𝝉 · 𝒗 = σ𝑖 𝜕𝑥 𝛕 · 𝑣 𝑖 = σ𝑖 σ𝑗 𝜕𝑥 𝜏𝑖𝑗 𝑣𝑗 (1)
𝑖 𝑖
𝜕
• 𝐯 · 𝛻 · 𝝉 = σ𝑗 𝑣𝑗 𝛻 · 𝝉 𝑗 = σ𝑗 σ𝑖 𝑣𝑗 𝜕𝑥 𝜏𝑖𝑗 (2)
𝑖
• Now LHS can be written in terms of components as:
𝜕 𝜕
• 𝛕: 𝛻𝐯 = σ𝑖 σ𝑗 𝜏𝑗𝑖 𝜕𝑥 𝑣𝑗 = σ𝑖 σ𝑗 𝜏𝑖𝑗 𝜕𝑥 𝑣𝑗 (3)
𝑖 𝑖
• 2nd form resulting from symmetry of τ
• Subtraction of Eq. (1) from Eq. (2) will give Eq. (3)
Vector and tensor algebra in curvilinear coordinates
𝜕 𝝏 𝝏 𝝏
• 𝛁𝐯 = σ𝑖 σ𝑗 𝛅𝑖 𝛅𝑗 𝜕𝑥 𝑣𝑗 = 𝜹𝒙 𝜹𝒙 𝝏𝒙 𝒗𝒙 + 𝜹𝒙 𝜹𝒚 𝝏𝒙 𝒗𝒚 + 𝜹𝒙 𝜹𝒛 𝝏𝒙 𝒗𝒛
𝑖
𝝏 𝝏 𝝏
+𝜹𝒚 𝜹𝒙 𝒗 + 𝜹𝒚 𝜹𝒚 𝒗 + 𝜹𝒚 𝜹𝒛 𝒗
𝝏𝒚 𝒙 𝝏𝒚 𝒚 𝝏𝒚 𝒛
𝝏 𝝏 𝝏
+𝜹𝒛 𝜹𝒙 𝒗 + 𝜹𝒛 𝜹𝒚 𝒗 + 𝜹𝒛 𝜹𝒛 𝒗
𝝏𝒛 𝒙 𝝏𝒛 𝒚 𝝏𝒛 𝒛
= 𝜹𝒙 𝜹𝒙 𝟎 + 𝜹𝒙 𝜹𝒚 𝒃 + 𝜹𝒙 𝜹𝒛 𝟎 + 𝜹𝒚 𝜹𝒙 𝒃 + 𝜹𝒚 𝜹𝒚 𝟎 + 𝜹𝒚 𝜹𝒛 𝟎 + 𝜹𝒛 𝜹𝒙 𝟎 + 𝜹𝒛 𝜹𝒚 𝟎 + 𝜹𝒛 𝜹𝒛 𝟎
= 𝜹𝒙 𝜹𝒚 𝒃 + 𝜹𝒚 𝜹𝒙 𝒃
• 𝛁 ∙ 𝐯𝐯 is of the order of Σ−2 = 1+1+1−2 = 1 so vector. Here 𝑣𝑥 = 𝑏𝑦 ; 𝑣𝑦 = 𝑏𝑥 ; 𝑣𝑧 = 0
𝝏 𝝏 𝝏
• 𝛁 ∙ 𝐯𝐯 = 𝜹𝒙 𝛁 ∙ 𝐯𝐯 𝒙 + 𝜹𝒚 𝛁 ∙ 𝐯𝐯 𝒚 + 𝜹𝒛 𝛁 ∙ 𝐯𝐯 𝒛 = 𝜹𝒙 𝝏𝒙 𝒗𝒙 𝒗𝒙 + 𝜹𝒙 𝝏𝒚 𝒗𝒚 𝒗𝒙 + 𝜹𝒙 𝝏𝒛 𝒗𝒛 𝒗𝒙
𝝏 𝝏 𝝏
+𝜹𝒚 𝝏𝒙 𝒗𝒙 𝒗𝒚 + 𝜹𝒚 𝝏𝒚 𝒗𝒚 𝒗𝒚 + 𝜹𝒚 𝝏𝒛 𝒗𝒛 𝒗𝒚
𝝏 𝝏 𝝏
+𝜹𝒛 𝝏𝒙 𝒗𝒙 𝒗𝒛 + 𝜹𝒛 𝝏𝒚 𝒗𝒚 𝒗𝒛 + 𝜹𝒛 𝝏𝒛 𝒗𝒛 𝒗𝒛
𝝏 𝟐 𝟐 𝝏 𝟐 𝝏
= 𝜹𝒙 𝒃 𝒚 + 𝜹𝒙 𝒃 𝒙𝒚 + 𝜹𝒙 𝟎
𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏 𝟐 𝝏 𝟐 𝟐 𝝏
+𝜹𝒚 𝒃 𝒙𝒚 + 𝜹𝒚 𝒃 𝒙 + 𝜹𝒚 𝟎
𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏 𝝏 𝝏
+𝜹𝒛 𝝏𝒙 𝟎 + 𝜹𝒛 𝝏𝒚 𝟎 + 𝜹𝒛 𝝏𝒛 𝟎 = 𝒃𝟐 𝒙𝜹𝒙 + 𝒃𝟐 𝒚𝜹𝒚
• Similarly do for second velocity field: (b). 𝑣𝑥 = −𝑏𝑦 ; 𝑣𝑦 = 𝑏𝑥 ; 𝑣𝑧 = 0.
𝜏𝑥𝑥 = 3 𝜏𝑥𝑦 = 2 𝜏𝑥𝑧 = −1
• (4). Components of a tensor 𝝉ധ are 𝜏𝑦𝑥 = 2 𝜏𝑦𝑦 = 2 𝜏𝑦𝑧 = 1 and the components of a
𝜏𝑧𝑥 = −1 𝜏𝑧𝑦 = 1 𝜏𝑧𝑧 = 4
vector 𝑣Ԧ are vx = 5, vy = 3 and vz = −2. Calculate the following: (a) 𝝉ധ ∙ 𝑣 , (b) 𝝉ധ : 𝝉ധ and (c)
𝑣 ∙ 𝝉ധ ∙ 𝑣
• Solution:
• (a) Order of 𝝉 ∙ 𝐯 is ∑−2 = 2+1−2 = 1, so vector
• 𝝉 ∙ 𝐯 = 𝜹𝒙 𝝉 ∙ 𝐯 𝒙 + 𝜹𝒚 𝝉 ∙ 𝐯 𝒚 + 𝜹𝒛 𝝉 ∙ 𝐯 𝒛 = 𝜹𝒙 𝝉𝒙𝒙 𝒗𝒙 + 𝜹𝒙 𝝉𝒙𝒚 𝒗𝒚 + 𝜹𝒙 𝝉𝒙𝒛 𝒗𝒛
+𝜹𝒚 𝝉𝒚𝒙 𝒗𝒙 + 𝜹𝒚 𝝉𝒚𝒚 𝒗𝒚 + 𝜹𝒚 𝝉𝒚𝒛 𝒗𝒛
+𝜹𝒛 𝝉𝒛𝒙 𝒗𝒙 + 𝜹𝒛 𝝉𝒛𝒚 𝒗𝒚 + 𝜹𝒛 𝝉𝒛𝒛 𝒗𝒛
= 𝜹𝒙 (𝟑)(𝟓) + 𝜹𝒙 (𝟐)(𝟑) + 𝜹𝒙 (−𝟏)(−𝟐)
+𝜹𝒚 (𝟐)(𝟓) + 𝜹𝒚 (𝟐)(𝟑) + 𝜹𝒚 (𝟏)(−𝟐)
+𝜹𝒛 −𝟏 𝟓 + 𝜹𝒛 𝟏 𝟑 + 𝜹𝒛 𝟒 −𝟐
= 𝟐𝟑𝜹𝒙 + 𝟏𝟒𝜹𝒚 − 𝟏𝟎𝜹𝒛
• (b). 𝝉ധ : 𝝉ധ is of order of ∑−4 = 2+2−4 = 0, so scalar
• 𝝉: 𝝉 = σ𝒊 σ𝒋 𝝉𝟐𝒊𝒋 = 𝟑𝟐 + 𝟐𝟐 + (−𝟏)𝟐 +𝟐𝟐 + 𝟐𝟐 + 𝟏𝟐 + (−𝟏)𝟐 +𝟏𝟐 + 𝟒𝟐 = 𝟒𝟏
• (c). 𝑣 ∙ 𝝉ധ ∙ 𝑣
• 𝝉ധ ∙ 𝑣 is of the order of ∑−2 = 2+1−2 = 1, so vector
• Then 𝑣 ∙ 𝝉ധ ∙ 𝑣 is of the order of ∑−2 = 1+1−2 = 0, so overall a scalar will be
resulted.
• 𝑣 ∙ 𝝉ധ ∙ 𝑣 = σ𝒊 𝒗𝒊 𝝉ധ ∙ 𝑣 𝒊 = 𝒗𝒙 𝝉ധ ∙ 𝑣 𝒙 + 𝒗𝒚 𝝉ധ ∙ 𝑣 𝒚 + 𝒗𝒛 𝝉ധ ∙ 𝑣 𝒛
= 𝒗𝒙 𝝉𝒙𝒙 𝒗𝒙 + 𝒗𝒙 𝝉𝒙𝒚 𝒗𝒚 + 𝒗𝒙 𝝉𝒙𝒛 𝒗𝒛
+𝒗𝒚 𝝉𝒚𝒙 𝒗𝒙 + 𝒗𝒚 𝝉𝒚𝒚 𝒗𝒚 + 𝒗𝒚 𝝉𝒚𝒛 𝒗𝒛
+𝒗𝒛 𝝉𝒛𝒙 𝒗𝒙 + 𝒗𝒛 𝝉𝒛𝒚 𝒗𝒚 + 𝒗𝒛 𝝉𝒛𝒛 𝒗𝒛
= 𝟓 𝟑 𝟓 + 𝟓 𝟐 𝟑 + 𝟓 −𝟏 −𝟐 + 𝟑 𝟐 𝟓 + 𝟑 𝟐 𝟑 + 𝟑 𝟏 −𝟐
+ −𝟐 −𝟏 𝟓 + −𝟐 𝟏 𝟑 + −𝟐 (𝟒) −𝟐 = 𝟏𝟕𝟕
References
• B. R. Bird, E. W. Stewart and N. E. Lightfoot, Transport
Phenomena, John Wiley & Sons, 2003 (Second Edition).
• W. M.Deen, Analysis of Transport Phenomena, Oxford University
Press, New York, 1998 (First Indian Edition).
• R.L. Panton, Incompressible Flow, Wiley India, 2011 (Third
Edition).
• R.P. Chhabra and J.F. Richardson, Non-Newtonian Flow and
Applied Rheology, 2nd Edition, Butterworth-Heinemann, Oxford,
UK, 2008.
• C.W. Macosko, Rheology: Principles, Measurements, and
Applications, Wiley-VCH, New York, 1994.
Thank you