77 research outputs found
A fractional B-spline collocation method for the numerical solution of fractional predator-prey models
We present a collocation method based on fractional B-splines for the solution of fractional differential problems. The key-idea is to use the space generated by the fractional B-splines, i.e., piecewise polynomials of noninteger degree, as approximating space. Then, in the collocation step the fractional derivative of the approximating function is approximated accurately and efficiently by an exact differentiation rule that involves the generalized finite difference operator. To show the effectiveness of the method for the solution of nonlinear dynamical systems of fractional order, we solved the fractional Lotka-Volterra model and a fractional predator-pray model with variable coefficients. The numerical tests show that the method we proposed is accurate while keeping a low computational cost
Efficient multistep methods for tempered fractional calculus: Algorithms and Simulations
In this work, we extend the fractional linear multistep methods in [C.
Lubich, SIAM J. Math. Anal., 17 (1986), pp.704--719] to the tempered fractional
integral and derivative operators in the sense that the tempered fractional
derivative operator is interpreted in terms of the Hadamard finite-part
integral. We develop two fast methods, Fast Method I and Fast Method II, with
linear complexity to calculate the discrete convolution for the approximation
of the (tempered) fractional operator. Fast Method I is based on a local
approximation for the contour integral that represents the convolution weight.
Fast Method II is based on a globally uniform approximation of the trapezoidal
rule for the integral on the real line. Both methods are efficient, but
numerical experimentation reveals that Fast Method II outperforms Fast Method I
in terms of accuracy, efficiency, and coding simplicity. The memory requirement
and computational cost of Fast Method II are and ,
respectively, where is the number of the final time steps and is the
number of quadrature points used in the trapezoidal rule. The effectiveness of
the fast methods is verified through a series of numerical examples for
long-time integration, including a numerical study of a fractional
reaction-diffusion model
Good (and Not So Good) practices in computational methods for fractional calculus
The solution of fractional-order differential problems requires in the majority of cases the use of some computational approach. In general, the numerical treatment of fractional differential equations is much more difficult than in the integer-order case, and very often non-specialist researchers are unaware of the specific difficulties. As a consequence, numerical methods are often applied in an incorrect way or unreliable methods are devised and proposed in the literature. In this paper we try to identify some common pitfalls in the use of numerical methods in fractional calculus, to explain their nature and to list some good practices that should be followed in order to obtain correct results
Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for
the study of astrophysical plasma dynamics. The code incorporates an
energy-conserving particle integrator into a stable, second-order--accurate,
three-stage predictor-predictor-corrector integration algorithm. The
constrained transport method is used to enforce the divergence-free constraint
on the magnetic field. A delta-f scheme is included to facilitate a
reduced-noise study of systems in which only small departures from an initial
distribution function are anticipated. The effects of rotation and shear are
implemented through the shearing-sheet formalism with orbital advection. These
algorithms are embedded within an architecture similar to that used in the
popular astrophysical magnetohydrodynamics code Athena, one that is modular,
well-documented, easy to use, and efficiently parallelized for use on thousands
of processors. We present a series of tests in one, two, and three spatial
dimensions that demonstrate the fidelity and versatility of the code.Comment: 27 pages, 12 figures, accepted for publication in Journal of
Computational Physic
Non-Linear Langevin and Fractional Fokker-Planck Equations for Anomalous Diffusion by Levy Stable Processes
The~numerical solutions to a non-linear Fractional Fokker--Planck (FFP)
equation are studied estimating the generalized diffusion coefficients. The~aim
is to model anomalous diffusion using an FFP description with fractional
velocity derivatives and Langevin dynamics where L\'{e}vy fluctuations are
introduced to model the effect of non-local transport due to fractional
diffusion in velocity space. Distribution functions are found using numerical
means for varying degrees of fractionality of the stable L\'{e}vy distribution
as solutions to the FFP equation. The~statistical properties of the
distribution functions are assessed by a generalized normalized expectation
measure and entropy and modified transport coefficient. The~transport
coefficient significantly increases with decreasing fractality which is
corroborated by analysis of experimental data.Comment: 20 pages 7 figure
A General Return-Mapping Framework for Fractional Visco-Elasto-Plasticity
We develop a fractional return-mapping framework for power-law
visco-elasto-plasticity. In our approach, the fractional viscoelasticity is
accounted through canonical combinations of Scott-Blair elements to construct a
series of well-known fractional linear viscoelastic models, such as
Kelvin-Voigt, Maxwell, Kelvin-Zener and Poynting-Thomson. We also consider a
fractional quasi-linear version of Fung's model to account for stress/strain
nonlinearity. The fractional viscoelastic models are combined with a fractional
visco-plastic device, coupled with fractional viscoelastic models involving
serial combinations of Scott-Blair elements. We then develop a general
return-mapping procedure, which is fully implicit for linear viscoelastic
models, and semi-implicit for the quasi-linear case. We find that, in the
correction phase, the discrete stress projection and plastic slip have the same
form for all the considered models, although with different property and
time-step dependent projection terms. A series of numerical experiments is
carried out with analytical and reference solutions to demonstrate the
convergence and computational cost of the proposed framework, which is shown to
be at least first-order accurate for general loading conditions. Our numerical
results demonstrate that the developed framework is more flexible, preserves
the numerical accuracy of existing approaches while being more computationally
tractable in the visco-plastic range due to a reduction of in CPU time.
Our formulation is especially suited for emerging applications of fractional
calculus in bio-tissues that present the hallmark of multiple viscoelastic
power-laws coupled with visco-plasticity
- …