Adaptive wavelet simulation of global ocean dynamics using a new Brinkman volume penalization

Abstract. In order to easily enforce solid-wall boundary conditions in the presence of complex coastlines, we propose a new mass and energy conserving Brinkman penalization for the rotating shallow water equations. This penalization does not lead to higher wave speeds in the solid region. The error estimates for the penalization are derived analytically and verified numerically for linearized one-dimensional equations. The penalization is implemented in a conservative dynamically adaptive wavelet method for the rotating shallow water equations on the sphere with bathymetry and coastline data from NOAA's ETOPO1 database. This code could form the dynamical core for a future global ocean model. The potential of the dynamically adaptive ocean model is illustrated by using it to simulate the 2004 Indonesian tsunami and wind-driven gyres.

[1]  Nicholas K.-R. Kevlahan,et al.  Computation of turbulent flow past an array of cylinders using a spectral method with Brinkman penalization , 2001 .

[2]  Jerrold E. Marsden,et al.  The Euler-Poincaré Equations in Geophysical Fluid Dynamics , 1999, chao-dyn/9903035.

[3]  C. Amante,et al.  ETOPO1 arc-minute global relief model : procedures, data sources and analysis , 2009 .

[4]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[5]  Stéphane Popinet,et al.  A tree-based solver for adaptive ocean modelling , 2007 .

[6]  T. Dubos,et al.  A conservative adaptive wavelet method for the shallow‐water equations on staggered grids , 2013 .

[7]  John B. Bell,et al.  A Cartesian Grid Projection Method for the Incompressible Euler Equations in Complex Geometries , 1994, SIAM J. Sci. Comput..

[8]  William C. Skamarock,et al.  A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C-grids , 2010, J. Comput. Phys..

[9]  Kenji Satake,et al.  Tsunami Source of the 2004 Sumatra–Andaman Earthquake Inferred from Tide Gauge and Satellite Data , 2007 .

[10]  Suppression of three-dimensional flow instabilities in tube bundles , 2005 .

[11]  Jacques Liandrat,et al.  A penalization method applied to the wave equation , 2005 .

[12]  H. Kreiss,et al.  Stability Theory of Difference Approximations for Mixed Initial Boundary Value Problems. II , 1972 .

[13]  Dmitry Kolomenskiy,et al.  A Fourier spectral method for the Navier-Stokes equations with volume penalization for moving solid obstacles , 2009, J. Comput. Phys..

[14]  Oleg V. Vasilyev,et al.  A Brinkman penalization method for compressible flows in complex geometries , 2007, J. Comput. Phys..

[15]  Vasily Titov,et al.  The Global Reach of the 26 December 2004 Sumatra Tsunami , 2005, Science.

[16]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[17]  Randall J. LeVeque,et al.  LeVeque adaptive refinement on the sphere Logically rectangular finite volume methods with , 2009 .

[18]  H. Kreiss.  Stability theory for difference approximations of mixed initial boundary value problems. I , 1968 .

[19]  Alistair Adcroft,et al.  How slippery are piecewise‐constant coastlines in numerical ocean models? , 1998 .

[20]  Oleg V. Vasilyev,et al.  An Adaptive Wavelet Collocation Method for Fluid-Structure Interaction at High Reynolds Numbers , 2005, SIAM J. Sci. Comput..

[21]  V. Guinot,et al.  Flux and source term discretization in two‐dimensional shallow water models with porosity on unstructured grids , 2006 .

[22]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[23]  Sol Hellerman,et al.  Normal Monthly Wind Stress Over the World Ocean with Error Estimates , 1983 .

[24]  Nicholas K.-R. Kevlahan,et al.  An adaptive wavelet method for turbulence in complex geometries , 2001 .

[25]  A. Adcroft,et al.  Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model , 1997 .

[26]  Matthias Aechtner,et al.  A conservative adaptive wavelet method for the shallow‐water equations on the sphere , 2013, 1404.0405.

[27]  O. Vasilyev,et al.  Hybrid wavelet collocation–Brinkman penalization method for complex geometry flows , 2002 .

[28]  Jörn Behrens,et al.  Tsunami simulations on several scales , 2008 .

[29]  D. Shirokoff,et al.  A Sharp-Interface Active Penalty Method for the Incompressible Navier–Stokes Equations , 2013, J. Sci. Comput..

[30]  Stéphane Popinet,et al.  Quadtree-adaptive tsunami modelling , 2011 .

[31]  Baylor Fox-Kemper,et al.  Adaptive volume penalization for ocean modeling , 2012, Ocean Dynamics.

[32]  J. Hervouet.  Hydrodynamics of Free Surface Flows: Modelling with the Finite Element Method , 2007 .

[33]  Kai Schneider,et al.  Adaptive Wavelet Simulation of a Flow around an Impulsively Started Cylinder Using Penalisation , 2002 .

[34]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[35]  Randall J. LeVeque,et al.  The GeoClaw software for depth-averaged flows with adaptive refinement , 2010, 1008.0455.