Approximate solution of time-fractional advection-dispersion equation via fractional variational iteration method.
ABSTRACT: This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs.
Project description:In this article, a new general integral identity involving generalized fractional integral operators is established. With the help of this identity new Hermite-Hadamard type inequalities are obtained for functions whose absolute values of derivatives are convex. As a consequence, the main results of this paper generalize the existing Hermite-Hadamard type inequalities involving the Riemann-Liouville fractional integral.
Project description:A high-order finite difference scheme is proposed for solving time fractional heat equations. The time fractional derivative is described in the Riemann-Liouville sense. In the proposed scheme a new second-order discretization, which is based on Crank-Nicholson method, is applied for the time fractional part and fourth-order accuracy compact approximation is applied for the second-order space derivative. The spectral stability and the Fourier stability analysis of the difference scheme are shown. Finally a detailed numerical analysis, including tables, figures, and error comparison, is given to demonstrate the theoretical results and high accuracy of the proposed scheme.
Project description:In this work, a novel fractional order Coronavirus (2019-nCov) mathematical model with modified parameters is presented. The new fractional operator can be written as a linear combination of a Riemann-Liouville integral and a Caputo derivative. The suggested system is ruled by eight fractional-order nonlinear differential equations. The optimal control of the suggested model is the main objective of this work. Three control variables are presented in this model to minimize the number of infected population. Necessary control conditions are derived. Two schemes are constructed to simulate the proposed optimal control system. Prove of the schemes- stability are given. In order to validate the theoretical results numerical simulations and comparative studies with Caputo derivative are given.
Project description:In this article, some practical software optimization methods for implementations of fractional order backward difference, sum, and differintegral operator based on Grünwald–Letnikov definition are presented. These numerical algorithms are of great interest in the context of the evaluation of fractional-order differential equations in embedded systems, due to their more convenient form compared to Caputo and Riemann–Liouville definitions or Laplace transforms, based on the discrete convolution operation. A well-known difficulty relates to the non-locality of the operator, implying continually increasing numbers of processed samples, which may reach the limits of available memory or lead to exceeding the desired computation time. In the study presented here, several promising software optimization techniques were analyzed and tested in the evaluation of the variable fractional-order backward difference and derivative on two different Arm® Cortex®-M architectures. Reductions in computation times of up to 75% and 87% were achieved compared to the initial implementation, depending on the type of Arm® core.
Project description:By using Krasnoselskii's fixed point theorem, we study the existence of at least one or two positive solutions to a system of fractional boundary value problems given by -D(0+)(ν1)y1(t) = λ1a1(t)f(y1(t), y2(t)), - D(0+)(ν2)y2(t) = λ2a2(t)g(y1(t), y2(t)), where D(0+)(ν) is the standard Riemann-Liouville fractional derivative, ν1, ν2 ∈ (n - 1, n] for n > 3 and n ∈ N, subject to the boundary conditions y1((i))(0) = 0 = y ((i))(0), for 0 ≤ i ≤ n - 2, and [D(0+)(α)y1(t)] t=1 = 0 = [D(0+ (α)y2(t)] t=1, for 1 ≤ α ≤ n - 2, or y1((i))(0) = 0 = y ((i))(0), for 0 ≤ i ≤ n - 2, and [D(0+)(α)y1(t)] t=1 = ϕ1(y1), [D(0+)(α)y2(t)] t=1 = ϕ2(y2), for 1 ≤ α ≤ n - 2, ϕ1, ϕ2 ∈ C([0,1], R). Our results are new and complement previously known results. As an application, we also give an example to demonstrate our result.
Project description:The data presented in this paper are related to the paper entitled "A Numerical Method for Solving Caputo's and Riemann-Liouville's Fractional Differential Equations which includes multi-order fractional derivatives and variable coefficients", available in the "Communications in Nonlinear Science and Numerical Simulation" journal. Here, data are included for three of the four examples of Fractional Differential Equation (FDE) reported in , the other data is already available in . Data for each example contain: the interval of the solution, the solution by using the proposed method, the analytic solution and the absolute error. Data were obtained through Octave 5.1.0 simulations. For a better comprehension of the data, a pseudo-code of three stages and nine steps is included.
Project description:We consider Riemann mappings from bounded Lipschitz domains in the plane to a triangle. We show that in this case the Riemann mapping has a linear variational principle: It is the minimizer of the Dirichlet energy over an appropriate affine space. By discretizing the variational principle in a natural way we obtain discrete conformal maps which can be computed by solving a sparse linear system. We show that these discrete conformal maps converge to the Riemann mapping in [Formula: see text], even for non-Delaunay triangulations. Additionally, for Delaunay triangulations the discrete conformal maps converge uniformly and are known to be bijective. As a consequence we show that the Riemann mapping between two bounded Lipschitz domains can be uniformly approximated by composing the discrete Riemann mappings between each Lipschitz domain and the triangle.
Project description:We study instantaneous, mixing-driven, bimolecular equilibrium reactions in a system where transport is governed by a multidimensional space fractional dispersion equation. The superdiffusive, nonlocal nature of the system causes the location and magnitude of reactions that take place to change significantly from a classical Fickian diffusion model. In particular, regions where reaction rates would be zero for the Fickian case become regions where the maximum reaction rate occurs when anomalous dispersion operates. We also study a global metric of mixing in the system, the scalar dissipation rate and compute its asymptotic scaling rates analytically. The scalar dissipation rate scales asymptotically as t(-()(d)(+)(?)()/)(?) , where d is the number of spatial dimensions and ? is the fractional derivative exponent.
Project description:Isotropic linear and nonlinear fractional derivative constitutive relations are formulated and examined in terms of many parameter generalized Kelvin models and are analytically extended to cover general anisotropic homogeneous or non-homogeneous as well as functionally graded viscoelastic material behavior. Equivalent integral constitutive relations, which are computationally more powerful, are derived from fractional differential ones and the associated anisotropic temperature-moisture-degree-of-cure shift functions and reduced times are established. Approximate Fourier transform inversions for fractional derivative relations are formulated and their accuracy is evaluated. The efficacy of integer and fractional derivative constitutive relations is compared and the preferential use of either characterization in analyzing isotropic and anisotropic real materials must be examined on a case-by-case basis. Approximate protocols for curve fitting analytical fractional derivative results to experimental data are formulated and evaluated.
Project description:The goal of this paper is to examine the possible extension of the singular perturbation differential equation to the concept of fractional order derivative. To achieve this, we presented a review of the concept of fractional calculus. We make use of the Laplace transform operator to derive exact solution of singular perturbation fractional linear differential equations. We make use of the methodology of three analytical methods to present exact and approximate solution of the singular perturbation fractional, nonlinear, nonhomogeneous differential equation. These methods are including the regular perturbation method, the new development of the variational iteration method, and the homotopy decomposition method.