Home / Math and Physics Files / Math / Curvilinear Systems / Tensor Doc and Support / Files related to May 2015 update
Lai version of DTijk v4
DOCX · 49.7 KB
Open DOCX file
Phil's derivation notes, dated 4.24.15, rewriting equations from Lai (around page 503) in his own curvilinear-coordinate notation. He expands a rank-2 tensor T in basis vectors, differentiates with the Christoffel-type coefficients G, and relates the rank-3 tensor M to the covariant derivative of T using scale factors h and the R matrices. The result is hk M_ijk = ∂kTij + Tqj Γqki + Tiq Γqkj, matching his equation (8A.26), later written up as Section J.9.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
Lai's Version of (T)'ijk v3 PhL 4.24.15
I have never understood how they do this on page 503, but today I feel confident that I can translate their equations into my notation.
1. Assume this Picture !!!!
The n vectors are in x-space, as usual. Assume the variations of n have this form
dn(x') = [Γ()]kjndx'j k(x')
I can rewrite this as
∂'jn(x') = [Γ()]kjnk(x')
To simplify, lets define
Gkjn ≡ [Γ()]kjn
Then I have
(∂'jn(x')) = Gkjn(x)k(x')
dn = Gkjn(x') dx'j k(x')
di = Gkji dx'j k dei = Γijk dxj ek Lai p 502 dev notation I think
2. Next consider ( Lai notation equations are in blue)
Lai Me
ei i
dxj dx'j
Γijk Gkji
M = Mijk eiejej M = Σijk [M()]ijk ijk [M()]ijk = M'ijk
T = Tij eiej T = Σij [T()]ij ij [T()]ij = T'ij
Let's now install arguments as chosen above
M(x') = Σijk [M()(x')]ijk i(x')j(x')k(x') M = rank 3
M(x') = Σijk M'ijk(x') i(x')j(x')k(x')
T(x') = Σij [T()(x')]ij i(x')j(x') T = rank 2
T(x') = Σij T'ij(x') i(x')j(x')
3. Consider now, for Picture B above,
T(x') = Σij T'ij(x') i(x')j(x')
Then apply ∂'k to get
∂'kT(x') = (∂'kT'ij(x')) i(x')j(x') + T'ij(x')(∂'k i(x'))j(x') + T'ij(x')i(x')(∂kj(x'))
Now use the relation above that
(∂'jn(x')) = Grjn(x')r(x')
so that
(∂'k i(x')) = Grki(x')r(x')
(∂'k j(x')) = Gskj(x')s(x')
and then
∂'kT(x') = (∂'kT'ij(x')) i(x')j(x') + T'ij(x')(∂'k i(x'))j(x') + T'ij(x')i(x')(∂'kj(x'))
= (∂'kT'ij(x')) i(x')j(x') + T'ij(x') Grki(x')r(x')j(x') + T'ij(x')i(x')Gskj(x')s(x')
Now some index shuffling. In the middle term do i↔r and in the last term j↔s:
= (∂'kT'ij(x')) i(x')j(x') + T'rj(x') Gikr(x')i(x')j(x') + T'is(x')i(x')Gjks(x')j(x')
= [ ∂'kT'ij(x') + T'rj(x') Gikr(x') + T'is(x')Gjks(x')] i(x')j(x')
and this verifies the result quoted above.
4. If we assume Picture B as above then
(∂kT)rs = ∂kTrs = Mrsk
and this is our first inkling of a connection between M and T. The index ordering above I take from this quote
(T)ijk ≡ ∂kTij (J.1.3)
with the idea that somehow M = T in the abstract.
What does this look like in x'-space? That is a sticky wicket in this notation because ∂kTrs is not a tensor! I will assume however that Mrsk is a true tensor and what I really mean is this
(T)rsk =Trs;k = Mrsk
and in this notation we have in x'-space
(T)'rsk =T'rs;k = M'rsk
5. Question 1: What tensor transformation rules apply to the above line going from x to x'-space? First, we know that
M'rsk(x') =Rrb Rsc Rka Mbca(x) T'rs(x') = RrdRseTde(x) ∂'k = Rkd∂d
6. Going down a blind path, I could then evaluate, as in earlier appendices
∂'kT'rs(x') = ( Rkd∂d)(RrdRseTde)
= Rkd [RrdRse(∂dTde) + (∂dRrd)RseTde + Rrd(∂dRse)Tde ]
and this at least makes a connection between [ ∂'kT'rs(x')] and [∂dTde(x)] .
7. I don't really see any connection between this last item and the expansions in 2. Well there is this;
Define
Q'ijk(x') ≡ [∂'kT'ij(x') + T'aj(x') Gika(x') + T'ia(x')Gjka(x') ]
I arbitrarily put Q in script with a prime but it means nothing right now, just a definition. Then from 3 above
(∂'kT(x')) = Σij Q'ijk(x') i(x')j(x')
dT(x') = (∂'kT(x'))dx'k = Σij Q'ijk(x')dx'k i(x')j(x')
So OK, the first line above does make some kind of connection between an expansion and ∂'kT. However, I don't really know what (∂'kT(x)) means having no indices. It is a vector in direct product space. Probably I could dot both sides with ua(x')ub(x') and use ua(x') i(x') = [i(x')]a = [i(x')]a . Then maybe this is true
∂'kTab(x') = Σij Q'ijk(x') [i(x')]a [j(x')]b
But the left side here is really unclear. Should there be a prime on the T? There is no prime when I just write the expansion T = expansion. It is really (∂'kT(x'))ab . This object is fuzzy and I cannot really continue to use it to mean anything. But OK, let's assume it has no prime somehow because we dotted with the u vectors, and let's see what happens. Then I have connected Q'ijk with something which at least has a meaning. Suppose I then multiply both sides of this by [r(x')]a [s(x')]b and sum on a and b. I end up then with things like
Σa [i(x')]a [r(x')]a = i(x') r(x')
These vectors are in x-space of Picture B. Therefore i(x') h'i = ei(x').
Σa [i(x')]a [r(x')]a = i(x') r(x') = (hi'hr;)-1 ei(x') er(x') = (h'ih'r)-1 g'ir(x)
But we know that g'ir(x) = (h'i)2 δi,r so then
(h'ih'r)-1 g'ir(x) = (h'ih'r)-1 (h'i)2 δi,r = δi,r
So carrying out this process we get
[r(x')]a [s(x')]b ∂'kTab(x') = Σij Q'ijk(x') δi,r δj,s = Q'krs(x')
so we have now an expression for Q'krs, for what it is worth. At least I now know this:
∂'kTab(x') = M'abk = Σij Q'ijk [i(x)]a [j(x)]b
This certainly looks fishy since we have some kind of mixed world. But carry on anyway.
8. What I am really interested in is M'abk = h'ah'bh'k M'abk. So I can write
M'abk = h'ah'bh'k M'abk // Picture A,B notation
M'abk = h'ah'bh'k RaARbB RkK MABK = h'ah'bh'k RaARbB RkK ∂KTAB
But x-space is Cartesian, so MABK = ∂KTAB = ∂KTAB. From above I know that
∂'kTab(x') = Σij Q'ijk(x') [i(x')]a [j(x')]b .
Now we know that ∂'k = RkK∂K and we can think of Tab(x') = Tab(x'(x) ) as a function of x. Then
RkK(∂KTab) = Σij Q'ijk [i]a [j]b x' = x'(x) everywhere
Now use the inversion rule to write this as
(∂KTab) = RkK Q'ijk [i]a [j]b
or
(∂KTAB) = RsK Q'ijs [i]A [j]B
so I then have
M'abk = h'ah'bh'k RaARbB RkK MABK = h'ah'bh'k RaARbB RkK ∂KTAB
= h'ah'bh'k RaARbB RkK {RsK Q'ijs [i]A [j]B }
= h'ah'bh'k RaARbB RkKRsK [i]A [j]B Q'ijs
Do I know anything about RkKRsK ?? Not an orthogonality rule. Could write as
RkKRsK = g'ka gKb RabRsK = g'ka RaKRsK = g'ka δas = g'ks (as I suspected, should be doc'd)
And for orthogs we get
g'ks = δk,s (h's)-2
Then we have showed that
RkKRsK = (h'k)-2δk,s
I might as well install this right now to get
M'abk = h'ah'bh'k RaARbB (h'k)-2 [i]A [j]B Q'ijk
or
M'abk = h'k-1 h'ah'bRaARbB [i]A [j]B Q'ijk
Now this last statement is valuable. It relates [M()]abk = M'abk to QijK . So now the question is this: Can this last line be simplified? I first quote
Mab ≡ h'a Rab (E.8.7)
and then I can write
M'abk = h'k-1 (h'aRaA)(h'bRbB) [i]A [j]B Q'ijk
or
M'abk = h'k-1 MaAMbB [i]A [j]B Q'ijk
I quote again.
Nin = h'n-1Rni = h'n-1(en)i = (n)i , // (7.13.1) (E.8.12)
Then I can further simplify
M'abk = h'k-1 MaAMbB [i]A [j]B Q'ijk
NAi NBj
M'abk = h'k-1 MaAMbBNAiNBj Q'ijk
= h'k-1 (MaANAi)(MbBNBj) Q'ijk
= h'k-1δaiδbj Q'ijk
= h'k-1Q'abk
which then says
h'kM'abk = Q'abk
or
h'kM'abk = [∂'kT'ab(x') + T'qb(x') Gakq(x') + T'aq(x')Gbkq(x') ]
9. Now how would I convert this to Lai notation? From above I already have
Mijk = M'ijk Tij = T'ij Γijk = Gkji hk = h'k
Lai's equations are for Picture C1 so I replace x' → x everywhere. So start with
h'kM'abk = [∂'kT'ab(x') + T'qb(x') Gakq(x') + T'aq(x')Gbkq(x') ]
and make the translations to get
hkMabk = [∂kTab + TqbGakq(x') + Taq Gbkq(x') ]
Γqka Γqkb
or
Mabkhk = [∂kTab + Tqb Γqka + Taq Γqkb ]
or
Mijmhm = [∂mTij + Tqj Γqmi + Tiq Γqmj ]
Finally, this agrees with (8A.26), this only took me two days!!
10. Now that I have finally gotten the right answer, I have to do it all again in some efficient manner. The secret was to get the Pictures right!!!
// This is now written up as Section J.9 !!