Abstract

I previously used Burgers' equation to introduce a new method of spatial discretisation. The analysis is based upon centre manifold theory so we are assured that the discretisation accurately models all the processes and their subgrid scale interactions and may be constructed systematically. Here I show how boundaries to the physical domain may be naturally incorporated into the numerical modelling of Burgers' equation. We investigate Neumann and Dirichlet boundary conditions. As well as modelling the nonlinear advection, the method naturally derives symmetric matrices with constant bandwidth to correspond to the self-adjoint diffusion operator. The techniques developed here may be used to accurately model the nonlinear evolution of quite general spatio-temporal dynamical systems on bounded domains.

Full Paper (Size: 212 KB)