High-Order Local Discontinuous Galerkin Algorithm with Time Second-Order Schemes for the Two-Dimensional Nonlinear Fractional Diffusion Equation

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 θ.

Chebyshev Pseudo-Spectral Method for Solving Fractional Advection-Dispersion Equation

Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.

A Local Discontinuous Galerkin Method for Two-Dimensional Time Fractional Diffusion Equations

For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.

Stability and Convergence of an Implicit Difference Approximation for the Space Riesz Fractional Reaction-Dispersion Equation

In this paper,we consider a Riesz space-fractional reaction-dispersion equation (RSFRDE).The RSFRDE is obtained from the classical reaction-dispersion equation by replacing the second-order space derivative with a Riesz derivative of orderβ∈(1,2]. We propose an implicit finite difference approximation for RSFRDE.The stability and convergence of the finite difference approximations are analyzed.Numerical results are found in good agreement with the theoretical analysis.

Modeling contaminant transport in homogeneous porous media with fractional advection-dispersion equation

The newly developed Fractional Advection-Dispersion Equation (FADE), which is FADE was extended and used in this paper for modelling adsorbing contaminant transport by adding an adsorbing term. A parameter estimation method and its corresponding FORTRAN based program named FADEMain were developed on the basis of Nonlinear Least Square Algorithm and the analytical solution for one-dimensional FADE under the conditions of step input and steady state flow. Data sets of adsorbing contaminants Cd and NH4+-N transport in short homogeneous soil columns and conservative solute NaCI transport in a long homogeneous soil column, respectively were used to estimate the transport parameters both by FADEMain and the advection-dispersion equation (ADE) based program CXTFIT2.1. Results indicated that the concentration simulated by FADE agreed well with the measured data. Compared to the ADE model, FADE can provide better simulation for the concentration in the initial lower concentration part and the late higher concentration part of the breakthrough curves for both adsorbing contaminants. The dispersion coefficients for ADE were from 0.13 to 7.06 cm2/min, while the dispersion coefficients for FADE ranged from 0.119 to 3.05 cm1.856/min for NaCI transport in the long homogeneous soil column. We found that the dispersion coefficient of FADE increased with the transport distance, and the relationship between them can be quantified with an exponential function. Less scale-dependent was also found for the dispersion coefficient of FADE with respect to ADE.

A compact scheme for two-dimensional nonlinear time fractional wave equations

Based on the equivalent integro-differential form of the considered problem, a numerical approach to solving the two-dimensional nonlinear time fractional wave equations(NTFWEs) is considered in this paper. To this end, an alternating direction implicit(ADI) numerical scheme is derived. The scheme is established by combining the secondorder convolution quadrature formula and Crank–Nicolson technique in time and afourth-order difference approach in space. The convergence and unconditional stability of the proposed compact ADI scheme are strictly discussed after a concise solvabilityanalysis. A numerical example is shown to demonstrate the theoretical analysis.

IMPLICIT DIFFERENCE APPROXIMATION FOR A TIME FRACTIONAL ADVECTION-DISPERSION EQUATION

In this paper, a time fractional advection-dispersion equation is considered. From the relationship between the Caputo derivative and the Griinwald derivative, the Caputo derivative is approximated by using the Griinwald derivative. An implicit difference approximation for this equation is proposed. We prove that this approximation is unconditionally stable and convergent. Finally, numerical examples are given.

An efficient computational technique for time-fractional modified Degasperis-Procesi equation arising in propagation of nonlinear dispersive waves

In this paper,an efficient hybrid numerical scheme which is based on a joint venture of the q-homotopy analysis method and Sumudu transform is applied to investigate the time-fractional modified Degasperis-Procesi(DP)equation.The present study considers the Caputo fractional derivative.The fractional order modified DP model is very important and plays a great role in study of ocean engineering and science.The proposed scheme provides a beautiful opportunity for proper selection of the auxiliary parameter h and the asymptotic parameterρ(≥1)to handle mainly the differential equations of nonlinear nature.The offered scheme produces the solution in the shape of a convergent series in a large admissible domain which is helpful to regulate the region of convergence of a series solution.The proposed work computes the approximate analytical solution of the fractional modified DP equation systematically and also presents graphically the variation of the obtained solution for diverse values of the fractional parameterβ.

Investigation of adequate closed form travelling wave solution to the space-time fractional non-linear evolution equations

This work aims to construct exact solutions for the space-time fractional(2+1)-dimensional dispersive longwave(DLW)equation and approximate long water wave equation(ALW)utilizing the twovariable(G′/G,1/G)-expansion method and the modified Riemann-Liouville fractional derivative.The recommended equations play a significant role to describe the travel of the shallow water wave.The fractional complex transform is used to convert fractional differential equations into ordinary differential equations.Several wave solutions have been successfully achieved using the proposed approach and the symbolic computer Maple package.The Maple package program was used to set up and validate all of the computations in this investigation.By choosing particular values of the embedded parameters,we pro-duce multiple periodic solutions,periodic wave solutions,single soliton solutions,kink wave solutions,and more forms of soliton solutions.The achieved solutions might be useful to comprehend nonlinear phenomena.It is worth noting that the implemented method for solving nonlinear fractional partial dif-ferential equations(NLFPDEs)is efficient,and simple to find further and new-fangled solutions in the arena of mathematical physics and coastal engineering.

一类分数阶对流弥散方程差分方法

讨论了带有分数阶初边值问题的分数阶对流弥散方程,分别利用标准和移位的Grünwald-Letnikov分数阶算子离散方程以及边界条件中的Riemann-Liouville分数阶导数,并构造了相应的隐式有限差分格式和矩阵格式,证明了该差分格式的稳定性和收敛性。最后通过数值算例验证了其有效性。

Fundamental calculus of the fractional derivative defined with Rabotnov exponential kernel and application to nonlinear dispersive wave model

Before going further with fractional derivative which is constructed by Rabotnov exponential kernel,there exist many questions that are not addressed.In this paper,we try to recapitulate all the fundamental calculus,which we can obtain with this new fractional operator.The problems in this paper are to determine the solutions of the fractional differential equations where the second members are constant functions,polynomial functions,exponential functions,trigonometric functions,or Mittag-Leffler functions.For all the fractional differential equations,the obtained solutions are represented graphically.The Laplace transform of the fractional derivative with Rabotnov exponential kernel is the primary tool in the investigations.Finally,we give the fundamental solution to the nonlinear time-fractional modified Degasperis-Procesi equation by considering the fractional operator with Rabotnov exponential kernel.

考虑时空相关的分数阶对流—弥散方程及其解

本文在考虑弥散过程的时空相关性的基础上,用非局域性的处理方法,将二阶对流—弥散方程进行推广得到了分数阶的对流—弥散方程,方程中弥散项和对时间的导数被分数阶导数所代替。此方程的柯西问题的格林函数解是一分数稳定分布密度函数。由方程的稳定分布密度函数解说明了局域等效弥散系数与弥散过程有关,得出了等效弥散系数与运移尺度有关,是运移距离的幂函数的结论。这一结论从理论上解释了弥散系数的尺度效应。最后,用一实验的实测数据对所得结果进行检验,检验结果很好地说明了弥散过程中的偏态特征和“拖尾”现象,而传统二阶对流—弥散方程的高斯分布解却不能解释。因此,用分数阶的对流—弥散方程比二阶对流—弥散方程能更好的描述溶质在多孔介质中的弥散行为。

一维均质与非均质土柱中溶质迁移的分数微分对流-弥散模拟

分数微分对流-弥散方程(FADE)是模拟溶质迁移问题的新理论,但应用FADE来模拟溶质迁移时能否克服弥散的尺度效应尚待验证。利用长土柱实验资料结合FADE的解析解拟合推求FADE的弥散系数,并分析其与尺度之间的相关关系。研究结果表明,FADE的弥散系数具有随尺度增大而增大的现象,且均质土柱中FADE的弥散系数尺度效应小于非均质土柱中弥散系数尺度效应。在均质土柱中,弥散系数与尺度之间成指数相关关系,在非均质土柱中,弥散系数与尺度之间成幂相关关系。考虑了弥散系数分别与迁移时间和迁移距离呈线性递增两种相关关系,进而分别构建了3种考虑弥散尺度效应的FADE模型,并提出了求解的差分方法。利用上述3种考虑弥散尺度效应的FADE来模拟和预测不同空间位置处的溶质迁移过程。结果表明,对均质土柱中的溶质迁移可得到较好的模拟结果;对于非均质土柱,其模拟结果与实测结果仍然存在一定的差异。

时间分数阶色散方程的有限差分方法

提出求解时间分数阶色散方程的一类隐式差分格式,并证明其无条件稳定性和收敛性,收敛阶为O(τ+h2)。该分数阶色散方程是将一般的色散方程中的时间一阶导数用α(0<α<1)阶导数代替所得到的。数值算例表明本方法是有效的。

二维分数阶对流-弥散方程的数值解

对二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程分别建立了差分格式,实现了对其的数值求解。针对理想算例进行计算求解,分析了时间和空间分数阶阶数取不同值时的扩散变化规律,验证了各自所描述的时间相关性与空间相关性。同时与传统的二维整数阶对流-弥散方程的求解结果作了对比。当时间和空间分数阶阶数α与γ分别取整数时,二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程都与传统二维整数阶对流-弥散方程的计算结果相同,说明提出的对二维分数阶对流-弥散方程的数值求解方法是可行的。其结果对地下水溶质运移的进一步研究提供了有效的手段。

分数阶对流——弥散方程的数值求解

对严格的时间分数阶对流——弥散方程和严格的空间分数阶对流——弥散方程分别建立了差分格式,并用所建立的两个差分格式对同一理想算例进行了求解.通过对分数阶导数取不同的参数值,得到一系列结果,分析了不同分数阶导数描述的反常扩散现象及其变化规律,并和传统的整数阶对流——弥散方程的求解结果进行了对比.当时间分数阶对流——弥散方程和空间分数阶对流——弥散方程的分数阶导数的参数分别取整数值时,时间分数阶对流——弥散方程、空间分数阶对流——弥散方程和传统整数阶对流——弥散方程的计算结果相同,表明本文提出的对时间分数阶对流——弥散方程和空间对流——弥散方程数值求解方法是可行的,且整数阶对流——弥散方程是分数阶对流——弥散方程的特殊情况.和正常扩散相比,时间分数阶对流——弥散方程中分数阶导数的参数值越小,溶质扩散得越慢,表现为拖尾分布:空间分数阶对流——弥散方程中分数阶导数的参数值越小,溶质扩散得越快,表明空间的非局域性相关性越强.

基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移

时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质.

改进时间分数阶模型模拟非Fick溶质运移

通过将经典时间分数阶对流-弥散方程的等待时间分布函数的尾部修改为指数型,推导出了改进时间分数阶对流-弥散方程,并提出有效的时空算子分裂数值求解方法。对两个理想算例和一个实际算例进行计算,结果表明,改进的时间分数阶对流-弥散方程继承了时间分数阶对流-弥散方程能模拟穿透曲线幂率型拖尾分布的优点,还可模拟穿透曲线尾部由幂率型转换到指数型的过程;特征时间λ、分数阶指数γ和两相容量比例系数β共同决定了运移行为。改进的新模型可以区分非均质介质中流动相和非流动相中的溶质浓度,更细微地模拟非Fick溶质运移行为。

Caputo分数阶反应-扩散方程的隐式差分逼近

分数阶微分方程在许多应用科学上比整数阶微分方程更能准确地模拟自然现象.本文考虑分数阶反应-扩散方程.将一阶的时间偏导数用Caputo分数阶导数替换,并给出了一个隐式的差分格式.利用能量方法给出此差分格式的稳定性与收敛性证明,最后用数值例子说明差分格式是有效的.

Riesz分数阶反应-扩散方程数值近似的稳定性与收敛性分析

分数阶微分方程可以用来模拟工程,物理,生物等科学领域中的许多现象,然而分数阶微分方程的数值方法与理论分析是一项困难的事,其理论分析与经典的数值方法之间有很大的差异.本文考虑一个Riesz分数阶反应-扩散方程.这个方程是将一般的反应-扩散方程的二阶导用Riesz导数来替换.利用Riemann-Liouville定义和Grünwald-Letnikov定义之间的关系,我们提出了一个显示的数值近似,同时讨论了稳定性与收敛性,并给出数值例子.