This paper presents a review of the complementary technique with the emphasis on computable and guaranteed upper bounds of the approximation error. For simplicity, the approach is described on a numerical solution of the Poisson problem. We derive the complementary error bounds, prove their fundamental properties, present the method of hypercircle, mention possible generalizations and show a couple of numerical examples.
The paper deals with fast solving of large saddle-point systems arising in wavelet-Galerkin discretizations of separable elliptic PDEs. The periodized orthonormal compactly supported wavelets of the tensor product type together with the fictitious domain method are used. A special structure of matrices makes it possible to utilize the fast Fourier transform that determines the complexity of the algorithm. Numerical experiments confirm theoretical results.
The general method of averaging for the superapproximation of an arbitrary partial derivative of a smooth function in a vertex of a simplicial triangulation of a bounded polytopic domain in for any is described and its complexity is analysed.
This paper concerns the composite grid finite element (FE) method for solving boundary value problems in the cases which require local grid refinement for enhancing the approximating properties of the corresponding FE space. A special interest is given to iterative methods based on natural decomposition of the space of unknowns and to the implementation of both the composite grid FEM and the iterative procedures for its solution. The implementation is important for gaining all benefits of the described...
This paper is concerned with the computation of 3D vertex singularities of anisotropic elastic fields with Dirichlet boundary conditions, focusing on the derivation of error estimates for a finite element method on graded meshes. The singularities are described by eigenpairs of a corresponding operator pencil on spherical polygonal domains. The main idea is to introduce a modified quadratic variational boundary eigenvalue problem which consists of two self-adjoint, positive definite sesquilinear...
This paper is concerned with the computation of 3D vertex singularities of anisotropic elastic fields with Dirichlet boundary conditions, focusing on the derivation of error estimates for a finite element method on graded meshes. The singularities are described by eigenpairs of a corresponding operator pencil on spherical polygonal domains. The main idea is to introduce a modified quadratic variational boundary eigenvalue problem which consists of two self-adjoint, positive definite sesquilinear...
A new method for computation of the fundamental solution of electrodynamics for general anisotropic nondispersive materials is suggested. It consists of several steps: equations for each column of the fundamental matrix are reduced to a symmetric hyperbolic system; using the Fourier transform with respect to space variables and matrix transformations, formulae for Fourier images of the fundamental matrix columns are obtained; finally, the fundamental solution is computed by the inverse Fourier transform....
The paper presents the results of numerical solution of the evolution law for the constrained mean-curvature flow. This law originates in the theory of phase transitions for crystalline materials and describes the evolution of closed embedded curves with constant enclosed area. It is reformulated by means of the direct method into the system of degenerate parabolic partial differential equations for the curve parametrization. This system is solved numerically and several computational studies are...
The paper presents the results of numerical solution of the Allen-Cahn equation with a non-local term. This equation originally mentioned by Rubinstein and Sternberg in 1992 is related to the mean-curvature flow with the constraint of constant volume enclosed by the evolving curve. We study this motion approximately by the mentioned PDE, generalize the problem by including anisotropy and discuss the computational results obtained.
We present a finite element method to compute guided modes in a stratified medium. The major difficulty to overcome is related to the unboundedness of the stratified medium. Our method is an alternative to the use of artificial boundary conditions and to the use of integral representation formulae. The domain is bounded in such a way we can write the solution on its lateral boundaries in terms of Fourier series. The series is then truncated for the computations over the bounded domain. The problem...
We present a finite element method to compute guided modes in a stratified medium. The major difficulty to overcome is related to the unboundedness of the stratified medium. Our method is an alternative to the use of artificial boundary conditions and to the use of integral representation formulae. The domain is bounded in such a way we can write the solution on its lateral boundaries in terms of Fourier series. The series is then truncated for the computations over the bounded domain. The problem...