Matlab implementation of a moving grid method based on the equidistribution principle

摘要

The objective of this paper is to report on the development of a method of lines (MOL) tool-box within MATLAB, and especially, on the implementation and test of a moving grid algorithm based on the equidistribution principle. This new implementation includes various spatial approximation schemes based on finite differences and slope limiters, the choice between several monitor functions, automatic grid adaptation to the initial condition, and provides a relatively easy tuning for the non-expert user. Several issues, including the sensitivity of the numerical results to the tuning parameters, are discussed. A few test problems characterized by solutions with steep moving fronts, including the Buckley-Leverett equation and an extended Fisher-Kolmogorov equation, are investigated so as to demonstrate the algorithm and software performance.