Retencja wody na przypadkowych powierzchniach
Retencja wody na przypadkowych powierzchniach to symulacja zatrzymywania wody w stawach na powierzchni komórek o różnej wysokości w regularnym układzie, takim jak kwadratowa siatka, gdzie woda spływa na każdą komórkę w systemie. Granice systemu są otwarte i umożliwiają wypływ wody. Woda zostanie uwięziona w stawach i ostatecznie wszystkie stawy wypełnią się do maksymalnej wysokości, a dodatkowa woda przepłynie przez przelewy i poza granice systemu. Problem polega na znalezieniu ilości wody uwięzionej lub zatrzymanej na danej powierzchni. Zostało to szeroko zbadane dla przypadkowych powierzchni.
Losowe powierzchnie
Jednym z systemów, w którym badano kwestię retencji, jest powierzchnia o losowych wysokościach. Tutaj można zmapować losową powierzchnię do perkolacji miejsca, a każda komórka jest odwzorowywana na miejsce na leżącym poniżej wykresie lub siatce reprezentującej system. Korzystając z teorii perkolacji , można wyjaśnić wiele właściwości tego systemu. Jest to przykład modelu inwazji perkolacji, w którym płyn jest wprowadzany do układu z dowolnego losowego miejsca.
W hydrologii zajmuje się odpływem i tworzeniem zlewni. Granica między różnymi zlewniami ( działami wodnymi w Ameryce Północnej) tworzy przepaść melioracyjną o wymiarze fraktalnym około 1,22.
Problem retencji można odwzorować na standardową perkolację. Na przykład dla systemu pięciu równie prawdopodobnych poziomów ilość zmagazynowanej wody R 5 jest po prostu sumą wody zmagazynowanej w systemach dwupoziomowych R 2 ( p ) z różnymi ułamkami poziomów p w najniższym stanie:
- R5 = R2 (1/5) + R2 ( 2/5 ) + R2 (3/5) + R2 ( 4/5 )
Typowe systemy dwupoziomowe 1,2 z p = 0,2, 0,4, 0,6, 0,8 pokazano po prawej stronie (niebieski: mokry, zielony: suchy, żółty: przelewy graniczące z terenami mokrymi). Retencja netto pięciopoziomowego systemu jest sumą tych wszystkich elementów. Najwyższy poziom nie zatrzymuje wody, ponieważ znajduje się znacznie powyżej progu perkolacji dla kwadratowej siatki, 0,592746.
Retencja systemu dwupoziomowego R 2 ( p ) to ilość wody podłączonej do stawów, które nie dotykają granicy systemu. Kiedy p jest powyżej krytycznego progu perkolacji p c , pojawi się gromada lub staw perkolacyjny, który odwiedzi cały system. Prawdopodobieństwo, że punkt należy do gromady perkolacyjnej lub „nieskończonej” jest zapisywane jako P ∞ w teorii perkolacji i jest powiązane z R 2 ( p ) przez R 2 ( p )/ L 2 = p − P ∞ gdzie L jest rozmiarem kwadratu. Zatem zachowanie systemu wielopoziomowego można powiązać z dobrze znaną wielkością w teorii perkolacji .
Do pomiaru retencji można zastosować algorytm zalewania , w którym woda jest wprowadzana od granic i zalewana najniższym przelewem w miarę podnoszenia się poziomu. Retencja to po prostu różnica w poziomie wody, w jakiej miejsce zostało zalane, pomniejszona o wysokość terenu pod nim.
Oprócz systemów dyskretnych poziomów opisanych powyżej, można uczynić zmienną terenową zmienną ciągłą, powiedzmy od 0 do 1. Podobnie można sprawić, aby sama wysokość powierzchni była ciągłą funkcją zmiennych przestrzennych. We wszystkich przypadkach pozostaje podstawowa koncepcja mapowania do odpowiedniego systemu perkolacji .
Ciekawym wynikiem jest to, że kwadratowy układ n dyskretnych poziomów może zatrzymać więcej wody niż układ n+1 poziomów, dla wystarczająco dużego rzędu L > L *. To zachowanie można zrozumieć za pomocą teorii perkolacji, której można również użyć do oszacowania L * ≈ ( p - p c ) - ν gdzie ν = 4/3, p = i */ n gdzie i * jest największą wartością i dla który ja / n < p c , a p c = 0,592746 jest progiem perkolacji miejsca dla sieci kwadratowej. Symulacje numeryczne dają następujące wartości L *, które ekstrapoluje się na wartości niecałkowite. Na przykład R 2 < R 3 dla L ≤ 51, ale R 2 > R 3 dla L ≥ 52:
N | n + 1 | L * | Retencja w L * |
---|---|---|---|
2 | 3 | 51.12 | 790 |
4 | 5 | 198,1 | 26000 |
7 | 8 | 440,3 | 246300 |
9 | 10 | 559.1 | 502000 |
12 | 13 | 1390,6 | 428850 |
14 | 15 | 1016.3 | 2607000 |
Gdy n staje się większe, przecinanie staje się coraz rzadsze, a wartość L *, w której występuje przecinanie, nie jest już monotoniczną funkcją n .
Retencja, gdy powierzchnia nie jest całkowicie przypadkowa, ale skorelowana z wykładnikiem Hursta H , jest omówiona w Schrenk i in.
Zobacz też
Dalsza lektura
- Pickover, Clifford (2002). Zen magicznych kwadratów, kół i gwiazd: wystawa zaskakujących struktur w różnych wymiarach . Princeton, NJ: Princeton University Press. ISBN 978-0-691-11597-9 .
- Stauffer, Dietrich; Aharony, A. (1994). Wprowadzenie do teorii perkolacji . Londyn Bristol, Pensylwania: Taylor i Francis. ISBN 978-0-7484-0253-3 .
Linki zewnętrzne
- https://commons.wikimedia.org/wiki/Category:Associative_magic_squares_of_order_4
- Hugo Pfoertnera. Sekwencja OEIS A201126 (Maksymalna retencja wody magicznego kwadratu rzędu n) , z linkami do obrazów magicznego kwadratu
- Hugo Pfoertnera. Sekwencja OEIS A201127 (maksymalna retencja wody półmagicznego kwadratu rzędu n)
- Strona dyskusyjna o konkursach programistycznych Ala Zimmermanna
- Przedmiot w daremnej szafie
- OEIS A261798 (maksymalna retencja wody asocjacyjnego magicznego kwadratu rzędu n)
- OEIS A268311 (Liczba wolnych poliomino, które tworzą ciągłą ścieżkę komórek połączonych krawędziami obejmujących n X n kwadrat w obu wymiarach) — wyliczenie poliomino i wzory jezior
- OEIS A275359 (Maksymalne uwięzienie liczb w kostkach liczbowych n X n X n z pełnymi objętościami uwięzionymi) — uaktualnienie modelu z 2D do 3D
- [1] Przyroda 2018
- [2] Histogram retencji wody jako problem obliczeniowy
- http://oeis.org/A331507/ Maksymalna liczba stawów