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 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 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.展开更多
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 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.展开更多
This study enhances the application of cross-sectional warping considered mixed finite element(WMFE)formulation to accurately determine natural vibration,static displacement response,and shear and normal stress evalua...This study enhances the application of cross-sectional warping considered mixed finite element(WMFE)formulation to accurately determine natural vibration,static displacement response,and shear and normal stress evaluation with very close to the precision of solid finite elements(FEs)in two-phase/multi-phase functionally graded(FG)laminated composite beams strength using carbon nanotubes(CNTs).The principles of three dimensional(3D)elasticity theory are used to derive constitutive equations.The mixed finite element(MFE)method is improved by accounting for warping effects by displacement-based FEs within the cross-sectional domain.The MFE with two nodes has a total of 24 degrees of freedom.The two-phase material consists of a polymer matrix reinforced with aligned CNTs that are FG throughout the beam thickness.The multi-phase FG beam is modeled as a three-component composite material,consisting of CNTs,a polymer matrix,and fibers.The polymer matrix is reinforced by longitudinally aligned fibers and randomly dispersed CNT particles.The fiber volume fractions are considered to change gradually through the thickness of the beam following a power-law variation.The W-MFE achieves satisfactory results with fewer degrees of freedom than 3D solid FEs.Benchmark examples examine the effects of ply orientation,configuration,and fiber gradation on FG beam behavior.展开更多
Beam is one of the common structures in engineering,with the development of technology,homogeneous beams no longer meet the needs of engineering structural design,for this reason,people have researched the non-homogen...Beam is one of the common structures in engineering,with the development of technology,homogeneous beams no longer meet the needs of engineering structural design,for this reason,people have researched the non-homogeneous beams.In this paper,we study the mixed finite element method for the vibration problem of non-homogeneous damped beams.The fourth-order differential equations are transformed into a system of low-order partial differential equations by introducing intermediate variables,constructing a semidiscrete extended mixed finite element format,proving the existence and uniqueness of the solution of the format,and utilizing the elliptic projection operator for the error estimation.The time derivative term is discretized by the central difference,and the fully discrete mixed element format is given to prove the stability and convergence of the format.The feasibility and effectiveness of the mixed method are verified by numerical examples,and the effects of different damping coefficientsμon beam vibration are investigated.展开更多
A weak Galerkin mixed finite element method is studied for linear elasticity problems without the requirement of symmetry.The key of numerical methods in mixed formulation is the symmetric constraint of numerical stre...A weak Galerkin mixed finite element method is studied for linear elasticity problems without the requirement of symmetry.The key of numerical methods in mixed formulation is the symmetric constraint of numerical stress.In this paper,we introduce the discrete symmetric weak divergence to ensure the symmetry of numerical stress.The corresponding stabilizer is presented to guarantee the weak continuity.This method does not need extra unknowns.The optimal error estimates in discrete H^(1) and L^(2) norms are established.The numerical examples in 2D and 3D are presented to demonstrate the efficiency and locking-free property.展开更多
In this paper,we propose mixed finite element methods for the Reissner-Mindlin Plate Problem by introducing the bending moment as an independent variable.We apply the finite element approximations of the stress field ...In this paper,we propose mixed finite element methods for the Reissner-Mindlin Plate Problem by introducing the bending moment as an independent variable.We apply the finite element approximations of the stress field and the displacement field constructed for the elasticity problem by Hu(J Comp Math 33:283–296,2015),Hu and Zhang(arXiv:1406.7457,2014)to solve the bending moment and the rotation for the Reissner-Mindlin Plate Problem.We propose two triples of finite element spaces to approximate the bending moment,the rotation,and the displacement.The feature of these methods is that they need neither reduction terms nor penalty terms.Then,we prove the well-posedness of the discrete problem and obtain the optimal estimates independent of the plate thickness.Finally,we present some numerical examples to demonstrate the theoretical results.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
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 paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bil...In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ｜｜·｜｜div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.展开更多
In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergenc...In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.展开更多
In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are def...In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are defined on three different meshes. This method is a conforming quadrilateral Q1 × Q1 - P0 element approximation for the Navier-Stokes equations. First-order error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to illustrate the effectiveness of the proposed method.展开更多
In this paper,we present an efficient method of two-grid scheme for the approximation of two-dimensional nonlinear parabolic equations using an expanded mixed finite element method.We use two Newton iterations on the ...In this paper,we present an efficient method of two-grid scheme for the approximation of two-dimensional nonlinear parabolic equations using an expanded mixed finite element method.We use two Newton iterations on the fine grid in our methods.Firstly,we solve an original nonlinear problem on the coarse nonlinear grid,then we use Newton iterations on the fine grid twice.The two-grid idea is from Xu's work[SIAM J.Numer.Anal.,33(1996),pp.1759–1777]on standard finite method.We also obtain the error estimates for the algorithms of the two-grid method.It is shown that the algorithm achieve asymptotically optimal approximation rate with the two-grid methods as long as the mesh sizes satisfy h=O(H^((4k+1)/(k+1))).展开更多
The two-level penalty mixed finite element method for the stationary Navier-Stokes equations based on Taylor-Hood element is considered in this paper. Two algorithms are proposed and analyzed. Moreover, the optimal st...The two-level penalty mixed finite element method for the stationary Navier-Stokes equations based on Taylor-Hood element is considered in this paper. Two algorithms are proposed and analyzed. Moreover, the optimal stability analysis and error estimate for these two algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms.展开更多
In this paper,we propose a new numerical scheme for the coupled Stokes-Darcy model with the Beavers-Joseph-Saffman interface condition.We use the weak Galerkin method to discretize the Stokes equation and the mixed fi...In this paper,we propose a new numerical scheme for the coupled Stokes-Darcy model with the Beavers-Joseph-Saffman interface condition.We use the weak Galerkin method to discretize the Stokes equation and the mixed finite element method to discretize the Darcy equation.A discrete inf-sup condition is proved and the optimal error estimates are also derived.Numerical experiments validate the theoretical analysis.展开更多
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.展开更多
基金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 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.
基金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.
文摘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 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.
基金funding provided by the Scientific and Technological Research Council of Türkiye(TÜBİTAK).
文摘This study enhances the application of cross-sectional warping considered mixed finite element(WMFE)formulation to accurately determine natural vibration,static displacement response,and shear and normal stress evaluation with very close to the precision of solid finite elements(FEs)in two-phase/multi-phase functionally graded(FG)laminated composite beams strength using carbon nanotubes(CNTs).The principles of three dimensional(3D)elasticity theory are used to derive constitutive equations.The mixed finite element(MFE)method is improved by accounting for warping effects by displacement-based FEs within the cross-sectional domain.The MFE with two nodes has a total of 24 degrees of freedom.The two-phase material consists of a polymer matrix reinforced with aligned CNTs that are FG throughout the beam thickness.The multi-phase FG beam is modeled as a three-component composite material,consisting of CNTs,a polymer matrix,and fibers.The polymer matrix is reinforced by longitudinally aligned fibers and randomly dispersed CNT particles.The fiber volume fractions are considered to change gradually through the thickness of the beam following a power-law variation.The W-MFE achieves satisfactory results with fewer degrees of freedom than 3D solid FEs.Benchmark examples examine the effects of ply orientation,configuration,and fiber gradation on FG beam behavior.
文摘Beam is one of the common structures in engineering,with the development of technology,homogeneous beams no longer meet the needs of engineering structural design,for this reason,people have researched the non-homogeneous beams.In this paper,we study the mixed finite element method for the vibration problem of non-homogeneous damped beams.The fourth-order differential equations are transformed into a system of low-order partial differential equations by introducing intermediate variables,constructing a semidiscrete extended mixed finite element format,proving the existence and uniqueness of the solution of the format,and utilizing the elliptic projection operator for the error estimation.The time derivative term is discretized by the central difference,and the fully discrete mixed element format is given to prove the stability and convergence of the format.The feasibility and effectiveness of the mixed method are verified by numerical examples,and the effects of different damping coefficientsμon beam vibration are investigated.
基金supported by the National Natural Science Foundation of China(Grant Nos.11871038,12131014).
文摘A weak Galerkin mixed finite element method is studied for linear elasticity problems without the requirement of symmetry.The key of numerical methods in mixed formulation is the symmetric constraint of numerical stress.In this paper,we introduce the discrete symmetric weak divergence to ensure the symmetry of numerical stress.The corresponding stabilizer is presented to guarantee the weak continuity.This method does not need extra unknowns.The optimal error estimates in discrete H^(1) and L^(2) norms are established.The numerical examples in 2D and 3D are presented to demonstrate the efficiency and locking-free property.
文摘In this paper,we propose mixed finite element methods for the Reissner-Mindlin Plate Problem by introducing the bending moment as an independent variable.We apply the finite element approximations of the stress field and the displacement field constructed for the elasticity problem by Hu(J Comp Math 33:283–296,2015),Hu and Zhang(arXiv:1406.7457,2014)to solve the bending moment and the rotation for the Reissner-Mindlin Plate Problem.We propose two triples of finite element spaces to approximate the bending moment,the rotation,and the displacement.The feature of these methods is that they need neither reduction terms nor penalty terms.Then,we prove the well-posedness of the discrete problem and obtain the optimal estimates independent of the plate thickness.Finally,we present some numerical examples to demonstrate the theoretical results.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金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 the National Natural Science Foundation of China(No.10971203,11271340,11101384)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20094101110006)
文摘In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ｜｜·｜｜div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.
基金supported by National Natural Science Foundation of China(Grant Nos.11001259,11031006,11071265,11201501 and 91230110)National Basic Research Program of China(973 Project)(Grant No. 2011CB309703)+3 种基金International S&T Cooperation Program of China(Grant No. 2010DFR00700)Croucher Foundation of Hong Kong Baptist Universitythe National Center for Mathematics and Interdisciplinary Science,CAS,the President Foundation of AMSS-CASthe Fundamental Research Funds for the CentralUniversities(Grant No. 2012121003)
文摘In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.
基金supported by National Nature Science Foundation of China (No.10471073)
文摘In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are defined on three different meshes. This method is a conforming quadrilateral Q1 × Q1 - P0 element approximation for the Navier-Stokes equations. First-order error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to illustrate the effectiveness of the proposed method.
基金supported by Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China 10971074+1 种基金the National Basic Research Program under the Grant 2005CB321703Hunan Provincial Innovation Foundation For Postgraduate CX2009B119。
文摘In this paper,we present an efficient method of two-grid scheme for the approximation of two-dimensional nonlinear parabolic equations using an expanded mixed finite element method.We use two Newton iterations on the fine grid in our methods.Firstly,we solve an original nonlinear problem on the coarse nonlinear grid,then we use Newton iterations on the fine grid twice.The two-grid idea is from Xu's work[SIAM J.Numer.Anal.,33(1996),pp.1759–1777]on standard finite method.We also obtain the error estimates for the algorithms of the two-grid method.It is shown that the algorithm achieve asymptotically optimal approximation rate with the two-grid methods as long as the mesh sizes satisfy h=O(H^((4k+1)/(k+1))).
基金This work was supported in part by the National Natural Science Foundation of China （Grant Nos. 11271298, 11271313, 61163027）, the Key Project of Chinese Ministry of Education （Grant No. 212197）, the Natural Science Foundation of Xinjiang Province （Grant No. 2013211B01）, and the Doctoral Foundation of Xinjiang University （Grant No. BS120102）.
文摘The two-level penalty mixed finite element method for the stationary Navier-Stokes equations based on Taylor-Hood element is considered in this paper. Two algorithms are proposed and analyzed. Moreover, the optimal stability analysis and error estimate for these two algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms.
基金National Natural Science Foundation of China(Grant Nos.11901015,11971198,91630201,11871245,11771179 and 11826101)the Program for Cheung Kong Scholars(Grant No.Q2016067)Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education,Jilin University。
文摘In this paper,we propose a new numerical scheme for the coupled Stokes-Darcy model with the Beavers-Joseph-Saffman interface condition.We use the weak Galerkin method to discretize the Stokes equation and the mixed finite element method to discretize the Darcy equation.A discrete inf-sup condition is proved and the optimal error estimates are also derived.Numerical experiments validate the theoretical analysis.
基金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.
