232 research outputs found
Discontinuous Petrov-Galerkin method based on the optimal test space norm for one-dimensional transport problems
We revisit the finite element analysis of convection dominated flow problems within the recently developed Discontinuous Petrov-Galerkin (DPG) variational framework. We demonstrate how test function spaces that guarantee numerical stability can be computed automatically with respect to the so called optimal test space norm by using an element subgrid discretization. This should make the DPG method not only stable but also robust, that is, uniformly stable with respect to the P'eclet number in the current application. The effectiveness of the algorithm is demonstrated on two problems for the linear advection-diffusion equation. © 2011 Published by Elsevier Ltd
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
Breaking spaces and forms for the DPG method and applications including Maxwell equations
Discontinuous Petrov Galerkin (DPG) methods are made easily implementable
using `broken' test spaces, i.e., spaces of functions with no continuity
constraints across mesh element interfaces. Broken spaces derivable from a
standard exact sequence of first order (unbroken) Sobolev spaces are of
particular interest. A characterization of interface spaces that connect the
broken spaces to their unbroken counterparts is provided. Stability of certain
formulations using the broken spaces can be derived from the stability of
analogues that use unbroken spaces. This technique is used to provide a
complete error analysis of DPG methods for Maxwell equations with perfect
electric boundary conditions. The technique also permits considerable
simplifications of previous analyses of DPG methods for other equations.
Reliability and efficiency estimates for an error indicator also follow.
Finally, the equivalence of stability for various formulations of the same
Maxwell problem is proved, including the strong form, the ultraweak form, and a
spectrum of forms in between
- …