Skip to content
Texas A&M University
Mathematics

Numerical Analysis Seminar

Spring 2022

 

Date:March 30, 2022
Time:3:00pm
Location:BLOC 628
Speaker:Jesse Chan, Rice University
Title:Recent advances in high order entropy stable discontinuous Galerkin methods
Abstract:High order methods are known to be unstable when applied to nonlinear conservation laws whose solutions exhibit shocks and turbulence. These methods have traditionally required additional filtering, limiting, or artificial viscosity to avoid solution blow up. Entropy stable schemes address this instability by ensuring that physically relevant solutions satisfy a semi-discrete entropy inequality independently of numerical resolution. In this talk, we will review the construction entropy stable discontinuous Galerkin methods and discuss recent work on positivity preservation, reduced order modeling, and robustness under-resolved simulations.

Date:May 4, 2022
Time:3:00pm
Location:BLOC 628
Speaker:Dr J.M. Melenk, TU Wien
Title:Stability and convergence of Galerkin discretizations of the Helmholtz equation in piecewise smooth media
Abstract:We consider the Helmholtz equation with variable coefficients at large wavenumber k. In order to understand how k affects the convergence properties of discretizations of such problems, we develop a regularity theory for the Helmholtz equation that is explicit in k. At the heart of our analysis is the decomposition of solutions into two components: the first component is a piecewise analytic, but highly oscillatory function and the second one has infinite regularity but features wavenumber-independent bounds. This decomposition generalizes earlier decompositions of Melenk and Sauter (2010 & 2022) which considered the Helmholtz equation with constant coefficients, to the case of piecewise analytic coefficients. This regularity theory for the Helmholtz equation with variable coefficients allows for the analysis of high order Galerkin discretizations of the Helmholtz equation that are explicit in the wavenumber k. We show that quasi-optimality is guaranteed under the following scale resolution condition: (a) the approximation order p is selected as p = O(log k) and (b) the mesh size h is such that kh/p is sufficiently small. This scale resolution condition ensures quasi-optimality for a variety of time-harmonic wave propagation problems including FEM-BEM coupling and Maxwell problems.