Earth
Tsunami erreicht die Küste
Heben Sie über tiefem Wasser einen drei Meter hohen Meeresbuckel und sehen Sie ihn sich teilen und laufen: der Kamm quert den tiefen Ozean mit sqrt(g h), fast zweihundert Meter in der Sekunde, erreicht den Pegel genau, wann die Formel es sagt, verlangsamt sich und wächst über dem Schelf und läuft zwischen den Häusern eines Fischerdorfs den Strand hinauf. Ein steilerer Strand wird weniger weit landeinwärts überflutet; schließen Sie die Küste zu einer Bucht, und die Welle bündelt sich zur doppelten Auflaufhöhe. Das Verfahren hält ein ruhendes Meer über jedem Grund in Ruhe, lässt nie eine Tiefe negativ werden und verbucht jeden Kubikmeter, der den offenen Seerand quert.
Was passiert
Die nichtlinearen Flachwassergleichungen entlang eines 60-km-Streifens in 1200 Zellen: 2000 m Ozean, ein Hang zu einem 100-m-Schelf, ein Strand von 1:50 bis zu 20 m hohem Land, die Breite wahlweise fünffach zu einer Bucht verengt. Die Seiten nutzen Audusses hydrostatische Rekonstruktion mit dem HLL-Fluss, so bleibt ein ruhender See über jedem Grund in Ruhe, trockenes Land eingeschlossen, und Tiefen bleiben positiv; die Manning-Reibung ist implizit; der Schritt folgt einer Courant-Zahl von 0,45. Ein Gauß-Buckel von 3 m Höhe und 5 km Breite ist der gegebene Start - der Bruch, der ihn hob, wird nicht berechnet. Die Auflaufhöhe ist die Höhe der Wasserlinie, die Volumenbilanz schließt ein, was über die offene See abfließt.
Zum Ausprobieren
Vergleichen Sie die Ankunft am Tiefseepegel mit der geschlossenen Langwellenzeit; sehen Sie die Welle auf dem Schelf wachsen; überfluten Sie das Dorf; machen Sie den Strand steiler; schließen Sie die Küste zu einer Bucht.
Das Modell und seine Grenzen
Eine horizontale Dimension, keine Dispersion, kein Brechungsmodell über das stoßerfassende Verfahren hinaus, eine glatte gegebene Quelle, Häuser ohne Einfluss auf die Strömung. Diese Lehrküste ist keine Gefahrenkarte.
Gleichungen des Moduls
Gleichungen ∂(b·h)/∂t + ∂(b·h·u)/∂x = 0 ∂(b·h·u)/∂t + ∂(b·(h·u² + g·h²/2))/∂x = g·h²/2·∂b/∂x − g·b·h·∂z/∂x − friction h* = max(0, h + z − max(z_L, z_R)), HLL flux c = √(g·h), A2/A1 = (h1/h2)^(1/4)·(b1/b2)^(1/2) (Green)
