Paralelní iterační řešič TERMO
Stáhnout paralelní řešič TERMO 1a (TAR.GZ), součást konečněprvkového systému GEM. Řešič TERMO je volný software distribuovaný podle ustanovení Obecné veřejné licence GNU, vydávané Free Software Foundation.
Popis
Kód TERMO se hodí k řešení diskrétního tvaru nestacionárních tepelných rovnic v 3D. Počátečně-okrajová úloha je diskretizována konečnými prvky v prostoru a konečnými diferencemi v čase. Odpovídající lineární systém algebraických rovnic v každém časovém kroku se řeší metodou sdružených gradientů s předpodmíněním založeným na aditivních Schwarzových metodách s překrytím.
Specifickým rysem kódu TERMO je využití pravidelných strukturovaných sítí, které zjednodušeně řečeno jsou sítěmi vzniklými jako výsledek deformace nějaké pravidelné pravoúhlé sítě (vzorová síť, šablona). Diskretizace je dána rozkladem oblasti na osmiuzlové cihly a následně do čtyřstěnných konečných prvků. Proto odpovídající matice mají všechny nenulové prvky uspořádáné ve 27 uzlové pravidelné šabloně. S předpokladem symetrie matice můžeme navíc uchovávat pouze horní trojúhelníkovou část matice. To lze realizovat řádek po řádku využitím 14 prvkové šablony pro ukládání nenulových prvků matice.
Předpokládáme, že materiál je anizotropní a heterogenní a není závislý na teplotních změnách. Tepelná radiace není zahrnuta. Tepelný zdroj se uvažuje exponenciálně klesající v čase (speciálním případem je konstantní zdroj). Čas se měří v letech. Podrobnější informace o programu lze nalézt v dokumentaci (PDF en).
Celý program je napsaný v jazyce Fortran 95 a pro paralelizaci se využívá přístupu OpenMP. Kód byl testován na symetrickém multiprocesoru Natan.