Obsah

Řídící parametry řešiče pro nestacionární úlohu

Řídící parametry řešiče XT2T/XT3T se zapisují do souboru name.iB:

IP KREST 0 KOUT INT3 NSAX NSTEPX 0 RP TIMS ERAL EDIF TOL DTRUN PIVOT PENAL

Význam parametrů

Pro každý zatěžovací stav, představující časový interval řešení, se musí zadat v rámci AS dávky:

AS číslo zatěžovacího stavu /$\dots$ /R TIMX STEP TSC /$\dots$

přičemž

Vhodné je též aktivovat AV sadu:

AV číslo sady T 6 N KAPPR KAUTO KPRED V 4*0

přiřazením v AS dávce

AS číslo zatěžovacího stavu /$\dots$ /A číslo sady /$\dots$

přičemž

V prvním zatěžovacím stavu je nutné zadat počáteční teplotního pole pomocí GV sady.

Každému zatěžovacímu stavu přísluší tolik záznamů v souboru name.STR a name.TEM, kolik bylo provedeno integračních kroků.

Integrační schema

Pro integraci diferenciálních rovnic v čase je použita zobecněná jednokroková metoda s parametrem $\mathtt{TSC}.$ Hodnota $\mathtt{TSC}=0$ odpovídá explicitní dopředné Eulerově metodě, která způsobí linearizaci rovnic, s řešením odpovídajícím posloupnosti po částech lineárních změn (dějů). Tato metoda sice nevyžaduje žádné iterace, může však dojít k nestabilitám, tj. k rozkmitání teplot v uzlových bodech. Nutná je pak pečlivá kontrola výsledků, takže se metodu nedoporučuje používat.

Hodnota $\mathtt{TSC}>0$ vede vždy k implicitní metodě s podmínkou nepodmíněné stability $\mathtt{TSC}\ge0{,}5.$ Případ $\mathtt{TSC}=1$ představuje plně implicitní zpětnou Eulerovu metodu. Implicitní integrace vyžaduje iterace, protože v každém časovém kroku je nutné řešit nelineární soustavu rovnic. Komentář týkající se parametrů $\mathtt{ERAL},$ $\mathtt{EDIF},$ $\mathtt{TOL}$ a $\mathtt{KAPPR}$ platí beze změny.

Hodnota $\mathtt{NSAX}$ určuje maximální počet iterací v každém časovém kroku. Doporučuje se přiřadit AV sadu s $\mathtt{KAPPR}=1$ v každém zatěžovacím stavu. Z důvodů stability je vhodné v každém časovém kroku iterovat až do konvergence, tj. zvolit hodnotu $\mathtt{NSAX}$ raději vyšší, aby počet iterací v každém kroku byl určen kritérii $\mathtt{EDIF}$ a $\mathtt{TOL}$ a nikoliv dosažením hodnoty $\mathtt{NSAX}.$

Časový krok

Řešení sleduje v čase posloupnost zatěžovacích stavů. Každému zatěžovacímu stavu musí být přiřazena hodnota $\mathtt{TIMX},$ což je čas, do kterého momentální předpis platí. Jakmile je dosaženo $t=\mathtt{TIMX},$ nahradí se dosavadní předpis novým zatěžovacím stavem.

Hodnota $\mathtt{STEP}$ je uživatelem předepsaná délka integračního kroku. Tato délka může být v každém zatěžovacím stavu (intervalu) různá. Je třeba mít ovšem na paměti, že změna integračního kroku vyžaduje nový přímý chod. Pokud je integrace řízena automaticky, je $\mathtt{STEP}$ výchozí délka kroku.

Automatické řízení je v každém zatěžovacím stavu (intervalu) inicializováno klíčem $\mathtt{KAUTO}=1$ zadávaným AV sadou (pro $\mathtt{TSC}=0{,}5$ nebo $\mathtt{TSC}=0$ se nedoporučuje používat). Při automatickém řízení se integrační krok prodlužuje, jestliže je $$\mathtt{PLTE}<0{,}25\,\mathtt{TOL},$$ kde $\mathtt{PLTE}$ je odhad maxima lokální chyby řešení [$^\circ\text{C}$]. Integrační krok se zkracuje pro $$\mathtt{TOL}<\mathtt{PLTE}.$$ Hodnota $\mathtt{TOL}$ se zadává na RP řádku. Vypočtená délka integračního kroku se zaokrouhluje na celistvý násobek $\mathtt{DTRUN}.$ Tento čas tak představuje elementární kvantum (s ohledem na výstupy).

V zatěžovacích stavech, v nichž je $\mathtt{KAUTO}=1,$ se hledá celkový počet časových kroků od začátku výpočtu. Překročí-li hodnotu $\mathtt{NSTEPX},$ výpočet se ukončí.

Restart úlohy

Pokud úloha skončila úspěšně, je možné pokračovat v řešení zadáním dalších zatěžovacích stavů s klíčem $\mathtt{KREST}=3.$ Na IP řádku je pak třeba zadat pořadové číslo kroku $\mathtt{INT3},$ od kterého se má začít (včetně). Pokud má být úloha restartována od skončení posledního řešení, zadá se $$\mathtt{INT3}=1+\sum\mathtt{NSTEP}_i,$$ kde $\mathtt{NSTEP}_i$ je počet integračních kroků v $i$-tém zatěžovacím stavu (intervalu) a sumace probíhá přes všechny zatěžovací stavy.

Predikce termofyzikálních vlastností

Před zahájením výpočtu v časovém intervalu je možné odhadnout hodnoty veličin extrapolací z minulého kroku, což může urychlit výpočet a zlepšit konvergenci. Predikce se spouští klíčem $\mathtt{KPRED}=1,$ který se zadává AV sadou a platí v rámci odpovídajícího zatěžovacího stavu.

1)
Parametr $\mathtt{TIMS}$ je jen formální veličinou. Posun časové osy slouží ke snadnějšímu popisu časových závislostí.