Šíření a reflexe elektromagnetického vlnění


Pokud se vrátíme ke vzorci pro komplexní amplitudu napětí na přenosovém vedení a spočítáme z ní skutečné hodnoty napětí např. pro případ zakončení charakteristickou impedancí, zjistíme, že imageS/01_12/image002.gif, kde imageS/01_12/image004.gif. Jedná se o  jednorozměrnou analogii šíření vln v prostoru. Tato analogie není náhodná, její opodstatnění vyplyne z řešení Maxwellových rovnic pro vlnovody.

Nejprve zopakujeme šíření elektromagnetických vln ve vakuu a v prostředí a zavedeme Poyntingův vektor.

Vyjdeme z rotačních Maxwellových rovnic:

imageS/01_12/image006.gif      ,      imageS/01_12/image007.gif.

První rovnici vynásobíme imageS/01_12/image008.gif druhou imageS/01_12/image010.gif a odečteme je od sebe. Dostaneme:

imageS/01_12/image012.gif.

Pokud bude závislost imageS/01_12/image016.gif na imageS/01_12/image008.gif a imageS/01_12/image020.gif na imageS/01_12/image010.gif lineární, tj. imageS/01_12/image024.gif(imageS/01_12/image025.gifje tenzor permitivity, ne Levi-Civitův tenzor) je možné napsat imageS/01_12/image026.gif.

Dále platí imageS/01_12/Image028.gifa pokud zavedeme Poyntingův vektor imageS/01_12/image030.gif(označení je z německého Strahlungs vektor - vektor záření - viz [3] str. 120 a dále) můžeme napsat rovnici imageS/01_12/image032.gif. Všechny členy mají rozměr J.s-1 . Jedná se o bilanční rovnici energie. Výraz imageS/01_12/image036.gifpředstavuje Jouleovské teplo vytvořené v objemové jednotce vodivého prostředí za jednotku času.

Hustota energie elektromagnetického pole je imageS/01_12/image034.gif.

Smysl Poyntingova vektoru můžeme ozřejmit, přepíšeme-li naši bilanční rovnici pro vakuum (σ = 0) tj. imageS/01_12/image038.gif. Integrujeme přes určitý objem V a použijeme Gaussovu větu:

imageS/01_12/image040.gif (první imageS/01_12/image042.gif je Poyntingův vektor, druhé imageS/01_12/image042.gif je element plochy).

Poyntingův vektor imageS/01_12/image042.gif je vektorem proudové hustoty elektromagnetické energie. Tento důkaz podává pouze informaci, že celkový úbytek energie z oblasti za jednotku času je dán integrálem imageS/01_12/image048.gif. Poyntingův vektor nás informuje o změnách a velikostech toku energie v každém místě prostoru. 

Mějme například válcový vodič, kterým teče proud Iz.

 

imageS/1121.gif  

Pokud označíme R odpor na jednotku délky, bude pak pole Ez dáno výrazem imageS/01_12/image050.gif (Ohmův zákon). Magnetické pole na povrchu vodiče (resp. v libovolné vzdálenosti r od  vodiče) je imageS/01_12/image052.gif. Poyntingův vektor směruje tedy radiálně do vodiče: imageS/01_12/image030.gif   ,   imageS/01_12/image057.gif.

Obr. 1.12.1 Válcový vodič

 

Spočítáme-li energii vnesenou tímto vektorem za 1s do vodiče integrací přes válec o poloměru vodiče a délce 1m, dostaneme výkon

imageS/01_12/image058.gif, který se rovná výkonu disipovanému Jouleovským teplem. Podle této představy baterie vytváří elektrické a magnetické pole. Energie teče tímto polem přes povrch vodiče do vodiče.

Pro stacionární střídavý případ se používá komplexní Poyntingův vektor, který udává střední hodnotu reálného Poyntingova vektoru za jednu periodu:

imageS/01_12/image060.gif              , kde * udává komplexní sdruženost.

V souvislosti s energií elektromagentického pole je vhodné si uvědomit, že elektromagnetické pole nese i hybnost a moment hybnosti, takže elektromagnetické pole může působit i mechanickými účinky - ponderomotivní účinky elektromagnetického pole. Například v tokamacích se v poslední době užívá ohřev na dolnohybridním kmitočtu. To je de facto využití ponderomotivních účinků pole pro urychlení částic podél křivek magnetického pole v tokamaku, čímž zpětně dochází ke zvyšování energie v jednotce objemu plazmatu.

Zatím jsme odvozovali rovnice pro imageS/01_12/image062.gif a φ v nevodivém prostředí a v kvazistacionárním případě ve vodiči, kdy jsme zanedbávali posuvný proud proti proudu vodivostnímu.

Nyní odvodíme rovnice pro imageS/01_12/image008.gif a imageS/01_12/image020.gifpro obecný případ, kdy budeme uvažovat vodivostní a posuvný proud (viz [3] str.135). Předpokládáme lineární matematické vztahy:

imageS/01_12/image066.gif      ,      imageS/01_12/image067.gif     ,     imageS/01_12/image068.gif

 a vyjdeme z rotačních rovnic:

imageS/01_12/image070.gif      ,     imageS/01_12/image072.gif.

Na první rovnici užijeme operátor rot a dosadíme ze druhé rovnice:

imageS/01_12/image074.gif.

V homogenním prostředí je imageS/01_12/image076.gif a tedy pravá strana rovnice je rovna 0.

Obdobně aplikací rotace na druhou rovnici a po úpravě:

imageS/01_12/image078.gif.

Kvazistacionární případ pak můžeme charakterizovat jako „pomalu se měnící pole“. Druhé derivace jsou mnohem menší než první derivace podle času:

imageS/01_12/image080.gif      ,      imageS/01_12/image081.gif.

Nerovnosti lze dosáhnout buď  „nízkým“ kmitočtem nebo velkou vodivostí.

V tomto případě lze “vypustit” členy obsahující druhé derivace z odvozených rovnic. V opačném případě, platí obrácení „silné“ nerovnosti, což je dáno většinou nevodivým prostředím σ → 0, dostáváme vlnové rovnice, které první derivace imageS/01_12/image008.gif a imageS/01_12/image020.gifpodle času neobsahují. Nás bude především zajímat případ harmonických časových závislostí:

imageS/01_12/image085.gif      ,      imageS/01_12/image086.gif,

což vede k Helmholzově rovnici:

imageS/01_12/image088.gif  , kde imageS/01_12/Image090.gif  ,  imageS/01_12/image092.gif

Pro ωt << 1 dostaneme kvazistacionární přiblížení (difuze pole k povrchu vodiče), pro ωt → ∞ (malé  σ) dostaneme řešení ve tvaru rovinné vlny šířící se rychlostí imageS/01_12/Image096.gif, kde imageS/01_12/image098.gif je relativní index lomu prostředí.

Obecně lze řešení vlnové rovnice v nevodivém prostředí psát ve tvaruimageS/01_12/image100.gif (viz [3] str.138 a dále) .imageS/01_12/image102.gif je jednotkový vektor ve směru šíření vlny, imageS/01_12/Image096.gif je fázová rychlost vlny, N je obecně funkcí kmitočtu N = N(ω).

Z analýzy řešení plyne, že každé nestatické řešení vyhovující Maxwellovým rovnicím se skládá pouze z příčných (transverzálních) komponent. Vektory imageS/01_12/image008.gif a imageS/01_12/image020.gif jsou rovnoběžné s rovinou kolmou na imageS/01_12/image102.gif. Stočíme-li souřadnou soustavou tak, že směr šíření bude totožný s osou z, soustava Maxwellových rovnic se nám podstatně zjednoduší. Předpokládejme dále, že derivace podle x a y jsou nulové. Rozebíráme tedy případ, kdy jediné změny imageS/01_12/image008.gif a imageS/01_12/image010.gifnastávají ve směru z. Pak dostaneme z rotačních Maxwellových rovnic

imageS/01_12/image113.gif      a       imageS/01_12/image115.gif

tyto rovnice:

imageS/01_12/image116.gif

imageS/01_12/image117.gif

imageS/01_12/image118.gif
imageS/01_12/image119.gif

imageS/01_12/image120.gif

imageS/01_12/image121.gif

Rovnice s nulou na levé straně potvrzují naše dřívější zjištění, že nestatické řešení nemá z-ovou složku. Řešení pro Ex ,Ey , Hz , Hy se dají napsat jako součet funkcí argumentu imageS/01_12/Image123.gif a imageS/01_12/Image125.gif . Vezmeme-li například pouze složku Ex postupující v kladném směru osy z:imageS/01_12/Image127.gif, pak imageS/01_12/Image129.gif.

Po integraci podle času a zanedbání konstantní komponenty:

 imageS/01_12/Image131.gif, kde imageS/01_12/Image133.gif .

Veličina η má rozměr Ω a nazývá se vlastní impedancí média. Pro vakuum má hodnotu imageS/01_12/Image135.gif .

Obdobnými úpravami bychom dostali

imageS/01_12/image137.gif

imageS/01_12/image138.gif.

Ze skalárního součinu imageS/01_12/image008.gif. imageS/01_12/image010.gif=0 plyne, že vektory imageS/01_12/image008.gif a imageS/01_12/image010.gif jsou na sebe kolmé. Z vektorového součinu imageS/01_12/image008.gifx imageS/01_12/image010.gif, dostaneme vektor v kladném směru osy z (opačně pro imageS/01_12/Image151.gifimageS/01_12/Image153.gif). Spočítáme-li energii elektrického pole imageS/01_12/Image155.gifa magnetického pole imageS/01_12/Image157.gif zjistíme, že U= UH , takže hustota energie v každém bodě a v každém okamžiku je rovnoměrně rozložena mezi elektrické a magnetické pole. Poyntingův vektor imageS/01_12/Image159.gif směřuje v kladném směru osy z. Časová střední hodnota Poyntingova vektoru musí být stejná ve všech rovinách podél vlny. Lišit se však mohou okamžité hodnoty podle toho, zda dochází ke vzrůstu nebo úbytku energie mezi těmito rovinami.

 

 Rovinná vlna je tedy charakterizována těmito vlastnostmi :

Z hlediska použití mají zejména význam vlny, jejichž argument se harmonicky mění. Zkoumejme takovou monochromatickou vlnu. Pokud směr šíření není totožný s jednou z os, je možné zavést proměnnou imageS/01_12/Image177.gif , kterou převedeme vlnovou rovnici na v podstatě jednorozměrný případ (viz [3] str. 139). Stačí proto uvažovat, že se vlna šíří podél osy z.

Řešení vlnové rovnice napíšeme ve tvaru:

imageS/01_12/Image179.gif           (1.12.1)

Případně ještě pomocí komplexních amplitud imageS/01_12/image181.gif , kde imageS/01_12/Image183.gif je vlnové číslo (vlnočet) imageS/01_12/Image185.gif . Na k je možné nahlížet jako na charakteristiku vlny, neboť imageS/01_12/image187.gif udává změnu fáze vlny na 1m šíření. Spíše se však nahlíží na k jako na konstantu prostředí, neboť je možné napsat imageS/01_12/image189.gif , kde N je index lomu prostředí.

Pro složku Hy je možné napsat:     imageS/01_12/image191.gif.

Okamžitou hodnotu udává reálná část vztahu (1.12.1). Postupné opakování se děje po vzdálenosti imageS/01_12/Image193.gif.

Střední hodnota Poyntingova vektoru: 

imageS/01_12/image195.gif.

 Pokud vlna pochází z jednoho zdroje případně z několika zdrojů s konstantním fázovým rozdílem, můžeme pro řešení vlnové rovnice napsat:

 imageS/01_12/image201.gif      ,      imageS/01_12/image203.gif.

Tyto vektory imageS/01_12/image008.gif a imageS/01_12/image020.gif budou obecně komplexní. Kvadráty imageS/01_12/image209.gifimageS/01_12/image210.gifbudou také komplexní. Například s fází -2αimageS/01_12/Image211.gif.

 Zavedeme-li vektor imageS/01_12/Image217.gif fázově posunutý vůči imageS/01_12/Image213.gif o úhel -α tj.imageS/01_12/image215.gif, bude imageS/01_12/image219.gifreálné a bude mít stejnou velikost jako imageS/01_12/Image213.gif. Pak imageS/01_12/image223.gif.

Vektor imageS/01_12/Image217.gif rozložíme na reálnou a imaginární složku, tj. imageS/01_12/Image226.gif , kde imageS/01_12/Image228.gif jsou reálné. Platí, že imageS/01_12/Image230.gif a oba vektory jsou kolmé ke směru šíření.

Pro imageS/01_12/image008.gif pak platí:

imageS/01_12/Image234.gif

kde

imageS/01_12/image236.gif

imageS/01_12/image237.gif.

Platí tedy imageS/01_12/Image238.gif . Vektor imageS/01_12/image008.gifrotuje v rovině kolmé na směr šíření vlny a jeho konec obecně opisuje elipsu. Taková vlna se nazývá elipticky polarizovaná. Je-li imageS/01_12/Image242.gif , je vlna kruhově polarizovaná. Je-li imageS/01_12/Image244.gif nebo imageS/01_12/image246.gif nulové, pak pole vlny je v každém okamžiku a místě paralelní jednomu směru, vlna je lineárně polarizovaná. Vektor imageS/01_12/image010.gifmagnetického pole se chová podobně. Vztahy mezi Ei a Hj lze získat pomocí impedance η.

Pokud má vlna stále konstantní amplitudu, nemůže pomocí ní přenášet informaci. Informaci je možné přenášet, měníme-li amplitudu, frekvenci nebo fázi vlny. Taková vlna pak bude superpozicí rovinných monochromatických vln o různé amplitudě, kmitočtu a fázi.

Předpokládejme opět jednorozměrné přiblížení imageS/01_12/Image248.gif. Jedním z řešení bude imageS/01_12/image250.gif , kde imageS/01_12/Image252.gif .

V obecném prostředí (nikoliv homogenním) může být imageS/01_12/Image254.gif , kde imageS/01_12/Image256.gif nemusí být lineární funkce. Obecné řešení lze pak napsat jako superpozici individuálních monochromatických vln s různou amplitudou imageS/01_12/Image258.gif. Změny fáze neuvažujeme. Lze ukázat, že složení vln s různou fází dá stejný výsledek jako složení vln s různou frekvencí: imageS/01_12/image260.gif. Budeme předpokládat, že vlna ψ je superpozicí vln, jejichž vlnočty leží v úzkém intervalu (viz [3] str.148) imageS/01_12/Image262.gif a že amplitudy se v tomto intervalu nemění, tj. imageS/01_12/Image264.gif, pak imageS/01_12/image266.gif.

Disperzní závislost linearizujeme imageS/01_12/Image268.gif a zavedeme imageS/01_12/Image270.gif .

Argument v integrandu je  imageS/01_12/Image272.gif.

Výsledné imageS/01_12/image274.gif , kde imageS/01_12/Image276.gif je amplituda vlnového balíku imageS/01_12/image278.gif.

Místa stejné amplitudy balíku se pohybují prostorem rychlostí  imageS/01_12/image280.gif, kde V se nazývá grupová rychlost.

Označíme-li imageS/01_12/Image282.gif , můžeme napsat, že amplituda balíku bude nenulová v úzké oblasti, kde imageS/01_12/Image284.gif . Zvyšujeme-li šířku pásma (imageS/01_12/Image286.gif), můžeme balík lépe prostorově lokalizovat. Vztah mezi v a V dostaneme ze vztahu  imageS/01_12/image189.gif, budeme-li derivovat podle ω:

 imageS/01_12/Image290.gif

neboť N je neklesající funkcí kmitočtu imageS/01_12/image292.gif .