Multiresolution schemes for the numerical solution of 2-D conservation laws I

A generalization of Harten's multiresolution algorithms to two-dimensional (2-D) hyperbolic conservation laws is presented. Given a Cartesian grid and a discretized function on it, the method computes the local-scale components of the function by recursive diadic coarsening of the grid. Since the function's regularity can be described in terms of its scale or multiresolution analysis, the numerical solution of conservation laws becomes more efficient by eliminating flux computations wherever the solution is smooth. Instead, in those locations, the divergence of the solution is interpolated from the next coarser grid level. First, the basic 2-D essentially nonoscillatory (ENO) scheme is presented, then the 2-D multiresolution analysis is developed, and finally the subsequent scheme is tested numerically. The computational results confirm that the efficiency of the numerical scheme can be considerably improved in two dimensions as well.