Error analysis of a mixed finite element approximation of the semilinear Sobolev equations

Based on a mixed finite element method, we construct semidiscrete approximations of the solution u and the flux term ?u+?u _t of the semilinear Sobolev equations. The existence and uniqueness of the semidiscrete approximations are demonstrated and the error estimates of optimal rate in L ² normed space are derived. And also we construct the fully discrete approximations of u and ?u+?u _t and analyze the convergence of optimal rate in L ² normed space.     

A COMPLETE ANALYSIS FOR SOME A POSTERIORI ERROR ESTIMATES WITH THE FINITE ELEMENT METHOD OF LINES FOR A NONLINEAR PARABOLIC EQUATION

ABSTRACT

A posteriori error estimates for semidiscrete finite element methods for a nonlinear parabolic initial-boundary value problem are considered. The error estimates are obtained by solving local parabolic or elliptic equations for corrections to the solution on each element. The convergence results improve previous results where unnecessary assumptions are imposed on the approximate solution and the elliptic projection of the exact solution.     

Error estimates for mixed finite element approximations to a linear Stefan problem

We introduce a semidiscrete mixed finite element approximation for the single-phase linear Stefan problem and show the unique existence of the approximation. And the optimal rate of convergence inL ² andH ¹ norms are derived.     

Finite element approximation to global stabilization of the Burgers’ equation by Neumann boundary feedback control law

In this article, we discuss global stabilization results for the Burgers’ equation using nonlinear Neumann boundary feedback control law. As a result of the nonlinear feedback control, a typical nonlinear problem is derived. Then, based on C ⁰-conforming finite element method, global stabilization results for the semidiscrete solution are analyzed. Further, introducing an auxiliary projection, optimal error estimates in \(L^{\infty }(L^{2})\), \(L^{\infty }(H^{1})\) and \(L^{\infty }(L^{\infty })\)-norms for the state variable are obtained. Moreover, superconvergence results are established for the first time for the feedback control laws, which preserve exponential stabilization property. Finally, some numerical experiments are conducted to confirm our theoretical findings.     

L ∞-convergence of collocation and galerkin approximations to linear two-point parabolic problems

Two semidiscrete collocation approximations using smooth cubic splines are developed as approximations to the solution of two-point linear parabolic boundary value problems.L _∞-convergence results are presented for these two approximations as well as the piecewise linear Galerkin approximation. Several computational examples are given to illustrate the convergence results and demonstrate the applicability of the method.     

A discrete finite element galerkin method for a unidimensional single-phase Stefan problem

Based on Landau-type transformation, a Stefan problem with nonlinear free boundary condition is transformed into a system consisting of parabolic equation and the ordinary differential equations. Semidiscrete approximations are constructed. Optimal orders of convergence of semidiscrete approximation inL ₂,H ¹ andH ² normed spaces are derived.     

Error estimates of H 1-Galerkin mixed finite element method for Schrödinger equation

An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.     

A finite element method for the Sivashinsky equation

The Sivashinsky equation is a nonlinear evolutionary equation of fourth order in space. In this paper we have analyzed a semidiscrete finite element method and completely discrete scheme based on the backward Euler method and Crank–Nicolson–Galerkin scheme. A linearized backward Euler method have been developed and error bounds are derived for an L² projection.     

A second order finite difference-spectral method for space fractional diffusion equations

A high order finite difference-spectral method is derived for solving space fractional diffusion equations,by combining the second order finite difference method in time and the spectral Galerkin method in space.The stability and error estimates of the temporal semidiscrete scheme are rigorously discussed,and the convergence order of the proposed method is proved to be O(τ2+Nα-m)in L2-norm,whereτ,N,αand m are the time step size,polynomial degree,fractional derivative index and regularity of the exact solution,respectively.Numerical experiments are carried out to demonstrate the theoretical analysis.     

Convergence of FEM with interpolated coefficients for semilinear hyperbolic equation

To solve spatially semidiscrete approximative solution of a class of semilinear hyperbolic equations, the finite element method (FEM) with interpolated coefficients is discussed. By use of semidiscrete finite element for linear problem as comparative function, the error estimate in L^∞$L^{\infty }$-norm is derived by the nonlinear argument in Chen [Structure theory of superconvergence of finite elements, Hunan Press of Science and Technology, Changsha, 2001 (in Chinese)]. This indicates that convergence of FEMs with interpolated coefficients for a semilinear equation is similar to that of classical FEMs.     

The relaxed inexact projection methods for the split feasibility problem

In this paper we present several relaxed inexact projection methods for the split feasibility problem (SFP). Each iteration of the first proposed algorithm consists of a projection onto a halfspace containing the given closed convex set. The algorithm can be implemented easily and its global convergence to the solution can be established under suitable conditions. Moreover,we present some modifications of the relaxed inexact projection method with constant stepsize by adopting Armijo-like search. We furthermore present a variable-step relaxed inexact projection method which does not require the computation of the matrix inverses and the largest eigenvalue of the matrix A^TA, and the objective function can decrease sufficiently at each iteration. We show convergence of these modified algorithms under mild conditions. Finally, we perform some numerical experiments, which show the behavior of the algorithms proposed.     

Asymptotic behavior and finite element error estimates of Kelvin-Voigt viscoelastic fluid flow model

In this article, the convergence of the solution of the Kelvin-Voigt viscoelastic fluid flow model to its steady state solution with exponential rate is established under the uniqueness assumption. Then, a semidiscrete Galerkin method for spatial direction keeping time variable continuous is considered and asymptotic behavior of the semidiscrete solution is derived. Moreover, optimal error estimates are achieved for large time using steady state error estimates. Based on linearized backward Euler method, asymptotic behavior for the fully discrete solution is studied and optimal error estimates are derived for large time. All the results are even valid for κ→0, that is, when the Kelvin-Voigt model converges to the Navier-Stokes system. Finally, some numerical experiments are conducted to confirm our theoretical findings.     

Local existence and non-relativistic limits of shock solutions to a multidimensional piston problem for the relativistic Euler equations

The multidimensional piston problem is a special initial-boundary value problem. The boundary conditions are given in two conical surfaces: one is the boundary of the piston, and the other is the shock whose location is to be determined later. In this paper, we are concerned with spherically symmetric piston problem for the relativistic Euler equations. A local shock front solution with the state equation p = a ² ρ, a is a constant and has been established by the Newton iteration. To overcome the difficulty caused by the free boundary, we introduce a coordinate transformation to fix it and employ the linear iteration scheme to establish a sequence of approximate solutions to the auxiliary problems by iteration. In each step, the value of the solution of the previous problem is taken as the data to determine the solution of the next problem. We obtain the existence of the original problem by establishing the convergence of these sequences. Meanwhile, we establish the convergence of the local solution as c → ∞ to the corresponding solution of the classical non-relativistic Euler equations.     

Contrôlabilité de l'équation des ondes à densité rapidement oscillante à une dimension d'espace

We consider the one-dimensional wave equation with periodic density of period ε → 0. By a counterexample due to Avellaneda, Bardos, and Rauch, we know that the boundary controllability property does not hold uniformly as ε → 0. We prove that the control remains uniformly bounded if we control the projection of the solution over the subspace generated by the eigenfunctions associated with the eigenvalues λ ≤ C_ε⁻², C > 0 being small enough. This result is sharp in the sense that the control diverges when the projection over the eigenfunctions such that λ ~ C_ε⁻², with C large, is controlled. We use the classical WKB asymptotic development that provides sharp results on the convergence of the spectrum and the theory of non-harmonic Fourier series.     

Error estimates for a finite volume element method for parabolic equations in convex polygonal domains

We analyze the spatially semidiscrete piecewise linear finite volume element method for parabolic equations in a convex polygonal domain in the plane. Our approach is based on the properties of the standard finite element Ritz projection and also of the elliptic projection defined by the bilinear form associated with the variational formulation of the finite volume element method. Because the domain is polygonal, special attention has to be paid to the limited regularity of the exact solution. We give sufficient conditions in terms of data that yield optimal order error estimates in L₂ and H 1 . The convergence rate in the L_∞ norm is suboptimal, the same as in the corresponding finite element method, and almost optimal away from the corners. We also briefly consider the lumped mass modification and the backward Euler fully discrete method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

A characteristics-mixed covolume method for a convection-dominated transport problem

In this paper, we propose a characteristics-mixed covolume method for approximating the solution to a convection dominated transport problem. The method is a combination of characteristic approximation to handle the convection term in time and mixed covolume method spatial approximation to deal with the diffusion term. The velocity and press are approximated by the lowest order Raviart-Thomas mixed finite element space on rectangles. The projection of a mixed covolume element is introduced. We prove its first order optimal rate of convergence for the approximate velocities in the L² norm as well as for the approximate pressures in the L² norm.     

Galerkin Approximations for the Linear Parabolic Equation with the Third Boundary Condition

We solve a linear parabolic equation in ^d , d 1, with the third nonhomogeneous boundary condition using the finite element method for discretization in space, and the -method for discretization in time. The convergence of both, the semidiscrete approximations and the fully discretized ones, is analysed. The proofs are based on a generalization of the idea of the elliptic projection. The rate of convergence is derived also for variable time step-sizes.     

A priori error estimates of finite volume element method for hyperbolic optimal control problems

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.     

An Augmented Lagrangian Method for a Class of Inverse Quadratic Programming Problems

We consider an inverse quadratic programming (QP) problem in which the parameters in the objective function of a given QP problem are adjusted as little as possible so that a known feasible solution becomes the optimal one. We formulate this problem as a minimization problem with a positive semidefinite cone constraint and its dual is a linearly constrained semismoothly differentiable (SC¹) convex programming problem with fewer variables than the original one. We demonstrate the global convergence of the augmented Lagrangian method for the dual problem and prove that the convergence rate of primal iterates, generated by the augmented Lagrange method, is proportional to 1/r, and the rate of multiplier iterates is proportional to  $1/\sqrt{r}$ , where r is the penalty parameter in the augmented Lagrangian. As the objective function of the dual problem is a SC¹ function involving the projection operator onto the cone of symmetrically semi-definite matrices, the analysis requires extensive tools such as the singular value decomposition of matrices, an implicit function theorem for semismooth functions, and properties of the projection operator in the symmetric-matrix space. Furthermore, the semismooth Newton method with Armijo line search is applied to solve the subproblems in the augmented Lagrange approach, which is proven to have global convergence and local quadratic rate. Finally numerical results, implemented by the augmented Lagrangian method, are reported.     

Higher order fully discrete scheme combined withH 1-Galerkin mixed finite element method for semilinear reaction-diffusion equations

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by anH ¹-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) ofindex one. Apriori error estimates for semidiscrete scheme are derived for both differential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.