Program FlexPDE w zastosowaniu do symulacji pól temperatury

COORDINATES
Cartesian3
SELECT
errlim = 1e-4
VARIABLES
Tp { temperatura }
DEFINITIONS
k { przewodność cieplna }
ro { gestość }
cp {cieplo wlasciwe }
r = 0.5 { promien dziury }
q = -k*grad(Tp) { strumien ciepla }
a1 = 0.3 { polowa grubosci zbiornika }
b1 = 0.4 { polowa szerokosci zbiornika }
c1 = 0.6 { polowa wysookosci zbiornika }
a2 = 0.25 { polowa grubosci wypelnienia }
b2 = 0.35 { polowa szerokosci zbiornika }
c2 = 0.55 { polowa wysokosci zbiornika }
T0 = 15 { temperatura poczatkowa }
Tpow = 25 { temperatura powietrza }
Tgleba = 2 { temperatura gleby }
hT = 10 { wspolczynnik wnikania ciepla }
V = VOL_INTEGRAL(1)
{ calkowita objetosc zbiornika wraz z
wypelnieniem }
Vz = VOL_INTEGRAL(1,'zbiornik') { objetosc zbiornika bez wypelnienia }
Vw = VOL_INTEGRAL(1,'wypelnienie') { objetosc wypelnienia }
Tsred = VOL_INTEGRAL(Tp)/V { srednia temperatura zbiornika z wypelnieniem
}
Tzsred = VOL_INTEGRAL(Tp,'zbiornik')/Vz { srednia temperatura zbiornika bez
wypelnienia }
Twsred = VOL_INTEGRAL(Tp,'wypelnienie')/Vw { srednia temperatura wypelnienia
}
Minuta = 60
Godzina = 60*Minuta