For some reason I am unable to find a source on the internet about this. I think I have an answer, but I want to be doubly sure about this. All I could find (here), is that for an incompressible fluid, the viscous stress term is written in conservative form as: $$ \nabla\cdot\left[\mu\left(\nabla\boldsymbol{u}+(\nabla\boldsymbol{u})^{\top}\right)\right]. $$ Knowing that this viscous stress acts as diffusion, I assumed that I could simply replace $\boldsymbol{u}$ by any vectorial quantity and $\mu$ by any diffusion coefficient. However, when working this out in spherical coordinates for the MHD induction equation, I find that I should get instead: $$ \nabla\cdot\left[\eta\left(\nabla\boldsymbol{B}-(\nabla\boldsymbol{B})^{\top}\right)\right]. $$ I'm confused. I can't seem to find any mistake I made which would explain the minus sign. Is one of them wrong or are they somehow both right?
At the same time, the advection term in the Navier-Stokes equation reads: $$ -\nabla\cdot (\rho\boldsymbol{u}\otimes\boldsymbol{u}), $$ while in the MHD induction equation, we have instead: $$ -\nabla\cdot (\boldsymbol{u}\otimes\boldsymbol{B}-\boldsymbol{B}\otimes\boldsymbol{u}). $$ I'm not sure I understand why they are both different, if they represent the same physics (advection), even though the flow field obviously advects itself.