A MATLAB method of lines template for transport equations

Abstract

Many environmental problems involve diffusion and convection processes, which can be described by partial differential equations (PDEs). This paper will describe the development of a MATLAB template that generates a numerical solution to PDEs using the method of lines. The template will be applied to various unsaturated flow problems within soil physics to demonstrate the versatility of the method. In particular, the template will generate solutions for three cases (1) one-dimensional Richards' equation for vertical infiltration; (2) coupled one-dimensional Richards' equation and solute transport equation for horizontal water and contaminant flow; and (3) two-dimensional Richard's equation for unsaturated flow over a complex geometry. Where possible, the results from the template will be compared against analytical solutions to determine the accuracy of the numerical solution. In addition, the paper will provide a discussion on possible extensions to the template and future directions.