Lecture 15 :: Vektoranalys

15

Lecture :: Vektoranalys
sammanfattning

Här introduceras detaljer kring nablaoperatorn. Divergensen \(\nabla\bullet \mathbf{F}\) och rotationen \(rot\ \mathbf{F}=curl\ \mathbf{F}=\nabla\times\mathbf{F}\) för ett vektorfält \(\mathbf{F}\) ges tolkningar. Vi introducerar divergenssatsen (Gauss sats) och Greens sats och ger exempel på hur de kan användas.





Läsanvisningar::
Avsnitt i kursboken
Detta material svarar mot Adams Kapitel 16.1-16.4
Nablaräkning
Nablaoperatorn
Nablaoperatorn definieras av \[ \nabla=(\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z}) \] Med \(\nabla\) så definierade vi tidigare gradienten till en funktion. Vi har tidigare sett att gradienten innehåller viktig information om funtionen. (pekar ut maximal tillväxtriktning och är ortogonal mot nivåkurvorna.) Gradienten bildar ett vektorfält och är modellen för konservativa vektorfält.
Nu introducerar vi två nya operationer på vektorfält: \[ \text{Divergensen :: }\quad div\ \mathbf{F}=\nabla\bullet\mathbf{F}= \frac{\partial F_1}{\partial x}+\frac{\partial F_2}{\partial y}+\frac{\partial F_3}{\partial z} \] som kan tolkas om information om hur fältet expanderar eller komprimeras.
\[ \begin{split} \text{Rotationen :: }\quad\nabla\times\mathbf{F}&= curl\ \mathbf{F}= rot\ \mathbf{F} =\det\left[\begin{array}{ccc}i & j & k \\\frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\F_1 & F_2 & F_3\end{array}\right]=\\ &=\left(\frac{\partial F_3}{\partial y}-\frac{\partial F_2}{\partial z},\frac{\partial F_1}{\partial z}-\frac{\partial F_3}{\partial x},\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\right) \end{split} \] Rotationen mäter hur stor rotationseffekten är i en punkt i ett flöde.

Övning 1::

Visa nablaidentiteten \(\nabla(\bullet\nabla\times\mathbf{F})=0\)

Lösning ::

Vi har att \[ \begin{split} &\nabla(\bullet\nabla\times\mathbf{F})=\\ &=\frac{\partial}{\partial x}\left(\frac{\partial F_3}{\partial y}-\frac{\partial F_2}{\partial z}\right)+\frac{\partial}{\partial y}\left(\frac{\partial F_1}{\partial z}-\frac{\partial F_3}{\partial x}\right)+ \frac{\partial}{\partial x}\left(\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\right) \end{split} \] Utnyttja att de blandade andraderivatorna är lika så följer det att uttrycket ovan blir noll.
Video :: Vektoranalys intro
Exempel på divergenssatsen

Övning 2 ::

Gör Adams 15.6.1 mha divergenssatsen.

Lösning ::

Vi har att \(\mathbf{F}=(x,z,0)\) och då får vi att \(\nabla\bullet\mathbf{F}=1\) vilket ger \[ \begin{split} \iint_S \mathbf{F}\bullet N dS&=\iiint_V \underbrace{\nabla\bullet\mathbf{F}}_{=1} dV=\\ &=\text{ volymen av tetraedern }= \frac{6\cdot 3\cdot 2}{6}=6 \end{split} se även videon på nästa sida! \]
Video :: Adams 15.6.1 igen.
Curl-räkning

Övning 3 ::

Låt \(\mathbf{F}=(-y,x,0)\). Beräkna integralen \[ \frac{1}{\pi\epsilon^2}\int_{C_\epsilon} \mathbf{F}\bullet dr, \] där \(C_\epsilon)\) är cirkeln centrerad i origo med radie epsilon. Vad blir integralen gränsvärde då \(\epsilon\to 0^+\)?
Vad blir \(N\bullet\nabla\times\mathbf{F}\), där \(N\) är enhetsnormal till cirkelskivan som ger positiv orientering av kuran

Lösning ::

Parametrisera cirkeln med \(r(t)=(\epsilon\cos t,\epsilon\sin t)\). Man får att \(dr=(-\epsilon\sin t,\epsilon\cos t)dt\). Med denna parametrisering får vi att \[ \begin{split} \frac{1}{\pi\epsilon^2}\int_{C_\epsilon} \mathbf{F}\bullet dr&= \frac{1}{\pi\epsilon^2}\int_0^{2\pi}(-\epsilon\sin t,\epsilon\cos t)\bullet(-\epsilon\sin t,\epsilon\cos t)dt=\\ &=\frac{1}{\pi\epsilon^2}\int_0^{2\pi}(\epsilon^2\sin^2 t+\epsilon^2\cos^2 t)dt=\\ &=\frac{1}{\pi\epsilon^2}\int_0^{2\pi} r^2dt=\frac{1}{\pi\epsilon^2}\cdot \epsilon^2\cdot2\cdot \pi=2, \end{split} \] Denna integral är alltså oberoende av \(\epsilon\) och gränsvärdet blir således \(2\). Vi får också att \[ \nabla\times\mathbf{F}=(0,0,2) \] Om kurvan har positiv orientering så får vi att cirkelskivans normal blir \(N=(0,0,1)\). Vi får därför att \[ N\bullet\nabla\times\mathbf{F}=2 \] Det är ingen slump att dessa två saker blir lika...
Video :: Tolkning av curl.
En prematur verifiering av Greens sats

Övning 4 ::

Fortsätt med situationen beskriven i övning 3. Beräkna nu
  1. \[\int_{C_\epsilon} \mathbf{F}\bullet dr\]
  2. \[\iint_{B_\epsilon} \nabla\times F\bullet N dS,\] där \(B_\epsilon\) är cirkelskivan \(x^2+y^2\leq\epsilon\)

Lösning ::

  1. Från övning 3 har vi att integralen blir \(2\pi\epsilon^2\)
  2. Övning 3 ger oss också att \((\nabla\times F)\bullet N=2\).
    I polära koordinater har vi att \(dS=rdrd\theta\) vilket ger att integralen blir \[ \iint_{B_\epsilon} \nabla\times F\bullet N dS=2\int_0^{2\pi}\int_0^\epsilon r dr d\theta= r^2\int_0^{2\pi}d\theta=2\pi \epsilon^2. \]
Att de två integralerna är lika är en konsekvens av Greens sats som vi pratar om på nästa video.
Man ska notera att ett tvådimensionellt vektorfält \(F=(F_1(x,y),F_2(x,y))\) kan betraktas som ett tredimensionellt vektorfält \(F=(F_1,F_2,0)\) och då får vi att \[ \nabla\times F=(0,0,\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}) \] vilket ger att \[ \nabla\times\mathbf{F}\bullet N=\nabla\times\mathbf{F}\bullet (0,0,1)=\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y} \] där vi utnyttjat att ett område i planet har vi normalriktningen \(N=(0,0,1)\).
Video :: Greens sats