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.

Error Estimates for Mixed Finite Element Methods for Sobolev Equation

The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).

Two new least-squares mixed finite element procedures for convection-dominated Sobolev equations

Two new convection-dominated are derived under the approximate solutions least-squares mixed finite element procedures are formulated for solving Sobolev equations. Optimal H（div;Ω）×H1（Ω） norms error estimates standard mixed finite spaces. Moreover, these two schemes provide the with first-order and second-order accuracy in time increment, respectively.

ANALYSIS OF AUGMENTED THREE-FIELD MACRO-HYBRID MIXED FINITE ELEMENT SCHEMES

On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, ＂mass＂ preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.

Mixed Finite Element Formats of any Order Based on Bubble Functions for Stationary Stokes Problem

Mixed element formats of any order based on bubble functions for the stationary Stokes problem are derived in triangular and tetrahedral meshes and the convergence of these formats are proved.

IMPLICIT-EXPLICIT MULTISTEP FINITE ELEMENT-MIXED FINITE ELEMENT METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equation is discretized by a mixed finite element method. The electron and hole density equations are treated by implicit-explicit multistep finite element methods. The schemes are very efficient. The optimal order error estimates both in time and space are derived.

Nonlinear simulation of arch dam cracking with mixed finite element method

This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.

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 METHOD OF HEXAHEDRAL ELEMENTS FOR NAVIER-STOKES PROBLEM

In this paper, we derive a new mixed element format of hexahedral elements for Navier-Stokes problem in three-dimensional space.

Discrete formulation of mixed finite element methods for vapor deposition chemical reaction equations

The vapor deposition chemical reaction processes, which are of extremely extensive applications, can be classified as a mathematical model by the following governing nonlinear partial differential equations containing velocity vector, temperature field, pressure field, and gas mass field. The mixed finite element （MFE） method is employed to study the system of equations for the vapor deposition chemical reaction processes. The semidiscrete and fully discrete MFE formulations are derived. And the existence and convergence （error estimate） of the semidiscrete and fully discrete MFE solutions are demonstrated. By employing MFE method to treat the system of equations for the vapor deposition chemical reaction processes, the numerical solutions of the velocity vector, the temperature field, the pressure field, and the gas mass field can be found out simultaneously. Thus, these researches are not only of important theoretical means, but also of extremely extensive applied vistas.

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.

MIXED FINITE ELEMENTS FOR PARABOLICINTEGRO-DIFFERENTIAL EQUATIONS

In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discrete and full-discrete mixed finite elements for the equations, and obtain the optimal L-2 error estimates.

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.

A STABILIZED MIXED FINITE ELEMENT FORMULATION FOR THE NON-STATIONARY INCOMPRESSIBLE BOUSSINESQ EQUATIONS

In this study, we employ mixed finite element （MFE） method, two local Gauss integrals, and parameter-free to establish a stabilized MFE formulation for the non-stationary incompressible Boussinesq equations. We also provide the theoretical analysis of the existence, uniqueness, stability, and convergence of the stabilized MFE solutions for the stabilized MFE formulation.

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper, a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure. The convergence analysis is presented and optimal error estimates of both broken H^1-norm and L^2-norm for velocity as well as the L^2-norm for the pressure are derived.

Direct Assessment of Interlaminar Stresses in Composite Multilayered Plates Using a Layer-Wise Mixed Finite Element Model

An accurate determination of intedaminar transversal stresses in composite multilayered plates, especially near free-edge, is of great importance in the study of inter-ply damage modes, mainly in the initiation and growth of delamination. In this paper, interlaminar stresses are determined by layer-wise mixed finite element model. Each layer is analyzed as an isolated one where the displacement continuity is ensured by means of Lagrange multipliers （which represent the statics variables）. This procedure allows the authors to work with any single plate model, obtaining the interlaminar stresses directly without loss of precision. The FSDT （first shear deformation theory） with transverse normal strain effects included is assumed in each layer, but Lagrange polynomials are used to describe the kinematic instead of Taylor＇s polynomial functions of the thickness coordinates, as is common. This expansion allows the authors to pose the interlaminar displacements compatibility simpler than the second one. The in-plane domain of the plate is discretized by four-node quadrilateral elements, both to the field of displacement and to the Lagrange multipliers. The mixed interpolation of tensorial components technique is applied to avoid the shear-locking in the finite element model. Several examples were carried out and the results have been satisfactorily compared with those available in the literature.

Determination of Natural Frequencies of Multilayered Plates by Mixed Finite Elements

Natural frequencies for multilayer plates are calculated by mixed finite element method. The main object of this paper is to use the mixed model for multilayer plates, analyzing each layer as an isolated plate, where the continuity of displacements is achieved by Lagrange multipliers （representing static variables）. This procedure allows us to work with any model for single plate （so as to ensure the proper behavior of each layer）, and the complexity of the multilayer system is avoided by ensuring the condition of displacements by the Lagrange multipliers （static variables）. The plate is discretized by finite element modeling based on a primary hybrid model, where the domain is divided by quadrilateral, both for the displacement field and static variables. This mixed element for plates was implemented and several examples of vibrations have been verified successfully by the results obtained by other methods in the literature.

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.

A MIXED FINITE ELEMENT APPROXIMATIONFOR COMPRESSIBLE FLOW OF CONTAMINATION FROM NUCLEAR WASTE IN POROUS MEDIA

A compressible nuclear waste disposal contamination in porous media. is modeled by a coupled system of partial differential equations. The approximation of this system using a finite element method for the brine, radionuclides, and heat combined with a mixed finite element method for the pressure and velocity are analyzed. Optimal order error estimates in H-1 and L-2 are derived. This paper improves upon previously derived estimates in two aspects. Firstly, the error analysis is given with no restriction on the diffusion tensors. That is, it has included the effects of molecular diffusion and dispersion. Secondly, the 'complete compressibility' case is considered.

UNIFORM SUPERCONVERGENCE ANALYSIS OF A TWO-GRID MIXED FINITE ELEMENT METHOD FOR THE TIME-DEPENDENT BI-WAVE PROBLEM MODELING D-WAVE SUPERCONDUCTORS

In this paper,a two-grid mixed finite element method(MFEM)of implicit Backward Euler(BE)formula is presented for the fourth order time-dependent singularly perturbed Bi-wave problem for d-wave superconductors by the nonconforming EQ_(1)^(rot) element.In this approach,the original nonlinear system is solved on the coarse mesh through the Newton iteration method,and then the linear system is computed on the fine mesh with Taylor’s expansion.Based on the high accuracy results of the chosen element,the uniform superclose and superconvergent estimates in the broken H^(1)-norm are derived,which are independent of the negative powers of the perturbation parameter appeared in the considered problem.Numerical results illustrate that the computing cost of the proposed two-grid method is much less than that of the conventional Galerkin MFEM without loss of accuracy.