For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid an...For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.展开更多
The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentr...The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.展开更多
A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlin...A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.展开更多
In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of fini...In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.展开更多
A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is g...A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is given, which shows that the same order of accuracy can be achieved as solving the system directly in the fine mesh when h = H2. Both theoretical analysis and numerical experiments illustrate the efficiency of the algorithm for solving the coupled problem.展开更多
Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level met...Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.展开更多
We design and analyze an iterative two-grid algorithm for the finite element discretizations of strongly nonlinear elliptic boundary value problems in this paper.We propose an iterative two-grid algorithm,in which a n...We design and analyze an iterative two-grid algorithm for the finite element discretizations of strongly nonlinear elliptic boundary value problems in this paper.We propose an iterative two-grid algorithm,in which a nonlinear problem is first solved on the coarse space,and then a symmetric positive definite problem is solved on the fine space.The main contribution in this paper is to establish a first convergence analysis,which requires dealing with four coupled error estimates,for the iterative two-grid methods.We also present some numerical experiments to confirm the efficiency of the proposed algorithm.展开更多
This paper aims to construct a two-grid scheme of fully discretized expanded mixed finite element methods for optimal control problems governed by parabolic integro-differential equations and discuss a priori error es...This paper aims to construct a two-grid scheme of fully discretized expanded mixed finite element methods for optimal control problems governed by parabolic integro-differential equations and discuss a priori error estimates.The state variables and co-state variables are discretized by the lowest order Raviart-Thomas mixed finite element,and the control variable is approximated by piecewise constant functions.The time derivative is discretized by the backward Euler method.Firstly,we define some new mixed elliptic projections and prove the corresponding error estimates which play an important role in subsequent convergence analysis.Secondly,we derive a priori error estimates for all variables.Thirdly,we present a two-grid scheme and analyze its convergence.In the two-grid scheme,the solution of the parabolic optimal control problem on a fine grid is reduced to the solution of the parabolic optimal control problem on a much coarser grid and the solution of a decoupled linear algebraic system on the fine grid and the resulting solution still maintains an asymptotically optimal accuracy.At last,a numerical example is presented to verify the theoretical results.展开更多
In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery...Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
In this paper,we propose an iterative two-grid method for the edge finite element discretizations(a saddle-point system)of Perfectly Matched Layer(PML)equations to the Maxwell scattering problem in two dimensions.Firs...In this paper,we propose an iterative two-grid method for the edge finite element discretizations(a saddle-point system)of Perfectly Matched Layer(PML)equations to the Maxwell scattering problem in two dimensions.Firstly,we use a fine space to solve a discrete saddle-point system of H(grad)variational problems,denoted by auxiliary system 1.Secondly,we use a coarse space to solve the original saddle-point system.Then,we use a fine space again to solve a discrete H(curl)-elliptic variational problems,denoted by auxiliary system 2.Furthermore,we develop a regularization diagonal block preconditioner for auxiliary system 1 and use H-X preconditioner for auxiliary system 2.Hence we essentially transform the original problem in a fine space to a corresponding(but much smaller)problem on a coarse space,due to the fact that the above two preconditioners are efficient and stable.Compared with some existing iterative methods for solving saddle-point systems,such as PMinres,numerical experiments show the competitive performance of our iterative two-grid method.展开更多
Expanded mixed finite element approximation of nonlinear reaction-diffusion equations is discussed. The equations considered here are used to model the hydrologic and bio-geochemical phenomena. To linearize the mixed-...Expanded mixed finite element approximation of nonlinear reaction-diffusion equations is discussed. The equations considered here are used to model the hydrologic and bio-geochemical phenomena. To linearize the mixed-method equations, we use a two-grid method involving a small nonlinear system on a coarse gird of size H and a linear system on a fine grid of size h. Error estimates are derived which demonstrate that the error is O（△t ＋ h k＋1 ＋ H 2k＋2 d/2 ） （k ≥ 1）, where k is the degree of the approximating space for the primary variable and d is the spatial dimension. The above estimates are useful for determining an appropriate H for the coarse grid problems.展开更多
A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire c...A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.展开更多
A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium.The concentration equation is treated by a mixed finite element method with characteristics(CMFEM)and the press...A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium.The concentration equation is treated by a mixed finite element method with characteristics(CMFEM)and the pressure equation is treated by a parabolic mixed finite element method(PMFEM).Two-grid algorithm is considered to linearize nonlinear coupled system of two parabolic partial differential equations.Moreover,the L q error estimates are conducted for the pressure,Darcy velocity and concentration variables in the two-grid solutions.Both theoretical analysis and numerical experiments are presented to show that the two-grid algorithm is very effective.展开更多
Poisson-Nernst-Planck equations are widely used to describe the electrodiffusion of ions in a solvated biomolecular system. Two kinds of two-grid finite element algorithms are proposed to decouple the steady-state Poi...Poisson-Nernst-Planck equations are widely used to describe the electrodiffusion of ions in a solvated biomolecular system. Two kinds of two-grid finite element algorithms are proposed to decouple the steady-state Poisson-Nernst-Planck equations by coarse grid finite element approximations. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the two-grid algorithms for solving Poisson-Nernst-Planck equations.展开更多
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))).展开更多
We discuss the inexact two-grid methods for solving eigenvalue problems, including both partial differential and integral equations. Instead of solving the linear system exactly in both traditional two-grid and accele...We discuss the inexact two-grid methods for solving eigenvalue problems, including both partial differential and integral equations. Instead of solving the linear system exactly in both traditional two-grid and accelerated two-grid method, we point out that it is enough to apply an inexact solver to the fine grid problems, which will cut down the computational cost. Different stopping criteria for both methods are developed for keeping the optimality of the resulting solution. Numerical examples are provided to verify our theoretical analyses.展开更多
A two grid technique for solving the steady incompressible Navier Stokes equations in a penalty method was presented and the convergence of numerical solutions was analyzed. If a coarse size H and a fine size ...A two grid technique for solving the steady incompressible Navier Stokes equations in a penalty method was presented and the convergence of numerical solutions was analyzed. If a coarse size H and a fine size h satisfy H=O(h 13-s )(s=0(n=2);s=12(n=3), where n is a space dimension), this method has the same convergence accuracy as the usual finite element method. But the two grid method can save a lot of computation time for its brief calculation. Moreover, a numerical test was couducted in order to verify the correctness of above theoretical analysis.展开更多
The purpose of this paper is to analyze an efficient method for the solution of the nonlinear system resulting from the discretization of the elliptic Monge-Ampere equation by a C0 interior penalty method with Lagrang...The purpose of this paper is to analyze an efficient method for the solution of the nonlinear system resulting from the discretization of the elliptic Monge-Ampere equation by a C0 interior penalty method with Lagrange finite elements.We consider the two-grid method for nonlinear equations which consists in solving the discrete nonlinear system on a coarse mesh and using that solution as initial guess for one iteration of Newton’s method on a finer mesh.Thus both steps are inexpensive.We give quasi-optimal W1,1 error estimates for the discretization and estimate the difference between the interior penalty solution and the two-grid numerical solution.Numerical experiments confirm the computational efficiency of the approach compared to Newton’s method on the fine mesh.展开更多
文摘For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.
基金Project supported by the State Key Program of National Natural Science Foundation of China(No.11931003)the National Natural Science Foundation of China(Nos.41974133,11671157,11971410)。
文摘The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.
文摘A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.
基金Project supported by the National Natural Science Foundation of China(Nos.11671157 and11826212)
文摘In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.
基金supported by National Foundation of Natural Science(11471092,11326231)Zhejiang Provincial Natural Science Foundation of China(LZ13A010003)
文摘A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is given, which shows that the same order of accuracy can be achieved as solving the system directly in the fine mesh when h = H2. Both theoretical analysis and numerical experiments illustrate the efficiency of the algorithm for solving the coupled problem.
文摘Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.
基金funded by the Science and Technology Development Fund,Macao SAR(Grant Nos.0070/2019/A2,0031/2022/A1)supported by the National Natural Science Foundation of China(Grant No.11901212)supported by the National Natural Science Foundation of China(Grant No.12071160).
文摘We design and analyze an iterative two-grid algorithm for the finite element discretizations of strongly nonlinear elliptic boundary value problems in this paper.We propose an iterative two-grid algorithm,in which a nonlinear problem is first solved on the coarse space,and then a symmetric positive definite problem is solved on the fine space.The main contribution in this paper is to establish a first convergence analysis,which requires dealing with four coupled error estimates,for the iterative two-grid methods.We also present some numerical experiments to confirm the efficiency of the proposed algorithm.
基金supported by National Natural Science Foundation of China(No.12571388)the Visiting scholar program of National Natural Science Foundation of China(No.12426616)Natural Science Research Start-up Foundation of Recruiting Talents of Nanjing University of Posts and Telecommunications(No.NY223127)。
文摘This paper aims to construct a two-grid scheme of fully discretized expanded mixed finite element methods for optimal control problems governed by parabolic integro-differential equations and discuss a priori error estimates.The state variables and co-state variables are discretized by the lowest order Raviart-Thomas mixed finite element,and the control variable is approximated by piecewise constant functions.The time derivative is discretized by the backward Euler method.Firstly,we define some new mixed elliptic projections and prove the corresponding error estimates which play an important role in subsequent convergence analysis.Secondly,we derive a priori error estimates for all variables.Thirdly,we present a two-grid scheme and analyze its convergence.In the two-grid scheme,the solution of the parabolic optimal control problem on a fine grid is reduced to the solution of the parabolic optimal control problem on a much coarser grid and the solution of a decoupled linear algebraic system on the fine grid and the resulting solution still maintains an asymptotically optimal accuracy.At last,a numerical example is presented to verify the theoretical results.
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金Supported by the Natural Science Foundation of Guangxi Province(Grant Nos.2023GXNSFAA026067,2024GXN SFAA010521)the National Natural Science Foundation of China(Nos.12361079,12201149,12261026).
文摘Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金This work was partially supported by NSFC Project(Grant No.11031006,91130002,11171281,10971059,11026091)the Key Project of Scientific Research Fund of Hunan Provincial Science and Technology Department(Grant No.2011FJ2011)Hunan Provincial Innovation Foundation for Postgraduate(CX2010B245,CX2010B246).
文摘In this paper,we propose an iterative two-grid method for the edge finite element discretizations(a saddle-point system)of Perfectly Matched Layer(PML)equations to the Maxwell scattering problem in two dimensions.Firstly,we use a fine space to solve a discrete saddle-point system of H(grad)variational problems,denoted by auxiliary system 1.Secondly,we use a coarse space to solve the original saddle-point system.Then,we use a fine space again to solve a discrete H(curl)-elliptic variational problems,denoted by auxiliary system 2.Furthermore,we develop a regularization diagonal block preconditioner for auxiliary system 1 and use H-X preconditioner for auxiliary system 2.Hence we essentially transform the original problem in a fine space to a corresponding(but much smaller)problem on a coarse space,due to the fact that the above two preconditioners are efficient and stable.Compared with some existing iterative methods for solving saddle-point systems,such as PMinres,numerical experiments show the competitive performance of our iterative two-grid method.
基金Supported by the National Natural Science Foundation of China Grant (No. 10771124)the Research Fund for Doctoral Program of High Education by State Education Ministry of China (No. 20060422006)+1 种基金the Program for Innovative Research Team in Ludong Universitythe Discipline Construction Fund of Ludong University
文摘Expanded mixed finite element approximation of nonlinear reaction-diffusion equations is discussed. The equations considered here are used to model the hydrologic and bio-geochemical phenomena. To linearize the mixed-method equations, we use a two-grid method involving a small nonlinear system on a coarse gird of size H and a linear system on a fine grid of size h. Error estimates are derived which demonstrate that the error is O（△t ＋ h k＋1 ＋ H 2k＋2 d/2 ） （k ≥ 1）, where k is the degree of the approximating space for the primary variable and d is the spatial dimension. The above estimates are useful for determining an appropriate H for the coarse grid problems.
基金Project supported by the National Natural Science Foundation of China(No.40074031)the Science Foundation of the Science and Technology Commission of Shanghai Municipalitythe Program for Young Excellent Talents in Tongji University(No.2007kj008)
文摘A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.
基金Natural Science Foundation of Guangdong province,China(2018A0303100016)Educational Commission of Guangdong Province,China(2019KTSCX174)+1 种基金The second author's work is supported by the State Key Program of National Natural Science Foundation of China(11931003)National Natural Science Foundation of China(41974133,11671157).
文摘A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium.The concentration equation is treated by a mixed finite element method with characteristics(CMFEM)and the pressure equation is treated by a parabolic mixed finite element method(PMFEM).Two-grid algorithm is considered to linearize nonlinear coupled system of two parabolic partial differential equations.Moreover,the L q error estimates are conducted for the pressure,Darcy velocity and concentration variables in the two-grid solutions.Both theoretical analysis and numerical experiments are presented to show that the two-grid algorithm is very effective.
文摘Poisson-Nernst-Planck equations are widely used to describe the electrodiffusion of ions in a solvated biomolecular system. Two kinds of two-grid finite element algorithms are proposed to decouple the steady-state Poisson-Nernst-Planck equations by coarse grid finite element approximations. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the two-grid algorithms for solving Poisson-Nernst-Planck equations.
基金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))).
基金The authors are grateful to Prof. Zhaojun Bai in University of Cali- fornia, Davis, and Prof. Carlos J. Garcia-Cervera in University of California, Santa Barbara for their helpful discussions. The authors are grateful to the editor and the referees for their valuable comments, which improves the quality of the paper greatly. Weiguo Gao is supported by the National Natural Science Foundation of China under grants 91330202, Special Funds for Major State Basic Research Projects of China （2015CB858560003）, and Shanghai Science and Technology Development Funds 13dz2260200 and 13511504300. Qun Gu acknowledges the financial support from China Scholarship Council （No. 2011610055）.
文摘We discuss the inexact two-grid methods for solving eigenvalue problems, including both partial differential and integral equations. Instead of solving the linear system exactly in both traditional two-grid and accelerated two-grid method, we point out that it is enough to apply an inexact solver to the fine grid problems, which will cut down the computational cost. Different stopping criteria for both methods are developed for keeping the optimality of the resulting solution. Numerical examples are provided to verify our theoretical analyses.
文摘A two grid technique for solving the steady incompressible Navier Stokes equations in a penalty method was presented and the convergence of numerical solutions was analyzed. If a coarse size H and a fine size h satisfy H=O(h 13-s )(s=0(n=2);s=12(n=3), where n is a space dimension), this method has the same convergence accuracy as the usual finite element method. But the two grid method can save a lot of computation time for its brief calculation. Moreover, a numerical test was couducted in order to verify the correctness of above theoretical analysis.
基金Gerard Awanou was partially supported by NSF DMS grant#1720276Hengguang Li by NSF DMS grant#1819041 and by NSF China Grant 11628104.
文摘The purpose of this paper is to analyze an efficient method for the solution of the nonlinear system resulting from the discretization of the elliptic Monge-Ampere equation by a C0 interior penalty method with Lagrange finite elements.We consider the two-grid method for nonlinear equations which consists in solving the discrete nonlinear system on a coarse mesh and using that solution as initial guess for one iteration of Newton’s method on a finer mesh.Thus both steps are inexpensive.We give quasi-optimal W1,1 error estimates for the discretization and estimate the difference between the interior penalty solution and the two-grid numerical solution.Numerical experiments confirm the computational efficiency of the approach compared to Newton’s method on the fine mesh.
