  A numerical study of the role of anomalous diffusion in front propagation in
reaction-diffusion systems is presented. Three models of anomalous diffusion
are considered: fractional diffusion, tempered fractional diffusion, and a
model that combines fractional diffusion and regular diffusion. The reaction
kinetics corresponds to a Fisher-Kolmogorov nonlinearity. The numerical method
is based on a finite-difference operator splitting algorithm with an explicit
Euler step for the time advance of the reaction kinetics, and a Crank-Nicholson
semi-implicit time step for the transport operator. The anomalous diffusion
operators are discretized using an upwind, flux-conserving, Grunwald-Letnikov
finite-difference scheme applied to the regularized fractional derivatives.
With fractional diffusion of order $\alpha$, fronts exhibit exponential
acceleration, $a_L(t) \sim e^{\gamma t/\alpha}$, and develop algebraic decaying
tails, $\phi \sim 1/x^{\alpha}$. In the case of tempered fractional diffusion,
this phenomenology prevails in the intermediate asymptotic regime
  $\left(\chi t \right)^{1/\alpha} \ll x \ll 1/\lambda$, where $1/\lambda$ is
the scale of the tempering. Outside this regime, i.e. for $x > 1/\lambda$, the
tail exhibits the tempered decay $\phi \sim e^{-\lambda x}/x^{\alpha+1}$, and
the front velocity approaches the terminal speed $v_*=
\left(\gamma-\lambda^\alpha \chi\right)/ \lambda$. Of particular interest is
the study of the interplay of regular and fractional diffusion. It is shown
that the main role of regular diffusion is to delay the onset of front
acceleration. In particular, the crossover time, $t_c$, to transition to the
accelerated fractional regime exhibits a logarithmic scaling of the form $t_c
\sim \log \left(\chi_d/\chi_f\right)$ where $\chi_d$ and $\chi_f$ are the
regular and fractional diffusivities.
