17 svar
283 visningar
naytte behöver inte mer hjälp
naytte 7419 – Moderator
Postad: 21 aug 14:16 Redigerad: 21 aug 14:20

Hur ska man ställa upp flödesintegralen över ytan S?

Halloj!

Jag sitter med frågan nedan från en gammal tentamen i flervarabielanalys:

Jag antar att ytan, låt oss kalla den SS, är ytan som utgörs av snittet av planet och cylindern som definieras i uppgiften. Jag har försökt göra en skiss för lite visualisering:

Tydligen är jag inte så bra på att rita men det är meningen att skivan ska "luta" på samma sätt som planet x+y+z=1x+y+z=1 gör. Jag tolkar det som att jag ska ställa upp någon integral på formen:

SFx,y,zn^dS\displaystyle \iint_S \mathbf{F}\left(x,y,z\right)\hat{\mathbf n} dS

Men hur ska man välja n^\hat{\mathbf n}? Jag antar att parantesen "(i riktning uppåt zz-axel)" ska ange hur ytan är orienterad men jag är inte helt med på hur man ska tolka det.

D4NIEL 3345
Postad: 21 aug 14:34 Redigerad: 21 aug 14:47

Riktning uppåt z-axel innebär att vi vill välja normalen till ytan så att z-elementet i uttrycket för normalen får ett positivt värde. (1,1,1) pekar i rätt riktning (-1,-1,-1) pekar i fel riktning. Tänk också på att få rätt normering för det du vill göra.

Eftersom ytan ligger i ett plan är det kanske inte jättesvårt att se vad normalen blir (den är dessutom konstant över hela ytan). Men för övning kan det vara nyttigt att

1. Parametrisera ytan

2. Bestäm dess normal (antingen med beräkning eller genom utantillkunskap om normalens utseende för din valda parametrisering)

Integralen du vill beräkna är alltså skalärprodukten av fältet och ytelementet

SF·dS\displaystyle \int_S \mathbf{F}\cdot \mathrm{d}\mathbf{S}

naytte 7419 – Moderator
Postad: 21 aug 15:22 Redigerad: 21 aug 15:23

Okej, jag ger en parametrisering ett försök. Jag tänker att vi vill beskriva ytan SS med en vektorvärd funktion r\mathbf{r} enligt r(s,t)=x(s,t)i+y(s,t)j+z(s,t)k\mathbf{r}(s,t)=x(s,t)\mathbf i+y(s,t)\mathbf j+z(s,t)\mathbf k. Eftersom "botten" till ytan är en cirkel med radie 22 verkar det rimligt att vi sveper över denna i xyxy-planet och inför parametriseringen:

r(r,θ)=rcosθi+rsinθj+(1-r(sinθ+cosθ))k\displaystyle \mathbf{r}(r,\theta)=r\cos\theta\mathbf i+r\sin\theta\mathbf j+(1-r(\sin\theta+\cos\theta))\mathbf k

där 0θ2π0\le \theta \le 2\pi och 0r20 \le r\le 2

Ser detta rätt ut? Jag är osäker på hur man bestämmer normalen när man väl har en parametrisering. Som du nämner råkar det ju vara lätt i just det här fallet eftersom det är "ett plan" men generellt vet jag inte hur man gör. Jag misstänker att det finns en kryssprodukt med någonstans i beräkningen eftersom vi vill ha en vinkelrät vektor.

Ytelementet dSdS avbildas vid denna transformation på:

dS=rr×rθdrdθ\displaystyle d{S}=\left|\frac{\partial\mathbf{r} }{\partial r} \times\frac{\partial \mathbf{r}}{\partial \theta} \right| drd\theta

Så det verkar rimligt (kanske?) att:

dS=rr×rθdrdθ\displaystyle d\mathbf{S}= \left(\frac{\partial\mathbf{r} }{\partial r} \times\frac{\partial \mathbf{r}}{\partial \theta}\right) drd\theta

Men nu är jag ute på djupt vatten.

D4NIEL 3345
Postad: 21 aug 15:41 Redigerad: 21 aug 15:51

Det är helt korrekt! Snyggt! :-)

Att det fungerar så beror alltså på att enhetsnormalen NN ges av

N=rr×rθ|rr×rθ|N=\frac{\frac{\partial \mathbf{r}}{\partial r}\times\frac{\partial \mathbf{r}}{\partial \theta}}{|\frac{\partial \mathbf{r}}{\partial r}\times\frac{\partial \mathbf{r}}{\partial \theta}|}

Sätter vi ihop det med ditt uttryck för dS\mathrm{d}S får vi

dS=NdS=rr×rθdrdθ  =r,r,rdrdθ\mathrm{d}\mathbf{S}=N\mathrm{d}S=\left(\frac{\partial \mathbf{r}}{\partial r}\times\frac{\partial \mathbf{r}}{\partial \theta}\right)\,\mathrm{d}r\mathrm{d}\theta \quad \left(=\left(r,r,r\right)\mathrm{d}r \mathrm{d}\theta\right)

Men i just den här uppgiften kan det som sagt vara enklare att välja parameterframställningen

r=(x,y,f(x,y))\mathbf{r}=(x,y,f(x,y)) och med tillägget (x2+y2<4)(x^2+y^2<4) och sedan bekymra sig om vad det betyder senare.

Med en sådan parameterframställning i (x,y)(x,y) brukar man efter några uppgifter lära sig utantill att 

rx×ry=-fx',-fy',1\frac{\partial \mathbf{r}}{\partial x}\times \frac{\partial\mathbf{r}}{\partial y}=\left(-f^\prime_x, -f^\prime_y, 1\right)

naytte 7419 – Moderator
Postad: 21 aug 18:16 Redigerad: 21 aug 18:27

Jag är lite förvirrad över vad vi ens har lyckats göra här. Vi kom fram till dS=(r,r,r)drdθd\mathrm{\mathbf S}=(r,r,r)drd\theta. Men vi behövde aldrig bekymra oss över vilken normal vi valde. Hur vet vi att vi räknade med "rätt" dSd\mathrm{\mathbf S}?

EDIT:

och en grej till. Vi valde ju ordningen rr×rθ\frac{\partial \mathbf{r}}{\partial r}\times \frac{\partial \mathbf{r}}{\partial \theta}, men principiellt hade vi väl lika gärna kunnat ta rθ×rr\frac{\partial \mathbf{r}}{\partial \theta}\times \frac{\partial \mathbf{r}}{\partial r}?

PATENTERAMERA 7132
Postad: 21 aug 18:31

Ja, men tänk på att du skall välja den med positiv z-komponent enligt problemtexten.

naytte 7419 – Moderator
Postad: 21 aug 18:44 Redigerad: 21 aug 18:45

Men hur vet man från början vilken ordning som blir rätt automatiskt?

PATENTERAMERA 7132
Postad: 21 aug 19:21

Det kan vara lite knepigt att på förhand inse vilken som är rätt. Enklast är nog att välja en och kolla om det blir rätt riktning när man räknat klart. Om inte, så byt bara tecken på din framräknade normal.

naytte 7419 – Moderator
Postad: 21 aug 21:45 Redigerad: 21 aug 21:55

Tolkar jag det rätt som att vi fick ut att enhetsnormalen är N=(1,1,1)N=(1,1,1)?


Tillägg: 21 aug 2025 22:04

EDIT: ser nu att det såklart inte blir en enhetsnormal, utan endast EN normal.

naytte 7419 – Moderator
Postad: 21 aug 22:08 Redigerad: 21 aug 23:58

Men hur som helst, om vi låter n^x,y,z\hat{\mathbf{n}}\left(x,y,z\right) vara en funktion som ger enhetsnormalen till SS i varje punkt (x,y,z)(x,y,z) borde vi alltså ha:

SFx,y,zn^x,y,zdS=EFrr,θn^rr,θJrdrdθ\displaystyle \iint_S \mathbf{F}\left(x,y,z\right)\hat{\mathbf{n}}\left(x,y,z\right)d{S}=\iint_{E}\mathbf{F}\left(\mathbf{r}\left(r,\theta\right)\right)\hat{\mathbf{n}}\left(\mathbf{r}\left(r,\theta\right)\right) J_\mathbf{r} drd\theta

för något område EE i (r,θ)(r,\theta)-planet, eller hur?

PATENTERAMERA 7132
Postad: 21 aug 22:10

naytte 7419 – Moderator
Postad: 21 aug 23:49 Redigerad: 22 aug 00:00

Okej, jag tror jag hänger med. Vi har det orienterade ytelementet dS=n^dSd\mathbf{S} = \hat{\mathbf{n}}dS. n^\hat{\mathbf{n}} är här enhetsnormalen till ytan SS i varje punkt (x,y,z)(x,y,z). När vi sedan inför vår parametrisering avbildas dSd\mathbf{S} på en annan vektor, och då slutar det vara relevant huruvida vektorn vi får är normal mot någon yta eller har längd 1, ty vi arbetar nu i ett enkelt tvådimensioellt område.

Är det rätt uppfattat?

D4NIEL 3345
Postad: 22 aug 02:05 Redigerad: 22 aug 03:04

Jag tror du har grepp om det nu, men jag blir lite stressad när du skriver att det slutar vara relevant huruvida vektorn vi får är en normal  eller hur lång den är :-) Så för säkerhets skull tänker jag försöka knyta parameterframställningen med Jakobianen (som jag tror du börjar bli van vid). 

Vektorn rr\displaystyle \frac{\partial \mathbf{r}}{\partial r} motsvarar första kolonnen i Jakobianen för den tvådimensionella avbildningen. Det är en basvektor för tangentplanet till ytan i punkten

Vektorn rθ\displaystyle \frac{\partial \mathbf{r}}{\partial \theta} motsvarar andra kolonnen i Jakobianen. Det är den andra basvektorn till tangentplanet till ytan i punkten.

När du kryssar två vektorer får du en ny vektor som är vinkelrät mot båda vektorerna. I det här fallet får du alltså en vektor som är vinkelrät mot tangentplanet till ytan i punkten. Den "onormerade" normalen

N=rr×rθ\displaystyle N=\frac{\partial \mathbf{r}}{\partial r}\times \frac{\partial \mathbf{r}}{\partial \theta}

har därmed exakt rätt längd för att kompensera avbildningen. Om du vill kan du istället se det som att du kryssar två differentialer för att se hur stor infinitesimal area de spänner upp

N^dS=(rrdr)×(rθdθ)\displaystyle \hat{N}\mathrm{d}S =(\frac{\partial \mathbf{r}}{\partial r} \mathrm{d}r)\times (\frac{\partial \mathbf{r}}{\partial \theta}\mathrm{d}\theta)

Om du ansätter din parameterframställning redan från början ska du alltså inte multiplicera med någon extra skalfaktor |JT||J_T| i din integral. Allt sköts av din parameterframställning. Om du däremot mellanlandar i r(x,y)r(x,y) får du först en skalfaktor i normalen (-fx',-fy',1)(-f^\prime_x, -f^\prime_y,1) och sedan eventuellt en "vanlig" Jakobian när/om du byter till polära (om du inte ser att det bara blir arean av en cirkel med radien 2 kvar) .

Slutligen, omvänd ordning av vektorerna i kryssprodukten innebär att normalen byter tecken. Det finns alltså en sorts automatisk orientering av ytan i den uppräkningsordning du väljer för dina parametrar.

naytte 7419 – Moderator
Postad: 22 aug 11:16 Redigerad: 22 aug 11:35

Okej, så för att vara tydlig har vi alltså:

SFx,y,z·dS=EFrr,θ·rr×rθdrdθ\displaystyle \iint_S \mathbf{F}\left(x,y,z\right)\cdot d\mathbf{S}=\iint_{E} \mathbf{F}\left(\mathbf{r}\left(r,\theta\right)\right)\cdot\left(\frac{\partial\mathbf{r}}{\partial r} \times \frac{\partial \mathbf{r}}{\partial \theta}\right)drd\theta

där EE är något område i (r,θ)(r,\theta)-planet.

Det jag menade med att det inte spelar någon roll om vektorn vi får är "normal" eller "normerad" är att vektorn:

rr×rθ=r,r,r\displaystyle \left(\frac{\partial\mathbf{r}}{\partial r} \times \frac{\partial \mathbf{r}}{\partial \theta}\right) = \left(r,r,r \right)

inte är normal mot något i (r,θ)(r,\theta)-planet eller normerad. Är det korrekt uppfattat?

EDIT:

Och angående orientering, går det att veta från början vilken parameterordning som ger rätt orientering? I facit pratade de inte ens om ordning utan de bara körde på enligt omskrivningen ovan i det här inlägget (fast med en annan parametrisering).

D4NIEL 3345
Postad: 22 aug 12:24 Redigerad: 22 aug 12:38

Du skapar ett nytt rθr\theta-plan för varje punkt i din parameterframställning. Du har kommit fram till att en normal med exakt rätt längd för just din valda parameterframställning är N=(r,r,r)N=(r,r,r) för varje värde på rr. Din normalvektor är alltså helt vinkelrät mot varje rθr\theta-plan du skapar.

Vektorn (1,1,1)(1,1,1) och vektorn (r,r,r)(r,r,r) är alltid helt parallella eftersom du valt att låta rr löpa mellan 0 och 2 (vi kan bortse från vad som händer då r=0r=0).

Vektorn (r,r,r)(r,r,r) har ett positivt z-element så länge du låter r ligga i ditt valda parameterintevall. Du har därmed orienterat ytan så att positiv riktning är "uppåt" z-axeln. Vidare har din vektor exakt rätt norm för vad du ska använda den till. Men den är inte normerad i bemärkelsen att ||(r,r,r)||=1||(r,r,r)||=1, tvärtom varierar normen av vektorn över ditt parameterspann.

Som PATENTERAMERA påpekar kan det ibland vara mer praktiskt att räkna på och kontrollera att ytans normal har rätt orientering i efterhand. Skulle man ha råkat göra fel är det bara att byta uppräkningsordning eller tecken. Men basvektorerna utgör tillsammans med normalen ett högerorienterat system. På samma sätt som x^×y^=z^\hat{x}\times\hat{y}=\hat{z} är r^×θ^=N^\hat{r}\times \hat{\theta}=\hat{N}. Det går att visualisera på förhand åt vilket håll normalen pekar.

Edit: Det slår mig nu i efterhand att du med rθr\theta-planet kanske menar parametermängden EE av rθr\theta-värden. Och i den betydelsen kan du se de det som att definitionsområdet eller punktmängden EE inte är "vinkelrät" mot något.

naytte 7419 – Moderator
Postad: 22 aug 15:17 Redigerad: 22 aug 17:53

Det jag menade var att när vi gör vår substitution integrerar vi istället över ett enkelt 2D-område i rθr\theta-världen där 0r20\le r \le 2 och 0θ2π0\le \theta \le 2\pi. Mot denna yta kommer (r,r,r)(r,r,r) aldrig vara vinkelrät, är vad jag menar. Visst stämmer det?

EDIT: men vad spelar det för roll att (r,r,r)(r,r,r) har positivt zz-element? Är det inte enhetsnormalen till ytan SS vi kräver ska ha positivt zz-element? 

D4NIEL 3345
Postad: 23 aug 02:06 Redigerad: 23 aug 02:18

Ja, vi kan inte jämföra geometriska objekt mellan de olika mängderna.

r\mathbf{r} är en kontinuerlig funktion from 2\mathbb{R}^2 till 3\mathbb{R}^3. Definitionsområdet EE utgör ett slutet, begränsat område. Värdemängden för r\mathbf{r} blir då en sammanhängande begränsad punktmängd SS i 3\mathbb{R}^3 (närmare bestämt en yta). Ungefär så här:

Det blir nonsens om vi ska prata om hur de olika mängderna EE och SS förhåller sig till varandra geometriskt eller hur vektorer i 3\mathbb{R}^3 förhåller sig till EE

Om vektorn N=(r,r,r)N=(r,r,r) har ett positivt zz-element måste även enhetsvektorn ha det, de pekar ju åt samma håll!

Vi kan normera vektorn N=(r,r,r)N=(r,r,r). Då blir det

N^=N||N||=(r,r,r)3r2=13(1,1,1)\hat{N}=\frac{N}{||N||}=\frac{(r,r,r)}{\sqrt{3r^2}}=\frac{1}{\sqrt{3}}(1,1,1)

eftersom r0r\geq 0. Är du med?

Vidare, vi säker alltså normalen som pekar "uppåt", eftersom det i uppgiften står att vi söker flödet "i riktning uppåt zz-axel" 

naytte 7419 – Moderator
Postad: 23 aug 02:21

Ah, nu är jag med! Jag förvirrade mig själv lite men nu är jag med till 100 %!

Tack för all hjälp! :D

Svara
Close