Lagrangian blocks on Eulerian mesh for shallow-water wave computations

Waves in shallow water are computed by moving blocks of water in the direction of the flow by the Lagrangian method. Mass and momentum in the blocks are re-distributed onto the fixed Eulerian mesh to minimize distortion at each increment of time. This Lagrangian blocks on Eulerian mesh (LBEM) method...

Full description

Bibliographic Details
Main Author: Tan, Lai Wai
Other Authors: Vincent H Chu (Supervisor)
Format: Others
Language:en
Published: McGill University 2010
Subjects:
Online Access:http://digitool.Library.McGill.CA:80/R/?func=dbin-jump-full&object_id=92363
Description
Summary:Waves in shallow water are computed by moving blocks of water in the direction of the flow by the Lagrangian method. Mass and momentum in the blocks are re-distributed onto the fixed Eulerian mesh to minimize distortion at each increment of time. This Lagrangian blocks on Eulerian mesh (LBEM) method is introduced as an alternative to the classical finite volume methods. The method is positive-depth definite and free of numerical oscillation. A large number of numerical computations are carried out to evaluate the accuracy and the computational stability of the LBEM method. The results of 6 publications in journals and 9 papers in conference proceedings are compiled to produce the total of 8 chapters and 4 appendices of this manuscript-based thesis. Computations using the LBEM method have been conducted for one-dimensional and two-dimensional waves in shallow water. These include the computations of the two-dimensional standing waves in a parabolic bowl, the propagation of oblique shock waves in a square basin, the shoaling of solitary waves and periodic waves over levee, the dam-break waves over dike, and the flood routing through an idealized city. These computational problems have been selected because they have either analytical solutions or available experimental data for comparison with the computations. When exact solutions are available, block refinements are conducted to show the convergence of the LBEM computations toward the exact solutions. The goals of these computations are to show how the LBEM method (i) can track the dry-and-wet water interface without using interface treatment, (ii) can capture shock wave without using any flux or slope limiter, and (iii) can calculate the intermittent flow such as the wave overtopping the levee. Although the tracking of interfaces, the capture of shock waves and the calculation of intermittent overflow are difficult problems for the classical methods, the block advection can compute these problems with absolute s === Des vagues en eau peu profonde sont calculées en déplaçant des blocs de l'eau dans la direction du flux par la méthode Lagrangienne. La masse et l'élan (momentum) dans les blocs sont redistribués sur la maille Eulérienne fixée pour minimiser la déformation à chaque incrément de temps. Cette méthode des blocs Lagrangiens sur la maille Eulérienne (LBEM) est présentée comme alternative aux méthodes classiques de volume finie. La méthode est positif-profondeur définie et exempte de l'oscillation numérique. Un grand nombre de calculs numériques sont effectués pour évaluer l'exactitude et la stabilité informatique de la méthode de LBEM. Les résultats de six publications en journaux et de neuf articles dans des actes de la conférence sont compilés pour produire le total de huit chapitres et de quatre annexes de ceci thèse manuscrit-basée. Des calculs suivre la méthode de LBEM avaient été conduits pour les vagues en eau peu profonde en une dimension et en deux dimensions. Ceux-ci incluent les calculs des vagues stationnaires bidimensionnelles dans une cuvette parabolique, de la propagation des vagues chocs obliques dans un bassin carré, du banc des vagues solitaires et les vagues périodiques par-dessus de la levée, des vagues de coupure de barrage par-dessus de digue, et du cheminement d'inondation à travers une ville idéalisée. Ces problèmes informatiques avaient été choisis parce qu'ils ont les solutions analytiques ou les données expérimentales disponibles. Quand les solutions exactes sont disponibles, des raffinements de bloc sont conduits pour montrer la convergence des calculs de LBEM vers les solutions exactes. Les buts de ces calculs sont de montrer comment la méthode de LBEM (i) peut tracer l'interface sèche-et-mouillée de l'eau sans employer le traitement d'interface, (ii) peut capturer le vague choc sans utiliser n'importe quel limiteur de flux ou de pente, et (iii) peut calculer le flux intermittent tel que la vag