For most mathematical equations,
analytical calculation of periodic travelling wave solutions is not possible, and therefore it is necessary to perform
numerical computations. For
partial differential equations, denote by
x and
t the (one-dimensional) space and time variables, respectively. Then periodic travelling waves are functions of the travelling wave variable
z=
x-
c t. Substituting this solution form into the partial differential equations gives a system of
ordinary differential equations known as the travelling wave equations. Periodic travelling waves correspond to
limit cycles of these equations, and this provides the basis for
numerical computations. The standard computational approach is
numerical continuation of the travelling wave equations. One first performs a continuation of a
steady state to locate a
Hopf bifurcation point. This is the starting point for a branch (family) of periodic travelling wave solutions, which one can follow by numerical continuation. In some (unusual) cases both end points of a branch (family) of periodic travelling wave solutions are
homoclinic solutions, in which case one must use an external starting point, such as a numerical solution of the partial differential equations. Periodic travelling wave
stability can also be calculated numerically, by computing the
spectrum. This is made easier by the fact that the spectrum of periodic travelling wave solutions of partial differential equations consists entirely of
essential spectrum. Possible numerical approaches include Hill's method and numerical continuation of the spectrum.
Software: The free,
open-source software package Wavetrain http://www.ma.hw.ac.uk/wavetrain is designed for the numerical study of periodic travelling waves. Using
numerical continuation, Wavetrain is able to calculate the form and stability of periodic travelling wave solutions of partial differential equations, and the regions of
parameter space in which waves exist and in which they are stable. == Applications ==