Earth
Le tsunami atteint le rivage
Soulevez une bosse de mer de trois mètres au-dessus d’eaux profondes et regardez-la se scinder et courir : la crête traverse l’océan profond à sqrt(g h), près de deux cents mètres par seconde, atteint le marégraphe exactement quand la formule le dit, ralentit et grandit sur le plateau continental, puis monte sur la plage parmi les maisons d’un village de pêcheurs. Une plage plus raide est inondée moins loin ; fermez la côte en baie et la vague s’y concentre jusqu’à doubler le run-up. Le schéma garde au repos une mer au repos sur tout fond, ne rend jamais une profondeur négative et compte chaque mètre cube qui franchit la limite de mer ouverte.
Ce qui se passe
Les équations non linéaires en eau peu profonde le long d’une bande de 60 km en 1200 cellules : 2000 m d’océan, un talus vers un plateau de 100 m, une plage de 1:50 jusqu’à une terre de 20 m, la largeur pouvant se rétrécir cinq fois en baie. Les faces utilisent la reconstruction hydrostatique d’Audusse avec le flux HLL, si bien qu’un lac au repos reste au repos sur tout fond, terre sèche comprise, et que les profondeurs restent positives ; le frottement de Manning est implicite ; le pas suit un nombre de Courant de 0,45. Une bosse gaussienne de 3 m de haut et 5 km de large est le départ donné - la rupture qui l’a soulevée n’est pas calculée. Le run-up est la hauteur de la ligne d’eau, le bilan de volume inclut ce qui sort par la mer ouverte.
À essayer
Comparez l’arrivée au marégraphe profond au temps d’onde longue de la formule ; regardez la vague grandir sur le plateau ; inondez le village ; rendez la plage plus raide ; fermez la côte en baie.
Le modèle et ses limites
Une dimension horizontale, sans dispersion, sans modèle de déferlement au-delà du schéma de capture des chocs, une source donnée lisse, des maisons sans effet sur l’écoulement. Cette côte pédagogique n’est pas une carte d’aléa.
Équations du module
Équations ∂(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)
