Galerkin approximations to non-linear pseudo-parabolic partial differential equations. (Short Communication).
Proper orthogonal decomposition (POD) is a powerful technique for model reduction of linear and non-linear systems. It is based on a Galerkin type discretization with basis elements created from the system itself. In this work, error estimates for Galerkin POD methods for linear elliptic, parameter-dependent systems are proved. The resulting error bounds depend on the number of POD basis functions and on the parameter grid that is used to generate the snapshots and to compute the POD basis. The...
The famous Zlámal’s minimum angle condition has been widely used for construction of a regular family of triangulations (containing nondegenerating triangles) as well as in convergence proofs for the finite element method in . In this paper we present and discuss its generalization to simplicial partitions in any space dimension.
This paper is concerned with the generalization of the finite element method via the use of non-polynomial enrichment functions. Several methods employ this general approach, e.g. the extended finite element method and the generalized finite element method. We review these approaches and interpret them in the more general framework of the partition of unity method. Here we focus on fundamental construction principles, approximation properties and stability of the respective numerical method. To...
In this paper we describe PDELab, an extensible C++ template library for finite element methods based on the Distributed and Unified Numerics Environment (Dune). PDELab considerably simplifies the implementation of discretization schemes for systems of partial differential equations by setting up global functions and operators from a simple element-local description. A general concept for incorporation of constraints eases the implementation of essential boundary conditions, hanging nodes and varying...
A general approach was proposed in this article to develop high-order exponentially fitted basis functions for finite element approximations of multi-dimensional drift-diffusion equations for modeling biomolecular electrodiffusion processes. Such methods are highly desirable for achieving numerical stability and efficiency. We found that by utilizing the one-to-one correspondence between the continuous piecewise polynomial space of degree k + 1 and the divergencefree vector space of degree k, one...
We deal with a posteriori error control of discontinuous Galerkin approximations for linear boundary value problems. The computational error is estimated in the framework of the Dual Weighted Residual method (DWR) for goal-oriented error estimation which requires to solve an additional (adjoint) problem. We focus on the control of the algebraic errors arising from iterative solutions of algebraic systems corresponding to both the primal and adjoint problems. Moreover, we present two different reconstruction...
The paper presents the solution to the geodetic boundary value problem by the finite element method in area of Slovak Republic. Generally, we have made two numerical experiments. In the first one, Neumann BC in the form of gravity disturbances generated from EGM-96 is used and the solution is verified by the quasigeoidal heights generated directly from EGM-96. In the second one, Neumann BC is computed from gravity measurements and the solution is compared to the quasigeoidal heights obtained by...
This contribution presents a general numerical method for computing lower and upper bound of the optimal constant in Friedrichs’ inequality. The standard Rayleigh-Ritz method is used for the lower bound and the method of is employed for the upper bound. Several numerical experiments show applicability and accuracy of this approach.
We derive a posteriori error estimates for singularly perturbed reaction–diffusion problems which yield a guaranteed upper bound on the discretization error and are fully and easily computable. Moreover, they are also locally efficient and robust in the sense that they represent local lower bounds for the actual error, up to a generic constant independent in particular of the reaction coefficient. We present our results in the framework of the vertex-centered finite volume method but their nature...
We deal with the a-posteriori estimation of the error for finite element solutions of nonstationary heat conduction problems with mixed boundary conditions on bounded polygonal domains. The a-posteriori error estimates are constucted by solving stationary “reconstruction” problems, obtained by replacing the time derivative of the exact solution by the time derivative of the finite element solution. The main result is that the reconstructed solutions, or reconstructions, are superconvergent approximations...
A method of characterizing all eigenvalues of a preconditioned discretized scalar diffusion operator with Dirichlet boundary conditions has been recently introduced in Gergelits, Mardal, Nielsen, and Strakoš (2019). Motivated by this paper, we offer a slightly different approach that extends the previous results in some directions. Namely, we provide bounds on all increasingly ordered eigenvalues of a general diffusion or elasticity operator with tensor data, discretized with the conforming finite...