Dr. Koncsos László egy. docens

1. KÖRNYEZETI RENDSZEREK MODELLEZÉSE Dr. Koncsos László egy. docens Matematikai modellezés

2. 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

3. 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

4. Bevezetés – modell fogalma 3. Példa – QUAL 2K vízmin. modell

5. Bevezetés – modell fogalma 3. Példa – QUAL 2K vízmin. modell

6. Bevezetés – modell fogalma 3. Példa – QUAL 2K vízmin. modell

7. 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

8. 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]

9. 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]

10. 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

11. Fizikai kisminta kísérletek

12. 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)

13. 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

14. 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:

15. 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ű

16. 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

17. 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

18. 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

19. 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ó

20. 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

21. 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

22. 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

23. 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

24. Módszertan – dekompozíció és aggregációPélda 1. • C - KONCENTRÁCIÓ VEKTOR • HIDRODINAMIKAI EGYENLETEK • KEZDETI- ÉS PEREMFELTÉTELEK

26. Módszertan – dekompozíció és aggregációPélda 2. • C - KONCENTRÁCIÓ VEKTOR • HIDRODINAMIKAI EGYENLETEK • KEZDETI- ÉS PEREMFELTÉTELEK

27. D(x,y)=f(v(x,y))

28. 3. 4. Helyszínrajz 1. 1. Hidegvíz-csatorna 2. 2. Melegvíz-csatorna 3. Sarkantyú 4. Keresztgát 5. 5. Uszodi sziget

30. Módszertan – dekompozíció és aggregációPélda 3. • C - KONCENTRÁCIÓ VEKTOR • HIDRODINAMIKAI EGYENLETEK • KEZDETI- ÉS PEREMFELTÉTELEK

31. 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

32. QUAL 2K vízmin. modell

33. QUAL 2K vízmin. modell

34. PAKS Morfológiai modell z(x,y) 2D Hidrodinamikai modell kst v(x,y) 2D Transzport modell Dy T(x,y)

35. Komplexitás

36. Elméleti megalapozottság