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.展开更多
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 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.展开更多
基金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 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.
基金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.