Treffer: A MATLAB method of lines template for transport equations

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. [Copyright &y& Elsevier]

Copyright of Environmental Modelling & Software is the property of Elsevier B.V. and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)