Übung 1 am 12.03.2010

Hiermit ist die Edyn-Beispieljagt eröffnet:
a01SS2010f.pdf (44.4 KB)
Hier sind meine Lösungsvorschläge :smiley: :

1.1. Tensorschreibweise

a) C_{i}=\epsilon_{ijk}\partial_{j}B_{k}=rot\vec B
b) l^2=a_j a_k \delta_{jk}=\vec a \vec a = \vec a^2
c) \rho=\partial_m \delta_{km} E_i \delta_{ik} = div \vec E
d) C_m=\partial_s \epsilon_{tsm} B_t = -rot \vec B
e) \rho=\delta_{km}\partial_i \delta_{mk} E_i=tr(\mathbb{1}) \cdot div \vec E
f) \epsilon_{abc} \partial_b \partial_c \Phi = \vec \nabla \times \vec \nabla \Phi = 0
g) \epsilon_{abc} \partial_l \epsilon_{klm} \partial_b a_m \delta_{kc} = \vec \nabla (\vec \nabla \vec a)-\Delta \vec a
h) q=\partial_a \partial_b \epsilon_{ace} \epsilon_{bcd} \Phi \delta_{ed} = (tr(\mathbb{1}) - 1) \Delta \Phi

1.2. Vektorfelder

a)

rot \vec B = \begin{pmatrix} 0\ 0\ 1 \end{pmatrix}

div \vec B = 0

rot \vec E = \vec 0

div \vec E = -1

b)

rot \vec E = 0

\vec E = -\vec \nabla \Phi \to \Phi = \frac{x^2+y^2}{4}-x-y

c)

div \vec B = 0


\vec B = \vec \nabla \times \vec A

A_1 = \begin{pmatrix} 0\ 0\ x+y-\frac{x^2+y^2}{4} \end{pmatrix}

A_2 = \begin{pmatrix} z(\frac{x}{2}-1)\ -z(1-\frac{y}{2})\ 0 \end{pmatrix}

rot( \vec A_1 - \vec A_2 ) = 0

ich krieg das Gleiche heraus, bis auf :wink:

ad 1.1:
statt tr(1) kann man auch einfach 3 schreiben, wenn man 3 Dimensionen annimmt.

ad 1.2:
a) div E = +1 (statt -1)
b) die Lösung stimmt für \vec{E}=\nabla \phi, andernfalls ist die Lösung -\phi

Du hast recht:
div \vec E = \partial_x (\frac{x}{2}-1)+\partial_y (\frac{y}{2}-1)=\frac{1}{2}+\frac{1}{2}=1

Stimmt ebenfalls. Ich hab zwar -\Delta \Phi hingeschrieben, aber mit \Delta \Phi gerechnet.

Die Lösung lautet: \Phi = x + y - \frac{x^2+y^2}{4}

Lösungsvorschläge für 1.3. a - c

a)

q \int_{-\infty}^\infty\int_{-\infty}^\infty\int_{-\infty}^\infty dx dy dz \delta(x) \delta(y) \delta(z-a) = q

b)

\Phi(x,y,z)=\frac{q k_1}{sqrt{x^2+y^2+(z-a)^2}}

c)

\vec E = \frac{q k_1}{\sqrt{x^2+y^2+(z-a)^2}^3}\begin{pmatrix}
x\
y\
z-a
\end{pmatrix}

d)

Komischerweise komme ich div \vec E = 0 wenn ich die Divergenz direkt berechne. Finde ich halt doch sehr eigenartig.

mit 1.3 a-c bin ich einverstanden :slight_smile:

ad d) stimmt schon, wenn du einfach die Divergenz des Ausdrucks für E bildest, ist diese 0 (das passt, weil der Bereich ja quellenfrei ist)
auf ganz R ohne den Punkt (0,0,a), wo die Punktladung sitzt (dort wird der Nenner unendlich).

Für dem Bereich in der Nähe der Punktladung könnte man sich an die Vorgehensweise im Skriptum halten,
wo die erste Maxwell-Gl. hergeleitet wird. Mit dem Satz von Gauß komme ich ja zu dem Ausdruck
div \vec{E} = k_1 4 \pi \rho (\vec{r})
Ich hab leider keine Ahnung, was ich mit den 4pi machen soll.

Zusätzliche DiskussionsteilnehmerInnen sind herzlich willkommen :smiley:

Ich vermute die 4pi sind das übliche Problem.
nämlich: Es heißt immer wir lernen nur noch das SI system, aber die theoretiker verwenden weiterhin das cgs system in den vorlesungen.

Geb ich dir recht.

Setzt man k_1=\frac{1}{4 \pi \epsilon_0}, erhält man den Ausdruck für das SI System: div \vec E=\frac{\rho}{\epsilon_0} oder div \vec D=\rho

Lösungsvorschläge für 1.3

e) \phi (r,\Theta ,\varphi ) = \frac{k_1 q}{\sqrt{r^2-2ar\cos \Theta +a^2}}
f) \vec{E} = \frac{k_1 q}{(r^2-2ar\cos \Theta +a^2)^{3/2}}\left [ (r-a\cos \Theta) \vec{e_r} + a \sin \theta \vec{e_\varphi } \right ]
Mit Hilfe der Koordinatentransformation, die im Plenum gezeigt wurde, kann man - unter großem Rechenaufwand - den Ausdruck für E in kartesischen Koordinaten auf Kugelkoordinaten transformieren. Für den Radialanteil stimmt es, weiter habe ich noch nicht gerechnet.

Weiß jemand einen einfacheren Weg, das Ergebnis zu überprüfen?

Hat jemand eine Ahnung, wie die allgemeine Lösung bei 1.2 c.) ausschauen könnt, wenn man die Eigenschaft mit der Rotation kennt. Ich habe keinen Schimmer, wie man da auf das Vektorpotenzial A rückschließen kann.

Meine Lösung sieht etwas anders aus:

\Phi(r,\theta,\phi) =q k_1 \frac{1}{r}

Auch mithilfe der Delta-Funktion in Kugelkoordinaten kommt man zu diesem Ergebnis.

Siehe Methoden-Skriptum Seite 217: \delta^3(\vec x -\vec x^')=\frac{1}{4 \pi r^2} \delta(r)

\vec E = \frac{q k_1}{r^2} \vec e_r wobei \vec e_r = \begin{pmatrix}
sin \theta cos \phi\
sin \theta sin \phi\
cos \theta
\end{pmatrix}

Meine Ergebnisse beruhen darauf, dass ich

x = r sin \theta \cos \phi

y = r sin \theta sin \phi

z=r cos \theta + a

gewählt habe.
Ist meines Erachtens auch sinnvoll, da ich hiermit die Kugelsymmetrie ausnutze.

Setzt man in \vec E die Kartesischen Koordinaten wieder ein, so zeigt sich das Ergebnis von 1.3 c.

Ich denke, dass hat etwas mit den Maxwell’schen Gleichungen zu tun. In den Präsentationsunterlagen zur Vorlesung findet man homogene und inhomogene Maxwell Gleichungen. Wenn man sich das ansieht merkt man, dass wir uns in den Übungen mit den homogenen Gleichungen beschäftigt haben. Für eine allgemeine Lösung benötigt man zum homogenen Teil auch eine Partikulärlösung. Ich glaube mit den inhomogenen Maxwell Gleichungen findet man einen partikulären Anteil.

Du hast - sinnvollerweise - den Ursprung der Kugelkoordinaten gleich der Position der Punktladung gesetzt.
Widerspricht das nicht der Angabe („die üblichen Kugelkoordinaten um den Ursprung“)?
Ich dachte, sie wollen es uns schwerer machen, indem der Ursprung des Kugelkoordinatensystems in (0,0,0) liegt,
aber rechentechnisch ist das natürlich sehr ungeschickt. Vielleicht haben sie aber auch das gemeint, was du gemacht hast.
:confused:

ist bei bsp. 1.1 wirklich nur gefragt, dass man das anders hinschreibt? also keine umrechnung oder sonstiges?

lg

Ich denke sie wollen, dass wir (falls das möglich ist) die Rechnungen vereinfachen und dann in die Vektorschreibweise zurückführen. (vgl. 1ten Beitrag)

So um jetzt alle Zweifel zu beseitigen war ich beim Herrn Ipp und er meinte, dass es vom Ursprung aus zu rechnen ist. Also nicht von der Kugelmitte.

Ich habe fast das gleiche Ergebnis nur ist bei mir \vec e_\varphi → \vec e_\vartheta

Danke für’s Nachfragen, das ist besser als spekulieren :wink:

Stimmt, da habe ich mich beim Eingeben geirrt.

Hey ich habe eine Möglichkeit gefunden wie man zeigt, dass das Ergebnis des E-Feldes in Kugelkoordinaten das selbe ist wie in Kartesischen :bulb: :

Man kann die Einheitsvektoren \vec e_r und \vec e_\vartheta ja als Linearkombination von \vec e_1, \vec e_2 und \vec e_3 darstellen:

\vec e_r = sin \vartheta cos\varphi \begin{pmatrix}
1\
0\
0
\end{pmatrix} + sin\vartheta sin \varphi \begin{pmatrix}
0\
1\
0
\end{pmatrix} - cos \vartheta \begin{pmatrix}
0\
0\
1
\end{pmatrix}

\vec e_\vartheta = cos \vartheta cos\varphi \begin{pmatrix}
1\
0\
0
\end{pmatrix} + cos\vartheta sin \varphi \begin{pmatrix}
0\
1\
0
\end{pmatrix} - sin \vartheta \begin{pmatrix}
0\
0\
1
\end{pmatrix}

Setzt man das für \vec e_r und \vec e_\vartheta ein kommt man dann auf die Form:
\begin{pmatrix}
r sin \vartheta cos \varphi\
r sin \vartheta sin \varphi\
r cos \vartheta - a
\end{pmatrix}

was ja das selbe ist wie

\begin{pmatrix}
x\
y\
z - a
\end{pmatrix}

Nun muss man nur noch den Wurzelausdruck in Kartesische Koordinaten umwandeln. Fertig :smiley:.