add separate functions for viscous dissipation