Allowing users to specify which multimaterial mean to use in 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.