next up previous
Next: Periodic wave conduction block Up: Excitation wave breaking Previous: Plane waves in linear

Solitary wave blocking

We assume the following properties of the unperturbed version of this equation ( tex2html_wrap_inline824 :

  1.   There exists a stable solution in the form of a solitary wave, tex2html_wrap_inline826 , with width tex2html_wrap_inline828 , defined, for example, as the separation between points in which some component of v has a chosen value.
  2.   This solitary wave can develop from initial conditions in the form of this same wave laterally squeezed, i.e. tex2html_wrap_inline832 , k;SPMgt;1, with width, tex2html_wrap_inline836 , shorter than normal but longer than some minimal width, tex2html_wrap_inline838 ( tex2html_wrap_inline840 ), but if tex2html_wrap_inline842 , the wave decays.
  3.   The typical time for development/establishment of the wave profile is tex2html_wrap_inline844 . If, as is often the case, the excitation waves are supported by processes of different time scales, then tex2html_wrap_inline844 will be the slowest of them (i.e. the time constant of the limiting stage). The meaning of this parameter will be seen more clearly from its use below.
If tex2html_wrap_inline848 but is any positive constant, then assumption 1 implies that (4) has a stationary solitary wave solution tex2html_wrap_inline850 , with width tex2html_wrap_inline852 . If K(t) is not constant but changes slowly, we may expect that the solution will have the same form as this wave, slowly adjusting its width accordingly. If K(t) changes too rapidly, only then may the wave solution collapse.

With the assumptions listed above, this can be formalised in a phenomenological model. Let us assume, for simplicity, that the dynamics of the wave width is linear with constant relaxation time tex2html_wrap_inline844 . As the instantaneous equilibrium state should be changing in accordance with K(t), this gives

  equation111

This equation is in terms of the spatial variable tex2html_wrap_inline862 used in (4), which physically corresponds to the width measured in the direction of the flow. The rescaled width tex2html_wrap_inline864 , corresponding to real width measured perpendicularly to the wavefront, then obeys

equation125

The starting and the final asymptotic value of tex2html_wrap_inline866 is tex2html_wrap_inline868 . In between, it decreases below tex2html_wrap_inline868 but always remains positive. The minimal value of tex2html_wrap_inline866 is achieved at tex2html_wrap_inline874 which gives

  equation142

The system of equations (8,10) determines, in principle, the critical shear tex2html_wrap_inline876 and the corresponding time tex2html_wrap_inline878 of the break; however, the exact solution is rather tedious.

To obtain a simple analytical estimate, let us consider the case of practical interest, when the ratio of the activation timescale to the inhibition timescale is

equation154

Propagation of excitation waves in this limit is described by the Fife approximation [7]. In particular, the limiting stage of wave formation is establishment of its width, which happens on the time scale of tex2html_wrap_inline880 , and minimal possible width from which the wave can recover is of the order of the activation front width, tex2html_wrap_inline882 , whereas, in contrast, the normal wave width is, obviously, tex2html_wrap_inline884 .

In this limit we should expect that the breakup occurs only if the effective diffusivity changes very quickly compared with tex2html_wrap_inline886 . Assuming, in a self-consistent way, that tex2html_wrap_inline888 and neglecting the change of L(t) during this time interval, tex2html_wrap_inline892 , to leading order we obtain a system of two equations:

   eqnarray173

Using (5) in (12), tex2html_wrap_inline894 , then (13) implies tex2html_wrap_inline896 , and so the solution is

  eqnarray198

in agreement with the original assumptions[8].

Thus, in this case, there is a critical shear, above which the wave is quenched, and the value of this shear depends on the difference between the time scales of activator and inhibitor, which is a measure of the excitability of the system; the more excitable the system, the bigger the shear needed to destroy the wave.


next up previous
Next: Periodic wave conduction block Up: Excitation wave breaking Previous: Plane waves in linear

Vadim Biktashev
Wed Aug 12 18:01:29 BST 1998