360 likes | 473 Views
KÖRNYEZETI RENDSZEREK MODELLEZÉSE. Dr. Koncsos László egy. docens. Matematikai modellezés. m 1. m 2. F. F. r. Bevezetés – modell fogalma. 1. Példa – Newton bolygómozgási modellje egyik első modern modell – úttörő eredmény két test probléma kiterjedés nélküli testek Nap + 1 bolygó
E N D
KÖRNYEZETI RENDSZEREK MODELLEZÉSE Dr. Koncsos László egy. docens Matematikai modellezés
m1 m2 F F r Bevezetés – modell fogalma 1. Példa – Newton bolygómozgási modellje • egyik első modern modell – úttörő eredmény • két test probléma • kiterjedés nélküli testek • Nap + 1 bolygó • Kepler törvényeinek „elegáns” leírása • több-test probléma – jóval bonyolultabb feladat
Reagáló anyagok Hűtővíz Hűtővíz Termék és melléktermékek Bevezetés – modell fogalma 2. Példa – kémai reakció kevert reaktorban Törvények, hipotézisek • anyag- és energiamegmaradás trv • Fourier hővezetési trv. • irreverzibilis reakció • rendszerparaméterek állandóak • tökéletes keverés – homogén részrendszerek … n. a köpeny reagálása azonnali Több modell is készíthető • egyszerűsítések • kezelhetőség
Bevezetés – modell fogalma 3. Példa – QUAL 2K vízmin. modell
Bevezetés – modell fogalma 3. Példa – QUAL 2K vízmin. modell
Bevezetés – modell fogalma 3. Példa – QUAL 2K vízmin. modell
Bevezetés – modell fogalma Modellalkotás • szabadságfoka gyakran igen nagy • gátló tényezők, lehetőségek figyelembe vétele • megfelelő stragtégia megválasztása Modellek alkalmazásával célunk • fizikai (biológiai,…) folyamatok előrejelzése, rendszerváltozók jövőbeli alakulásának meghatározása • további kísérletek, megfigyelések módosítása • fogalmaink fejlesztése megértése • tervezési célok
Bevezetés – modell fogalma Modell – nehéz általánosan definiálni „Egy matematikai modell matematikai egyenletek (struktúrák) tetszőleges teljes és konzisztens halmaza, amelyet arra terveztek, hogy más tulajdonságok összességét, azok prototípusát írja le. A prototípus lehet fizikai, biológiai, társadalmi, pszichológiai vagy konceptuális (vázlatos) tulajdonság, vagy esetleg éppen egy másik matematikai modell.” [Aris, 1978]
Bevezetés – modell fogalma Modell – nehéz általánosan definiálni „…Tehát a modellek jellemzően a valóság egyszerűsített megfelelői, amelyek alkalmasak a vizsgálatra. … Ezért úgy definiálhatjuk a matematikai modellt, mint egy idealizált szabályrendszer, ami a fizikai rendszer külső behatásra adott válaszát adja meg.” [Chapra, 1997]
Bevezetés – modell fogalma Modell – nehéz általánosan definiálni „A modellek a modern tudomány alapvető eszközei közé tartoznak.” [Frigg, Hartmann 2006] x időleges, célszerű, kényelmes megközelítés örök igazság kifejeződése x túl egyszerű / túl bonyolult hasznos / haszontalan igaz / hamis
Modellek elemei Változók (valós rendszer jellemzői) Függvénykapcsolatok (kölcsönhatások) Paraméterek (ráták, együtthatók, állandók)
Modellek besorolása …származtatás alapján: • Determinisztikus • fizika (más tudományág) alaptörvényeire, igazolt összefüggésere alapul • ált. differenciál- és integrálegyenletek rendszereként írható le Empirikus • megfigyelések alapján felállított, de nem általános érvényű (~empirikus) összefüggésekre alapul • független megfigyelések segítségével • Sztochasztikus • Valószínűség-elméleti törvényszerűségekre alapul • Hibrid • determinisztikus + empirikus elemek + sztochasztikus • biológiai, meteorológiai, … folyamatok leírása során hasznos
Dinamikus Statikus Input Input xout,1 xin,1 xout,1 xout,1 xin,2 xout,1 xin,2 xin,3 xin,3 xin,1 Output Output time time time time time time time time time time Modellek besorolása …időbeliség alapján:
Modellek besorolása … paraméterek alapján: • Időben • állandó paraméterű • változó paraméterű • Térben • „halmozott” (lumped), állandó paraméterű • osztott paraméterű
Modellalkotás folyamata • Identifikáció • alaptörvények, feltevések számbavétele • elégséges részrendszer kiválasztása • matematikai realizáció (papír-ceruza / szoftver) • Kalibráció • modellállandók beállítása • számítások és megfigyelések összevetésén alapul • Validáció • kalibrált modell igazolása • független megfigyelések segítségével • Érzékenységvizsgálat
t=0 h d=2R x(t) t talaj rés Identifikáció – determinisztikus Szennyzőanyag szivárgási probléma Kérdés: x(0)=h x(t)=? Feltételezések • R=áll. • rfolyadék=áll. • talaj homogén, befogadóképessége állandó • Rések rendszere állandó Felírható modell: kifolyási sebesség ~ folyadék szint Szivárgási modell Térfogat t időpontban
http://www.nanowerk.com/ p=1/2 p=1/2 Identifikáció – sztochasztikus Galton deszka Milyen lesz a golyók eloszlása a gyűjtőcellákban? • azonos golyók • n sor • n-1 ütközés • n+1 gyűjtőcella n x=0,1,2,…,n –tartály indexe A Moivre-Laplace tétel szerint: (a normális eloszlással közelíthető) n+1
t=0 h d=2R x(t) t talaj rés Kalibráció – szivárgási probléma Levezetett szivárgási modell • input (ftlen) t • output (függő) x(t) • paraméterek k,h, R Modellparaméter • matematikai modell közvetlenül nem mérhető paramétere • fizikai tartalommal, gyakran dimenzióval rendelkezik • meghatározása közvetett módon zajlik kalibráció
x [m] h t [s] 0 T Kalibráció – szivárgási probléma Kalibrációs feladat • x(t)n-szeri megfigyelése a {0;T} idő-intervallumban x=(xmért,1,xmért,2,…….,xmért,n) – mért idősor • x(t) számítása a {0;T} idő-intervallumban, különböző k értékekre x(k)=(xszámított,1,xszámított,2,…….,xszámított,n) – számított idősor • célfüggvény k-tól függő minimuma? • Modellparaméterek beállítása úgy, hogy a számított eredmények a legjobban közelítsék a valóságban megfigyelt értékeket • manuális • algoritmikus Mért értékek k=0.01 k=kopt k=100 Kiértékelés: r(k1) > r(k2) > … > r(kn)kopt=kn
Kalibráció – szivárgási probléma Távolság - hiba • absztrakt matematikai fogalom • jelentősége – célfüggvény felírása • általános összefüggések: • módosított alakú összefüggések: diszkrét függvényre folytonos függvényre folytonos függvényre diszkretizált súlyozott hibafüggvény
x [m] h t [s] 0 T Validáció – szivárgási probléma Validációs feladat • igazolni a modell és a kalibrált paraméter-vektor megfelelő voltát • a kalibrált modell optimalizációs időszakon kívüli egyeztetése a valósággal Mért értékek k=kopt
Módszertan – dekompozíció és aggregáció 2000 1980 PROBLÉMA & RENDSZER MEGOLDÁS MEGOLDÁS ? MEGOLDÁS AGGREGÁCIÓ AGGREGÁLT MODELL DEKOMPOZÍCIÓ DEKOMPOZÍCIÓ DEKOMPONÁLT RENDSZER
Módszertan – dekompozíció és aggregációPélda 1. • C - KONCENTRÁCIÓ VEKTOR • HIDRODINAMIKAI EGYENLETEK • KEZDETI- ÉS PEREMFELTÉTELEK
Módszertan – dekompozíció és aggregációPélda 2. • C - KONCENTRÁCIÓ VEKTOR • HIDRODINAMIKAI EGYENLETEK • KEZDETI- ÉS PEREMFELTÉTELEK
3. 4. Helyszínrajz 1. 1. Hidegvíz-csatorna 2. 2. Melegvíz-csatorna 3. Sarkantyú 4. Keresztgát 5. 5. Uszodi sziget
Módszertan – dekompozíció és aggregációPélda 3. • C - KONCENTRÁCIÓ VEKTOR • HIDRODINAMIKAI EGYENLETEK • KEZDETI- ÉS PEREMFELTÉTELEK
Módszertan – dekompozíció és aggregációPélda 3. • R(C, P) - REAKCIÓ TAG (félempírikus) • P - PARAMÉTER VEKTOR • IDENTIFIKÁCIÓ SZÜKSÉGES • HIPOTÉZISEK • KALIBRÁLÁS ÉS IGAZOLÁS • ÉRZÉKENYSÉGI ÉSBIZONYTALANSÁGI ELEMZÉSEK
PAKS Morfológiai modell z(x,y) 2D Hidrodinamikai modell kst v(x,y) 2D Transzport modell Dy T(x,y)