We deal with the goal-oriented error estimates and mesh adaptation for nonlinear partial differential equations. The setting of the adjoint problem and the resulting estimates are not based on a differentiation of the primal problem but on a suitable linearization which guarantees the adjoint consistency of the numerical scheme. Furthermore, we develop an efficient adaptive algorithm which balances the errors arising from the discretization and the use of nonlinear as well as linear iterative solvers. Several numerical examples demonstrate the efficiency of this algorithm.

This paper is concerned with the analysis of the finite element method for the numerical solution of an elliptic boundary value problem with a nonlinear Newton boundary condition in a two-dimensional polygonal domain. The weak solution loses regularity in a neighbourhood of boundary singularities, which may be at corners or at roots of the weak solution on edges. The main attention is paid to the study of error estimates. It turns out that the order of convergence is not dampened by the nonlinearity if the weak solution is nonzero on a large part of the boundary. If the weak solution is zero on the whole boundary, the nonlinearity only slows down the convergence of the function values but not the convergence of the gradient. The same analysis is carried out for approximate solutions obtained by numerical integration. The theoretical results are verified by numerical experiments.

We deal with the numerical solution of a linear convection-diffusion-reaction equation using the discontinuous Galerkin method of arbitrary polynomial approximation degree on anisotropic triangular grids. We derive a posteriori goal-oriented error estimates taking into account the anisotropy of mesh elements. The resulting error estimates are employed for the construction of an anisotropic mesh adaptation algorithm which locally optimizes the size and shape of mesh elements. The computational performance is demonstrated by several numerical experiments.

The paper is concerned with the numerical analysis of an elliptic equation in a polygon with a nonlinear Newton boundary condition, discretized by the finite element or discontinuous Galerkin methods. Using the monotone operator theory, it is possible to prove the existence and uniqueness of the exact weak solution and the approximate solution. The main attention is paid to the study of error estimates. To this end, the regularity of the weak solution is investigated and it is shown that due to the boundary corner points, the solution looses regularity in a vicinity of these points. It comes out that the error estimation depends essentially on the opening angle of the corner points and on the parameter defining the nonlinear behaviour of the Newton boundary condition. Theoretical results are compared with numerical experiments confirming a nonstandard behaviour of error estimates.