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...
This paper addresses the derivation of new second-kind Fredholm combined field integral equations for the Krylov iterative solution of tridimensional acoustic scattering problems by a smooth closed surface. These integral equations need the introduction of suitable tangential square-root operators to regularize the formulations. Existence and uniqueness occur for these formulations. They can be interpreted as generalizations of the well-known Brakhage-Werner [A. Brakhage and P. Werner, Arch....
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 purpose of our work is to develop an automatic shape optimization tool for runner wheel blades in reaction water turbines, especially in Kaplan turbines. The fluid flow is simulated using an in-house incompressible turbulent flow solver based on recently introduced isogeometric analysis (see e.g. J. A. Cotrell et al.: Isogeometric Analysis: Toward Integration of CAD and FEA, Wiley, 2009). The proposed automatic shape optimization approach is based on a so-called hybrid optimization which combines...
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...
Making use of a line integral defined without use of the partition of unity, Green’s theorem is proved in the case of two-dimensional domains with a Lipschitz-continuous boundary for functions belonging to the Sobolev spaces
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...