High-Resolution Methods and Adaptive Refinement for Tsunami Propagation and Inundation

Abstract

Summary. We describe the extension of high resolution finite volume methods and adaptive refinement for the shallow water equations in the context of tsunami modeling. Godunov-type methods have been used extensively for modeling the shallow water equations in many contexts, however, tsunami modeling presents some unique challenges that must be overcome. We describe some of the specific difficulties associated with tsunami modeling, and summarize some numerical approaches that we have used to overcome these challenges. For instance, we have developed a well-balanced Riemann solver that is appropriate in the deep ocean regime as well as robust in near-shore and dry regions. Additionally, we have extended adaptive refinement algorithms to this application. We briefly describe some of the modifications necessary for using these adaptive methods for tsunami modeling.

Clawpack: building an open source ecosystem for solving hyperbolic PDEs

“…The AMRClaw code formed the starting point but it was necessary to make many modifications to support the requirements of this application, as described briefly below. This code originated with the work of George (George, 2004;George, 2006;George, 2008) and was initially called TsunamiClaw. Later it became clear that many other geophysical flow applications have similar requirements and the code was generalized as…”

Adaptive modelling of long-distance wave propagation and fine-scale flooding during the Tohoku tsunami

“…The quadtree-adaptive tsunami model used in this article, implemented within the Gerris flow solver framework is one of the few adaptive tsunami models published (George and LeVeque, 2008;Harig et al, 2008;Popinet, 2011c). A summary of the model is given in the next section.…”

Lagrangian particle path formulation of multilayer shallow-water flows dynamically coupled to vessel motion

“…Equation (2.11) implies that 15) while adding (2.11) to (2.13) and using (2.1) lead to the mass flux condition …”

“…The studies most relevant to the one in this paper examine the two-layer sloshing problem using a numerical scheme based upon a class of high resolution wave-propagating finite volume methods known as f-wave methods for both the forced [12] and the coupled problem [13]. This f-wave approach is very effective and can be readily be extended to multilayer systems [14] and systems with bottom topography [15], but [13] find the scheme is limited to layer density ratios of ρ 2 /ρ 1 0.7, where ρ 1 and ρ 2 are the fluid densities in the lower and upper layers, respectively, due to a linear growth in the system constraint error. Therefore this approach is not able to model the interior workings of the OWEL WEC, where the air/water density ratio is ρ 2 /ρ 1 ≈ 10 −3 .…”