Program FlexPDE w zastosowaniu do symulacji pól temperatury

Godzina = 60*Minuta
Liczba_godzin = 3
CzasSym = Godzina*Liczba_godzin { czas symulacji w godzinach }
DeltaCzas = 60 { interwal czasowy w sekundach }
INITIAL VALUES
Tp = T0
EQUATIONS
ro*cp*dt(Tp)+div(q) = 0
EXTRUSION
SURFACE 'denko zbiornika' z = -c1
LAYER 'dolna warstwa zbiornika'
SURFACE 'denko wypelnienia' z = - c2
LAYER 'wypelnienie'
SURFACE 'wieczko wypelnienia' z = c2
LAYER 'gorna warstwa zbiornika'
SURFACE 'wieczko zbiornika' z = c1
BOUNDARIES
SURFACE 'denko zbiornika' VALUE(Tp) = Tgleba
SURFACE 'wieczko zbiornika' NATURAL(Tp) = hT*(Tp - Tpow)
REGION 1 'zbiornik' { plastik }
k = 0.22 { przewodność cieplna scianek }
ro = 1200 { gestość scianek }
cp = 1250 {cieplo wlasciwe scianek }
START(-a1,-b1)
NATURAL(Tp) = hT*(Tp - Tpow) LINE TO (a1,-b1)
LINE TO (a1,b1)
LINE TO (-a1,b1)
LINE TO CLOSE
LIMITED REGION 2 'wypelnienie' { pustka }
LAYER 2 VOID { proznia w zbiorniku }
START(-a2,-b2)
LINE TO (a2,-b2)
LINE TO (a2,b2)
LINE TO (-a2,b2)
LINE TO CLOSE