Allow the specification of the mean type when computing the bulk_viscosity

Allow users to specify which mean they want to use when applying the bulk_viscosity algorithm.

This has been implemented more generally so that all multimaterial averaging routines are capable of it regardless of rank but nothing other than bulk_viscosity actually exposes the functionality.

A (very) short test has been added.