Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadrati...Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.展开更多
In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular famil...In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.展开更多
Based on an asymptotic expansion of finite element, an extrapolation cascadic multigrid method （EXCMG） is proposed, in which the new extrapolation and quadratic interpolation are used to provide a better initial val...Based on an asymptotic expansion of finite element, an extrapolation cascadic multigrid method （EXCMG） is proposed, in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid. In the case of multiple grids, both superconvergence error in H^1-norm and the optimal error in l2-norm are analyzed. The numerical experiment shows the advantage of EXCMG in comparison with CMG.展开更多
In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-d...In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.展开更多
In this article,we study a 2D nonlinear time-fractional Rayleigh-Stokes problem,which has an anomalous subdiffusion term,on triangular meshes by quadratic finite volume element schemes.Time-fractional derivative,defin...In this article,we study a 2D nonlinear time-fractional Rayleigh-Stokes problem,which has an anomalous subdiffusion term,on triangular meshes by quadratic finite volume element schemes.Time-fractional derivative,defined by Caputo fractional derivative,is discretized through L2−1σformula,and a two step scheme is used to approximate the time first-order derivative at time tn−α/2,where the nonlinear term is approximated by using a matching linearized difference scheme.A family of quadratic finite volume element schemes with two parameters are proposed for the spatial discretization,where the range of values for two parameters areβ1∈(0,1/2),β2∈(0,2/3).For testing the precision of numerical algorithms,we calculate some numerical examples which have known exact solution or unknown exact solution by several kinds of quadratic finite volume element schemes,and contrast with the results of an existing quadratic finite element scheme by drawing diversified comparison plots and showing the detailed data of L2 error results and convergence orders.Numerical results indicate that,L2 error estimate of one scheme with parameters β_(1)=(3−√3)/6,β2=(6+√3−√21+6√3)/9 is O(h^(3)+△t^(2)),and L^(2) error estimates of other schemes are O(h^(2)+△t^(2)),where h and △t denote the spatial and temporal discretization parameters,respectively.展开更多
In this paper,we construct an H1-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates.The element has optimal approximation rates.Using this quadratic element,tw...In this paper,we construct an H1-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates.The element has optimal approximation rates.Using this quadratic element,two stable discretizations for the Stokes equations are developed,which can be viewed as the extensions of the P2-P0 and the Q2-(discontinuous)P1 elements,respectively,to polygonal meshes.Numerical results are presented,which support our theoretical claims.展开更多
基金supported by the National Natural Science Foundation of China(Nos.11871009,12271055)the Foundation of LCP and the Foundation of CAEP(CX20210044).
文摘Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.
基金The first author is partially supported by the U.S. Department of Energy, Office of Science, Office of Advanced Scientific Computing Research as part of the Collaboratory on Mathematics for Mesoscopic Modeling of Materials under Award Number DE-SC-0009249, and the Key Program of National Natural Science Foundation of China with Grant No. 91430215. The second author is supported by State Key Laboratory of Scientific and Engineering Computing （LSEC）, National Center for Mathematics and Interdisciplinary Sciences of Chinese Academy of Sciences （NCMIS）, and National Natural Science Foundation of China with Grant No. 11471026 he is thankful to the Center for Computational Mathematics and Applications, the Pennsylvania State University, where he worked on this manuscript as a visiting scholar. The authors are grateful to Professor Jinchao Xu, Dr. Yuanming Xiao and Dr. Maximilian Metti for their valuable suggestions and discussions, to Professor Haijun Wu for his valuable help on preparing the numerical example, and to the anonymous referee for the valuable comments and suggestion which lead to improvements of the paper.
文摘In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.
基金Supported by National Natural Science Foundation of China (10771063)the Doctor Programme of the National Education Committee (20050542006)
文摘Based on an asymptotic expansion of finite element, an extrapolation cascadic multigrid method （EXCMG） is proposed, in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid. In the case of multiple grids, both superconvergence error in H^1-norm and the optimal error in l2-norm are analyzed. The numerical experiment shows the advantage of EXCMG in comparison with CMG.
基金the National Natural Science Fund(11661058,11761053)Natural Science Fund of Inner Mongolia Autonomous Region(2016MS0102,2017MS0107)+1 种基金Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07)National Undergraduate Innovative Training Project of Inner Mongolia University(201710126026).
文摘In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.
基金This work was partially supported by the National Natural Science Foundation of China(No.11871009).
文摘In this article,we study a 2D nonlinear time-fractional Rayleigh-Stokes problem,which has an anomalous subdiffusion term,on triangular meshes by quadratic finite volume element schemes.Time-fractional derivative,defined by Caputo fractional derivative,is discretized through L2−1σformula,and a two step scheme is used to approximate the time first-order derivative at time tn−α/2,where the nonlinear term is approximated by using a matching linearized difference scheme.A family of quadratic finite volume element schemes with two parameters are proposed for the spatial discretization,where the range of values for two parameters areβ1∈(0,1/2),β2∈(0,2/3).For testing the precision of numerical algorithms,we calculate some numerical examples which have known exact solution or unknown exact solution by several kinds of quadratic finite volume element schemes,and contrast with the results of an existing quadratic finite element scheme by drawing diversified comparison plots and showing the detailed data of L2 error results and convergence orders.Numerical results indicate that,L2 error estimate of one scheme with parameters β_(1)=(3−√3)/6,β2=(6+√3−√21+6√3)/9 is O(h^(3)+△t^(2)),and L^(2) error estimates of other schemes are O(h^(2)+△t^(2)),where h and △t denote the spatial and temporal discretization parameters,respectively.
基金supported by the NSFC grant 11671210 and 12171244.
文摘In this paper,we construct an H1-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates.The element has optimal approximation rates.Using this quadratic element,two stable discretizations for the Stokes equations are developed,which can be viewed as the extensions of the P2-P0 and the Q2-(discontinuous)P1 elements,respectively,to polygonal meshes.Numerical results are presented,which support our theoretical claims.
