Numerické modelování elektromagnetických polí se s rozvojem výpočetní techniky a neustále rostoucí výkonností počítačů stalo spolu s optimalizačními technikami nepostradatelnou složkou návrhu konstrukcí nových elektrotechnických a elektronickýchzařízení i zařízení z ostatních oblastí technické praxe. Numerické modelování je také bezesporu nedílnou součástí komplexních analýz chování časoprostorových polí, které jsou důležité pro posouzení nových požadavků na kvalitu zařízení jako je elektromagnetická kompatibilita. Složité problémy řešené v současné technické praxi nelze zvládnout ve většině případů jinými prostředky než pomocí vhodných numerických metod za použití výkonných počítačů.
Je vidět, matice koeficientů získaná aplikací MKD řídká; jednom řádku je
nejvýše nenulových koeficientů při libovolném počtu rovnic. Jednou nevýhod
MKD je, způsobuje vážné komplikace rozhraních složitou geometrií rozhraních
se skokovou změnou permitivity.
. Proto lze spolehlivě řešit na
současných 500 000 rovnic výkonných počítačích několik milionů rovnic
tohoto typu.
Výhodou těchto metod malý nárok paměť počítače, neboť uschovávají jen pole
s nenulovými koeficienty dále vysoká stabilita řešení. instruktivní učebnici [4]. Soustavy uvedenými vlastnostmi jsou dobře
řešitelné některou moderních iteračních metod, např. intenzitu mezi uzly 17
a směru osy +y. něj pak lze určit kapacitu elektrod.
Metodu lze snadno rozšířit úlohy viz [1]. Metodu lze použít pouze pro řešení úloh
s uzavřenou hranicí (známá Dirichletova nebo Neumannova podmínka). Hledané řešení získáme vždy uzlech
sítě, ostatních bodech prostoru interpolací.
Po nalezení uzlových potenciálů lze stanovit další veličiny, např.
Při výpočtu statických polí dnes MKD pro uvedené komplikace téměř vytlačena metodou
konečných prvků, kterou seznámíme podrobně následujícím textu. rovnici pro uzel lze napsat tvaru
nebo01004 4798 =−−−− φφφφ
1004 4798 =−−− φφφφ
Všechny sestavené rovnice uspořádáme přehledného maticového zápisu
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
−−
−−−
−−−
−−
−−
−−−
−−−
−−−
−
100
100
100
0
0
0
0
0
0
420010000
141001000
014100100
002400020
100042000
010014100
001001410
000100141
000000024
9
8
7
6
5
4
3
2
1
φ
φ
φ
φ
φ
φ
φ
φ
φ
Soustavu řešíme některou numerických metod. Podrobný postup popsán např.01
y E
l
φ φ
φ −
−
−
− =
∆
Z okrajové podmínky elektrodě lze stanovit hustotu náboje její
numerickou integrací celkový náboj. Pro tento úsek platí
4 17
4 17
4 17
29. metodou konjugovaných gradientů.9
, odkud 2990 V/m
0.FEKT Vysokého učení technického Brně
převedeme tuto hodnotu druhou stranu. Např. Koeficient diagonále je
větší nebo roven součtu absolutních hodnot koeficientů mimo diagonálu, soustava tedy
diagonálně dominantní positivně definitní. Hledaný vektor potenciálů vyřešení
φ [9,1 18,1 26,2 29,9 31,0 37,3 56,5 62,6 64,1]Τ