Skin Effekt Belag und Impedanz: Update



In meinem letzten Post [1] habe ich mich gefragt, warum meine Implementierung von Skin Effekt Belag (Drähte mit einem endlichen Widerstand) nicht mit den Resultaten von NEC-2 übereinstimmten. Es stellt sich heraus dass ich falsch lag in der Annahme dass Skin-Effekt Belag nur eine rein ohmsche Komponente hat, durch den Skin Effekt ergibt sich in Wirklichkeit eine relativ große reaktive Komponente.

Hinweis: Dieser Artikel enthält ziemlich viel Mathematik, vielleicht möchten Sie ja ans Ende scrollen und checken wie die neue Implementierung in pymininec sich verglichen mit meinem letzten Artikel verhält wo nur ein rein ohmscher Widerstand angenommen wurde.

Der englische Wikipedia Artikel zum skin effect enthält die korrekte Formel:

Zint=kρ2πrJ0(kr)J1(kr)

wobei

k=jωμρ

und r ist der Draht-Radius, Jv sind die Bessel-Funktionen der ersten Gattung mit der Ordnung v. Zint ist die Impedanz pro Längeneinheit des Drahtes.

NEC-2 verwendet eine andere Formel [2] (S.75):

Zi=jΔiaiωμ2πσ[ber(q)+jbei(q)ber(q)+jbei(q)]

mit

q=(ωμσ)1/2ai

und ai ist der Draht-Radius, Δi die Länge des Drahtsegments, σ die Elektrische Leitfähigkeit des Drahtes und ber und bei die Kelvin-Funktionen. Im folgenden werde ich wieder r für den Draht-Radius, und nicht ai verwenden. Nachdem ber die Ableitung von ber ist, können wir ber als ber0 (Ordnung 0) und ber als ber1 (Ordnung 1) schreiben und analog für bei.

Im folgenden bin ich etwas ungenau in der Verwendung der magnetischen Permeabilität im Freiraum μ0 vs. der Permeabilität eines Drahtes μ, für nicht ferromagnetische Drähte sind diese (fast) identisch: μ=μrμ0 mit μr1 für nicht ferromagnetische Materialien.

Der NEC-2 code verwendet die Fortran Funktion ZINT mit zwei Parametern, der Leitfähigkeit multipliziert mit der Wellenlänge λ (Parameter SIGL) und der Radius geteilt durch die Wellenlänge λ (Parameter ROLAM). Neben der Approximation der Kelvin-Funktionen (der Term mit dem Bruch aus Kelvin-Funktionen ist BR1), berechnet ZINT:

X = SQRT (TPCMU * SIGL) * ROLAM
ZINT = FJ * SQRT (CMOTP / SIGL) * BR1 / ROLAM

mit TPCMU = 2368.705, CMOTP=60.0, und FJ ist j. Das berechnete X wird der Berechnung der Kelvin-Funktionen übergeben mit BR1 als Resultat des Bruch-Terms. Das berechnete ZINT enthält nicht die Segment-Länge, ist also äquivalent mit der Wikipedia Definition in pro Längeneinheit. Mit den Reverse-Engineerten magischen Zahlen

CMOTP μ0c2π

TPCMU 2πcμ0

(c ist die Lichtgeschwindigkeit) bekommen wir für X aus dem Code (was wir jetzt q nennen um mit der obigen Formel aus [2] konsistent zu sein):

q=rλ2368.705σλ=r2368.705σλr2πcμ0σλ=r2πcμ0σfc=r2πμ0σω2π=rωμ0σ

Was mit der obigen Definition von q in [2] übereinstimmt.

Für Zint der Fortran-Implementierung bekommen wir:

Zint=jrμ0c2πσλ[ber0(q)+jbei0(q)ber1(q)+jbei1(q)]=jrμ0cf2πσc[ber0(q)+jbei0(q)ber1(q)+jbei1(q)]=jrμ0f2πσ[ber0(q)+jbei0(q)ber1(q)+jbei1(q)]=jrμ0ω2π2πσ[ber0(q)+jbei0(q)ber1(q)+jbei1(q)]=j2πrμ0ωσ[ber0(q)+jbei0(q)ber1(q)+jbei1(q)]

Diese Formel stimmt NICHT mit der Formel für Zi aus [2] weiter oben überein. Es sieht so aus als ob die Formel aus [2] (S. 75) die Frequenz f statt der Kreisfrequenz ω unter der Wurzel haben sollte. Die aus der NEC-2 Implementierung abgeleitete Formel ist aber konsistent mit der obigen Formel aus Wikipedia.

Der Term vor dem Bruch mit den Kelvin-Funktionen unterscheidet sich um einen Faktor von der Wikipedia Formel (Ich lasse Δi weg und habe so eine Formel pro Längeneinheit wie in Wikipedia). Ausserdem gilt ρ=1/σ:

Aus dem Code abgeleitet:

j2πrωμσ

Aus Wikipedia:

kρ2πr=jωμρρ2πr=12πrjωμρ=12πrjωμσ

Das ist ein Faktor von

jj=1+j2

weil j=1j2

Für die Kelvin-Funktionen gilt die folgende Äquivalenz zu den Bessel Funktionen [3] (10.61, Kelvin Functions):

bervx+jbeivx=Jv(xe3πj/4)=evπjJv(xeπj/4)

Mit den Spezialfällen der Ordnung 0 und 1:

ber0x+jbei0x=e0J0(xeπj/4)=J0(1j2x)ber1x+jbei1x=ejπJ1(xeπj/4)=J1(1j2x)

Die Wikipedia Formel übergibt kr an die Bessel-Funktionen:

kr=rjωμρ=rjωμσ

während NEC-2 folgendes übergibt:

q=rωμσ

Diese unterscheiden sich um den Faktor j. Das bedeutet wir übergeben einen Faktor um 2 zu groß an die Bessel-Funktionen. Daraus ergibt sich wegen der Symmetrie der Bessel Funktionen und er Tatsache dass die Funktionen im Nenner die Ableitung der Funktionen im Zähler sind, dass das Ergebnis um einen Faktor von 2 zu groß ist was den Faktor im ersten Term ausgleicht.

Nach der Implementierung der neuen Formel in pymininec verhält sich der reaktive Teil der Impedanz ähnlich wie bei NEC-2 und der ohmsche Widerstand ist noch näher an den Ergebnissen von NEC als vorher.

Auch wenn man den Imaginärteil (die Reaktanz) plottet sind die Resultate jetzt wie erwartet. Die Werte von MMANA sind immer noch überraschend, eine mögliche Ursache ist dass ich MMANA unter WINE in Linux verwendet habe, die ohmschen Werte der Impedanz die von MMANA berechnet wurden schauen allerdings korrekt aus.