A multigrid method for Reissner-Mindlin plates.
The author studies the behaviour of a multi-level method that combines the Jacobi iterations and the correction by aggragation of unknowns. Our considerations are restricted to a simple one-dimensional example, which allows us to employ the technique of the Fourier analysis. Despite of this restriction we are able to demonstrate differences between the behaviour of the algorithm considered and of multigrid methods employing interpolation instead of aggregation.
We propose a new type of multilevel method for solving eigenvalue problems based on Newton's method. With the proposed iteration method, solving an eigenvalue problem on the finest finite element space is replaced by solving a small scale eigenvalue problem in a coarse space and a sequence of augmented linear problems, derived by Newton step in the corresponding sequence of finite element spaces. This iteration scheme improves overall efficiency of the finite element method for solving eigenvalue...
A multilevel preconditioner based on the abstract framework of the auxiliary space method, is developed for the mortar method for the nonconforming P1 finite element or the lowest order Crouzeix-Raviart finite element on nonmatching grids. It is shown that the proposed preconditioner is quasi-optimal in the sense that the condition number of the preconditioned system is independent of the mesh size, and depends only quadratically on the number of refinement levels. Some numerical results confirming...
In this work, we consider singular perturbations of the boundary of a smooth domain. We describe the asymptotic behavior of the solution uE of a second order elliptic equation posed in the perturbed domain with respect to the size parameter ε of the deformation. We are also interested in the variations of the energy functional. We propose a numerical method for the approximation of uE based on a multiscale superposition of the unperturbed solution u0 and a profile defined in a model domain. We...
In this paper, multiscale finite element methods (MsFEMs) and domain decomposition techniques are developed for a class of nonlinear elliptic problems with high-contrast coefficients. In the process, existing work on linear problems [Y. Efendiev, J. Galvis, R. Lazarov, S. Margenov and J. Ren, Robust two-level domain decomposition preconditioners for high-contrast anisotropic flows in multiscale media. Submitted.; Y. Efendiev, J. Galvis and X. Wu, J. Comput. Phys. 230 (2011) 937–955; J. Galvis and...
We propose a multiscale model reduction method for partial differential equations. The main purpose of this method is to derive an effective equation for multiscale problems without scale separation. An essential ingredient of our method is to decompose the harmonic coordinates into a smooth part and a highly oscillatory part so that the smooth part is invertible and the highly oscillatory part is small. Such a decomposition plays a key role in our construction of the effective equation. We show...
In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite element method that reduces to cell-centered finite differences on irregular grids. The subdomain grids do not have to match across the interfaces. Continuity of flux between coarse elements is imposed via a mortar finite element space on a coarse grid...
In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite element method that reduces to cell-centered finite differences on irregular grids. The subdomain grids do not have to match across the interfaces. Continuity of flux between coarse elements is imposed via a mortar finite element space on a coarse grid...
In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite element method that reduces to cell-centered finite differences on irregular grids. The subdomain grids do not have to match across the interfaces. Continuity of flux between coarse elements is imposed via a mortar finite element space on a coarse grid...
We consider the numerical approximation of a first order stationary hyperbolic equation by the method of characteristics with pseudo time step k using discontinuous finite elements on a mesh . For this method, we exhibit a “natural” norm || ||h,k for which we show that the discrete variational problem is well posed and we obtain an error estimate. We show that when k goes to zero problem (resp. the || ||h,k norm) has as a limit problem (Ph) (resp. the || ||h norm) associated to the...
A new error correction method for the stationary Navier-Stokes equations based on two local Gauss integrations is presented. Applying the orthogonal projection technique, we introduce two local Gauss integrations as a stabilizing term in the error correction method, and derive a new error correction method. In both the coarse solution computation step and the error computation step, a locally stabilizing term based on two local Gauss integrations is introduced. The stability and convergence of the...
Finite element methods with piecewise polynomial spaces in space for solving the nonstationary heat equation, as a model for parabolic equations are considered. The discretization in time is performed using the Crank-Nicolson method. A new a priori estimate is proved. Thanks to this new a priori estimate, a new error estimate in the discrete norm of is proved. An -error estimate is also shown. These error estimates are useful since they allow us to get second order time accurate approximations...
In this paper we present a novel exponentially fitted finite element method with triangular elements for the decoupled continuity equations in the drift-diffusion model of semiconductor devices. The continuous problem is first formulated as a variational problem using a weighted inner product. A Bubnov-Galerkin finite element method with a set of piecewise exponential basis functions is then proposed. The method is shown to be stable and can be regarded as an extension to two dimensions of the...