Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Math / Curvilinear Systems / Tensor Doc and Support / tensor doc support 3_9_12 thru 3_28_12

vec lap the other way A

DOCX · 23.9 KB
Open DOCX file

Phil's note dated 3.24.12 tries to derive the vector Laplacian as u_{i;j;j} and show it equals grad(div u) minus curl curl u. He checks the identity in Cartesian coordinates using epsilon-epsilon contractions, then expands covariant derivatives with Christoffel symbols and compares with derivatives of sqrt(g). He ends by noting a resemblance to the Riemann curvature tensor. He says at the start that the approach did not work out.

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
Another approach to vector Laplacian PhL 3.24.12 I was unable to make this fly, but will keep it around for some future day. Well, consider Qi ≡ 2(ui) = div grad ui. If ui is a vector, grad ui we know is not a rank-2 tensor, and then div [ grad ui] is also not a tensor. Think about the covariant approach. Start with Qi = 2(ui) = ∂j∂j ui = ∂j ui,j = ui,j,j This thing is not a tensor until we put in the semicolons Qi = ui;j;j Only now is Qi a true vector. Maybe this is a simpler approach compared to what I did in Section 15. How does this involve that curl I wonder? Let's use Appendix XXX to write ui;j;k = gkα ui;j;α = gkα [ ∂α ui;j + Γiαn un;j – Γnjα ui;n ] Then you go install this: ui;j = ∂j ui + Γimjum to get ui;j;k = gkα [ ∂α {∂j ui + Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum } ] = [ ∂k {∂j ui + Γinjun} + gkα Γiαn {∂j un + Γnmjum } – gkα Γnjα {∂n ui + Γimnum } ] Now contract k with j to get ui;j;j = [ ∂j {∂j ui + Γinjun} + gjα Γiαn {∂j un + Γnmjum } – gjα Γnjα {∂n ui + Γimnum } ] And now we wonder if this equals (u)i ??? Maybe consider at an earlier point: Qi = ui;j;j (u)i = (uj;j);i – g-1/2εiab(g-1/2εbdeue;d);a If these really were equal, that would imply that ( I do the two semicolon replacements on the right) uj;j;i – ui;j;j = g-1/2εiab(g-1/2εbdeue,d),a That is pretty non-obvious to me! Is it even true in Cartesians? Then it would say uj,j,i – ui,j,j = εiabεbde(ue,d,a) But in this case we know that εiabεbde = εbiaεbde = δidδae – δieδad so RHS = εiabεbde(ue,d,a) = [δidδae – δieδad](ue,d,a) = (ua,i,a) – (ui,a,a) = (uj,i,j) – (ui,j,j) = uj,j,i – ui,j,j and yes, it works.! So armed with this info, let's try to do it more generally. Start with ui;j;k = gkα [ ∂α {∂j ui + Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum }] First, replace i→I uI;j;k = gkα [ ∂α {∂j uI + ΓInjun} + ΓIαn {∂j un + Γnmjum } – Γnjα {∂n uI + ΓImnum }] Now set I = j and k = i to get uj;j;i = giα [ ∂α {∂j uj + Γjnjun} + Γjαn {∂j un + Γnmjum } – Γnjα {∂n uj + Γjmnum }] Now start again with uI;j;k = gkα [ ∂α {∂j uI + ΓInjun} + ΓIαn {∂j un + Γnmjum } – Γnjα {∂n uI + ΓImnum }] This time, replace I = i and k = j ui;j;j = gjα [ ∂α {∂j ui + Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum }] We can then write uj;j;i – ui;j;j = giα [ ∂α {∂j uj + Γjnjun} + Γjαn {∂j un + Γnmjum } – Γnjα {∂n uj + Γjmnum }] – gjα [ ∂α {∂j ui + Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum }] Somehow I expect the entire RHS to reduce to g-1/2εiab(g-1/2εbdeue,d),a There are 12 terms on the RHS now. Here are the 6 double Γ terms (only i is a free index) giα Γjαn Γnmjum – giα Γnjα Γjmnum – gjα Γiαn Γnmjum + gjα Γnjα Γimnum In the first term do j↔n, = giα Γnαj Γjmnum – giα Γnjα Γjmnum – gjα Γiαn Γnmjum + gjα Γnjα Γimnum The first and fourth term can then be combined like this: giα Γnαj Γjmnum + STOP. Back up to an earlier point, uI;j;k = gkα [ ∂α uI;j + ΓIαn un;j – Γnjα uI;n ] First, replace i→j then k→i : uj;j;i = giα [ ∂α uj;j + Γjαn un;j – Γnjα uj;n ] Now start again uI;j;k = gkα [ ∂α uI;j + ΓIαn un;j – Γnjα uI;n ] and this time take I → i and k→ j ui;j;j = gjα [ ∂α ui;j + Γiαn un;j – Γnjα ui;n ] Now we can say uj;j;i – ui;j;j = giα [ ∂α uj;j + Γjαn un;j – Γnjα uj;n ] – gjα [ ∂α ui;j + Γiαn un;j – Γnjα ui;n ] I would like now to show that RHS = g-1/2εiab(g-1/2εbdeue,d),a . Maybe start with this g-1/2εiab(g-1/2εbdeue,d),a = g-1/2εiab(g1/2εbdeue,d),a = g-1/2εiabεbde (g1/2ue,d),a = g-1/2εbiaεbde (g1/2ue,d),a = g-1/2[ δid δae - δie δad ] (g1/2ue,d),a = g-1/2 { (g1/2ua,i),a – (g1/2ui,a),a } = g-1/2 { ∂a(g1/2ua,i) – ∂a (g1/2ui,a) } = g-1/2 { g1/2∂a(ua,i) +∂a(g1/2) ua,i – g1/2∂a (ui,a) – ∂a (g1/2) ui,a)} = [∂a(ua,i) – ∂a (ui,a) ] + ∂a(g1/2)[ ua,i – ui,a] = [∂j(uj,i) – ∂j (ui,j) ] + ∂j(g1/2)[ uj,i – ui,j] = [ uj,i,j –ui,j,j ] + ∂j(g1/2)[ uj,i – ui,j] = [ uj,j,i –ui,j,j ] + ∂j(g1/2)[ giαuj,α – gjαui,α] So I am done if I can show this: uj;j;i – ui;j;j = [ uj,j,i –ui,j,j ] + ∂j(g1/2)[ giαuj,α – gjαui,α] (*) Earlier I showed that uj;j;i – ui;j;j = giα [ ∂α {∂j uj + Γjnjun} + Γjαn {∂j un + Γnmjum } – Γnjα {∂n uj + Γjmnum }] – gjα [ ∂α {∂j ui + Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum }] The first terms on each line may be combined to get giα ∂α {∂j uj) – gjα ∂α (∂j ui) = ∂i {∂j uj) –∂j (∂j ui) = uj,j,i – ui,j,j This matches the first term on the RHS of (*), so we would then have to show that = giα [ ∂α {Γjnjun} + Γjαn {∂j un + Γnmjum } – Γnjα {∂n uj + Γjmnum }] – gjα [ ∂α { Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum }] = ∂j(g1/2)[ giαuj,α – gjαui,α] = g1/2 Γaaj[ giαuj,α – gjαui,α] I do know from Appendix E (d) that ∂j(g1/2) = g1/2 Γaaj to give the second RHS shown above. How on earth could this complicated thing ever be true? What action can generate g1/2 ? The square root of a determinant of gij. The only helpful theorems I know of are these: (∂cgab) = – [gan Γ bcn + gbn Γacn] 1 (∂cgab) = + [gan Γncb + gbn Γnca] 2 Γdab = (1/2) gdc [ ∂agbc + ∂bgca – ∂cgab] Γaan = (1/2) gad ∂ngad = (1/2)(1/g)∂ng = (1/) ∂n() Go back now to what we want to show is true: = giα [ ∂α {Γjnjun} + Γjαn {∂j un + Γnmjum } – Γnjα {∂n uj + Γjmnum }] – gjα [ ∂α { Γinjun} + Γiαn {∂j un + Γnmjum } – Γnjα {∂n ui + Γimnum }] = ∂j(g1/2)[ giαuj,α – gjαui,α] This would have to be true for any vector ui. It seems then that all the non-derivative terms would then have to add up to 0: giα [ ∂α {Γjnj}un + Γjαn {Γnmjum } – Γnjα {Γjmnum }] – gjα [ ∂α { Γinj}un + Γiαn { Γnmjum } – Γnjα {Γimnum }] = 0 ? giα [ ∂α {Γjmj}um + Γjαn {Γnmjum } – Γnjα {Γjmnum }] – gjα [ ∂α { Γimj}um + Γiαn { Γnmjum } – Γnjα {Γimnum }] = 0 ? giα [ ∂α {Γjmj} + ΓjαnΓnmj – ΓnjαΓjmn ] – gjα [ ∂α { Γimj} + Γiαn Γnmj – ΓnjαΓimn] = 0 ? Looking at Weinberg p 133, the above structure looks a lot like the curvature tensor. He says Rλμνκ = ∂κΓλμν – ∂νΓλμκ + ΓnμνΓλκn – ΓnμκΓλνn Rjmνα = ∂αΓjmν – ∂νΓjmα + ΓnmνΓjαn – ΓnmαΓjνn Rjmjα = ∂αΓjmj – ∂jΓjmα + ΓnmjΓjαn – ΓnmαΓjjn