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

*(English)*Zbl 1179.65120Summary: The objective of this paper is to report on the development of a method of lines (MOL) toolbox 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.

##### MSC:

65M20 | Method of lines for initial value and initial-boundary value problems involving PDEs |

##### Keywords:

method of lines; moving grids; Burgers equation; Buckley-Leverett equation; Fisher-Kolmogorov equation; algorithm; finite differences; numerical results
PDF
BibTeX
XML
Cite

\textit{P. Saucez} et al., Appl. Math. Comput. 215, No. 5, 1821--1829 (2009; Zbl 1179.65120)

Full Text:
DOI

##### References:

[1] | Blom, J.G.; Zegeling, P.A., Algorithm 731: a moving-grid interface for systems of one-dimensional time-dependent partial differential equations, ACM transactions on mathematical software, 20, 194-214, (1994) · Zbl 0889.65099 |

[2] | Krysl, P., A pragmatic introduction to the finite element method for thermal and stress analysis with the Matlab toolkit sofea, (2006), World Scientific Publishing · Zbl 1128.80001 |

[3] | Kurganov, A.; Tadmor, E., New high-resolution central schemes for nonlinear conservation laws and convection-diffusion equations, Journal of computational physics, 160, 241-282, (2000) · Zbl 0987.65085 |

[4] | Madsen, N.R.; Sincovec, R.F., General software for partial differential equations, (), 236 |

[5] | Saucez, P.; Schiesser, W.E.; Vande Wouwer, A., Upwinding in the method of lines, Mathematics and computers in simulation, 56, 171-185, (2001) · Zbl 0977.68022 |

[6] | Shampine, L.F.; Gladwell, I.; Thompson, S., Solving ODEs with MATLAB, (2003), Cambridge University Press · Zbl 1040.65058 |

[7] | Trefethen, L.N., Spectral methods in MATLAB, (2000), SIAM Philadephia · Zbl 0953.68643 |

[8] | Verwer, J.G.; Blom, J.G.; Furzeland, R.M.; Zegeling, P.A., A moving-grid method for one-dimensional pdes based on the method of lines, (), 160-175 · Zbl 0697.65086 |

[9] | Weideman, J.A.C.; Reddy, S.C., A Matlab differentiation matrix suite, ACM transactions on mathematical software, 26, 465-519, (2000) |

[10] | Vande Wouwer, A.; Saucez, P.; Schiesser, W.E., Adaptive method of lines, (2001), Chapman and Hall/CRC · Zbl 0986.65083 |

[11] | Vande Wouwer, A.; Saucez, Ph.; Schiesser, W.E., Simulation of distributed parameter systems using a Matlab-based method of lines toolbox – chemical engineering applications, Industrial engineering and chemistry research, 43, 3469-3477, (2004) |

[12] | Vande Wouwer, A.; Saucez, P.; Schiesser, W.E.; Thompson, S., A Matlab implementation of upwind finite differences and adaptive grids in the method of lines, Journal of computational and applied mathematics, 183, 245-258, (2005) · Zbl 1071.65544 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.