Simulations are carried out for varying mesh sizes, and the numerical solutions are compared by computing errors in two ways. In the case of an analytic solution being available, the errors in the numerical solutions are computed directly from the analytic solution. In the case of no availability of an analytic solution, an approximation to the error is computed using a very fine mesh numerical solution as the reference solution.
MSC:65M50, 35R37