Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn efficient differential quadrature method for fractional advection-diffusion equation
96
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This article studies a direct numerical approach for fractional advection-diffusion equations (ADEs). Using a set of cubic trigonometric B-splines as test functions, a differential quadrature (DQ) method is firstly proposed for the 1D and 2D time-fractional ADEs of order $(0,1]$. The weighted coefficients are determined, and with them, the original equation is transformed into a group of general ordinary differential equations (ODEs), which are discretized by an effective difference scheme or Runge-Kutta method. The stability is investigated under a mild theoretical condition. Secondly, based on a set of cubic B-splines, we develop a new Crank-Nicolson type DQ method for the 2D space-fractional ADEs without advection. The DQ approximations to fractional derivatives are introduced and the values of the fractional derivatives of B-splines are computed by deriving explicit formulas. The presented DQ methods are evaluated on five benchmark problems and the concrete simulations of the unsteady propagation of solitons and Gaussian pulse. In comparison with the existing algorithms in the open literature, numerical results finally illustrate the validity and accuracy.
rate research
Read More
In mathematical physics, the space-fractional diffusion equations are of particular interest in the studies of physical phenomena modelled by L{e}vy processes, which are sometimes called super-diffusion equations. In this article, we develop the differential quadrature (DQ) methods for solving the 2D space-fractional diffusion equations on irregular domains. The methods in presence reduce the original equation into a set of ordinary differential equations (ODEs) by introducing valid DQ formulations to fractional directional derivatives based on the functional values at scattered nodal points on problem domain. The required weighted coefficients are calculated by using radial basis functions (RBFs) as trial functions, and the resultant ODEs are discretized by the Crank-Nicolson scheme. The main advantages of our methods lie in their flexibility and applicability to arbitrary domains. A series of illustrated examples are finally provided to support these points.
In this article, we propose an exponential B-spline collocation method to approximate the solution of the fractional sub-diffusion equation of Caputo type. The present method is generated by use of the Gorenflo-Mainardi-Moretti-Paradisi (GMMP) scheme in time and an efficient exponential B-spline based method in space. The unique solvability is rigorously discussed. Its stability is well illustrated via a procedure closely resembling the classic von Neumann approach. The resulting algebraic system is tri-diagonal that can rapidly be solved by the known algebraic solver with low cost and storage. A series of numerical examples are finally carried out and by contrast to the other algorithms available in the literature, numerical results confirm the validity and superiority of our method.
In this paper, we provide a framework of designing the local discontinuous Galerkin scheme for integral fractional Laplacian $(-Delta)^{s}$ with $sin(0,1)$ in two dimensions. We theoretically prove and numerically verify the numerical stability and convergence of the scheme with the convergence rate no worse than $mathcal{O}(h^{k+frac{1}{2}})$.
In this paper, we present an embedded staggered discontinuous Galerkin method for the convection-diffusion equation. The new method combines the advantages of staggered discontinuous Galerkin (SDG) and embedded discontinuous Galerkin (EDG) method, and results in many good properties, namely local and global conservations, free of carefully designed stabilization terms or flux conditions and high computational efficiency. In applying the new method to convection-dominated problems, the method provides optimal convergence in potential and suboptimal convergence in flux, which is comparable to other existing DG methods, and achieves $L^2$ stability by making use of a skew-symmetric discretization of the convection term, irrespective of diffusivity. We will present numerical results to show the performance of the method.
In this paper, we focus on designing a well-conditioned Glarkin spectral methods for solving a two-sided fractional diffusion equations with drift, in which the fractional operators are defined neither in Riemann-Liouville nor Caputo sense, and its physical meaning is clear. Based on the image spaces of Riemann-Liouville fractional integral operators on $L_p([a,b])$ space discussed in our previous work, after a step by step deduction, three kinds of Galerkin spectral formulations are proposed, the final obtained corresponding scheme of which shows to be well-conditioned---the condition number of the stiff matrix can be reduced from $O(N^{2alpha})$ to $O(N^{alpha})$, where $N$ is the degree of the polynomials used in the approximation. Another point is that the obtained schemes can also be applied successfully to approximate fractional Laplacian with generalized homogeneous boundary conditions, whose fractional order $alphain(0,2)$, not only having to be limited to $alphain(1,2)$. Several numerical experiments demonstrate the effectiveness of the derived schemes. Besides, based on the numerical results, we can observe the behavior of mean first exit time, an interesting quantity that can provide us with a further understanding about the mechanism of abnormal diffusion.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
