Nonlinear Galerkin methods for a system of PDEs with Turing instabilities

We address and discuss the application of nonlinear Galerkin methods for the model reduction and numerical solution of partial differential equations (PDE) with Turing instabilities in comparison with standard (linear) Galerkin methods. The model considered is a system of PDEs modelling the pattern formation in vegetation dynamics. In particular, by constructing the approximate inertial manifold on the basis of the spectral decomposition of the solution, we implement the so-called Euler-Galerkin method and we compare its efficiency and accuracy versus the linear Galerkin methods. We compare the efficiency of the methods by (a) the accuracy of the computed bifurcation points, and, (b) by the computation of the Hausdorff distance between the limit sets obtained by the Galerkin methods and the ones obtained with a reference finite difference scheme. The efficiency with respect to the required CPU time is also accessed. For our illustrations we used three different ODE time integrators, from the Matlab ODE suite. Our results indicate that the performance of the Euler-Galerkin method is superior compared to the linear Galerkin method when either explicit or linearly implicit time integration scheme are adopted. For the particular problem considered, we found that the dimension of approximate inertial manifold is strongly affected by the lenght of the spatial domain. Indeeed, we show that the number of modes required to accurately describe the long time Turing pattern forming solutions increases as the domain increases.