|
7. Calcule globale utilizand ecuatia de transport diferentiala
În acest capitol se realizeaza trecerea de la formularile in limbaj general matematic, la scrierea si testarea algoritmilor numerici de rezolvare a ecuatiei de transport. Calculul numeric este supus unor constrangeri care sunt subliniate in paragraful 6.1. Sunt schitate in acest paragraf principiile programarii in FORTRAN si MATLAB, presupunand ca cititorul are cunostinte de baza in informatica. În paragraful 6.2, pornind de la metoda proiectorilor propunem construirea primului program de transport in geometrie placa. De asemenea prezentam doua probleme test pentru verificarea algoritmului de calcul. Pe baza metodei proiectorilor pentru geometrie placa prezentata in paragraful 4.5 facem generalizarea la geometrii 2D si 3D. În paragrafele 7.3 si 7.4 prezentam pas cu pas etapele scrierii unui algoritm de transport in geometrie carteziana 3D. Aceleasi etape si aceleasi probleme le intalnim si in cazul gemetriei x-y. Algoritmii de calcul rezultati au o forma foarte simpla daca utilizam limbajul automatelor celulare. Scrierea unui algoritm este un proces foarte dificil daca tinem seama de problemele care apar in practica. Foarte multe dintre aceste probleme conduc la nedeterminari numerice. Pentru eliminarea acestora este bine ca fiecare detaliu al unui algoritm sa fie inteles pe deplin. În cazul ecuatiei de transport problemele provin din reprezentarea functiei exponentiale in calculator. Prezentam algoritmul Lentz de reprezentare a acestei functii in paragraful 7.5. In paragraful 7.6 este descrisa structura generala a programului LEGENTR care este utilizat intr-o serie de lucrari de laborator in Capitolul 10.
Acest capitol se incheie cu o discutie comparativa intre schemele cu diferente discrete utilizate in agoritmii de difuzie si metoda proiectorilor.
|