Nauka #25 – Elektromagnetyzm: równania Maxwella, fale EM, anteny i linie transmisyjne

Synteza elektryczności i magnetyzmu w jeden zestaw czterech równań pola, dokonana przez Jamesa Clerka Maxwella w 1865 roku, jest jednym z największych osiągnięć teoretycznych w fizyce. Równania przewidują prędkość światła wyłącznie na podstawie stałych elektromagnetycznych, zjednoczyły optykę z elektrycznością i położyły podwaliny pod komunikację bezprzewodową, systemy energetyczne oraz elektrodynamikę kwantową. Każda antena, transformator, światłowód i kuchenka mikrofalowa działają na zasadach zapisanych w równaniach Maxwella.

1. Równania Maxwella — pełny obraz teorii

Cztery równania Maxwella wiążą ze sobą pole elektryczne E, pole magnetyczne B, gęstość ładunku ρ oraz gęstość prądu J. W postaci różniczkowej obowiązują w każdym punkcie przestrzeni i czasu; w postaci całkowej wiążą strumień i krążenie pola po powierzchniach i pętlach. Razem stanowią kompletny klasyczny opis wszystkich zjawisk elektromagnetycznych.

Postać różniczkowa (układ SI, w próżni lub ośrodku liniowym):

  ∇·E = ρ/ε_0                            (prawo Gaussa: pole E rozbiega się od ładunków)
  ∇·B = 0                                   (brak monopoli magnetycznych: B bez dywergencji)
  ∇×E = −∂B/∂t                        (Faraday: zmienne B indukuje krążące E)
  ∇×B = μ_0(J + ε_0 ∂E/∂t)         (Ampère + Maxwell: prądy i zmienne E generują B)

Postać całkowa:
  &oiint; E·dA = Q_enc/ε_0                   (całkowity strumień przez zamkniętą powierzchnię = ładunek wewnątrz)
  &oiint; B·dA = 0                               (brak izolowanych biegunów magnetycznych)
  ∮ E·dl = −dΦ_B/dt                     (SEM wokół pętli = −szybkość zmiany strumienia magnetycznego)
  ∮ B·dl = μ_0(I_enc + ε_0 dΦ_E/dt)   (krążenie magnetyczne = prąd wewnątrz + prąd przesunięcia)

Dodatek Maxwella: prąd przesunięcia ε_0 ∂E/∂t
  Bez niego ∇·J = −∂ρ/∂t (zachowanie ładunku) jest niezgodne z pierwotnym prawem Ampère'a.
  Sens fizyczny: zmienne pole E w szczelinie kondensatora działa jak prąd — mimo że żaden ładunek jej nie przekracza.

Stałe:
  ε_0 = 8,854 × 10^−12 F/m  (przenikalność elektryczna próżni)
  μ_0 = 4π × 10^−7 H/m      (przenikalność magnetyczna próżni)
  c = 1/√(μ_0ε_0) = 2,998 × 10^8 m/s  (prędkość światła, przewidziana przez Maxwella w 1865 r.)

2. Fale elektromagnetyczne i wektor Poyntinga

W przestrzeni wolnej od źródeł równania Maxwella można połączyć w równanie falowe: każda składowa E i B spełnia (∇² − (1/c²)∂²/∂t²) f = 0. Rozwiązaniami są poprzeczne fale elektromagnetyczne (TEM): E i B oscylują prostopadle do siebie i do kierunku propagacji, w fazie, przy czym |E|/|B| = c.

Monochromatyczna fala płaska rozchodząca się w kierunku +z:
  E(z,t) = E_0 cos(kz − ωt) x̂
  B(z,t) = (E_0/c) cos(kz − ωt) ŷ
  k = ω/c = 2π/λ   (liczba falowa)

Gęstość energii i pędu:
  u_E = (1/2)ε_0 E²    (gęstość energii elektrycznej)
  u_B = (1/2μ_0) B²         (gęstość energii magnetycznej)
  u_E = u_B dla fali płaskiej (równy podział)
  Razem: u = ε_0 E² = B²/μ_0

Wektor Poyntinga (gęstość strumienia energii):
  S = (1/μ_0) E × B   [W/m²]
  |S| = u·c = E²/(μ_0 c)   dla fali płaskiej

Natężenie uśrednione w czasie (napromienienie):
  ⟨S⟩ = E_0² / (2μ_0 c)  [W/m²]
  Na powierzchni Ziemi: napromienienie słoneczne ~ 1361 W/m²

Ciśnienie promieniowania:
  P_rad = ⟨S⟩/c   (dla powierzchni doskonale pochłaniającej)
  P_rad = 2⟨S⟩/c  (dla powierzchni doskonale odbijającej)
  Ciśnienie promieniowania słonecznego w odległości 1 AU: ~ 9 × 10^−6 Pa (małe, ale istotne dla statków kosmicznych)

Widmo EM (zakresy częstotliwości):
  Fale radiowe:     f < 300 MHz   (λ > 1 m)
  Mikrofale:        300 MHz – 300 GHz  (λ 1 mm – 1 m)
  Podczerwień:      300 GHz – 430 THz  (λ 700 nm – 1 mm)
  Światło widzialne: 430 – 770 THz  (λ 390 – 700 nm)
  Ultrafiolet:      770 THz – 30 PHz  (λ 10 – 390 nm)
  Promieniowanie X: 30 PHz – 30 EHz  (λ 0,01 – 10 nm)
  Promieniowanie gamma: f > 30 EHz    (λ < 0,01 nm)

3. Indukcja Faradaya i prawo Lenza

Prawo indukcji Faradaya mówi, że zmienny strumień magnetyczny przechodzący przez obwód indukuje siłę elektromotoryczną (SEM). Prawo Lenza określa jej znak: prąd indukowany zawsze przeciwdziała zmianie, która go wywołała (to zasada zachowania energii w innej postaci). Razem indukcja Faradaya i prawo Lenza są zasadą działania generatorów, transformatorów, silników indukcyjnych, ładowania bezprzewodowego i wykrywaczy metali.

Prawo Faradaya:
  &EMcycl; = −dΦ_B / dt   gdzie Φ_B = &iint; B·dA (strumień magnetyczny przez powierzchnię)
  Dla cewki o N zwojach:  &EMcycl; = −N dΦ_B/dt

Prawo Lenza: znak minus zapewnia, że indukowana SEM wywołuje prąd przeciwdziałający zmianie strumienia.

SEM ruchu:
  Przewodnik o długości L poruszający się z prędkością v prostopadle do B:
  &EMcycl; = BLv   (siła Lorentza działająca na swobodne ładunki w przewodniku)
  Generator Faradaya: obracająca się pętla (pole A, N zwojów) z prędkością kątową ω:
  &EMcycl;(t) = NBAω sin(ωt)   (sinusoidalne napięcie wyjściowe AC)

Indukcyjność własna L:
  Φ_B = LI   (strumień przez cewkę proporcjonalny do jej własnego prądu)
  V_L = −L dI/dt   (napięcie na cewce przeciwdziała zmianie prądu)
  Energia zmagazynowana: U_L = (1/2)LI²

Indukcyjność własna solenoidu:
  L = μ_0 n² V   (n = liczba zwojów na metr, V = objętość)
  Dla solenoidu o 1000 zwojach, promieniu 10 cm i długości 50 cm: L ≈ 3,9 mH

Transformator (idealny):
  V_1/V_2 = N_1/N_2   (stosunek napięć = stosunkowi zwojów)
  I_1/I_2 = N_2/N_1   (stosunek prądów odwrotny)
  P_1 = P_2            (zachowanie energii, przy sprawności 100%)
  Rzeczywiste transformatory: straty rezystancyjne (I²R), straty na prądy wirowe, straty histerezowe

4. Promieniowanie anten — od prądu do pola dalekiego

Antena zamienia zmienne w czasie prądy elektryczne na rozchodzące się fale elektromagnetyczne (tryb nadawania) albo odwrotnie (tryb odbioru). Charakterystyka promieniowania — kątowy rozkład mocy promieniowanej — wynika z potencjałów opóźnionych i przyspieszenia ładunków. Nawet najprostsza antena, dipol Hertza, zawiera w sobie istotę fizyki każdego promieniowania.

Dipol Hertza (nieskończenie krótki element prądu I dℓ):
  Dalekie pole elektryczne:
    E_θ = (I dℓ / 4π) · (μ_0 ω²/c) · (sinθ/r) · cos(ωt − kr)
  Dalekie pole magnetyczne:
    H_φ = E_θ / η_0   (η_0 = √(μ_0/ε_0) = 377 Ω impedancja falowa)
  Oba pola maleją jak 1/r (pole dalekie); pole bliskie ma człony 1/r² i 1/r³.

Charakterystyka promieniowania:
  Moc na jednostkę kąta bryłowego: dP/dΩ ∝ sin²θ  (kształt torusa, zero na osi dipola)

Całkowita moc promieniowana (wzór Larmora dla przyspieszanego ładunku):
  P = q²a² / (6πε_0 c³)   (a = przyspieszenie, q = ładunek)
  Dla elementu prądu: P = (I² dℓ² ω&sup4; μ_0) / (12π c³)  ∝ f&sup4;
  Wyższe częstotliwości promieniują znacznie wydajniej (dlatego AM potrzebuje ogromnego nadajnika, a FM znacznie mniejszego).

Kluczowe parametry anten:
  Kierunkowość D: stosunek natężenia głównej wiązki do średniej izotropowej
    Dipol Hertza: D = 1,5 (1,76 dBi)
    Dipol półfalowy: D = 1,64 (2,15 dBi)
    10-elementowa antena Yagi: D ~ 10–14 dBi
  Zysk G = D · η (η = sprawność promieniowania, < 1 z powodu strat omowych)
  Szerokość pasma: zakres f dla SWR < 2 (często 3–5% dla prostych dipoli)
  Impedancja: dipol półfalowy Z_in ≈ 73 + j42,5 Ω przy rezonansie

Wzór transmisyjny Friisa (bilans łącza):
  P_received = P_transmitted · G_T · G_R · (λ/4πd)²
  Tłumienie w przestrzeni swobodnej: FSPL = (4πd/λ)²
  Przy 2,4 GHz na dystansie 100 m: FSPL ~ 80 dB (krytyczne dla bilansu Wi-Fi)

5. Efekt naskórkowości — dlaczego prądy wysokiej częstotliwości unikają wnętrza

Gdy prąd zmienny płynie w przewodniku, indukowane prądy wirowe przeciwdziałają przenikaniu strumienia, wypychając prąd ku powierzchni. Gęstość prądu maleje wykładniczo wraz z głębokością, z charakterystyczną głębokością zwaną głębokością wnikania δ. Przy wysokich częstotliwościach δ może wynosić ułamek milimetra, co oznacza, że przewód miedziany przy 10 GHz przenosi prąd tylko w cienkiej powłoce — jego rezystancja objętościowa staje się nieistotna. Efekt ten decyduje o stratach na ścieżkach PCB, doborze rdzenia transformatora oraz konstrukcji ekranów RF.

Równanie falowe w przewodniku (zespolona przenikalność elektryczna):
  ∇²E = jωμσE   (σ = przewodność)
  Rozwiązanie: E(x) = E_0 exp(−x/δ) exp(−jx/δ)
  (zanik i przesunięcie fazy wraz z głębokością x)

Głębokość wnikania:
  δ = √(2/ωμσ) = 1/√(π f μ σ)
  f ↑ → δ ↓ (prąd skupia się w cieńszej warstwie)

Wartości dla miedzi (σ = 5,96 × 10^7 S/m):
  50 Hz:     δ = 9,3 mm   (pełne wykorzystanie przekroju przy częstotliwości sieciowej dla typowych kabli)
  1 kHz:     δ = 2,1 mm
  1 MHz:     δ = 66 µm
  100 MHz:   δ = 6,6 µm
  10 GHz:    δ = 0,66 µm  (obwody mikrofalowe)

Rezystancja powierzchniowa (efektywna rezystancja na kwadrat):
  R_s = 1/(σδ) = √(π f μ / σ)   [Ω/kwadrat]
  Miedziana ścieżka przy 1 GHz: R_s = 0,026 Ω/kwadrat — istotne dla ścieżek o długości milimetrów.

W konsekwencji:
  R_AC / R_DC = a/(2δ) dla promienia przewodu a ≫ δ
  Drut litzendraht: wiele cienkich żył o średnicy < 2δ skręconych razem
  → zmniejsza efektywną rezystancję przy częstotliwościach akustycznych i RF.

6. Linie transmisyjne — fale prowadzone i współczynnik fali stojącej

Linia transmisyjna — kabel koncentryczny, mikropasek, falowód — prowadzi fale elektromagnetyczne od źródła do obciążenia. Gdy impedancja obciążenia nie jest dopasowana do impedancji charakterystycznej linii Z_0, część sygnału zostaje odbita, tworząc falę stojącą. Współczynnik fali stojącej (SWR) mierzy stopień niedopasowania; maksymalny transfer mocy wymaga SWR = 1 (idealne dopasowanie). Teoria linii transmisyjnych jest niezbędna w projektowaniu RF, rozmieszczaniu ścieżek PCB powyżej częstotliwości gigahercowych oraz w pracach nad integralnością sygnału.

Równania telegrafistów (model obwodu rozłożonego):
  ∂V/∂x = −(R + jωL) I
  ∂I/∂x = −(G + jωC) V
  Stała propagacji: γ = α + jβ = √((R+jωL)(G+jωC))
  Impedancja charakterystyczna: Z_0 = √((R+jωL)/(G+jωC))

Linia transmisyjna bezstratna (R=0, G=0):
  Z_0 = √(L/C),  β = ω√(LC),  v_p = 1/√(LC)
  Kabel koncentryczny: Z_0 = (60/√ε_r) ln(b/a) Ω   (b/a = stosunek promienia zewnętrznego/wewnętrznego)
  Mikropasek:    Z_0 zależy od stosunku w/h i przenikalności podłoża ε_r (obliczane numerycznie)

Współczynnik odbicia przy obciążeniu Z_L:
  Γ = (Z_L − Z_0) / (Z_L + Z_0)   (zespolony, |Γ| ≤ 1)
  Zwarcie: Z_L = 0 → Γ = −1  (całkowite odbicie, z inwersją fazy)
  Rozwarcie:  Z_L = ∞ → Γ = +1  (całkowite odbicie, bez inwersji)
  Dopasowanie:       Z_L = Z_0 → Γ = 0  (brak odbicia)

Współczynnik fali stojącej napięcia:
  SWR = (1 + |Γ|) / (1 − |Γ|) ∈ [1, ∞)
  SWR = 1: idealne dopasowanie. SWR = ∞: całkowite odbicie (zwarcie lub rozwarcie).
  Moc odbita: |Γ|² mocy padającej.

Transformator ćwierćfalowy (dopasowanie impedancji):
  Odcinek λ/4 o Z_1 = √(Z_0 Z_L) przekształca Z_L w Z_0.
  Szerokość pasma ~ 20% wokół częstotliwości projektowej.

Wykres Smitha:
  Zespolona płaszczyzna Γ naniesiona tak, by okręgi stałej rezystancji dawały się łatwo wyznaczyć.
  Używany do dopasowania impedancji, analizy stabilności i projektowania filtrów.
  Narzędzia programowe: ADS, CST, SPICE / AWR Microwave Office.