Improved Difference Schemes for the Dirichlet Problem of Poisson's Equation in the Neighbourhood of Corners.
We present families of scalar nonconforming finite elements of arbitrary order with optimal approximation properties on quadrilaterals and hexahedra. Their vector-valued versions together with a discontinuous pressure approximation of order form inf-sup stable finite element pairs of order r for the Stokes problem. The well-known elements by Rannacher and Turek are recovered in the case r=1. A numerical comparison between conforming and nonconforming discretisations will be given. Since higher order...
In this paper, we propose a new numerical method for solving elliptic equations in unbounded regions of . The method is based on the mapping of a part of the domain into a bounded region. An appropriate family of weighted spaces is used for describing the growth or the decay of functions at large distances. After exposing the main ideas of the method, we analyse carefully its convergence. Some 3D computational results are displayed to demonstrate its efficiency and its high performance.
In this paper, we propose a new numerical method for solving elliptic equations in unbounded regions of . The method is based on the mapping of a part of the domain into a bounded region. An appropriate family of weighted spaces is used for describing the growth or the decay of functions at large distances. After exposing the main ideas of the method, we analyse carefully its convergence. Some 3D computational results are displayed to demonstrate its efficiency and its high performance.
We consider high order finite difference approximations to the Helmholtz equation in an exterior domain. We include a simplified absorbing boundary condition to approximate the Sommerfeld radiation condition. This yields a large, but sparse, complex system, which is not self-adjoint and not positive definite. We discretize the equation with a compact fourth or sixth order accurate scheme. We solve this large system of linear equations with a Krylov subspace iterative method. Since the method converges...
We consider high order finite difference approximations to the Helmholtz equation in an exterior domain. We include a simplified absorbing boundary condition to approximate the Sommerfeld radiation condition. This yields a large, but sparse, complex system, which is not self-adjoint and not positive definite. We discretize the equation with a compact fourth or sixth order accurate scheme. We solve this large system of linear equations with a Krylov subspace iterative method. Since the method converges...
The computation of nonlinear quasistationary two-dimensional magnetic fields leads to a nonlinear second order parabolic-elliptic initial-boundary value problem. Such a problem with a nonhomogeneous Dirichlet boundary condition on a part of the boundary is studied in this paper. The problem is discretized in space by the finite element method with linear functions on triangular elements and in time by the implicit-explicit method (the left-hand side by the implicit Euler method and the right-hand...
We developed a mimetic finite difference method for solving elliptic equations with tensor coefficients on polyhedral meshes. The first-order convergence estimates in a mesh-dependent norm are derived.
We developed a mimetic finite difference method for solving elliptic equations with tensor coefficients on polyhedral meshes. The first-order convergence estimates in a mesh-dependent H1 norm are derived.
We state a superconvergence result for the lowest order Raviart-Thomas approximation of eigenvalue problems. It is known that a similar superconvergence result holds for the mixed approximation of Laplace problem; here we introduce a new proof, since the one given for the source problem cannot be generalized in a straightforward way to the eigenvalue problem. Numerical experiments confirm the superconvergence property and suggest that it also holds for the lowest order Brezzi-Douglas-Marini...
This contribution deals with a mixed variational formulation of 3D contact problems with the simplest model involving friction. This formulation is based on a dualization of the set of admissible displacements and the regularization of the non-differentiable term. Displacements are approximated by piecewise linear elements while the respective dual variables by piecewise constant functions on a dual partition of the contact zone. The rate of convergence is established provided that the solution...
This contribution deals with a mixed variational formulation of 3D contact problems with the simplest model involving friction. This formulation is based on a dualization of the set of admissible displacements and the regularization of the non-differentiable term. Displacements are approximated by piecewise linear elements while the respective dual variables by piecewise constant functions on a dual partition of the contact zone. The rate of convergence is established provided that the solution...