Krumlinjete koordinater

Fra testwiki
Hopp til navigering Hopp til søk
Krumlinjete koordinater (q1,q2,q3)  og tilsvarende basisvektorer (e1, e2, e3) lar en vilkårlig vektor v uttrykkes til venstre ved sine kontravariante komponenter. I den duale basisen (e1, e2, e3) vises den samme vektoren med sine kovariante komponenter.

Krumlinjete koordinater er koordinatsystem hvor en eller flere av koordinatlinjene er krumme. På samme måte vil også en eller flere av koordinatflatene ha krumning. For å angi punkter på Jordens overflate benyttes for eksempel lengdegrad og breddegrad, som er slike koordinater. I astronomien bruker man på tilsvarende vis forskjellige himmelkoordinater.

Krumlinjete koordinater kan benyttes i euklidske rom hvor bestemte symmetriforhold gjør det naturlig og dermed tillater en enklere, matematisk beskrivelse. Et eksempel er utledningen av Keplers lover for planetbevegelse under påvirkning av tyngdekraften fra Solen. Dette gjøres enklest i kulekoordinater, da tyngdekraften bare avhenger av avstanden og derfor er symmetrisk om alle sentrale akser. På samme måte benyttes disse koordinatene i atomfysikken, hvor elektronene beveger seg i Coulomb-potensialet fra atomkjernen. For et molekyl som består av to atomer vil derimot polarkoordinater være naturlige, da molekylet kun har symmetri om aksen som forbinder atomene.

I ikke-euklidske rom beskrevet ved riemannsk geometri, er slike koordinater påkrevet. Dette gjelder på krumme flater og i Einsteins generelle relativitetsteori. Den forklarer hvordan det firedimensjonale tidrommet er krummet av den masse og energi som det inneholder. Moderne kosmologi er derfor beskrevet i slike koordinatsystem.

Basisvektorer

Krumlinjete koordinater kan defineres i et euklidsk rom ved å uttrykke de tre kartesiske koordinater xm = (x,y,z) som ikke-lineære funksjoner av tre nye parametre xμ = (x1,x2,x3). Skrives dette som Mal:Nowrap må disse koordinattransformasjonene være «invertible» i et lokalt område av rommet, slik at man også kan beregne de inverse transformasjonene Mal:Nowrap For eksempel, når man benytter polarkoordinater

x=rcosθ
y=rsinθ,

i et todimensjonalt plan, så er den inverse transformasjonen

r=x2+y2
θ=arcsinyr

gyldig i halvplanet x > 0. I det andre halfplanet finnes en tilsvarende, invers transformasjon.

Hvis de kartesiske basisvektorene kalles for em , hvor den latinske indeksen m tar like mange verdier som dimensjonen på rommet, kan et vilkårlig punkt i rommet angis ved «posisjonsvektoren» Mal:Nowrap. Vektorrommet har en euklidsk metrikk, og basisvektoren har per definisjon et indreprodukt

𝐞m𝐞n=δmn

uttrykt ved Kronecker-deltaet på høyre side. Ved bruk av Einsteins summekonvensjon, som sier at man alltid skal summere over to like indekser i et matematisk uttrykk, kan posisjonsvektoren mer kompakt skrives som Mal:Nowrap. På den måten er dens form uavhengig av hvor mange dimensjoner rommet har. Indeksen til den kartesiske koordinaten xm kunne like godt stått nede, men det er hensiktsmessig å la den stå oppe.

Krumlinjet basis

Når man varierer bare én av koordinatene xμ  i ligningene xm = xm(xμ ), beskriver de koordinatlinjen for denne koordinaten. Da koordinattransformasjonene er ikke-lineære, vil dette i alminnelighet være en krum kurve. I hvert punkt kan man nå velge å bruke tangentvektorene Mal:Nowrap til disse koordinatlinjene som nye, lokale basisvektorer. Brukes det kartesiske uttrykket for posisjonsvektoren r, blir

𝐞μ=xmxμ𝐞m

når man benytter Einsteins summekonvensjon og summerer over like indekser. Denne krumlinjete basisen vil derfor være forskjellig fra punkt til punkt, gitt ved matrisen Mal:Nowrap som er Jacobi-matrisen for transformasjonen. Lokalt i hvert endelig område av rommet kan den inverse matrisen Mal:Nowrap finnest.

I en liten omegn om et punkt kan man nå uttrykke enhver vektor V  ved sine komponenter på denne lokale basisen. Det betyr at Mal:Nowrap hvor

Vμ=xμxmVm

Disse nye komponentene sies å være kontravariante, da de transformerer på motsatt vis sammenlignet med basisvektorene Mal:Nowrap. Alle vektorer skal ha komponenter som transformerer på denne måten. Men dette betyr også at i dette koordinatsystemet er ikke lenger posisjonsvektoren r noen vektor, ettersom den ikke innfrir denne betingelsen. For separasjonen mellom to nære punkt gjelder derimot Mal:Nowrap = Mal:Nowrap  hvor

dxμ=xμxmdxm

Derfor er den differensielle separasjonen Mal:Nowrap en ekte vektor.

Avstanden ds  mellom to nærliggende punkt separert med vektoren Mal:Nowrap følger fra Mal:Nowrap = Mal:Nowrap Dette «differensielle linjeelementet» kan derfor skrives som

ds2=gμνdxμdxν

hvor Mal:Nowrap er komponentene til den metriske tensoren som må benyttes for slike krumlinjete koordinater. Tilsammen utgjør de en matrise med elementer som følger direkte fra transformasjonsligningene,

gμν=xmxμxnxνδmn

Endelige avstander mellom forskjellige punkter i rommet må beregnes fra linjeelementet ved integrasjon. Den korteste avstanden mellom to punkt er en geodetisk kurve. Den tilsvarer en rett linje i det euklidske rommet.

Volumelement

Metrikken kan også benyttes til å beregne volumer i krumlinjete koordinater. Har rommet tre dimensjoner, betraktes først et lite parallellepiped definert ved de infinitesimale vektorene Mal:Nowrap, Mal:Nowrap  og Mal:Nowrap  langs koordinatlinjene gjennom et punkt. Volumet av dette er da gitt ved trippelproduktet

dV=(𝐞1×𝐞2)𝐞3dx1dx2dx3

Men fra transformasjonsligningene for basisvektorene følger at trippelproduktet her kan uttrykkes ved determinanten til Jacobi-matrisen, Mal:Nowrap = Mal:Nowrap. Dette igjen er forbundet til metrikken via Mal:Nowrap ved å bruke regelen for beregning av determinanten til produktet av matriser. På den måten blir det differensielle volumelementet

dV=gdx1dx2dx3

ved å bruke den mer kompakte notasjonen Mal:Nowrap Dette resultatet kan generaliseres til å være gyldig i rom med et vilkårlig antall dimensjoner.

Eksempel

Polarkoordinater (r,θ ) med ortogonale basisvektorer (er, eθ).

I polarkoordinater r = r cosθex + r sinθey  finnes de transformerte basisvektorene ved direkte derivasjon,

𝐞r=𝐫r=cosθ𝐞x+sinθ𝐞y𝐞θ=𝐫θ=rsinθ𝐞x+rcosθ𝐞y

Indreproduktene gir de metriske komponentene Mal:Nowrap Mal:Nowrap og Mal:Nowrap Basisvektorene står derfor vinkelrett på hverandre slik at den metriske matrisen er diagonal,

gμν=[100r2]

På denne måten blir Mal:Nowrap det differensielle linjeelementet i disse koordinatene. I andre koordinatsystem kan alltid den samme fremgangsmåten benyttes for å konstruere den nye koordinatbasisen.

I to dimensjoner gir det differensielle volumelementet arealet av tilsvarende lite flateelement Mal:Nowrap I polarkoordinater er determinanten Mal:Nowrap Skal man for eksempel i disse koordinatene beregne arealet av en sirkel med radius R, blir det

A=0Rrdr02πdθ=12R22π=πR2

Mer generelt kan alle lignendene integral regnes ut på tilsvarende måte når integrandene er uavhengige av polarvinkelen θ.

Dual basis og kovariante komponenter

Når koordinaten Mal:Nowrap holdes fast, beskriver denne ene ligningen en flate i det euklidske rommet. På samme måte beskriver de andre transformasjonsligningene tilsvarende flater. Hvert punkt i rommet ligger på skjæringspunktet mellom slike flater. I nærmest omegn kan man benytte normalvektorene Mal:Nowrap til disse flatene som et alternativt set med basisvektorer. De danner en dual basis som også blir kalt for en « kobasis» og blir skrevet med den greske indeksen i hevet posisjon. Hver slik vektor står vinkelrett på tangentvektorene Mal:Nowrap i de andre retningene. Det følger fra direkte fra definisjonene,

𝐞μ𝐞ν=xμ𝐫xν=xμxν=δνμ

De kontravariante komponentene til vektoren Mal:Nowrap  kan derfor bestemmes direkte fra indreproduktet Mal:Nowrap Alternativt kan man nå skrive vektoren uttrykt ved dens komponenter i den duale basisen, det vil si at Mal:Nowrap. Dette er vektorens kovariante komponenter Mal:Nowrap Men dette er også lik med Mal:Nowrap  som betyr at

Vμ=gμνVν

Metrikken Mal:Nowrap kan derfor brukes til å «senke» en kontravariant indeks slik at man får en kovariant komponent. På tilsvarende måte er en kontravariant komponent gitt som Mal:Nowrap hvor Mal:Nowrap utgjør elementene i den metriske matrisen

gμν=xμxmxνxnδmn

At begge indeksene står oppe, betyr at denne inneholder de kontravariante komponentene av metrikken. Denne er igjen den inverse av den kovariante metrikken Mal:Nowrap Det er en direkte konsekvens av definisjonene

gμλgλν=xμxmxλxmxnxλxnxν=xμxmxmxν=xμxν=δνμ

da Kronecker-deltaet på høyre side utgjør elementene i enhetsmatrisen.

Med disse to metriske matrisene kan man skifte mellom kontravariante og kovariante komponenter. For eksempel, indreproduktet mellom vektorene V og U som er Mal:Nowrap kan derfor skrives på de ekvivalente formene Mal:Nowrap

Koordinattransformasjoner

Hvis man velger å benytte et annet, krumlinjet koordinatsystem Mal:Nowrap i det samme rommet, vil man også ha en sammenheng Mal:Nowrap i de områdene hvor begge koordinatsystemene overlapper. Sammenhengene mellom deres basisvektorer er dermed

𝐞ν=𝐫xν=𝐫xμxμxν=xμxν𝐞μ

Det betyr at for en vektor Mal:Nowrap vil de kontravariante komponentene transformere som

Vν=xνxμVμ,

mens komponentene til metrikken Mal:Nowrap transformerer som

gμν=xλxμxμxνgλμ,

det vil si som de kovariante komponentene til en tensor av rang to. I krumlinjete koordinater blir derfor metrikken omtalt som en «metrisk tensor».

Gradient

I det euklidske rommet E3 er den vanlige gradienten til en skalar funksjon Φ(x) definert som

Φ=Φx𝐢+Φy𝐣+Φz𝐤=Φxm𝐞m

Det er her naturlig å skrive basisvektorene som em selv om det vanligvis ikke er av betydning. Men ved en transformasjon til de krumlinjete koordinatene oppstår forbindelsen

𝐞m=xmxμ𝐞μ

med den duale basisen. Gradienten tar da formen

Φ=Φxmxmxμ𝐞μ=Φxμ𝐞μ

i dette mer generelle koordinatsystemet når den defineres på denne måten. På dette viset er gradienten automatisk definert med kovariant komponenter i motsetning til en vanlig vektor som har kontravariante komponenter. Gradienten får dermed samme form i alle koordinatsystem. Det gir den en særstilling i vektoranalysen hvor den også kalles for en differensial 1-form.

Christoffel-symbol

Hvis man betrakter et vektorfelt Mal:Nowrap uttrykt i den kartesiske basisen, så forandrer det seg fra sted til sted i rommet gjennom varisjonene av komponentene. Hvis man i stedet benytter krumlinjete koordinater slik at Mal:Nowrap, så vil også basisvektorene Mal:Nowrap variere med posisjonen i rommet. Den totale forandringen av vektorfeltet får dermed to bidrag som kan uttrykkes ved en utvidet form for derivasjon. Den kalles den «kovariante deriverte» av vektorfeltet.

Hvordan komponentene varierer, kan uttykkes ved de partielle deriverte Mal:Nowrap hvor operatoren Mal:Nowrap På samme måte er forandringene til basisvektoren gitt ved de partielle deriverte Mal:Nowrap Denne vektoren må kunne uttrykkes som en lineærkombinasjon av de andre basisvektorene slik at man kan skrive

α𝐞μ=𝐞νΓμαν

når man igjen benytter Einsteins summekonvensjon og summerer over like indekser. De nye og foreløbig ukjente koeffisientene Mal:Nowrap kalles Christoffel-symbol etter den tyske matematiker Elwin Christoffel som innførte disse størrelsene for omtrent 150 år siden. Da Mal:Nowrap = Mal:Nowrap så er dette symbolet symmetrisk i sine to nedre indekser. Det vil si at man alltid har at Mal:Nowrap

Størrelsen på Christoffelsymbolene finnes fra metrikken Mal:Nowrap Deriveres denne i retning Mal:Nowrap får man

αgμν=(α𝐞μ)𝐞ν+𝐞μ(α𝐞ν)=gμλΓναλ+gνλΓμαλ

Dette forenkles ved å definere

Γμνα=gμλΓναλ

som kalles et «Christoffel-symbol av første sort», mens det tidligere symbolet Mal:Nowrap da blir et «Christoffel-symbol av andre sort». Den deriverte av metrikken kan da skrives som  Mal:Nowrap. To lignende uttrykk kan herav finnes ved å bytte om på indeksene og bruke de symmetriene som de har. Ved å innføre den alternative skrivemåten ∂αf = f,α for partiell derivasjon som derfor kan kalles «komma-derivasjon», gir disse tre relasjonene resultatet

Γλμν=12(gλμ,ν+gλν,μgμν,λ)

for Christoffel-symbolet. Det er også symmetrisk i de to siste indeksene. Selv om det ser ut som en tensor med tre indekser, er det bare tilsynelatende. En direkte utregning viser at det transformerer på en mer komplisert måte under en koordinattransformasjon.

Ved direkte utregning i todimensjonale polarkoordinater finner man Mal:Nowrap, mens Mal:Nowrap. Likedan er Mal:Nowrap, og Mal:Nowrap Vanligvis er flertallet av symbolene lik med null, men det er ikke lett uten videre å vite hvilke.

Kovariant derivasjon

Når variasjonen av basisvektorene fra sted til sted nå er etablert, kan man også finne hvordan et vektorfelt Mal:Nowrap  varierer når det beskrives i krumlinjete koordinater. I retning Mal:Nowrap skrives denne forandringen ved å benytte en generalisert nabla-operator. Når den virker på en vektor, gir den en ny vektor

α𝐀(x)=(αAμ)𝐞μ+Aμ(α𝐞μ)=(αAμ+AλΓλαμ)𝐞μ

når man i siste ledd bytter om på indeksene som det summeres over. De to leddene i parentesen kalles for den kovariante deriverte av denne vektorkomponenten. Det er vanlig å angi denne ved et semikolon eller ved bruk av det samme nabla-symbolet,

αAμA;αμ=αAμ+AλΓλαμ

Navnet «semikolon-derivert» blir derfor også benyttet. Mens den partielt deriverte av en vektorkomponent ikke har noen tensoregenskaper, er den kovariante deriverte av en vektorkomponent en komponent av en tensor med rang to.

Den kovariant deriverte av en tensor kan finnes ut fra dette. For eksempel, en tensor av andre rang transformerer som tensorproduktet av to vektorer A  og B. Som for all annen derivasjon gjelder også da

α(𝐀𝐁)=(α𝐀)𝐁+𝐀(α𝐁)

Dermed blir den kovariante deriverte av produktet av de kontravariante komponentene

α(AμBν)=α(AμBν)+AλBνΓλαμ+AμBλΓλαν

Mer generelt vil det samme mønster gjelde for den kovariante deriverte av en annenrangs tensor T med kontravariante komponenter Tμν. Kovariant derivasjon øker tensorens rang med en.

For å finne den kovariante deriverte av en kovariant vektorkomponent Aμ, kan man benytte at indreproduktet Mal:Nowrap er en skalar størrelse. Derfor må

α(AμAμ)=(αAμ)Aμ+Aμ(αAμ)=(αAμ)Aμ+Aμ(αAμ)

Benytter man her det etablerte resultatet for Mal:Nowrap følger at

αAμAμ;α=αAμAλΓμαλ

Herav kan kovariante deriverte av kovariante komponenter av tensorer av høyere rang finnes. Et viktig eksempel er den kovariante deriverte av den metriske tensoren som blir null,

αgμν=0.

Denne egenskapen til den kovariante deriverte kan utnyttes i uttrykk som ∇α Aμ som dermed kan beregnes fra Mal:Nowrap I motsetning til partiell derivasjon, opptrer metrikken som en konstant under kovariant derivasjon.

Retningsderivasjon

Derivasjonsoperatoren ∂α  benyttes til beregning av forandring i retning av basisvektoren Mal:Nowrap Bruken kan utvides til å gjelde i en vilkårlig retning gitt ved vektoren Mal:Nowrap Det kan gjøres ved å kombinere forandringene langs hver basisvektor. Det gir den mer generelle operatoren

𝐮uαα

som kalles for den retningsderiverte i denne retningen. Når den virker på en skalar funksjon f(x), er resultatet som ventet Mal:Nowrap Den gir også forandringen av et vektorfelt i samme retning ved at

𝐮𝐀uαα𝐀=uαA;αμ𝐞μ

som er en ny vektor. Samme operator kan også benyttes til å gi forandringen av vektorfeltet langs en kurve Mal:Nowrap i rommet. Da tangenten til kurven er Mal:Nowrap = Mal:Nowrap blir variasjonen av feltet langs denne gitt ved den totalderiverte

d𝐀dλ=𝐮𝐀

hvor retningsvektoren u nå har komponenter Mal:Nowrap som er gitt ved kurvens form. Denne sammenhengen involverer ingen indekser og er derfor gyldig i alle koordinatsystem. Så lenge det underliggende rommet er euklidsk, kan operatoren Mal:Nowrap defineres slik som her ved vanlige partialderiverte. Derimot i krumme rom beskrevet ved Riemanns differensialgeometri og i generell relativitetsteori har denne operatoren en mer fundamental betydning.

Fysiske komponenter

Fysiske komponenter av en vektor eller tensor angis vanligvis i et kartesisk koordinatsystem hvor alle basisvektorene har samme lengde. I krumlinjete koordinater kan man finne de samme komponentene ved å opprette et lite, lokalt koordinatsystem i hvert punkt med nye basisvektorer med samme lengde og som står vinkelrett på hverandre. I dette lokale aksekorset kan de fysiske komponentene avleses.

Dette er mest aktuelt å gjøre i rom med D = 3  dimensjoner og krumlinjete koordinater som har ortogonale basisvektorer. Da vil den metriske tensoren ha formen

gμν=[h12000h22000h32]

hvor 𝐞1𝐞1=h12 og tilsvarende for de to andre basisvektorene. Normerte basisvektorer med samme lengde eller «enhetsvektorer» er da 𝐞^1=𝐞1/h1, 𝐞^2=𝐞2/h2 og 𝐞^3=𝐞3/h3 slik at

𝐞^μ𝐞^ν=δμν

Videre gjelder de fundamentale vektorproduktene 𝐞^1×𝐞^2=𝐞^2×𝐞^1=𝐞^3 og de andre som følger fra syklisk ombytte av indekser. Disse sammenhengene holder bare i tredimensjonale rom.

Alternativt kan man benytte de duale basisvektorene slik at 𝐞^1=h1𝐞1=h1x1 og tilsvarende for de to andre. En vektor som 𝐀=A1𝐞1+A2𝐞2+A3𝐞3 kan da skrives i denne nye basisen som 𝐀=h1A1𝐞^1+h2A2𝐞^2+h3A3𝐞^3. Sammenhengen mellom de kontravariante og de fysiske kompontene til vektoren er dermed A1^=h1A1, A2^=h2A2 og A3^=h3A3. På samme måte gjelder for de kovariante komponentene at A1^=A1/h1, A2^=A2/h2 og A3^=A3/h3.

Denne formalismen kan benyttes i mange praktiske sammenhenger, både i relativistisk og ikke-relativistisk fysikk. For eksempel kan Navier-Stokes-ligningene skrives i andre koordinater enn kartesiske når det er mer praktisk.

Eksempel med partikkelbevegelse

Posisjonsvektoren 𝐫 til en partikkel som beveger seg i et plan, forandrer seg med tiden t. Den beskriver derfor en kurve 𝐫=𝐫(t) hvis tangentvektor 𝐯=d𝐫/dt er lik med hastigheten til partikkelen. Foregår bevegelsen i et plan, kan det være en fordel å beskrive bevegelsen ved bruk av polarkoordinater x=rcosθ,y=rsinθ. Kvadratet av linjeelementet er dermed ds2=dr2+r2dθ2 slik at den diagonale metrikken er gitt ved hr=1,hθ=r.

I dette systemet kan hastighetsvektoren skrives 𝐯=r˙𝐞r+θ˙𝐞θ hvor de kontravariante komponentene til hastigheten er vr=dr/dtr˙ og vθ=dθ/dtθ˙. Akselerasjonen til partikkelen er definert som den retningsderiverte 𝐚=d𝐯/dt i retning av selve hastigheten. Denne følger fra den kovariante deriverte og blir

𝐚=𝐯˙=(v˙λ+Γμνλv˙μv˙ν)𝐞λ=(r¨+Γθθrθ˙θ˙)𝐞r+(θ¨+2Γrθθθ˙r˙)𝐞θ

når man setter inn for Christoffelsymbolene Γθθr=r,Γrθθ=1/r som bidrar. De kontravariante komponentene til akselerasjonen kommer dermed frem som ar=r¨rθ˙2 og aθ=θ¨+2θ˙r˙/r.

Ifølge Newtons 2. lov er akselerasjonen proporsjonal med kraften som virker på partikkelen. I radiell retning er den fysiske komponenten av akselerasjon ar^=hrar=ar=r¨rθ˙2. Det siste leddet her gir opphav til den fiktive sentrifugalkraften som virker på partikkelen. Tilsvarende er den fysiske komponenten av akselerasjonen i angulær retning aθ^=hθaθ=rθ¨+2θ˙r˙ hvor det siste leddet representerer den fiktive Coriolis-kraften.

Hvis den fysiske kraften er gitt ved Newtons gravitasjonslov, virker denne i radiell retning. I dette viktige tilfellet er derfor den fysiske komponenten aθ^=0. Da det kan skrives som

ddtr2θ˙=0,

har det som konsekvens at r2θ˙ er konstant under bevegelsen. Dette gjelder for en planets bevegelse rundt Solen og uttrykker Keplers 2. lov om flatehastighetens konstans. Den radielle komponenten til akselerasjonen kan så benyttes til å utlede de to andre Kepler-lovene.

Vektorderivasjon

I kartesiske koordinater blir vektorderivasjon utført ved hjelp av nabla-operatoren. Den enkleste operasjonen er gradienten av en skalar funksjon Φ=Φ(𝐫). De partialderiverte μΦ er de kovariante komponentene av gradienten til funksjonen. Den første, fysiske komponenten av gradienten er derfor 1Φ/h1 og likedan for de andre. I denne normerte basis er derfor gradienten gitt som

Φ=1h1Φx1𝐞^1+1h2Φx2𝐞^2+1h3Φx3𝐞^3

Mer komplisert er beregningen av divergensen til en vektor 𝐀=A1^𝐞^1+A2^𝐞^2+A3^𝐞^3. Den første termen gir et bidrag (A1^𝐞^1) som vil involvere den kovariante deriverte av basisvektoren 𝐞^1. Men spesielt i D = 3  kan utregningen forenkles ved å benytte at 𝐞^1=h2h3x2×x3 som gir

(A1^𝐞^1)=(A1^h2h3x2×x3)=(A1^h2h3)x2×x3+A1^h2h3(x2×x3)

Det siste leddet her gir null da curl til en gradient er alltid lik null. I det første leddet kan gradienten nå skrives ut og gir

(A1^𝐞^1)=(1h1A1^h2h3x1𝐞^1+1h2A1^h2h3x2𝐞^2+1h3A1^h2h3x3𝐞^3)𝐞^1h2h3=1h1h2h3A1^h2h3x1

Ved å ta med de tilsvarende bidragene fra de to andre komponentene til vektoren, er da dens divergens uttrykt ved deriverte av de fysiske komponentene, gitt som

𝐀=1h1h2h3(A1^h2h3x1+A2^h3h1x2+A3^h1h2x3)

På samme måte kan curl til vektoren regnes ut. Den første komponenten gir nå

×(A1^𝐞^1)=×(A1^h1x1=(A1^h1)×x1+A1^h1×x1

Igjen er det siste leddet en curl av en gradient og er dermed null. Ved å skrive gradienten i det første leddet ut, blir det

×(A1^𝐞^1)=(1h1A1^h1x1𝐞^1+1h2A1^h1x2𝐞^2+1h3A1^h1x3𝐞^3)×𝐞^1h1=1h3h1A1^h1x3𝐞^21h1h2A1^h1x2𝐞^3

Bidragene fra de to andre leddene i vektoren vil ha en lignende form. Legges de sammen, kan resultatet uttrykkes som en determinant og skrives på den kompakte formen

×𝐀=1h1h2h3|h1𝐞^1h2𝐞^2h3𝐞^3x1x2x3A1^h1A2^h2A3^h3|

på samme måte som for curl i kartesiske koordinater.

Det er nå ganske rett frem å finn Laplace-operatoren i slike kurvelineære koordinater. Den er definert som divergensen av gradienten til en funksjon og er derfor

2Φ=Φ=1h1h2h3[x1(h2h3h1Φx1)+x2(h3h1h2Φx2)+x3(h1h2h3Φx3)]

Eksempel med kulekoordinater

I kulekoordinater x1=r,x2=θ,x3=ϕ er den diagonale metrikken gitt ved h1=hr=1,h2=hθ=r,h3=hϕ=rsinθ slik at gradienten til en skalar funksjon blir

Φ=Φr𝐞^r+1rΦθ𝐞^θ+1rsinθΦϕ𝐞^ϕ

Likedan blir divergensen til et vektorfelt

𝐀=1r2Ar^r2r+1rsinθAθ^sinθθ+1rsinθAϕ^ϕ,

mens ×𝐀 følger fra å skrive ut den tilsvarende determinanten i det generelle resultatet. Det fullstendige uttrykket for Laplace-operatoren blir likedan

2Φ=1r2r(r2Φr)+1r2sinθθ(sinθΦθ)+1r2sin2θ2Φϕ2.

Se også

Litteratur

  • Ø. Grøn and A. Næss, Einstein's Theory, Springer, New York (2011). ISBN 978-1-4614-0706-5.
  • M.L. Boas, Mathematical Methods in the Physical Sciences, John Wiley & Sons, New York (1983). ISBN 0-471-04409-1.
  • M. R. Spiegel, Vector Analysis, Schaum's Outline Series, New York, (1959).
  • E. Kreyzig, Differential Geometry, Dover Publications, New York (1991). ISBN 0-486-66721-9.

Eksterne lenker

Mal:Autoritetsdata