Side 1 av 1

Beregning av vektorfluks gjennom en kule

Lagt inn: 20/04-2009 13:53
av chrtsta
Heisann.
Jeg skal regne ut vektorfluksen gjennom en kule plassert i origo med radius a. Vektorfeltet er gitt ved
[tex]F = x\vec{i} + (2y+z)\vec{j} + (z+x^2)\vec{k}[/tex].
Vektorfluksen er gjerne gitt ved
[tex]Q = \int_{\sigma} \vec{F} \cdot \vec{n} \, d\sigma[/tex].
Noen som kan hjelpe meg på vei med dette?

Lagt inn: 20/04-2009 14:12
av Gustav
Gjør om til kulekoordinater.

Har at

[tex]dA=a^2\sin(\phi)d\theta d\phi[/tex] så

[tex]\iint_{x^2+y^2+z^2=a^2}\vec{F}\cdot \vec{n} a^2\sin(\phi)d\theta d\phi[/tex]

Sett [tex]x=a\sin(\phi)\cos(\theta) \\ y=a\sin(\phi)\sin(\theta) \\ z=a\cos(\phi)[/tex] i uttrykkene for vektorene F og n og beregn integralet.

Lagt inn: 20/04-2009 14:19
av chrtsta
Hm, ok, men hva mener du med dA? Og er egentlig normalvektoren som kanskje er det største problemet mitt. Hva er uttrykket for denne?

Lagt inn: 20/04-2009 14:30
av Gustav
[tex]dA=d\sigma[/tex], altså infinitesimalt flateareal på kuleflata.

Normalvektoren blir vel

[tex]\frac{1}{a}<x,y,\sqrt{a^2-x^2-y^2}>[/tex] på øvre halvkule.

Lagt inn: 20/04-2009 15:07
av Gustav
[tex]\vec{F}\cdot \vec{n}=\frac{1}{a}<x,2y+z,z+x^2><x,y,\sqrt{a^2-x^2-y^2}>[/tex].

Har at [tex]x^2+y^2=a^2\sin^2(\phi)[/tex] så da får vi

[tex]\frac{1}{a}(x^2+2y^2+yz+z\sqrt{a^2-x^2-y^2}+x^2\sqrt{a^2-x^2-y^2})\\=\frac{1}{a}(a^2\sin^2(\phi)+y^2+yz+a(z+x^2)\cos(\phi))[/tex]

[tex]az\cos(\phi)=a^2\cos^2(\phi)[/tex] så

uttrykket forenkles til

[tex]\frac{1}{a}(a^2\sin^2(\phi)+y^2+yz+a(z+x^2)\cos(\phi))=\frac{1}{a}(a^2(\sin^2(\phi)+\cos^2(\phi))+y^2+yz+ax^2\cos(\phi))\\=\frac{1}{a}(a^2+y^2+yz+ax^2\cos(\phi))[/tex].

[tex]=a+a\sin^2(\phi)\sin^2(\theta)+a\sin(\phi)\cos(\phi)\sin(\theta)+a^2\sin^2(\phi)\cos^2(\theta)\cos(\phi)[/tex]

Huff, voldsomt mye og lett og gjøre feil :)

Lagt inn: 20/04-2009 17:38
av chrtsta
Takk for det!
Jeg prøvde å regne ut integralet, men jeg fikk feil svar.
Det som derimot ser ut til å bli riktig er en mye enklere metode: å bruke Gauss sats (divergensteoremet). Da har vi
[tex]\int_\sigma \vec{F} \cdot \vec{n} \, d\sigma = \int_\tau \nabla \cdot \vec{F} \, d\tau[/tex], hvor [tex]\int_\tau 1 \, d\tau = V_\tau = V_{kule}[/tex],
så jeg får [tex]Q = 4 V_{kule} = \frac{16 \pi a^3}{3}[/tex]

Lagt inn: 20/04-2009 20:48
av Gustav
chrtsta skrev:Takk for det!
Jeg prøvde å regne ut integralet, men jeg fikk feil svar.
Det som derimot ser ut til å bli riktig er en mye enklere metode: å bruke Gauss sats (divergensteoremet). Da har vi
[tex]\int_\sigma \vec{F} \cdot \vec{n} \, d\sigma = \int_\tau \nabla \cdot \vec{F} \, d\tau[/tex], hvor [tex]\int_\tau 1 \, d\tau = V_\tau = V_{kule}[/tex],
så jeg får [tex]Q = 4 V_{kule} = \frac{16 \pi a^3}{3}[/tex]
Smart!

Hadde vært artig å fått rett svar med den tungvinte metoden også da;) Men opplagt lettest å bruke satsen.

Husket du på å bruke rett enhetsnormalvektor på den nedre halvkula?

Lagt inn: 26/10-2009 14:24
av tool-nes
Låner denne tråden litt, i og med at jeg sitter med nøyaktig samme oppgave nå.
chrtsta skrev:Takk for det!
Jeg prøvde å regne ut integralet, men jeg fikk feil svar.
Det som derimot ser ut til å bli riktig er en mye enklere metode: å bruke Gauss sats (divergensteoremet). Da har vi
[tex]\int_\sigma \vec{F} \cdot \vec{n} \, d\sigma = \int_\tau \nabla \cdot \vec{F} \, d\tau[/tex], hvor [tex]\int_\tau 1 \, d\tau = V_\tau = V_{kule}[/tex],
så jeg får [tex]Q = 4 V_{kule} = \frac{16 \pi a^3}{3}[/tex]
Noen som kan gi en utfyllende forklaring på dette her?
Jeg er ikke helt sikker på hva som skjer fra ledd til ledd..

Og hvordan kan jeg finne Q, når sentrum ikke går gjennom Origo lenger, men gjennom et annet punkt i vektorfeltet?
Kan jeg gjøre dette på samme måte?

Håper på svar her :)

Lagt inn: 26/10-2009 14:51
av Karl_Erik
Det som skjer er at vi bruker Gauss' sats, som sier at vektorfluksen av et felt gjennom en flate er lik integralet av divergensen til feltet over volumet omsluttet av flaten. Her betyr det at integralet vi vil regne ut (vektorfluksen av feltet over en kuleflate) er lik integralet av divergensen til feltet (som er konstant lik 4) over volumet av kulen omsluttet av kuleflata. Vi vet at integralet av en konstant funksjon over et volum er lik volumet ganger konstanten, så dette integralet blir altså lik fire ganger kulas volum.

Er dette MEK1100-obligen, forresten?

EDIT: Ser ut som posten ble redigert mens jeg svarte, så for å svare på det andre spørsmålet ditt: Ja, du kan gjøre det på samme måte. Om du bruker Gauss' sats blir utregningene akkurat det samme, da kulas volum ikke endres av at vi flytter den rundt i feltet. Du kan selvfølgelig også bruke kulekoordinater med sentrum et annet sted enn i origo selv om dette ikke blir spesielt pent.

Lagt inn: 26/10-2009 14:55
av tool-nes
Tror jeg skjønte litt mere nå ja :)

Ja, det var det jeg trodde, at den ikke ville forandre seg noe, og svaret vil bli samme på oppgave a) og b). Og ville blitt veldig mye utregning med kulekoordinater osv..

Dette er vel måten å finne divergensen til feltet?
[tex]\bigtriangledown \cdot \vec{F} = \frac{\delta \vec{F}_x}{\delta_x} + \frac{\delta \vec{F}_y}{\delta_y} + \frac{\delta \vec{F}_z}{\delta_z} = 4[/tex]

Og deretter integralet av dette
[tex]\int_{\sigma} \bigtriangledown \cdot \vec{F} d \sigma= \int_{\tau} \vec{F} d \tau [/tex]
Vi vet at [tex] \int 1 d \tau = V_{kule}[/tex]
Dermed har vi funnet ut at [tex]Q = 4V_{kule} = \frac{16 \pi a^3}{3}[/tex]
:o :o