Wavy fronts in a hyperbolic FitzHugh-Nagumo system and the effects of cross diffusion

Oscillatory periodic pattern dynamics in hyperbolic reaction-advection-diffusion models

In this work we consider a quite general class of two-species hyperbolic reaction-advection-diffusion system with the main aim of elucidating the role played by inertial effects in the dynamics of oscillatory periodic patterns. To this aim, first, we use linear stability analysis techniques to deduce the conditions under which wave (or oscillatory Turing) instability takes place. Then, we apply multiple-scale weakly nonlinear analysis to determine the equation which rules the spatiotemporal evolution of pattern amplitude close to criticality. This investigation leads to a cubic complex Ginzburg-Landau (CCGL) equation which, owing to the functional dependence of the coefficients here involved on the inertial times, reveals some intriguing consequences. To show in detail the richness of such a scenario, we present, as an illustrative example, the pattern dynamics occurring in the hyperbolic generalization of the extended Klausmeier model. This is a simple two-species model used to describe the migration of vegetation stripes along the hillslope of semiarid environments. By means of a thorough comparison between analytical predictions and numerical simulations, we show that inertia, apart from enlarging the region of the parameter plane where wave instability occurs, may also modulate the key features of the coherent structures, solution of the CCGL equation. In particular, it is proven that inertial effects play a role, not only during transient regime from the spatially-homogeneous steady state toward the patterned state, but also in altering the amplitude, the wavelength, the angular frequency, and even the stability of the phase-winding solutions.

Diffusive instability in hyperbolic reaction-diffusion equation with different inertia

This work considers a two-dimensional hyperbolic reaction-diffusion system with different inertia and explores criteria for various instabilities, like a wave, Turing, and Hopf, both theoretically and numerically. It is proven that wave instability may occur in a two-species hyperbolic reaction-diffusion system with identical inertia if the diffusion coefficients of the species are nonidentical but cannot occur if diffusion coefficients are identical. Wave instability may also arise in a two-dimensional hyperbolic reaction-diffusion system if the diffusivities of the species are equal, which is never possible in a parabolic reaction-diffusion system, provided the inertias are different. Interestingly, Turing instability is independent of inertia, but the stability of the corresponding local system depends on the inertia. Theoretical results are demonstrated with an example where the local interaction is represented by the Schnakenberg system.

Multifront regime of a piecewise-linear FitzHugh-Nagumo model with cross diffusion

Oscillatory reaction-diffusion fronts are described analytically in a piecewise-linear approximation of the FitzHugh-Nagumo equations with linear cross-diffusion terms, which correspond to a pursuit-evasion situation. Fundamental dynamical regimes of front propagation into a stable and into an unstable state are studied, and the shape of the waves for both regimes is explored in detail. We find that oscillations in the wave profile may either be negligible due to rapid attenuation or noticeable if the damping is slow or vanishes. In the first case, we find fronts that display a monotonic profile of the kink type, whereas in the second case the oscillations give rise to fronts with wavy tails. Further, the oscillations may be damped with exponential decay or undamped so that a saw-shaped pattern forms. Finally, we observe an unexpected feature in the behavior of both types of the oscillatory waves: the coexistence of several fronts with different profile shapes and propagation speeds for the same parameter values of the model, i.e., a multifront regime of wave propagation.

Diffusive instabilities in a hyperbolic activator-inhibitor system with superdiffusion

We investigate analytically and numerically the conditions for wave instabilities in a hyperbolic activator-inhibitor system with species undergoing anomalous superdiffusion. In the present work, anomalous superdiffusion is modeled using the two-dimensional Weyl fractional operator, with derivative orders α∈ [1,2]. We perform a linear stability analysis and derive the conditions for diffusion-driven wave instabilities. Emphasis is placed on the effect of the superdiffusion exponent α, the diffusion ratio d, and the inertial time τ. As the superdiffusive exponent increases, so does the wave number of the Turing instability. Opposite to the requirement for Turing instability, the activator needs to diffuse sufficiently faster than the inhibitor in order for the wave instability to occur. The critical wave number for wave instability decreases with the superdiffusive exponent and increases with the inertial time. The maximum value of the inertial time for a wave instability to occur in the system is τ_{max}=3.6. As one of the main results of this work, we conclude that both anomalous diffusion and inertial time influence strongly the conditions for wave instabilities in hyperbolic fractional reaction-diffusion systems. Some numerical simulations are conducted as evidence of the analytical predictions derived in this work.

Incorporating inductances in tissue-scale models of cardiac electrophysiology

In standard models of cardiac electrophysiology, including the bidomain and monodomain models, local perturbations can propagate at infinite speed. We address this unrealistic property by developing a hyperbolic bidomain model that is based on a generalization of Ohm's law with a Cattaneo-type model for the fluxes. Further, we obtain a hyperbolic monodomain model in the case that the intracellular and extracellular conductivity tensors have the same anisotropy ratio. In one spatial dimension, the hyperbolic monodomain model is equivalent to a cable model that includes axial inductances, and the relaxation times of the Cattaneo fluxes are strictly related to these inductances. A purely linear analysis shows that the inductances are negligible, but models of cardiac electrophysiology are highly nonlinear, and linear predictions may not capture the fully nonlinear dynamics. In fact, contrary to the linear analysis, we show that for simple nonlinear ionic models, an increase in conduction velocity is obtained for small and moderate values of the relaxation time. A similar behavior is also demonstrated with biophysically detailed ionic models. Using the Fenton-Karma model along with a low-order finite element spatial discretization, we numerically analyze differences between the standard monodomain model and the hyperbolic monodomain model. In a simple benchmark test, we show that the propagation of the action potential is strongly influenced by the alignment of the fibers with respect to the mesh in both the parabolic and hyperbolic models when using relatively coarse spatial discretizations. Accurate predictions of the conduction velocity require computational mesh spacings on the order of a single cardiac cell. We also compare the two formulations in the case of spiral break up and atrial fibrillation in an anatomically detailed model of the left atrium, and we examine the effect of intracellular and extracellular inductances on the virtual electrode phenomenon.

Oscillatory pulses and wave trains in a bistable reaction-diffusion system with cross diffusion

We study waves with exponentially decaying oscillatory tails in a reaction-diffusion system with linear cross diffusion. To be specific, we consider a piecewise linear approximation of the FitzHugh-Nagumo model, also known as the Bonhoeffer-van der Pol model. We focus on two types of traveling waves, namely solitary pulses that correspond to a homoclinic solution, and sequences of pulses or wave trains, i.e., a periodic solution. The effect of cross diffusion on wave profiles and speed of propagation is analyzed. We find the intriguing result that both pulses and wave trains occur in the bistable cross-diffusive FitzHugh-Nagumo system, whereas only fronts exist in the standard bistable system without cross diffusion.

Diffusive instabilities in hyperbolic reaction-diffusion equations

We investigate two-variable reaction-diffusion systems of the hyperbolic type. A linear stability analysis is performed, and the conditions for diffusion-driven instabilities are derived. Two basic types of eigenvalues, real and complex, are described. Dispersion curves for both types of eigenvalues are plotted and their behavior is analyzed. The real case is related to the Turing instability, and the complex one corresponds to the wave instability. We emphasize the interesting feature that the wave instability in the hyperbolic equations occurs in two-variable systems, whereas in the parabolic case one needs three reaction-diffusion equations.