Ich habe ein Problem bei der Lösung des folgenden magnetostatischen Problems. Über Hilfe und Anleitung würde ich mich sehr freuen.
So wird das Problem beschrieben:
Eine Kugelschale mit RadiusR
, die eine gleichmäßige Oberflächenladung tragenσ
, wird mit einer Winkelgeschwindigkeit in Rotation versetztω
, siehe Abbildung 5.45 . Berechnen Sie das Vektorpotential, indem Sie das entsprechende partielle Differential lösen∇2A =μ0J
.
Ich habe versucht, das Problem zu lösen, indem ich zuerst den Oberflächenstrom unter Verwendung eines sphärischen Koordinatensystems aufgeschrieben habe.
J (R') = σv =σω × ρ =σωR'Sündeθδ(R'−R ) _e^ϕ
Wo
ρ
ist der Abstand von der z-Achse, wenn die z-Achse als Rotationsachse genommen wird und
θ
ist der von der z-Achse gemessene Winkel. Daraus schließen wir das
A
hängt nur von der ab
Aϕ
Komponente, da der Oberflächenstrom nur von der Phi-Komponente abhängt, und somit haben wir
∇2Aϕ=μ0Jϕ
. Und wir müssen nur noch eine "einfache"
Poission-Gleichung lösen . Unser Problem hat eindeutig
Azimutsymmetrie .
Wir fahren dann fort, die Laplace-Gleichung zu lösen
∇2Aϕ( r , θ ) = 0
. Wir lösen diese partielle Differentialgleichung mit Hilfe der Separation der Variablen. Am Ende stellen wir fest, dass die allgemeine Lösung eine Linearkombination separierbarer Lösungen ist.
Aϕ( r , θ ) =∑l = 0∞(AlRl+BlRl + 1)Pl( weilθ )
Für den Innenbereich wo
( r ≤ R ) ,Bl= 0
sonst würde das Potential am Ursprung explodieren.
Aich n( r , θ ) =∑l = 0∞AlRlPl( weilθ )
Im Außenbereich
( r ≥ R )
wir bekommen das
Al= 0
, weil sie im Unendlichen nicht auf Null gehen.
Adu t _( r , θ ) =∑l = 0∞BlRl + 1Pl( weilθ )
Diese müssen wir nun an der Grenze zu Funktionen zusammenfügen
r = R
. Das Vektorpotential sollte am Rand stetig sein. Wir haben also folgende Randbedingungen
Adu t _=Aich n
(∂Adu t _∂R−∂Aich n∂R)∣∣r = R= −μ0Jϕ= −μ0σω R sinθ
Aus der ersten Randbedingung finden wir das
Bl=AlR2 l + 1
Wenn wir es in die zweite Randbedingungsgleichung einsetzen, erhalten wir
∑l = 0∞( 2 Liter + 1 )AlRl - 1Pl( weilθ ) =μ0σω R sinθ
Verwenden Sie den Fourier-Trick, indem Sie beide Seiten mit multiplizieren
Pl'( weilθ ) Sündeθ
und integrieren. Wir können die Koeffizienten bestimmen
Al
da die Legendre-Polynome orthogonale Funktionen sind,
∫10Pl'( x )Pl( x ) d x =∫π0Pl'( weilθ )Pl( weilθ ) Sündeθd θ=22 l + 1δll'
Al=μ0σωR _2Rl - 1∫π0SündeθPl( weilθ ) Sündeθdθ _
Hier bin ich auf eine Mauer gestoßen. Ich habe versucht, das Problem durch Ausdrücken zu lösen
Sündeθ
in Bezug auf Legendre-Polynome und verwenden Sie dann die orthogonale Bedingung, um die Koeffizienten zu bestimmen
Al
. ich weiß, dass
Sünde2θ
kann unter Verwendung von Legendre-Polynomen auf folgende Weise ausgedrückt werden.
Sünde2θ =23(P0( weilθ ) −P2( weilθ ) )
Aber wenn ich nur haben willSündeθ
Ich muss aus allem die Quadratwurzel ziehen, und dann wird es wirklich chaotisch. Wo habe ich bei meiner Lösung einen Fehler gemacht?
![Bild des Problems](https://i.stack.imgur.com/o8JA8.png)
Kyle Kanos