Kimatická simulace s jednou hladinou na Icosahedrální síti

Místo pro debatu o věcech, které nelze zařadit do předchozích kategorií, přesto s počasím souvisejících.

Moderátor: Moderátoři

Odpovědět
TommyAst
Příspěvky: 1556
Registrován: 05 bře 2013, 22:44
Meteostanice: V 850 hPa pod -40 °C a zvoleni Babise v CR nechcem

Kimatická simulace s jednou hladinou na Icosahedrální síti

Příspěvek od TommyAst »

První dobře vypadacící výsledek - klimatické modelování 2D s využitím šestiúhelníkové sítě, stupeň rozlišení 5 - 40 962 bodů. Výpočet trval poměrně dlouho - přibližně 60 dní. Nejvíce výpočetního času sebrala interpolace (z 40962 bodů na 2000x1000, jeden obrázek trval okolo 23 s, celkem 32777 map pro 6 veličin - 196 662 map).

Pro exoplanetu Tau Cet e (ještě nebyla existence zcela prokázána) zvoleny parametry: Tlak 2 bar, tepelná kapacita 1007 J/K/kg, rovnovážná teplota 289 K (Země 255 K), Sklon rotační osy 37° (Země 23,45°), délka siderického roku 14525568 s/ 168.12 našich dní, perioda rotace 158 400 s, délka dne 160 127 s, eccentricita 0,05, rozdíl perihelium-rovnodenost 15°.

Teplota
http://uloz.to/xHpgWTPK/qx-000000040962 ... d-temp-avi

Forcing teploty
http://uloz.to/xGAhCtrQ/qx-000000040962 ... rcingt-avi

Změna teploty
http://uloz.to/xk2jqKEA/qx-000000040962 ... rtialt-avi

Insolace
http://uloz.to/xdrVAtK1/qx-000000040962 ... ated-s-avi

Radiační bilance
http://uloz.to/xMEpUhxC/qx-000000040962 ... ated-q-avi

K - Sinus zenitového úhlu
http://uloz.to/xA3LH2Fm/qx-000000040962 ... ated-k-avi

Jiný běh modelu - s atmosférou jen 0.2 bar a stejnými ostatními hodnotami.
Mapy jsou bez interpolace, neboť interpolace trvá doluho. Video vytvořeno z 182 345 map.
Teplota http://uloz.to/xvRQGhuS/map-classic-mola-temp-test-avi
Ro rychlejší download https://mega.co.nz/#!VxJgFQqT!LyKmgxe6j ... KGVN6sEnik

Parametry programu :
SigmaST=5.670373E-8 ! St.-B. konstanta
Mplanet=2.565E+25 Hmotnost planety
a=11582000.0 ! Rovníkový poloměr
asp=21000.0 ! Rozdíl poloměru rovník a póly
g=(6.6732E-11)*Mplanet/(a*a) ! povrchová gravitace, hmotnost 2.565E25 kg/4.25 Me, poloměr 11 582 km, hustota 3941.4
Efektivní tíhové zrychlení:
!gm=((6.6732E-11)*Mplanet/((a-(asp*cos(souradnice(ID,2))))**2)-((2.0*Pi*(1.0/DelkaDne))**2)*((a-(asp*cos(souradnice(ID,2))))*cos(souradnice(ID,2)))-(U0*U0+V0*V0+W0*W0)/((a-(asp*cos(souradnice(ID,2))))**2)-(KorF*U0*(-sin(souradnice(ID,3)))))
Sklon=0.64577182323790194346176558434079 ! Sklon rotační osy, 37 deg
e=0.05 ! ekcentricita
DelkaRoku=14525568.0 ! Sidercká, 168.12 Earth days
DelkaDne=158400.0 ! Siderická délka dne - perioda rotace, 44 hodin
DelkaDneSyn=160127.34 ! Sinodická délka dne, delší o 1727.34 s
EquinoxShift=0.26179939 ! Úhlový rozdíl v dráze mezi periheliem a rovnodenností, perihelium 15°před rovnodenností
HlAxis=82579200000.0 ! 8.25792E+10 m or 0.552 AU Hlavní polosa
KonstT=44.682256021 ! Konstanta teplota vs. výška H/T
Tsl=5344.0 ! Povrchová teplota hvězdy
aSL=551928000.0 ! POloměr hězdy, 0.793 Slunečního
Time=0.0 ! Čas start
dT=4.0 ! Časový krok
TimeW=14525568.0 ! První zápis hodnot (po 1 roce)
TimeWG=14525568.0 ! Čas pro první zápis globálních průměrů
TimeWSU=14525568.0 ! Čas pro start výpoču dlouhodobých průměrů
TimeWO=1329.5 ! Čas mezi zápisem hodnot pro vytvoření mapy
TimeWOG=41.55 !Čas mezi zápisi globálních průměrů
TimeMax=58102272.0 ! Délka běhu modelu - 4 místní roky
CasSU=0 ! Proměnná času pro dlouhodobé průměry, 0 na začátku, integer
PHIpsu=159686.38 ! průměrný geopotenciál na začátku, tlak 2 bar, hustota 1.252455, g 12.76 ms-2
Hsu=PHIpsu/g ! Průměrná výška na začaátku (tloušťka atmosféry)
EPSm=0.55 ! ckoeficien emisivitz vyzařování (může být 0-1)
Albedo=0.303 ! Albedo - průměr. Pro albedo 303 rovnovážná teplota 280.08 K (rovnovážná subsellární teplota 396.09 K) - hustota atmosféry 1.3954 kg/m3.
AlbedoAtm=0.434 ! Albedo atmosféry (pro nekonečnou vzduchovou hmotu, pro nulovou vzuchovou hmotu je to 0.303)
OpAtm=0.785 ! Propustnost atmosféry (vzduchová hmota 1, v zenitu)
Rho=1.252455 ! Hustota atmosféry homogenní, odpovídající 1 bar a 280.08 K, plynová konstanta R 288.05 pro suchý vzduch (R = 288.844 for mixing ratio 10 g/kg)
KpMerna=1007.0 ! Měrná tepelná kapacita cp - J/(kg.K).
V 850 hPa pod -40 °C a Babise prezidentem v CR nechceme
TommyAst
Příspěvky: 1556
Registrován: 05 bře 2013, 22:44
Meteostanice: V 850 hPa pod -40 °C a zvoleni Babise v CR nechcem

Re: Kimatická simulace s jednou hladinou na Icosahedrální sí

Příspěvek od TommyAst »

Nový kodek nalezen a soubory jsou nyní podstatně menší. Jedná se o h265, kodování nejpomalejší, quality 17. Rozlišení 2000x1000 30 fps. Velikost je nesrovnatelně nižší, kodování videa je výpočetně extrémně náročné. Staré velké soubory budou smazány.

Teplota (K)
http://uloz.to/xfCFumK1/qx-000000040962 ... d-temp-mp4

Forcing teploty dT/dt (K/s)
http://uloz.to/xc11MtP6/qx-000000040962 ... rcingt-mp4

Insolace S (W/m2)
http://uloz.to/xnd3crkh/qx-000000040962 ... ated-s-mp4

Radiační bilance Q (W/m2)
http://uloz.to/xGpp1gny/qx-000000040962 ... ated-q-mp4

Sinus výšky Slunce nad obzorem K
http://uloz.to/xwvvYUo4/qx-000000040962 ... ated-k-mp4
V 850 hPa pod -40 °C a Babise prezidentem v CR nechceme
Odpovědět

Zpět na „Témata související s počasím“