Create residual viscosity package +rv and generalize the ResidualViscosity class - Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update. - Create separate RungekuttaRV specifically using interior RV updates - Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class

classdef Definition
    methods (Abstract)

        % Returns a multiblock.Grid given some parameters
        g = getGrid(obj, varargin)

        % label is the type of label used for plotting,
        % default is block name, 'id' show the index for each block.
        show(obj, label, gridLines, varargin)
    end
end