In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.
In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discrete optimal system is obtained.We prove the existence and uniqueness of the solution to the semidiscrete optimal system and obtain the optimal order error estimates in L ∞(J;L 2)-and L ∞(J;H 1)-norm.Numerical experiments are presented to test these theoretical results.
A cascadic multigrid method is proposed for eigenvalue problems based on the multilevel correction scheme. With this new scheme, an eigenvalue problem on the finest space can be solved by linear smoothing steps on a series of multilevel finite element spaces and nonlinear correcting steps on special coarsest spaces. Once the sequence of finite element spaces and the number of smoothing steps are appropriately chosen, the optimal convergence rate with the optimal computational work can be obtained. Some numerical experiments are presented to validate our theoretical analysis.
In this paper,we develop a correction operator for the canonical interpolation operator of the Adini element.We use this new correction operator to analyze the discrete eigenvalues of the Adini element method for the fourth order elliptic eigenvalue problem in the three dimensions.We prove that the discrete eigenvalues are smaller than the exact ones.
In this paper,we discuss the a posteriori error estimates of the mixed finite element method for quadratic optimal control problems governed by linear parabolic equations.The state and the co-state are discretized by the high order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions.We derive a posteriori error estimates for both the state and the control approximation.Such estimates,which are apparently not available in the literature,are an important step towards developing reliable adaptive mixed finite element approximation schemes for the control problem.
The ordered patterns formed bymicrophase-separated block copolymer systems demonstrate periodic symmetry,and all periodic structures belong to one of 230 space groups.Based on this fact,a strategy of estimating the initial values of self-consistent field theory to discover ordered patterns of block copolymers is developed.In particular,the initial period of the computational box is estimated by the Landau-Brazovskii model as well.By planting the strategy into the whole-space discrete method,several new metastable patterns are discovered in diblock copolymers.
We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations.Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We prove that the computational work of this new scheme is truly optimal,the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.
In this paper,the Crank-Nicolson linear finite volume element method is applied to solve the distributed optimal control problems governed by a parabolic equation.The optimal convergent order O(h^(2)+k^(2))is obtained for the numerical solution in a discrete L^(2)-norm.A numerical experiment is presented to test the theoretical result.
