Do you want to publish a course? Click here

Space-Fractional Diffusion with Variable Order and Diffusivity: Discretization and Direct Solution Strategies

195   0   0.0 ( 0 )
 Added by George Turkiyyah
 Publication date 2021
and research's language is English




Ask ChatGPT about the research

We consider the multidimensional space-fractional diffusion equations with spatially varying diffusivity and fractional order. Significant computational challenges are encountered when solving these equations due both to the kernel singularity in the fractional integral operator and to the resulting dense discretized operators, which quickly become prohibitively expensive to handle because of their memory and arithmetic complexities. In this work, we present a singularity-aware discretization scheme that regularizes the singular integrals through a singularity subtraction technique adapted to the spatial variability of diffusivity and fractional order. This regularization strategy is conveniently formulated as a sparse matrix correction that is added to the dense operator, and is applicable to different formulations of fractional diffusion equations. We also present a block low rank representation to handle the dense matrix representations, by exploiting the ability to approximate blocks of the resulting formally dense matrix by low rank factorizations. A Cholesky factorization solver operates directly on this representation using the low rank blocks as its atomic computational tiles, and achieves high performance on multicore hardware. Numerical results show that the singularity treatment is robust, substantially reduces discretization errors, and attains the first-order convergence rate allowed by the regularity of the solutions. They also show that considerable savings are obtained in storage ($O(N^{1.5})$) and computational cost ($O(N^2)$) compared to dense factorizations. This translates to orders-of-magnitude savings in memory and time on multi-dimensional problems, and shows that the proposed methods offer practical tools for tackling large nonlocal fractional diffusion simulations.



rate research

Read More

Time fractional PDEs have been used in many applications for modeling and simulations. Many of these applications are multiscale and contain high contrast variations in the media properties. It requires very small time step size to perform detailed computations. On the other hand, in the presence of small spatial grids, very small time step size is required for explicit methods. Explicit methods have many advantages as we discuss in the paper. In this paper, we propose a partial explicit method for time fractional PDEs. The approach solves the forward problem on a coarse computational grid, which is much larger than spatial heterogeneities, and requires only a few degrees of freedom to be treated implicitly. Via the construction of appropriate spaces and careful stability analysis, we can show that the time step can be chosen not to depend on the contrast or scale as the coarse mesh size. Thus, one can use larger time step size in an explicit approach. We present stability theory for our proposed method and our numerical results confirm the stability findings and demonstrate the performance of the approach.
In this paper, we develop a robust fast method for mobile-immobile variable-order (VO) time-fractional diffusion equations (tFDEs), superiorly handling the cases of small or vanishing lower bound of the VO function. The valid fast approximation of the VO Caputo fractional derivative is obtained using integration by parts and the exponential-sum-approximation method. Compared with the general direct method, the proposed algorithm ($RF$-$L1$ formula) reduces the acting memory from $mathcal{O}(n)$ to $mathcal{O}(log^2 n)$ and computational cost from $mathcal{O}(n^2)$ to $mathcal{O}(n log^2 n)$, respectively, where $n$ is the number of time levels. Then $RF$-$L1$ formula is applied to construct the fast finite difference scheme for the VO tFDEs, which sharp decreases the memory requirement and computational complexity. The error estimate for the proposed scheme is studied only under some assumptions of the VO function, coefficients, and the source term, but without any regularity assumption of the true solutions. Numerical experiments are presented to verify the effectiveness of the proposed method.
We investigate diffusion equations with time-fractional derivatives of space-dependent variable order. We examine the well-posedness issue and prove that the space-dependent variable order coefficient is uniquely determined among other coefficients of these equations, by the knowledge of a suitable time-sequence of partial Dirichlet-to-Neumann maps.
112 - Dongdong Hu , Xuenian Cao 2019
In this paper, a second-order backward difference formula (abbr. BDF2) is used to approximate first-order time partial derivative, the Riesz fractional derivatives are approximated by fourth-order compact operators, a class of new alternating-direction implicit difference scheme (abbr. ADI) is constructed for two-dimensional Riesz space fractional nonlinear reaction-diffusion equation. Stability and convergence of the numerical method are analyzed. Numerical experiments demonstrate that the proposed method is effective.
76 - Lei Li , Jian-Guo Liu 2019
We consider a discretization of Caputo derivatives resulted from deconvolving a scheme for the corresponding Volterra integral. Properties of this discretization, including signs of the coefficients, comparison principles, and stability of the corresponding implicit schemes, are proved by its linkage to Volterra integrals with completely monotone kernels. We then apply the backward scheme corresponding to this discretization to two time fractional dissipative problems, and these implicit schemes are helpful for the analysis of the corresponding problems. In particular, we show that the overdamped generalized Langevin equation with fractional noise has a unique limiting measure for strongly convex potentials and establish the convergence of numerical solutions to the strong solutions of time fractional gradient flows. The proposed scheme and schemes derived using the same philosophy can be useful for many other applications as well.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا