Matrix polynomials play an important role in the theory of matrix differential equations. We develop a fixed point method to compute solutions of matrix polynomials equations, where the matricial elements of the matrix polynomial are considered separately as complex polynomials. Numerical examples illustrate the method presented.
A standard method for proving the inf-sup condition implying stability of finite element approximations for the stationary Stokes equations is to construct a Fortin operator. In this paper, we show how this can be done for two-dimensional triangular and rectangular Taylor-Hood methods, which use continuous piecewise polynomial approximations for both velocity and pressure.
We study a two-grid scheme fully discrete in time and space for solving the Navier-Stokes system. In the first step, the fully non-linear problem is discretized in space on a coarse grid with mesh-size H and time step k. In the second step, the problem is discretized in space on a fine grid with mesh-size h and the same time step, and linearized around the velocity uH computed in the first step. The two-grid strategy is motivated by the fact that under suitable assumptions, the contribution of uH...
A full multigrid finite element method is proposed for semilinear elliptic equations. The main idea is to transform the solution of the semilinear problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and semilinear problems on a very low dimensional space. The linearized boundary value problems are solved by some multigrid iterations. Besides the multigrid iteration, all other efficient numerical methods can also serve as...