Project van het vak Modelleren-B - AM2050-B voor Technische Wiskunde, aan de Technische Universiteit Delft
J. Bus
J. Imhof
B. Verhage
R. Roovers
Het model kan gestart worden door het bestand Execution.py te runnen.
Model dat gebruik maakt van classes zodat variabelen handig over te dragen zijn naar verschillende files. Naast het basismodel, zijn stroming en en diffusie meegenomen. Moleculair diffusie vindt niet plaats op dezelfde manier voor fytoplankton, maar op grotere schaal kunnen we de uitweiding van concencreerde benaderen met dezelfde soort vergelijkingen (met weliswaar diffusieconstanten met andere orde van grootte). De parameters die hier een rol in spelen moeten nog worden gecalibreerd/gefit of in ieder geval onderbouwd.
Er zijn drie verschillende randvoorwaarden geimplementeerd:
- Dirichelet
- Neumann
- Periodiek
Hier leggen we een constante rand op aan het systeem. Deze constante rand leidt tot gekke compensatie respons in het systeem waardoor het zeker niet fysisch oogt.
Een idee geopperd door prof. Heemink is een periodiek systeem introduceren (dus eind gaat helemaal terug naar begin), nu ook geïmplementeerd. Hierin gaan we ervanuit dat de linker rand en de rechter rand aan elkaar geplakt zijn, vergelijk een lijn over de oceaan die rond de aarde gaat, of een circulair aquarium. Geen gekke respons, maar toepassing is uiteraard beperkt.
Daarnaast is er nog een Neumann boundary gemaakt met zero flux, deze leidt tot het meest fysische (ogende) resultaat. We gaan er hier bij uit de concentraties net buiten het reservoir gelijk is aan die op de randen van het reservoir.