Zadržování vody na matematických površích - Water retention on mathematical surfaces

Voda se nalije na povrch lega.
Ilustrace zadržování vody na povrchu.

Zadržování vody na matematických površích je zachycování vody v rybnících na povrchu buněk různých výšek na pravidelném poli, jako je čtvercová mříž, kde voda prší na každou buňku v systému. Hranice systému jsou otevřené a umožňují odtékání vody. Voda bude zachycena v rybnících a nakonec se všechny rybníky naplní do své maximální výšky, přičemž jakákoli další voda bude proudit přes přelivy a za hranice systému. Problém je najít množství vody zachycené nebo zadržené pro daný povrch. Toto bylo rozsáhle studováno pro dva matematické povrchy: magické čtverce a náhodné povrchy. Model lze také použít na trojúhelníkovou mřížku.[1]

Kouzelné čtverce

Udržení na magickém čtverci 5x5.
Kouzelný čtverec 5 × 5 s maximální retencí.

Kouzelné čtverce byly studovány více než 2000 let. V roce 2007 byla navržena myšlenka studia zadržování vody na magickém čtverci.[2] V roce 2010 proběhla soutěž na programovacích soutěžích Al Zimmermanna[3] vytvořila v současnosti známé maximální retenční hodnoty pro magické čtverce v pořadí 4 až 28.[4] Zde najdete výpočetní nástroje používané k prozkoumání a ilustraci tohoto problému.[5][6][7][8]


Existuje 4 211 744 různých retenčních vzorů pro čtverec 7 × 7. K dosažení maximální retence je nejlepší kombinace jezera a rybníků. Žádné známé vzory pro maximální retenci nemají ostrov v rybníku nebo jezeře.[2]

Níže jsou zobrazeny magické čtverce s maximální retencí pro objednávky 7-9:[4]

Níže uvedené obrázky ukazují magický čtverec 10x10. Je možné podívat se na výše uvedené vzory a předpovědět, jaký bude vzorec pro maximální zadržení čtverce 10x10? Nebyla vyvinuta žádná teorie, která by dokázala předpovědět správnou kombinaci jezera a rybníků pro všechny objednávky, nicméně platí některé principy. První barevně kódovaná čísla ukazují konstrukční princip toho, jak je kolem jezera a rybníků rozmístěno největší možné množství. Druhá a třetí figurka ukazují slibné vzorce, které byly vyzkoušeny, ale nedosáhly maximální retence.[2]

Několik objednávek má více než jeden vzor pro maximální uchování. Obrázek níže ukazuje dva vzory magického čtverce 11x11 se zřejmou maximální retencí 3 492 jednotek:[4]

The nejdokonalejší kouzelné čtverce vyžadují, aby všechny (n-1) ^ 2 nebo v tomto případě všech 121 2x2 rovinných podmnožin měly stejný součet. (několik příkladů označených žlutým pozadím, červeným písmem). Oblasti zcela obklopené větším počtem jsou zobrazeny s modrým pozadím.[9]

Nejdokonalejší magický square.jpg

Před rokem 2010, pokud jste chtěli příklad kouzelného čtverce většího než 5 × 5, bylo nutné dodržovat chytrá konstrukční pravidla, která poskytla velmi izolované příklady. 13x13 pandiagonální magický čtverec níže je takový příklad. Utility Harryho Whitea CompleteSquare [5] umožňuje komukoli použít magický čtverec, protože hrnčíř by použil kus hlíny. Druhý obrázek ukazuje kouzelný čtverec 14x14, který byl zformován do podoby rybníků, které píší data 1514 - 2014. Animace zaznamenává, jak byl povrch vytesán tak, aby naplnil všechny rybníky, než voda odtéká z náměstí. Toto náměstí ctí 500. výročí slavného kouzelnického náměstí Durer v Melencolia I..

13 x 13 Pandiagonal Magic Sqauare.png
Magic square water retention.gif

Tento obrázek také poskytuje příklad čtverce a jeho doplňku, které mají stejný vzor retence. K dispozici jsou magické čtverce 137 řádu 4 a 3 254 798 řádu 5, které nezadržují vodu.[2]

16 x 16 asociativní magický čtverec zadržující 17840 jednotek. Jezero na prvním obrázku vypadá trochu ošklivěji než obyčejně. Jarek Wroblewski poznamenává, že dobré vzory pro maximální retenci budou mít stejný nebo téměř stejný počet zadržovacích buněk na každé okrajové hraně (v tomto případě 7 buněk na každé hraně) [3] Druhý obrázek je upraven a zastíněn v horních a dolních 37 hodnotách.

Magický čtvercový vzor.jpeg
Asociativní magický čtverec 2013.jpeg

Na obrázku níže je kouzelný čtverec formátu Luo-Shu 17x17.[10]Zdá se, že konstrukční metoda formátu Luo-Shu produkuje maximální počet rybníků. Odtoková cesta pro buňku zeleně je dlouhá a nakonec se vylije z náměstí u žluté přepadové buňky.

Obrázek vpravo ukazuje, jaké informace lze odvodit z pohledu na skutečný obsah vody v každé buňce. Pouze 144 hodnot je zvýrazněno, aby čtverec nevypadal příliš zaneprázdněn. Zaměření na zelenou buňku se základní hodnotou 7, nejvyšší překážkou na cestě ven je její sousední buňka s hodnotou 151 (ponecháno 151-7 = 144 jednotek). Voda pršela do této buňky a opouští čtverec u žluté buňky 10.


Mario Mamzeris vynalezl vlastní metodu pro konstrukci magických čtverců zvláštního řádu. Níže je zobrazen jeho asociativní magický čtverec řádu 19.[11]



Objednat 19 Magic Square.png


Na magickém čtverci 21 x 21 pod všemi sudými čísly tvoří přehrady a rybníky a všechna lichá čísla poskytují únikové cesty.[12]

Magic Square Sudé číslo Zadržování vody.png



Počítačový věk nyní umožňuje zkoumat fyzikální vlastnosti magických čtverců libovolného řádu. Obrázek níže ukazuje největší magický čtverec studovaný v soutěži. U L> 20 se počet proměnných / rovnic zvyšuje do bodu, kdy je předvídatelný vzor pro maximální retenci.

L = 28: 219822 zadržovacích jednotek
152596592577137122822562836572556562847267252542856542532865567367447164573
96406642571726277167157142866829744303173074350681680679515267820664611
265665355722496618714849512172140077441813017629354174947910617538914823068243644
66314724356313513751891984492137758747853913932660451750461566141442638477677276
2667201645723542264911715121177762472445034358562940614463475159246212513451444672
711155651161003575791126377771084694335468055952546852622714675236855732821246671
710161531742221193536277786429745654447417847341056351533140338775340256930445670
70917703251685094457791663664018392482129338408492585529369298424754582519676275
26771915610345553178039135853776142367309522245320437632386545497224123755161675277
264718444600508781196553653524883446241042165519861637029423310141649010975647652
66218723417782310564606420483359518548246475586283855716914922333523586113733274
2636692187831274295817739913688351602538636635371220745709963354349850217348727
66119784407179184195609393495203567360576394384388137625154523229489485219314738279
26874859730750561544131558356219454244635037258831644312016289102560317110329737272
7292052117723234012841115212233424160538336141257820261973611549589587432568736278
26274668580242187558183398601594182373296460349332556205419614323547586207114735273
269745458131111783376105326126223825936555444836261382574172493466126145630734280
2617471584655982214592145241673746085334093193305953481814283054535841996176533651
6602177353656194345165204381621528447211500135452342363301396527185225764306666281
270694517772392431312240375190617151913243335202312155113475402389776341370749650
2606931054057715502953803023363116202341334271975161509060736442576248667530703271
536923001636317703761911575524144155554226265903395077918814776143030843613270254
6832239742353537976915542149432245439021751062310720059118676034134659323711524696
684232281183775753037683275344875734384724575994644391437596041381607239512432697
480691209378440504140501767812011594042104675775716975819342647093596639180701499
64825869529919220848132131876646396635068423623975734370845024317043460370662641
10649386903937689688687407323674274174073973141667357054234137046664312
266472876866852882522516582897282502492902482916552926535669869970047664284
Jarek Wroblewski 24. března 2010

Toto je panmagický čtverec 32x32. Dwane Campbell pomocí binárních konstrukčních metod vytvořil tento zajímavý příklad zadržování vody.[13] Obslužný program GET TYPE aplikovaný na tento čtverec ukazuje, že má následující vlastnosti: 1) normální magie 2) pandiagonální 3) ohnutá úhlopříčka obousměrně 4) samoplnění.[Citace je zapotřebí ]

32 bimagic square.png

Náhodné povrchy

Zadržování vody na náhodném povrchu.
Zadržování vody na náhodném povrchu 10 úrovní.
Pět úrovní

Dalším systémem, ve kterém byla studována retenční otázka, je povrch náhodných výšek. Zde je možné namapovat náhodnou plochu na perkolaci místa a každá buňka je mapována na místo na podkladovém grafu nebo mřížce, která představuje systém. Použitím teorie perkolace, lze vysvětlit mnoho vlastností tohoto systému. Je to příklad modelu perkolační invaze, ve kterém je tekutina zaváděna do systému z libovolného náhodného místa.[14][15][16]

v hydrologie, jeden se zabývá odtokem a tvorbou povodí.[17] Hranice mezi různými povodí (povodí v Severní Americe) tvoří a drenážní předěl s fraktální dimenze asi 1,22.[18][19][20]

Problém zadržení lze mapovat na standardní perkolaci.[21][22][23] Například pro systém pěti stejně pravděpodobných úrovní je množství skladované vody R5 je jen součet vody uložené ve dvouúrovňových systémech R2p) s různými zlomky úrovní p v nejnižším stavu:

R5 = R2(1/5) + R2(2/5) + R2(3/5) + R2(4/5)

Vpravo jsou zobrazeny typické dvouúrovňové systémy 1,2 s p = 0,2, 0,4, 0,6, 0,8 (modrá: mokrá, zelená: suchá, žlutá: přelivy hraničící s mokrými místy). Čistá retence pětistupňového systému je součtem všech těchto. Nejvyšší úroveň nezachytí žádnou vodu, protože je vysoko nad práh perkolace pro čtvercovou mřížku, 0,592746.

Zachování dvouúrovňového systému R.2(p) je množství vody připojené k rybníkům, které se nedotýkají hranice systému. Když p je nad kritickým prahem perkolace p C, bude existovat prosakující shluk nebo rybník, který navštíví celý systém. Pravděpodobnost, že bod patří do prosakujícího nebo „nekonečného“ klastru, je zapsána jako P v teorii perkolace a souvisí s R2(p) od R2(p)/L2p − P kde L je velikost čtverce. Retence víceúrovňového systému tedy může souviset se známým množstvím v teorie perkolace.

K měření retence lze použít a zaplavovací algoritmus ve kterém je voda přiváděna z hranic a zaplavuje nejnižším přepadem, jak je hladina zvýšena. Zadržení je pouze rozdíl ve vodní hladině, kterou místo zaplavilo minus výška terénu pod ním.

Kromě výše popsaných systémů diskrétních úrovní lze z terénní proměnné udělat spojitou proměnnou, řekněme od 0 do 1. Podobně lze z výšky povrchu udělat spojitou funkci prostorových proměnných. Ve všech případech je základní koncept mapování vhodný perkolace systém zůstává.

Zvláštní výsledek je, že čtvercový systém n diskrétních úrovní může zadržet více vody než systém n + 1 úrovní pro dostatečně velký řád L> L *. Tomuto chování lze porozumět prostřednictvím teorie perkolace, kterou lze také použít k odhadu L * ≈ (p - pC)−ν kde ν = 4/3, p = i * / n kde i * je největší hodnota i taková, že i / n Ca strC = 0,592746 je prahová hodnota perkolace stránek pro čtvercovou mříž. Numerické simulace poskytují následující hodnoty L *, které jsou extrapolovány na jiné než celočíselné hodnoty. Například, R2 < R3 pro L ≤ 51, ale R2 > R3 pro L ≥ 52:[21]

nn + 1L *Retence ve společnosti L *
2351.12790
45198.126000
78440.3246300
910559.1502000
12131390.6428850
14151016.32607000

Jak se n zvětšuje, křížení je čím dál méně časté a hodnota L *, kde dochází ke křížení, již není monotónní funkcí n.

Zadržení, když povrch není zcela náhodný, ale souvisí s a Hurstův exponent H je popsána v.[23]

Algoritmy

Následující časová osa ukazuje použití různých algoritmů, které rozšířily velikost čtverce, který lze vyhodnotit pro uchování

2007 Definujte všechny procházky, které se vyhýbají sousedům, z každé vnitřní buňky do vnější a poté seřaďte všechny tyto cesty, abyste měli co nejmenší překážku nebo hodnotu buňky. Nejmenší hodnota překážky minus hodnota vnitřní buňky zajišťuje zadržování vody pro tuto vnitřní buňku (záporné hodnoty jsou nastaveny na hodnotu zadržení 0). Počet procházek vyhýbajících se sousedům, které mají být vyhodnoceny, roste exponenciálně s velikostí čtverce a omezuje tak tuto metodiku na L <6.[2]

2009 Algoritmus zaplavení - voda je přiváděna z hranic a zaplavuje nejnižším přelivem při zvyšování hladiny. Zadržení je pouze rozdíl ve vodní hladině, kterou místo zaplavilo minus výška terénu pod ním. Algoritmus zaplavení umožňuje vyhodnocení zadržování vody až do L <10 000.[21] Tento algoritmus je podobný Meyerův zaplavovací algoritmus který byl použit při analýze topografických povrchů.

2011 S vědomím, že systém na úrovni n lze rozdělit na kolekci systémů na dvou úrovních s různou pravděpodobností, lze použít standardní perkolační algoritmy k nalezení retence jako jednoduše celkového počtu webů na nižší úrovni minus odtokové oblasti (shluky webů nízké úrovně dotýkajících se hranice). Nová aplikace Algoritmus Hoshen-Kopelman ve kterém jsou přidány oba řádky a sloupce jeden po druhém, umožňuje L být velmi velký (až 109), ale úvahy o výpočetním čase omezují L řádově na 107.[24]

Cesty, které odvádějí vodu z náměstí, používané v algoritmu procházky vyhýbajícím se sousedům

Panel dole zleva doprava ukazuje: 1) tři jedinečné vnitřní polohy pro čtverec 5 × 5; 2 a 4) opravte cesty od čtverce šedě pro vnitřní rohovou buňku červeně; 3) nesprávná šedá cesta, protože voda nemůže cestovat po úhlopříčkách; 5) tato cesta je správná, ale mezi šedými buňkami je možný zkrat. Procházky vyhýbající se sousedům definují jedinečné nebo neredundantní cesty, které odvádějí vodu z náměstí.

Viz také

Reference

  1. ^ https://oeis.org/A303295 OEIS A303295
  2. ^ A b C d E Craig Knecht, http://www.knechtmagicsquare.paulscomputing.com
  3. ^ A b Al Zimmermann http://www.azspcs.net/Contest/MagicWater/FinalReport
  4. ^ A b C Harvey Heinz, http://www.magic-squares.net/square-update-2.htm#Knecht
  5. ^ A b Harry White, http://budshaw.ca/Download.html
  6. ^ Walter Trump http://www.trump.de/magic-squares/
  7. ^ Johan Ofverstedt,http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-176018
  8. ^ Hasan M., Masbaul Alam Polash M. (2020) Efektivní místní vyhledávání založené na omezeních pro maximalizaci zadržování vody na magických čtvercích. In: Hitendra Sarma T., Sankar V., Shaik R. (eds) Emerging Trends in Electrical, Communications, and Information Technologies. Přednášky v elektrotechnice, sv. 569. Springer, Singapur
  9. ^ Sloane, N. J. A. (vyd.). "Pořadí A270205 (počet 2 X 2 planárních podmnožin v krychli n X n X n)". The On-line encyklopedie celočíselných sekvencí. Nadace OEIS.
  10. ^ Harvey Heinz,http://www.magic-squares.net/square-update.htm
  11. ^ https://www.oddmagicsquares.com
  12. ^ „Mapování oblasti“.
  13. ^ http://magictesseract.com
  14. ^ Chayes, J. T .; L. Chayes; C. M. Newman (1985). "Stochastická geometrie perkolace invaze". Komunikace v matematické fyzice. 101 (3): 383–407. Bibcode:1985CMaPh.101..383C. doi:10.1007 / BF01216096.
  15. ^ Damron, Michael; Artëm Sapozhnikov; Bálint Vágvölgyi (2009). "Vztahy mezi perkolací invaze a kritickou perkolací ve dvou dimenzích". Annals of Probability. 37 (6): 2297–2331. arXiv:0806.2425. doi:10.1214 / 09-AOP462.
  16. ^ van den Berg, Jacob; Antal Járai; Bálint Vágvölgyi (2007). Msgstr "Velikost rybníka ve 2D prosakování". Pravděpodobnost elektronických komunikací. 12: 411–420. arXiv:0708.4369. Bibcode:2007arXiv0708.4369V. doi:10.1214 / ECP.v12-1327.
  17. ^ Tetzlaff, D .; McDonnell, J. J .; Uhlenbrook, S .; McGuire, K. J .; Bogaart, P. W .; Naef, F .; Baird, A. J .; Dunn, S. M .; Soulsby, C. (2011). „Konceptualizace spádových procesů: jednoduše příliš složité?“. Hydrologické procesy. 22 (11): 1727–1730. Bibcode:2008HyPr ... 22.1727T. doi:10,1002 / hyp.7069.
  18. ^ Fehr, E .; D. Kadau; N. A. M. Araújo; J. S. Andrade ml .; H. J. Herrmann (2011). "Škálování vztahů pro povodí". Fyzický přehled E. 84 (3): 036116. arXiv:1106.6200. Bibcode:2011PhRvE..84c6116F. doi:10.1103 / PhysRevE.84.036116. PMID  22060465.
  19. ^ Schrenk, K. J .; N. A. M. Araújo; J. S. Andrade ml .; H. J. Herrmann (2012). „Zlomení hodnocených povrchů“. Vědecké zprávy. 2: 348. arXiv:1103.3256. Bibcode:2012NatSR ... 2E.348S. doi:10.1038 / srep00348. PMC  3317236. PMID  22470841.
  20. ^ Fehr, E .; D. Kadau; J. S. Andrade Jr; H. J. Herrmann (2011). "Dopad poruch na povodí". Dopisy o fyzické kontrole. 106 (4): 048501. arXiv:1101.5890. Bibcode:2011PhRvL.106d8501F. doi:10.1103 / PhysRevLett.106.048501. PMID  21405368.
  21. ^ A b C Knecht, Craig; Walter Trump; Daniel ben-Avraham; Robert M. Ziff (2012). „Retenční kapacita náhodných povrchů“. Dopisy o fyzické kontrole. 108 (4): 045703. arXiv:1110.6166. Bibcode:2012PhRvL.108d5703K. doi:10.1103 / PhysRevLett.108.045703. PMID  22400865.
  22. ^ Baek, Seung Ki; Beom Jun Kim (2012). „Kritický stav modelu zadržování vody“. Fyzický přehled E. 85 (3): 032103. arXiv:1111.0425. Bibcode:2012PhRvE..85c2103B. doi:10.1103 / PhysRevE.85.032103. PMID  22587136.
  23. ^ A b Schrenk, K. J .; N. A. M Araújo; R. M. Ziff; H. J. Herrmann (2014). "Retenční kapacita korelovaných povrchů". Fyzický přehled E. 89 (6): 062141. arXiv:1403.2082. Bibcode:2014PhRvE..89f2141S. doi:10.1103 / PhysRevE.89.062141. PMID  25019758.
  24. ^ Hoshen, Joseph (1998). "O aplikaci vylepšeného Hoshen-Kopelmanova algoritmu pro analýzu obrazu". Písmena pro rozpoznávání vzorů. 19 (7): 575–584. doi:10.1016 / S0167-8655 (98) 00018-x.

Další čtení

  • Pickover, Clifford (2002). Zen magických čtverců, kruhů a hvězd: Výstava překvapivých struktur napříč dimenzemi. Princeton, NJ: Princeton University Press. ISBN  978-0-691-11597-9.
  • Stauffer, Dietrich; Aharony, A. (1994). Úvod do teorie perkolace. London Bristol, PA: Taylor & Francis. ISBN  978-0-7484-0253-3.

externí odkazy