Ceramics-Silikáty 33, (4) 311 - 323 (1989)

Je uvedena účinná výpočetní procedura pro třírozměrné stacionární rychlostní a teplotní pole ve sklovině. Algoritmus je založen na svázaném diferenčním řešení rovnic hybnosti, kontinuity a energie v primitivních proměnných pomocí koncepce multigrid. Modelové parciální diferenciální rovnice jsou diskretizovány v třikrát přeložené síti pomocí hybridního diferenčního schématu, které využívá metodu dílčích oblastí pro zajištění numerické stability výpočtu. Pro interační proceduru je použita symetrická svázaná Gauss-Seidelova technika, která vykazuje velmi rychlou konvergenci. Testovací výpočet proudění skloviny v modelové peci se svislými elektrodami byl proveden pomocí dvou úrovní sítí o 10 x 10 x 12 a 18 x 18 x 22 uzlech. Byla zjištěna dobrá stabilita výpočtu a značné, řádové zvýšení konvergenční rychlosti oproti dříve užívaným segregačním a jednosíťovým metodám.