A Posteriori Error Estimate of Two Grid Mixed Finite Element Methods for Semilinear Elliptic Equations

In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.

Hermite Finite Element Method for Vibration Problem of Euler-Bernoulli Beam on Viscoelastic Pasternak Foundation

Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.

A Priori and A Posteriori Error Estimates of Streamline Diffusion Finite Element Method for Optimal Control Problem Governed by Convection Dominated Diffusion Equation

In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.

Error estimates of H^1-Galerkin mixed finite element method for Schrdinger 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.

Adaptive Finite Element Method Based on Optimal Error Estimates for Linear Elliptic Problems on Nonconvex Polygonal Domains

The subject of this work is to propose adaptive finite element methods based on an optimal maximum norm error control estimate.Using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions,the proposed procedures are analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that they generate the correct type of refinement and lead to the desired control under consideration.

ADAPTIVE FINITE ELEMENT METHOD BASED ON OPTIMAL ERROR ESTIMATES FOR LINEAR ELLIPTIC PROBLEMS ON CONCAVE CORNER DOMAINS

An adaptive finite element procedure designed for specific computational goals is presented,using mesh refinement strategies based on optimal or nearly optimal a priori error estimates for the finite element method and using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions.The proposed procedure is analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that the method can generate the correct type of refinements and lead to the desired control under consideration.

A Posteriori Error Estimates for Finite Element Methods for Systems of Nonlinear,Dispersive Equations

The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.

A lumped mass nonconforming finite element method for nonlinear parabolic integro-differential equations on anisotropic meshes

A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.

COARSE-MESH-ACCURACY IMPROVEMENT OF BILINEAR Q_4-PLANE ELEMENT BY THE COMBINED HYBRID FINITE ELEMENT METHOD

The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads to enhancement of accuracy at coarse meshes, and that the combination parameter plays an important role in the enhancement. As an improvement of conforming bilinear Q(4)-plane element, the combined hybrid method adopted the most convenient quadrilateral displacements-stress mode, i.e.,the mode of compatible isoparametric bilinear displacements and pure constant stresses. By adjusting the combined parameter, the optimized version of the combined hybrid element was obtained and numerical tests indicated that this parameter-adjusted version behaves much better than Q(4)-element and is of high accuracy at coarse meshes. Due to elimination of stress parameters at the elemental level, this combined hybrid version is of the same computational cost as that of Q(4)-element.

Goal-oriented error estimation applied to direct solution of steady-state analysis with frequency-domain finite element method

Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.

THE SPACE-TIME FINITE ELEMENT METHOD FOR PARABOLIC PROBLEMS

Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.

P_1-nonconforming triangular finite element method for elliptic and parabolic interface problems

The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.

A streamline diffusion nonconforming finite element method for the time-dependent linearized Navier-Stokes equations

A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method （FEM） based on a proper orthogo- nal decomposition （POD） theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save mem- ory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element （FE） scheme is shown to be un- conditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations （FDEs）.

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION (Ⅱ)——THE DISCRETE-TIME CASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTA-TION (Ⅰ)-THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations,silt transport equation, the equation of bottom topography change,and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element(MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived.The error estimates are optimal.

LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR A CLASS OF STOKES EQUATION

A least-squares mixed finite element method was formulated for a class of Stokes equations in two dimensional domains. The steady state and the time-dependent Stokes' equations were considered. For the stationary equation, optimal H-t and L-2-error estimates are derived under the standard regularity assumption on the finite element partition ( the LBB-condition is not required). Far the evolutionary equation, optimal L-2 estimates are derived under the conventional Raviart-Thomas spaces.

DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR A FORWARD-BACKWARD HEAT EQUATION

A space-time finite element method,discontinuous in time but continuous in space,is studied to solve the nonlinear forward-backward heat equation.A linearized technique is introduced in order to obtain the error estimates of the approximate solutions.And the numerical simulations are given.

Anisotropic adaptive finite element method for magnetohydrodynamic flow at high Hartmann numbers

This paper presents an anisotropic adaptive finite element method （FEM） to solve the governing equations of steady magnetohydrodynamic （MHD） duct flow. A resid- ual error estimator is presented for the standard FEM, and two-sided bounds on the error independent of the aspect ratio of meshes are provided. Based on the Zienkiewicz-Zhu es- timates, a computable anisotropic error indicator and an implement anisotropic adaptive refinement for the MHD problem are derived at different values of the Hartmann number. The most distinguishing feature of the method is that the layer information from some directions is captured well such that the number of mesh vertices is dramatically reduced for a given level of accuracy. Thus, this approach is more suitable for approximating the layer problem at high Hartmann numbers. Numerical results show efficiency of the algorithm.

Adaptive Finite Element Method for Steady Convection-Diffusion Equation

This paper examines the numerical solution of the convection-diffusion equation in 2-D. The solution of this equation possesses singularities in the form of boundary or interior layers due to non-smooth boundary conditions. To overcome such singularities arising from these critical regions, the adaptive finite element method is employed. This scheme is based on the streamline diffusion method combined with Neumann-type posteriori estimator. The effectiveness of this approach is illustrated by different examples with several numerical experiments.