Discontinuous Galerkin unsteady discrete adjoint method for real-time efficient tsunami simulations

Authors:
SéBastien Blaise;Amik St-Cyr;Dimitri Mavriplis;Brian Lockwood
Affiliations:
Institute for Mathematics Applied to Geosciences, National Center for Atmospheric Research, Boulder, CO, USA and Institute of Mechanics, Materials and Civil Engineering, Université catholique ...;Institute for Mathematics Applied to Geosciences, National Center for Atmospheric Research, Boulder, CO, USA and Shell Technical Center Houston, 3333 Highway 6 South, Houston, Texas 77082, USA;Department of Mechanical Engineering, University of Wyoming, Laramie, WY, USA;Department of Mechanical Engineering, University of Wyoming, Laramie, WY, USA
Venue:
Journal of Computational Physics
Year:
2013

Citing 11
Cited 0

A limited memory algorithm for bound constrained optimization

SIAM Journal on Scientific Computing
Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization

ACM Transactions on Mathematical Software (TOMS)
Grid adaptation for functional outputs: application to two-dimensional inviscid flows

Journal of Computational Physics
A steepest descent method for vector optimization

Journal of Computational and Applied Mathematics
A study of spectral element and discontinuous Galerkin methods for the Navier-Stokes equations in nonhydrostatic mesoscale atmospheric modeling: Equation sets and test cases

Journal of Computational Physics
A Fully Implicit Jacobian-Free High-Order Discontinuous Galerkin Mesoscale Flow Solver

ICCS 2009 Proceedings of the 9th International Conference on Computational Science
Adjoint-based h-p adaptive discontinuous Galerkin methods for the 2D compressible Euler equations

Journal of Computational Physics
A standard test set for numerical approximations to the shallow water equations in spherical geometry

Journal of Computational Physics
Computational aspects of chemical data assimilation into atmospheric models

ICCS'03 Proceedings of the 2003 international conference on Computational science
Discrete and continuous adjoint approaches to estimate boundary heat fluxes in falling films

Optimization Methods & Software
Extreme-Scale AMR

Proceedings of the 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis

Quantified Score

Hi-index	31.45

Visualization

Abstract

An unsteady discrete adjoint implementation for a discontinuous Galerkin model solving the shallow water wave equations on the sphere is presented. Its use for tsunami simulations is introduced to reconstruct the initial condition automatically from buoy measurements. Based on this feature, a real-time tsunami model is developed, using several numerical tools such as a high-order discretization, hp-refinement, parallel dynamic load balancing and adjoint-based data assimilation. The model is able to reconstruct the tsunami source and accurately forecast its far-field propagation (e.g. from Japan to Chile, at a distance of about 17000km) in a computational time 20 times faster than the physical propagation time, to which the data collecting time needs to be added. The work presented constitutes a step towards an efficient nonlinear tsunami warning model. Additional features could be added for more complete realistic forecasts.