898 research outputs found
A Plane Wave Virtual Element Method for the Helmholtz Problem
We introduce and analyze a virtual element method (VEM) for the Helmholtz
problem with approximating spaces made of products of low order VEM functions
and plane waves. We restrict ourselves to the 2D Helmholtz equation with
impedance boundary conditions on the whole domain boundary. The main
ingredients of the plane wave VEM scheme are: i) a low frequency space made of
VEM functions, whose basis functions are not explicitly computed in the element
interiors; ii) a proper local projection operator onto the high-frequency
space, made of plane waves; iii) an approximate stabilization term. A
convergence result for the h-version of the method is proved, and numerical
results testing its performance on general polygonal meshes are presented
First order least squares method with weakly imposed boundary condition for convection dominated diffusion problems
We present and analyze a first order least squares method for convection
dominated diffusion problems, which provides robust L2 a priori error estimate
for the scalar variable even if the given data f in L2 space. The novel
theoretical approach is to rewrite the method in the framework of discontinuous
Petrov - Galerkin (DPG) method, and then show numerical stability by using a
key equation discovered by J. Gopalakrishnan and W. Qiu [Math. Comp. 83(2014),
pp. 537-552]. This new approach gives an alternative way to do numerical
analysis for least squares methods for a large class of differential equations.
We also show that the condition number of the global matrix is independent of
the diffusion coefficient. A key feature of the method is that there is no
stabilization parameter chosen empirically. In addition, Dirichlet boundary
condition is weakly imposed. Numerical experiments verify our theoretical
results and, in particular, show our way of weakly imposing Dirichlet boundary
condition is essential to the design of least squares methods - numerical
solutions on subdomains away from interior layers or boundary layers have
remarkable accuracy even on coarse meshes, which are unstructured
quasi-uniform
- …