Energy stable numerical methods for hyperbolic partial differential equations using overlapping domain decomposition

  • Authors:

  • Adam Reichert;Michael T. Heath;Daniel J. Bodony

  • Affiliations:

  • Department of Computer Science, University of Illinois at Urbana-Champaign, United States;Department of Computer Science, University of Illinois at Urbana-Champaign, United States;Department of Aerospace Engineering, University of Illinois at Urbana-Champaign, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2012

Quantified Score

Hi-index 31.45

Visualization

Abstract

Overlapping domain decomposition methods, otherwise known as overset grid or chimera methods, are useful for simplifying the discretization of partial differential equations in or around complex geometries. Though in wide use, such methods are prone to numerical instability unless numerical diffusion or some other form of regularization is used, especially for higher-order methods. To address this shortcoming, high-order, provably energy stable, overlapping domain decomposition methods are derived for hyperbolic initial boundary value problems. The overlap is treated by splitting the domain into pieces and using generalized summation-by-parts derivative operators and polynomial interpolation. New implicit and explicit operators are derived that do not require regularization for stability in the linear limit. Applications to linear and nonlinear problems in one and two dimensions are presented, where it is found the explicit operators are preferred to the implicit ones.