Assume that L is a non-negative self-adjoint operator on L^(2)(ℝ^(n))with its heat kernels satisfying the so-called Gaussian upper bound estimate and that X is a ball quasi-Banach function space onℝ^(n) satisfying som...Assume that L is a non-negative self-adjoint operator on L^(2)(ℝ^(n))with its heat kernels satisfying the so-called Gaussian upper bound estimate and that X is a ball quasi-Banach function space onℝ^(n) satisfying some mild assumptions.Let HX,L(ℝ^(n))be the Hardy space associated with both X and L,which is defined by the Lusin area function related to the semigroup generated by L.In this article,the authors establish various maximal function characterizations of the Hardy space HX,L(ℝ^(n))and then apply these characterizations to obtain the solvability of the related Cauchy problem.These results have a wide range of generality and,in particular,the specific spaces X to which these results can be applied include the weighted space,the variable space,the mixed-norm space,the Orlicz space,the Orlicz-slice space,and the Morrey space.Moreover,the obtained maximal function characterizations of the mixed-norm Hardy space,the Orlicz-slice Hardy space,and the Morrey-Hardy space associated with L are completely new.展开更多
Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted...Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted Kato square root problem for L.More precisely,we prove that the square root L^(1/2)satisfies the weighted L^(p)estimates||L^(1/2)(f)||L_(ω)^p(R^(n))≤C||■f||L_(ω)^p(R^(n);R^(n))for any p∈(1,∞)andω∈Ap(ℝ^(n))(the class of Muckenhoupt weights),and that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,2+ε)andω∈Ap(ℝ^(n))∩RH_(2+ε/p),(R^(n))(the class of reverse Hölder weights),whereε∈(0,∞)is a constant depending only on n and the operator L,and where(2+ε/p)'denotes the Hölder conjugate exponent of 2+ε/p.Moreover,for any given q∈(2,∞),we give a sufficient condition to obtain that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,q)andω∈A_(p)(R^(n))∩pRH_(q/p),(R^(n)).As an application,we prove that when the coefficient matrix A that appears in L satisfies the small BMO condition,the Riesz transform∇L^(−1/2)is bounded on L_(ω)^(p)(ℝ^(n))for any given p∈(1,∞)andω∈Ap(ℝ^(n)).Furthermore,applications to the weighted L^(2)-regularity problem with the Dirichlet or the Neumann boundary condition are also given.展开更多
In this article, we consider positive steady state solutions and dynamics for a spatially heterogeneous predator-prey system with modified Leslie-Gower and Holling-Type II schemes. The heterogeneity here is created by...In this article, we consider positive steady state solutions and dynamics for a spatially heterogeneous predator-prey system with modified Leslie-Gower and Holling-Type II schemes. The heterogeneity here is created by the degeneracy of the intra-specific pressures for the prey. By the bifurcation method, the degree theory, and a priori estimates, we discuss the existence and multiplicity of positive steady states. Moreover, by the comparison argument, we also discuss the dynamical behavior for the diffusive predator-prey system.展开更多
This paper provides a finite-difference discretization for the one-and two-dimensional tempered fractional Laplacian and solves the tempered fractional Poisson equation with homogeneous Dirichlet boundary conditions.T...This paper provides a finite-difference discretization for the one-and two-dimensional tempered fractional Laplacian and solves the tempered fractional Poisson equation with homogeneous Dirichlet boundary conditions.The main ideas are to,respectively,use linear and quadratic interpolations to approximate the singularity and non-singularity of the one-dimensional tempered fractional Laplacian and bilinear and biquadratic interpolations to the two-dimensional tempered fractional Laplacian.Then,we give the truncation errors and prove the convergence.Numerical experiments verify the convergence rates of the order O(h^2−2s).展开更多
Regularization methods have been substantially applied in image restoration due to the ill-posedness of the image restoration problem.Different assumptions or priors on images are applied in the construction of image ...Regularization methods have been substantially applied in image restoration due to the ill-posedness of the image restoration problem.Different assumptions or priors on images are applied in the construction of image regularization methods.In recent years,matrix low-rank approximation has been successfully introduced in the image denoising problem and significant denoising effects have been achieved.Low-rank matrix minimization is an NP-hard problem and it is often replaced with the matrix’s weighted nuclear norm minimization(WNNM).The assumption that an image contains an extensive amount of self-similarity is the basis for the construction of the matrix low-rank approximation-based image denoising method.In this paper,we develop a model for image restoration using the sum of block matching matrices’weighted nuclear norm to be the regularization term in the cost function.An alternating iterative algorithm is designed to solve the proposed model and the convergence analyses of the algorithm are also presented.Numerical experiments show that the proposed method can recover the images much better than the existing regularization methods in terms of both recovered quantities and visual qualities.展开更多
Separable multi-block convex optimization problem appears in many mathematical and engineering fields.In the first part of this paper,we propose an inertial proximal ADMM to solve a linearly constrained separable mult...Separable multi-block convex optimization problem appears in many mathematical and engineering fields.In the first part of this paper,we propose an inertial proximal ADMM to solve a linearly constrained separable multi-block convex optimization problem,and we show that the proposed inertial proximal ADMM has global convergence under mild assumptions on the regularization matrices.Affine phase retrieval arises in holography,data separation and phaseless sampling,and it is also considered as a nonhomogeneous version of phase retrieval,which has received considerable attention in recent years.Inspired by convex relaxation of vector sparsity and matrix rank in compressive sensing and by phase lifting in phase retrieval,in the second part of this paper,we introduce a compressive affine phase retrieval via lifting approach to connect affine phase retrieval with multi-block convex optimization,and then based on the proposed inertial proximal ADMM for 3-block convex optimization,we propose an algorithm to recover sparse real signals from their(noisy)affine quadratic measurements.Our numerical simulations show that the proposed algorithm has satisfactory performance for affine phase retrieval of sparse real signals.展开更多
In this paper,we consider the strong convergence of the time-space fractional diffusion equation driven by fractional Gaussian noise with Hurst index H∈(1/2,1).A sharp regularity estimate of the mild solution and the...In this paper,we consider the strong convergence of the time-space fractional diffusion equation driven by fractional Gaussian noise with Hurst index H∈(1/2,1).A sharp regularity estimate of the mild solution and the numerical scheme constructed by finite element method for integral fractional Laplacian and backward Euler convolution quadrature for Riemann-Liouville time fractional derivative are proposed.With the help of inverse Laplace transform and fractional Ritz projection,we obtain the accurate error estimates in time and space.Finally,our theoretical results are accompanied by numerical experiments.展开更多
In this work,we analyze the three-step backward differentiation formula(BDF3)method for solving the Allen-Cahn equation on variable grids.For BDF2 method,the discrete orthogonal convolution(DOC)kernels are positive,th...In this work,we analyze the three-step backward differentiation formula(BDF3)method for solving the Allen-Cahn equation on variable grids.For BDF2 method,the discrete orthogonal convolution(DOC)kernels are positive,the stability and convergence analysis are well established in[Liao and Zhang,Math.Comp.,90(2021),1207–1226]and[Chen,Yu,and Zhang,arXiv:2108.02910,2021].However,the numerical analysis for BDF3 method with variable steps seems to be highly nontrivial due to the additional degrees of freedom and the non-positivity of DOC kernels.By developing a novel spectral norm inequality,the unconditional stability and convergence are rigorously proved under the updated step ratio restriction rk:=τk/τk−1≤1.405 for BDF3 method.Finally,numerical experiments are performed to illustrate the theoretical results.To the best of our knowledge,this is the first theoretical analysis of variable steps BDF3 method for the Allen-Cahn equation.展开更多
The generalized fractional elastic models govern the stochastic motion of several many-body systems,e.g., polymers, membranes, and growing interfaces. This paper focuses on the exact formulations and their asymptotic ...The generalized fractional elastic models govern the stochastic motion of several many-body systems,e.g., polymers, membranes, and growing interfaces. This paper focuses on the exact formulations and their asymptotic behaviors of the average of the solutions of the generalized fractional elastic models. So we directly analyze the Cauchy problem of the averaged generalized elastic model involving time fractional derivative and the convolution integral of a radially symmetric friction kernel with space fractional Laplacian.展开更多
This paper focuses on the adaptive discontinuous Galerkin(DG)methods for the tempered fractional(convection)diffusion equations.The DG schemes with interior penalty for the diffusion term and numerical flux for the co...This paper focuses on the adaptive discontinuous Galerkin(DG)methods for the tempered fractional(convection)diffusion equations.The DG schemes with interior penalty for the diffusion term and numerical flux for the convection term are used to solve the equations,and the detailed stability and convergence analyses are provided.Based on the derived posteriori error estimates,the local error indicator is designed.The theoretical results and the effectiveness of the adaptive DG methods are,respectively,verified and displayed by the extensive numerical experiments.The strategy of designing adaptive schemes presented in this paper works for the general PDEs with fractional operators.展开更多
Following the fractional cable equation established in the letter [B.I. Henry, T.A.M. Langlands, and S.L.Wearne, Phys. Rev. Lett. 100(2008) 128103], we present the time-space fractional cable equation which describes ...Following the fractional cable equation established in the letter [B.I. Henry, T.A.M. Langlands, and S.L.Wearne, Phys. Rev. Lett. 100(2008) 128103], we present the time-space fractional cable equation which describes the anomalous transport of electrodiffusion in nerve cells. The derivation is based on the generalized fractional Ohm's law;and the temporal memory effects and spatial-nonlocality are involved in the time-space fractional model. With the help of integral transform method we derive the analytical solutions expressed by the Green's function; the corresponding fractional moments are calculated; and their asymptotic behaviors are discussed. In addition, the explicit solutions of the considered model with two different external current injections are also presented.展开更多
To model wave propagation in inhomogeneous media with frequency dependent power-law attenuation,it is needed to use the fractional powers of symmetric coercive elliptic operators in space and the Caputo tempered fract...To model wave propagation in inhomogeneous media with frequency dependent power-law attenuation,it is needed to use the fractional powers of symmetric coercive elliptic operators in space and the Caputo tempered fractional derivative in time.The model studied in this paper is semilinear stochastic space-time fractional wave equations driven by infinite dimensional multiplicative Gaussian noise and additive fractional Gaussian noise,because of the potential fluctuations of the external sources.The purpose of this work is to discuss the Galerkin finite element approximation for the semilinear stochastic fractional wave equation.First,the space-time multiplicative Gaussian noise and additive fractional Gaussian noise are discretized,which results in a regularized stochastic fractional wave equation while introducing a modeling error in the mean-square sense.We further present a complete regularity theory for the regularized equation.A standard finite element approximation is used for the spatial operator,and a mean-square priori estimates for the modeling error and the approximation error to the solution of the regularized problem are established.Finally,numerical experiments are performed to confirm the theoretical analysis.展开更多
Fractional operators are widely used in mathematical models describing abnormal and nonlocal phenomena.Although there are extensive numerical methods for solving the corresponding model problems,theoretical analysis s...Fractional operators are widely used in mathematical models describing abnormal and nonlocal phenomena.Although there are extensive numerical methods for solving the corresponding model problems,theoretical analysis such as the regularity result,or the relationship between the left-side and right-side fractional operators is seldom mentioned.Instead of considering the fractional derivative spaces,this paper starts from discussing the image spaces of Riemann-Liouville fractional integrals of L_(p)(Ω) functions,since the fractional derivative operators that are often used are all pseudo-differential.Then the high regularity situation-the image spaces of Riemann-Liouville fractional integral operators on the W^(m,p)(Ω) space is considered.Equivalent characterizations of the defined spaces,as well as those of the intersection of the left-side and right-side spaces are given.The behavior of the functions in the defined spaces at both the nearby boundary point/points and the points in the domain is demonstrated in a clear way.Besides,tempered fractional operators are shown to be reciprocal to the corresponding Riemann-Liouville fractional operators,which is expected to contribute some theoretical support for relevant numerical methods.Last,we also provide some instructions on how to take advantage of the introduced spaces when numerically solving fractional equations.展开更多
We make the split of the integral fractional Laplacian as(−△)^(s)u=(−△)(−△)^(s−1)u,where s∈(0,1/2)∪(1/2,1).Based on this splitting,we respectively discretize the oneand two-dimensional integral fractional Laplaci...We make the split of the integral fractional Laplacian as(−△)^(s)u=(−△)(−△)^(s−1)u,where s∈(0,1/2)∪(1/2,1).Based on this splitting,we respectively discretize the oneand two-dimensional integral fractional Laplacian with the inhomogeneous Dirichlet boundary condition and give the corresponding truncation errors with the help of the interpolation estimate.Moreover,the suitable corrections are proposed to guarantee the convergence in solving the inhomogeneous fractional Dirichlet problem and an O(h^(1+α)2s))convergence rate is obtained when the solution u∈C^(1,α)(Ω_(n)^(δ)),where n is the dimension of the space,∈(max(0,2s−1),1],δis a fixed positive constant,and h denotes mesh size.Finally,the performed numerical experiments confirm the theoretical results.展开更多
The weak Galerkin （WG） finite element method was first introduced by Wang and Ye for solving second order elliptic equations, with the use of weak functions and their weak gradients. The basis function spaces depend...The weak Galerkin （WG） finite element method was first introduced by Wang and Ye for solving second order elliptic equations, with the use of weak functions and their weak gradients. The basis function spaces depend on different combinations of polynomial spaces in the interior subdomains and edges of elements, which makes the WG methods flexible and robust in many applications. Different from the definition of jump in discontinuous Galerkin （DG） methods, we can define a new weaker jump from weak functions defined on edges. Those functions have double values on the interior edges shared by two elements rather than a limit of functions defined in an element tending to its edge. Naturally, the weak jump comes from the difference between two weak flmctions defined on the same edge. We introduce an over-penalized weak Galerkin （OPWG） method, which has two sets of edge-wise and element-wise shape functions, and adds a penalty term to control weak jumps on the interior edges. Furthermore, optimal a priori error estimates in H1 and L2 norms are established for the finite element （Pk（K）, Pk（e）, RTk（K））. In addition, some numerical experiments are given to validate theoretical results, and an incomplete LU decomposition has been used as a preconditioner to reduce iterations from the GMRES, CO, and BICGSTAB iterative methods.展开更多
By combining the characteristicmethod and the local discontinuous Galerkin method with carefully constructing numerical fluxes,variational formulations are established for time-dependent incompressible Navier-Stokes e...By combining the characteristicmethod and the local discontinuous Galerkin method with carefully constructing numerical fluxes,variational formulations are established for time-dependent incompressible Navier-Stokes equations in R^(2).The nonlinear stability is proved for the proposed symmetric variational formulation.Moreover,for general triangulations the priori estimates for the L^(2)−norm of the errors in both velocity and pressure are derived.Some numerical experiments are performed to verify theoretical results.展开更多
基金supported by the National Key Research and Development Program of China(2020YFA0712900)the National Natural Science Foundation of China(12371093,12071197,12122102 and 12071431)+2 种基金the Key Project of Gansu Provincial National Science Foundation(23JRRA1022)the Fundamental Research Funds for the Central Universities(2233300008 and lzujbky-2021-ey18)the Innovative Groups of Basic Research in Gansu Province(22JR5RA391).
文摘Assume that L is a non-negative self-adjoint operator on L^(2)(ℝ^(n))with its heat kernels satisfying the so-called Gaussian upper bound estimate and that X is a ball quasi-Banach function space onℝ^(n) satisfying some mild assumptions.Let HX,L(ℝ^(n))be the Hardy space associated with both X and L,which is defined by the Lusin area function related to the semigroup generated by L.In this article,the authors establish various maximal function characterizations of the Hardy space HX,L(ℝ^(n))and then apply these characterizations to obtain the solvability of the related Cauchy problem.These results have a wide range of generality and,in particular,the specific spaces X to which these results can be applied include the weighted space,the variable space,the mixed-norm space,the Orlicz space,the Orlicz-slice space,and the Morrey space.Moreover,the obtained maximal function characterizations of the mixed-norm Hardy space,the Orlicz-slice Hardy space,and the Morrey-Hardy space associated with L are completely new.
基金supported by the Key Project of Gansu Provincial National Science Foundation(23JRRA1022)the National Natural Science Foundation of China(12071431)+1 种基金the Fundamental Research Funds for the Central Universities(lzujbky-2021-ey18)the Innovative Groups of Basic Research in Gansu Province(22JR5RA391).
文摘Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted Kato square root problem for L.More precisely,we prove that the square root L^(1/2)satisfies the weighted L^(p)estimates||L^(1/2)(f)||L_(ω)^p(R^(n))≤C||■f||L_(ω)^p(R^(n);R^(n))for any p∈(1,∞)andω∈Ap(ℝ^(n))(the class of Muckenhoupt weights),and that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,2+ε)andω∈Ap(ℝ^(n))∩RH_(2+ε/p),(R^(n))(the class of reverse Hölder weights),whereε∈(0,∞)is a constant depending only on n and the operator L,and where(2+ε/p)'denotes the Hölder conjugate exponent of 2+ε/p.Moreover,for any given q∈(2,∞),we give a sufficient condition to obtain that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,q)andω∈A_(p)(R^(n))∩pRH_(q/p),(R^(n)).As an application,we prove that when the coefficient matrix A that appears in L satisfies the small BMO condition,the Riesz transform∇L^(−1/2)is bounded on L_(ω)^(p)(ℝ^(n))for any given p∈(1,∞)andω∈Ap(ℝ^(n)).Furthermore,applications to the weighted L^(2)-regularity problem with the Dirichlet or the Neumann boundary condition are also given.
基金supported by the National Natural Science Foundation of China(11361053,11201204,11471148,11471330,145RJZA112)
文摘In this article, we consider positive steady state solutions and dynamics for a spatially heterogeneous predator-prey system with modified Leslie-Gower and Holling-Type II schemes. The heterogeneity here is created by the degeneracy of the intra-specific pressures for the prey. By the bifurcation method, the degree theory, and a priori estimates, we discuss the existence and multiplicity of positive steady states. Moreover, by the comparison argument, we also discuss the dynamical behavior for the diffusive predator-prey system.
基金the National Natural Science Foundation of China under Grant No.11671182the Fundamental Research Funds for the Central Universities under Grant No.lzujbky-2018-ot03.
文摘This paper provides a finite-difference discretization for the one-and two-dimensional tempered fractional Laplacian and solves the tempered fractional Poisson equation with homogeneous Dirichlet boundary conditions.The main ideas are to,respectively,use linear and quadratic interpolations to approximate the singularity and non-singularity of the one-dimensional tempered fractional Laplacian and bilinear and biquadratic interpolations to the two-dimensional tempered fractional Laplacian.Then,we give the truncation errors and prove the convergence.Numerical experiments verify the convergence rates of the order O(h^2−2s).
基金This work is supported by the National Natural Science Foundation of China nos.11971215 and 11571156,MOE-LCSMSchool of Mathematics and Statistics,Hunan Normal University,Changsha,Hunan 410081,China.
文摘Regularization methods have been substantially applied in image restoration due to the ill-posedness of the image restoration problem.Different assumptions or priors on images are applied in the construction of image regularization methods.In recent years,matrix low-rank approximation has been successfully introduced in the image denoising problem and significant denoising effects have been achieved.Low-rank matrix minimization is an NP-hard problem and it is often replaced with the matrix’s weighted nuclear norm minimization(WNNM).The assumption that an image contains an extensive amount of self-similarity is the basis for the construction of the matrix low-rank approximation-based image denoising method.In this paper,we develop a model for image restoration using the sum of block matching matrices’weighted nuclear norm to be the regularization term in the cost function.An alternating iterative algorithm is designed to solve the proposed model and the convergence analyses of the algorithm are also presented.Numerical experiments show that the proposed method can recover the images much better than the existing regularization methods in terms of both recovered quantities and visual qualities.
基金Supported by the Natural Science Foundation of China(Grant Nos.12271050,12201268)CAEP Foundation(Grant No.CX20200027)+2 种基金Key Laboratory of Computational Physics Foundation(Grant No.6142A05210502)Science and Technology Program of Gansu Province of China(Grant No.21JR7RA511)the National Science Foundation(DMS 1816313)。
文摘Separable multi-block convex optimization problem appears in many mathematical and engineering fields.In the first part of this paper,we propose an inertial proximal ADMM to solve a linearly constrained separable multi-block convex optimization problem,and we show that the proposed inertial proximal ADMM has global convergence under mild assumptions on the regularization matrices.Affine phase retrieval arises in holography,data separation and phaseless sampling,and it is also considered as a nonhomogeneous version of phase retrieval,which has received considerable attention in recent years.Inspired by convex relaxation of vector sparsity and matrix rank in compressive sensing and by phase lifting in phase retrieval,in the second part of this paper,we introduce a compressive affine phase retrieval via lifting approach to connect affine phase retrieval with multi-block convex optimization,and then based on the proposed inertial proximal ADMM for 3-block convex optimization,we propose an algorithm to recover sparse real signals from their(noisy)affine quadratic measurements.Our numerical simulations show that the proposed algorithm has satisfactory performance for affine phase retrieval of sparse real signals.
基金supported by the National Natural Science Foundation of China(Grant Nos.12071195,12301509,12225107)by the Innovative Groups of Basic Research in Gansu Province(Grant No.22JR5RA391)+3 种基金by the Major Science and Technology Projects in Gansu Province-Leading Talents in Science and Technology(Grant No.23ZDKA0005)by the Science and Technology Plan of Gansu Province(Grant No.22JR5RA535)by the Fundamental Research Funds for the Central Universities(Grant No.lzujbky-2023-pd04)by the China Postdoctoral Science Foundation(Grant No.2023M731466).
文摘In this paper,we consider the strong convergence of the time-space fractional diffusion equation driven by fractional Gaussian noise with Hurst index H∈(1/2,1).A sharp regularity estimate of the mild solution and the numerical scheme constructed by finite element method for integral fractional Laplacian and backward Euler convolution quadrature for Riemann-Liouville time fractional derivative are proposed.With the help of inverse Laplace transform and fractional Ritz projection,we obtain the accurate error estimates in time and space.Finally,our theoretical results are accompanied by numerical experiments.
基金supported by the Science Fund for Distinguished Young Scholars of Gansu Province(Grant No.23JRRA1020)the Fundamental Research Funds for the Central Universities(Grant No.lzujbky-2023-06).
文摘In this work,we analyze the three-step backward differentiation formula(BDF3)method for solving the Allen-Cahn equation on variable grids.For BDF2 method,the discrete orthogonal convolution(DOC)kernels are positive,the stability and convergence analysis are well established in[Liao and Zhang,Math.Comp.,90(2021),1207–1226]and[Chen,Yu,and Zhang,arXiv:2108.02910,2021].However,the numerical analysis for BDF3 method with variable steps seems to be highly nontrivial due to the additional degrees of freedom and the non-positivity of DOC kernels.By developing a novel spectral norm inequality,the unconditional stability and convergence are rigorously proved under the updated step ratio restriction rk:=τk/τk−1≤1.405 for BDF3 method.Finally,numerical experiments are performed to illustrate the theoretical results.To the best of our knowledge,this is the first theoretical analysis of variable steps BDF3 method for the Allen-Cahn equation.
基金Supported by the Program for New Century Excellent Talents in University under Grant No.NCET-09-0438the National Natural Science Foundation of China under Grant Nos.11271173 and 11101330+1 种基金the Starting Research Fund from the Xi’an University of Technology under Grant No.108-211206the Scientific Research Program Funded by Shaanxi Provincial Education Department under Grant No.2013JK0581
文摘The generalized fractional elastic models govern the stochastic motion of several many-body systems,e.g., polymers, membranes, and growing interfaces. This paper focuses on the exact formulations and their asymptotic behaviors of the average of the solutions of the generalized fractional elastic models. So we directly analyze the Cauchy problem of the averaged generalized elastic model involving time fractional derivative and the convolution integral of a radially symmetric friction kernel with space fractional Laplacian.
基金the National Natural Science Foundation of China under grant no.11671182the Fundamental Research Funds for the Central Universities under grants no.lzujbky-2018-ot03 and no.lzujbky 2019-it17.
文摘This paper focuses on the adaptive discontinuous Galerkin(DG)methods for the tempered fractional(convection)diffusion equations.The DG schemes with interior penalty for the diffusion term and numerical flux for the convection term are used to solve the equations,and the detailed stability and convergence analyses are provided.Based on the derived posteriori error estimates,the local error indicator is designed.The theoretical results and the effectiveness of the adaptive DG methods are,respectively,verified and displayed by the extensive numerical experiments.The strategy of designing adaptive schemes presented in this paper works for the general PDEs with fractional operators.
基金Supported by the Program for New Century Excellent Talents in University under Grant No.NCET-09-0438the National Natural Science Foundation of China under Grant No.11271173+2 种基金the training Program of the Major Research Plan of the National Natural Science Foundation of China under Grant No.91120014the Starting Research Foundation from the Xi’an University of Technology under GrantNo.108-211206the Scientific Research Foundation of the Education Department of Shaanxi Province under Grant No.2013JK0581
文摘Following the fractional cable equation established in the letter [B.I. Henry, T.A.M. Langlands, and S.L.Wearne, Phys. Rev. Lett. 100(2008) 128103], we present the time-space fractional cable equation which describes the anomalous transport of electrodiffusion in nerve cells. The derivation is based on the generalized fractional Ohm's law;and the temporal memory effects and spatial-nonlocality are involved in the time-space fractional model. With the help of integral transform method we derive the analytical solutions expressed by the Green's function; the corresponding fractional moments are calculated; and their asymptotic behaviors are discussed. In addition, the explicit solutions of the considered model with two different external current injections are also presented.
基金supported by the National Natural Science Foundation of China(Grants No.41875084,11801452,12071195,12225107)the AI and Big Data Funds(Grant No.2019620005000775)+1 种基金the Innovative Groups of Basic Research in Gansu Province(Grant No.22JR5RA391)NSF of Gansu(Grant No.21JR7RA537).
文摘To model wave propagation in inhomogeneous media with frequency dependent power-law attenuation,it is needed to use the fractional powers of symmetric coercive elliptic operators in space and the Caputo tempered fractional derivative in time.The model studied in this paper is semilinear stochastic space-time fractional wave equations driven by infinite dimensional multiplicative Gaussian noise and additive fractional Gaussian noise,because of the potential fluctuations of the external sources.The purpose of this work is to discuss the Galerkin finite element approximation for the semilinear stochastic fractional wave equation.First,the space-time multiplicative Gaussian noise and additive fractional Gaussian noise are discretized,which results in a regularized stochastic fractional wave equation while introducing a modeling error in the mean-square sense.We further present a complete regularity theory for the regularized equation.A standard finite element approximation is used for the spatial operator,and a mean-square priori estimates for the modeling error and the approximation error to the solution of the regularized problem are established.Finally,numerical experiments are performed to confirm the theoretical analysis.
基金supported by National Natural Science Foundation of China(Grant No.11801448)the Natural Science Basic Research Plan in Shaanxi Province of China(Grant No.2018JQ1022).supported by National Natural Science Foundation of China(Grant No.11271173).
文摘Fractional operators are widely used in mathematical models describing abnormal and nonlocal phenomena.Although there are extensive numerical methods for solving the corresponding model problems,theoretical analysis such as the regularity result,or the relationship between the left-side and right-side fractional operators is seldom mentioned.Instead of considering the fractional derivative spaces,this paper starts from discussing the image spaces of Riemann-Liouville fractional integrals of L_(p)(Ω) functions,since the fractional derivative operators that are often used are all pseudo-differential.Then the high regularity situation-the image spaces of Riemann-Liouville fractional integral operators on the W^(m,p)(Ω) space is considered.Equivalent characterizations of the defined spaces,as well as those of the intersection of the left-side and right-side spaces are given.The behavior of the functions in the defined spaces at both the nearby boundary point/points and the points in the domain is demonstrated in a clear way.Besides,tempered fractional operators are shown to be reciprocal to the corresponding Riemann-Liouville fractional operators,which is expected to contribute some theoretical support for relevant numerical methods.Last,we also provide some instructions on how to take advantage of the introduced spaces when numerically solving fractional equations.
基金supported by the National Natural Science Foundation of China(Grant No.12071195)the AI and Big Data Funds(Grant No.2019620005000775)+1 种基金by the Fundamental Research Funds for the Central Universities(Grant Nos.lzujbky-2021-it26,lzujbky-2021-kb15)NSF of Gansu(Grant No.21JR7RA537).
文摘We make the split of the integral fractional Laplacian as(−△)^(s)u=(−△)(−△)^(s−1)u,where s∈(0,1/2)∪(1/2,1).Based on this splitting,we respectively discretize the oneand two-dimensional integral fractional Laplacian with the inhomogeneous Dirichlet boundary condition and give the corresponding truncation errors with the help of the interpolation estimate.Moreover,the suitable corrections are proposed to guarantee the convergence in solving the inhomogeneous fractional Dirichlet problem and an O(h^(1+α)2s))convergence rate is obtained when the solution u∈C^(1,α)(Ω_(n)^(δ)),where n is the dimension of the space,∈(max(0,2s−1),1],δis a fixed positive constant,and h denotes mesh size.Finally,the performed numerical experiments confirm the theoretical results.
基金The research of the second author is supported in part by the Natural Science Foundation of Gansu Province, China （Grant 18JR3RA290）, and Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund （the second phase）.
文摘The weak Galerkin （WG） finite element method was first introduced by Wang and Ye for solving second order elliptic equations, with the use of weak functions and their weak gradients. The basis function spaces depend on different combinations of polynomial spaces in the interior subdomains and edges of elements, which makes the WG methods flexible and robust in many applications. Different from the definition of jump in discontinuous Galerkin （DG） methods, we can define a new weaker jump from weak functions defined on edges. Those functions have double values on the interior edges shared by two elements rather than a limit of functions defined in an element tending to its edge. Naturally, the weak jump comes from the difference between two weak flmctions defined on the same edge. We introduce an over-penalized weak Galerkin （OPWG） method, which has two sets of edge-wise and element-wise shape functions, and adds a penalty term to control weak jumps on the interior edges. Furthermore, optimal a priori error estimates in H1 and L2 norms are established for the finite element （Pk（K）, Pk（e）, RTk（K））. In addition, some numerical experiments are given to validate theoretical results, and an incomplete LU decomposition has been used as a preconditioner to reduce iterations from the GMRES, CO, and BICGSTAB iterative methods.
基金Thisworkwas partially supported by theNational Basic Research(973)Programof China under Grant 2011CB706903the National Natural Science Foundation of China under Grant 11271173,Grant 11471150,Grant 11671182,and the CAPES and CNPq in Brazil.
文摘By combining the characteristicmethod and the local discontinuous Galerkin method with carefully constructing numerical fluxes,variational formulations are established for time-dependent incompressible Navier-Stokes equations in R^(2).The nonlinear stability is proved for the proposed symmetric variational formulation.Moreover,for general triangulations the priori estimates for the L^(2)−norm of the errors in both velocity and pressure are derived.Some numerical experiments are performed to verify theoretical results.
