A High Order Numerical Method for the Solution of the Advection Equation

Abstract

<p> This report presents a numerical method which can be used to solve the advection equation </p> <p> (∂ɸ/∂t) + (∂[u(x,t)ɸ]/∂x) = S(x,t) </p> where: </p> <p> ɸ ≣ concentration field </p> <p> u(x,t) ≣ velocity field </p> <p> S(x,t) ≣ source term </p> <p> Central to this method are the concept of particle path and the Eulerian interpretation of the time rate of change of the concentration field ɸ. </p> In actual comparison tests for particular cases with known solutions this method proved to be at least two orders of magnitude more accurate than the usual one sided upwind finite difference method. </p>