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 m...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.展开更多
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 t...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.展开更多
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 ...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 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 estima...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.展开更多
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 dualiza...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 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.
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 equ...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.展开更多
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...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.展开更多
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 s...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.展开更多
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...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.展开更多
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 stu...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.展开更多
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-discr...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.展开更多
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 ...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.展开更多
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...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 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, radi...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.展开更多
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 n...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.展开更多
In this article,we propose a new finite element spaceΛh for the expanded mixed finite element method(EMFEM)for second-order elliptic problems to guarantee its computing capability and reduce the computation cost.The ...In this article,we propose a new finite element spaceΛh for the expanded mixed finite element method(EMFEM)for second-order elliptic problems to guarantee its computing capability and reduce the computation cost.The new finite element spaceΛh is designed in such a way that the strong requirement V h⊂Λh in[9]is weakened to{v h∈V h;d i v v h=0}⊂Λh so that it needs fewer degrees of freedom than its classical counterpart.Furthermore,the newΛh coupled with the Raviart-Thomas space satisfies the inf-sup condition,which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix,and thus the existence,uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in R d,d=2,3 and for triangular partitions in R 2.Also,the solvability of the EMFEM for triangular partition in R 3 can be directly proved without the inf-sup condition.Numerical experiments are conducted to confirm these theoretical findings.展开更多
Mixed finite element methods are considered for a ferrofluid flow model with magnetization paralleled to the magnetic field.The ferrofluid model is a coupled system of the Maxwell equations and the incompressible Navi...Mixed finite element methods are considered for a ferrofluid flow model with magnetization paralleled to the magnetic field.The ferrofluid model is a coupled system of the Maxwell equations and the incompressible Navier-Stokes equations.By skillfully introducing some new variables,the model is rewritten as several decoupled subsystems that can be solved independently.Mixed finite element formulations are given to discretize the decoupled systems with proper finite element spaces.Existence and uniqueness of the mixed finite element solutions are shown,and optimal order error estimates are obtained under some reasonable assumptions.Numerical experiments confirm the theoretical results.展开更多
Superconvergence of the mixed finite element methods for 2-d Maxwell equations is studied in this paper. Two order of superconvergent factor can be obtained for the k-th Nedelec elements on the rectangular meshes.
文摘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.
基金Supported by the National Natural Science Foundation of China (10601022)Natural Science Foundation of Inner Mongolia Autonomous Region (200607010106)Youth Science Foundation of Inner Mongolia University(ND0702)
文摘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.
文摘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).
基金Supported by by the National Science Foundation for Young Scholars of China(11101431)the Fundamental Research Funds for the Central Universities (12CX04082A,10CX04041A)Shandong Province Natural Science Foundation of China(ZR2010AL020)
文摘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.
文摘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.
基金Supported by National Natural Science Foundation of China(11371331)Supported by the Natural Science Foundation of Education Department of Henan Province(14B110018)
文摘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.
文摘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.
基金supported by the National Nature Science Foundation of China (Grant No 90510017)
文摘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.
文摘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.
基金Project supported by the National Natural Science Foundation of China (Nos.10471100 and 40437017)the Science and Technology Foundation of Beijing Jiaotong University
文摘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.
文摘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.
文摘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.
文摘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.
基金supported by the National Science Foundation of China(11271127)Science Research Project of Guizhou Province Education Department(QJHKYZ[2013]207)
文摘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 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.
基金supported by the National Natural Science Foundation of China(Grant Nos.12201640,12071443).
文摘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.
基金supported by NSF of China grant 11971276H.Chen was supported by NSF of China grants 12171287,10971254 and 11471196+1 种基金H.Wang was supported by the ARO MURI Grant W911NF-15-1-0562by the National Science Foundation under Grant DMS-2012291.
文摘In this article,we propose a new finite element spaceΛh for the expanded mixed finite element method(EMFEM)for second-order elliptic problems to guarantee its computing capability and reduce the computation cost.The new finite element spaceΛh is designed in such a way that the strong requirement V h⊂Λh in[9]is weakened to{v h∈V h;d i v v h=0}⊂Λh so that it needs fewer degrees of freedom than its classical counterpart.Furthermore,the newΛh coupled with the Raviart-Thomas space satisfies the inf-sup condition,which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix,and thus the existence,uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in R d,d=2,3 and for triangular partitions in R 2.Also,the solvability of the EMFEM for triangular partition in R 3 can be directly proved without the inf-sup condition.Numerical experiments are conducted to confirm these theoretical findings.
基金supported by the National Natural Science Foundation of China(Grant Nos.11971094,12171340).
文摘Mixed finite element methods are considered for a ferrofluid flow model with magnetization paralleled to the magnetic field.The ferrofluid model is a coupled system of the Maxwell equations and the incompressible Navier-Stokes equations.By skillfully introducing some new variables,the model is rewritten as several decoupled subsystems that can be solved independently.Mixed finite element formulations are given to discretize the decoupled systems with proper finite element spaces.Existence and uniqueness of the mixed finite element solutions are shown,and optimal order error estimates are obtained under some reasonable assumptions.Numerical experiments confirm the theoretical results.
基金This work is subsidized by the special funds for major state basic research projects (No. 1999032800).
文摘Superconvergence of the mixed finite element methods for 2-d Maxwell equations is studied in this paper. Two order of superconvergent factor can be obtained for the k-th Nedelec elements on the rectangular meshes.