In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.T...In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.展开更多
In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-depe...In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.展开更多
基金supported by Natural Science Foundation of Jiangsu Province(No.BK20170374)Nature Science Research Program for Colleges and Universities of Jiangsu Province(No.17KJB110016)Scientific Research Project for University Students of Suzhou University of Science and Technology in 2017-2018
基金This work is supported by the National Natural Science Foundation of China(11661058,11761053)the Natural Science Foundation of Inner Mongolia(2017MS0107)the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07).
文摘In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.
文摘In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.