Od pewnego czasu kombinuję z programem do liczenia rozkładu pola elektrostatycznego w próżni metodą sukcesywnej relaksacji (Liebmanna). Opierałem się na książkach do optyki elektronowej (są tam poruszane tematy obliczania rozkładów pól elektrycznych i magnetycznych w soczewkach elektronowych). Pisałem w języku ANSI C (kompilator gcc). Problem elektrostatyczny jest reprezentowany przez dwie tablice. Demonstrator DEMO_XY wykorzystuje tablice dwuwymiarowe, Pierwsza typu unsigned char przechowuje warunek ograniczający danego węzła (0 to próżnia, 1 - elektroda typu 1, 2 - elektroda typu 2 itd.). Druga tablica typu double przechowuje wartości potencjału pola elektrostatycznego w danych węzłach. Rozwiązywanie problemu polega na trawersowaniu siatki i wykonywania procedur relaksacji dla poszczególnych węzłów (jeżeli dany węzeł reprezentuje próżnię). Jeżeli kolejne wykonanie procedury relaksacji dla węzłów typu próżnia prawie nic nie zmienia, wówczas obliczania są zakończone, a dane wyjściowe wypisywane do plików teksowych (dat) i graficznych (PPM P3 i BMP).
Program wgrałem na swoją stronę internetową. Jest tam więcej szczegółów.
Podaję linki:
https://marcinkulbaka.prv.pl/Liebmann/index.htmlhttps://marcinkulbaka.prv.pl/Liebmann/galleries/index.htmlpozdrawiam
Marcin Kulbaka