Oppgave med Stokes' teorem

[krje1980]

Hei.

Jeg har benyttet sommerferien til å repetere litt av multivariabel kalkulus pensumet jeg hadde for halvannet år siden. Forsøkte meg på en oppgave om Stokes' Teorem, men kommer ikke helt i mål:


Use Stokes's Theorem to show that

[tex]\oint_{C} y ~dx + z ~dy + x ~dz = \sqrt{3} \pi a^2[/tex]

where [tex]C[/tex] is the suitably oriented intersection of the surfaces [tex]x^2 + y^2 + z^2 = a^2[/tex] and [tex]x + y + z = 0[/tex].

OK, så Stokes' teorem sier:

[tex]\oint_{C}\vec{F} \cdot d\vec{r} = \iint_{S}\operatorname{curl} \vec{F} \cdot \vec{N} ~dS[/tex]

Jeg har regnet ut at:

[tex]\operatorname{curl} \vec{F} = -\vec{i} - \vec{j} - \vec{k}.[/tex]

Videre har jeg regnet ut at på overflaten [tex]S[/tex] har vi:

[tex]\vec{N}dS = \frac{\nabla G(x,y,z)}{G_{3}(x,y,z)}dxdy = \vec{i} + \vec{j} + \vec{k}dxdy[/tex]

(dette følger av ligningen for det oppgitte planet).

Da får jeg imidlertid:

[tex]\operatorname{curl} \vec{F} \cdot \vec{N}dS = -1 -1 -1 = -3[/tex]

Og da får jeg, når jeg projiserer dette ned på [tex]xy[/tex]-planet:

[tex]\iint_{S} \operatorname{curl} \vec{F} \cdot \vec{N} ~dS = -3 \iint_{A} dA = -3 \pi a^2[/tex]

som jo ikke stemmer.

Dersom noen kan se hva jeg gjør feil her, så vil jeg sette stor pris på det .

[Gustav]

Retorisk spørsmål: Har du normalisert enhetsnormalvektoren til planet?

Fortegnet har vel å gjøre med retningen man integrerer, og dette er jo avgjørende for retningen til enhetsnormalen. Du må se litt på hva som blir riktig der.

[krje1980]

Hei.

Jeg ser jo at ved å normalisere normalvektoren, så vil jeg få riktig svar (jeg ser jo også nå at jeg selvsagt kan skifte fortegn).

Men i følge min tekstbok så er:

[tex]\vec{N} = \pm \frac{\nabla G(x,y,z)}{|\nabla G(x,y,z)|}[/tex]

og:

[tex]dS = |\frac{\nabla G(x,y,z)}{G_{3}(x,y,z)}|[/tex]

Da får vi:

[tex]\vec{N} \bullet dS = \frac{1}{sqrt{3}}(\vec{i} + \vec{j} + \vec{k}) \cdot \sqrt{3} = \vec{i} + \vec{j} + \vec{k}[/tex]

Så dersom jeg skifter fortegn så får jeg svaret [tex]3 \pi a^2[/tex], men det er jo fremdeles ikke rett. Som sagt, jeg ser at normalisering av normalvektoren gir riktig svar, men vi må vel også ta hensyn til [tex]dS[/tex]?

[Gustav]

Ah, var litt kjapp der ja. Hjernen er åpenbart i feriemodus. For trøtt til å tenke også. Får se på det i morgen.

[krje1980]

OK. Takk skal du ha . Jeg skal på lang fjelltur i morgen, så du får unnskylde hvis jeg ikke får svart deg før i morgen kveld i så fall.

[krje1980]

Har du fått sett noe mer på dette, plutarco?

[Gustav]

Nå vil vel curl F stå normalt på planet og i motsatt retning av positiv enhetsnormalvektor, altså blir [tex]\iint curl F \cdot d\vec{s}=\int_{-a}^a\int_{-\sqrt{a^2-v^2}}^{\sqrt{a^2-v^2}} -|(-1,-1,-1)|dudv[/tex].

(Fortegnet avhenger jo av hvilken retning kurven C går i.)

Du må nesten skifte til koordinater u og v gitt ved å rotere xy-planet 45 grader om z-aksen til koordinater x´y´og deretter 45 grader om x´.

Projeksjonen av snittet av kula og planet ned på xy-planet blir vel ikke helt enkelt å finne? Kanskje feilen ligger der. Hvilke integrasjonsgrenser har du brukt?

[Masamune]

Jeg tror det er enklest å ikke projisere ned i xy-planet. Integralet [tex] \int_S dS [/tex] er jo lik [tex]\pi a^2[/tex].

[krje1980]

Takk til dere begge.

Dere har nok et poeng i at det jeg har gjort når jeg har forsøkt å projisere overflaten til xy-planet er feil (jeg satt bare [tex]z = 0[/tex] i ligningen for kulen, og fikk da [tex]x^2 + y^2 = a^2[/tex]). Men dette blir jo selvsagt ikke rett ettersom planet ikke skjærer kulen med så stor radius.

Nå er som sagt dette noe som jeg ikke har syslet med på nesten to år, så jeg sparker meg selv litt bak nå for at jeg ikke har holdt dette mer aktivt ved like. Jeg skal forsøke å se videre på oppgaven.

[Gustav]

Takk til dere begge.

Dere har nok et poeng i at det jeg har gjort når jeg har forsøkt å projisere overflaten til xy-planet er feil (jeg satt bare [tex]z = 0[/tex] i ligningen for kulen, og fikk da [tex]x^2 + y^2 = a^2[/tex]). Men dette blir jo selvsagt ikke rett ettersom planet ikke skjærer kulen med så stor radius.

Nettopp. Projeksjonen av snittkurven blir jo [tex]x^2+y^2+(-x-y)^2=a^2[/tex].

Det er jo mye enklere å bare lage seg nye koordinater u og v som "utspenner" planet x+y+z=0. Integrasjonen blir da slik jeg skrev over.

[krje1980]

Takk skal du ha! Nå er jeg med . Er jo også riktig som Masamune påpeker, at jeg kan bare regne ut integralet av dS uten å projisere ned i xy-planet. Normaliserer vi normalvektoren får vi da det ønskede svaret.