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.展开更多
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.展开更多
Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R...Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R) variational principle. The main features of the GME formulations are that the common C0-continuous polynomial shape functions for displacement methods are used to express both displacement and stress variables, and the coefficient matrix of these formulations is not only automatically symmetric but also invertible. Hence, the numerical results of the generalized mixed methods based on the GME formulations are stable. Displacement as well as stress results can be obtained directly from the algebraic system for finite element analysis after introducing stress and displacement boundary conditions simultaneously. Numerical examples show that displacement and stress results retain the same accuracy. The results of the noncompatible generalized mixed method proposed herein are more accurate than those of the standard noncompatible displacement method. The noncompatible generalized mixed element is less sensitive to element geometric distortions.展开更多
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.展开更多
Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditi...Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.展开更多
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).展开更多
The flows of rigid polyvinyl chloride (R-PVC) in co-rotating twin screw extruders with screw mixing elements and regular screw elements were simulated by using the finite element method. The three-dimensional, non-i...The flows of rigid polyvinyl chloride (R-PVC) in co-rotating twin screw extruders with screw mixing elements and regular screw elements were simulated by using the finite element method. The three-dimensional, non-isothermal flow fields of R-PVC in the two kinds of screw elements were calculated. The mixing performance of each type of element was studied by the particle tracking analysis method. The results show that the temperature distribution and shear-rate distribution are more uniform in the flow channel with screw mixing elements than in the flow channel with regular screw elements. Screw mixing elements provide better distributive and dispersive mixing performance but worse conveying capacity than regular screw elements.展开更多
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.展开更多
A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order...A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.展开更多
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.
Two simplifled and stabilized mixed element formats for the Stokes problem are derived by bubble function, and their convergence, i.e., error analysis, are proved. These formats can save more freedom degrees than othe...Two simplifled and stabilized mixed element formats for the Stokes problem are derived by bubble function, and their convergence, i.e., error analysis, are proved. These formats can save more freedom degrees than other usual formats.展开更多
The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least sq...The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least squares mixed finite element format for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of its solution are derived. Through this method, the combination among the mixed finite element spaces does not demand the discrete Babuska-Brezzi stability conditions so that the mixed finite element spaces could be chosen arbitrartily and the error estimates with optimal order could be obtained.展开更多
A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution...A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution are first discussed, and then a posteriori error estimator based on the NGME method is derived.展开更多
A nonlinear Galerkin/Petrov-least squares mixed element (NGPLSME) method for the stationary Navier-Stokes equations is presented and analyzed. The scheme is that Petrov-least squares forms of residuals are added to th...A nonlinear Galerkin/Petrov-least squares mixed element (NGPLSME) method for the stationary Navier-Stokes equations is presented and analyzed. The scheme is that Petrov-least squares forms of residuals are added to the nonlinear Galerkin mixed element method so that it is stable for any combination of discrete velocity and pressure spaces without requiring the Babu*lka-Brezzi stability condition. The existence, uniqueness and convergence (at optimal rate) of the NGPLSME solution is proved in the case of sufficient viscosity (or small data).展开更多
This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established with...This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.展开更多
The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted...The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.展开更多
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.展开更多
A nonlinear Galerkin mixed element (NGME) method for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of the NGME solution are derived.
文摘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.
文摘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 (Grant 11502286)
文摘Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R) variational principle. The main features of the GME formulations are that the common C0-continuous polynomial shape functions for displacement methods are used to express both displacement and stress variables, and the coefficient matrix of these formulations is not only automatically symmetric but also invertible. Hence, the numerical results of the generalized mixed methods based on the GME formulations are stable. Displacement as well as stress results can be obtained directly from the algebraic system for finite element analysis after introducing stress and displacement boundary conditions simultaneously. Numerical examples show that displacement and stress results retain the same accuracy. The results of the noncompatible generalized mixed method proposed herein are more accurate than those of the standard noncompatible displacement method. The noncompatible generalized mixed element is less sensitive to element geometric distortions.
基金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.
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.
文摘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).
文摘The flows of rigid polyvinyl chloride (R-PVC) in co-rotating twin screw extruders with screw mixing elements and regular screw elements were simulated by using the finite element method. The three-dimensional, non-isothermal flow fields of R-PVC in the two kinds of screw elements were calculated. The mixing performance of each type of element was studied by the particle tracking analysis method. The results show that the temperature distribution and shear-rate distribution are more uniform in the flow channel with screw mixing elements than in the flow channel with regular screw elements. Screw mixing elements provide better distributive and dispersive mixing performance but worse conveying capacity than regular screw elements.
基金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.
基金supported by the National Natural Science Foundation of China (No. 10601022)NSF ofInner Mongolia Autonomous Region of China (No. 200607010106)513 and Science Fund of InnerMongolia University for Distinguished Young Scholars (No. ND0702)
文摘A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.
文摘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.
文摘Two simplifled and stabilized mixed element formats for the Stokes problem are derived by bubble function, and their convergence, i.e., error analysis, are proved. These formats can save more freedom degrees than other usual formats.
基金Project supported by the National Natural Science Foundation of China (Nos.10471100 and 40437017)the Science and Technology Foundation of Beijing Jiaotong University
文摘The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least squares mixed finite element format for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of its solution are derived. Through this method, the combination among the mixed finite element spaces does not demand the discrete Babuska-Brezzi stability conditions so that the mixed finite element spaces could be chosen arbitrartily and the error estimates with optimal order could be obtained.
文摘A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution are first discussed, and then a posteriori error estimator based on the NGME method is derived.
文摘A nonlinear Galerkin/Petrov-least squares mixed element (NGPLSME) method for the stationary Navier-Stokes equations is presented and analyzed. The scheme is that Petrov-least squares forms of residuals are added to the nonlinear Galerkin mixed element method so that it is stable for any combination of discrete velocity and pressure spaces without requiring the Babu*lka-Brezzi stability condition. The existence, uniqueness and convergence (at optimal rate) of the NGPLSME solution is proved in the case of sufficient viscosity (or small data).
文摘This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.
文摘The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.
基金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.
基金Project supported by the National Natural Science Foundation of China (Nos.10471100 and 40437017)
文摘A nonlinear Galerkin mixed element (NGME) method for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of the NGME solution are derived.