... | ... | @@ -24,13 +24,13 @@ We indicate with $`\kappa_p`$ and $`\kappa_r`$ respectively the Planck and t |
|
|
|
|
|
In order to compute the viscous heating the stress tensor is given by:
|
|
|
|
|
|
$`\bar{ \bar \tau} = 2\rho \nu(\bar{ \bar D} -{1\over 3}(\nabla \vec v)\bar {\bar I})`$
|
|
|
$`\bar{ \bar \tau} = 2\rho \nu(\bar{ \bar D} -{1\over 3}(\nabla \cdot \vec v)\bar {\bar I})`$
|
|
|
where $`\bar {\bar D}`$ is the strain tensor an $`\nu`$ the shear viscosity.
|
|
|
In spherical coordinates the components of the sress tensor
|
|
|
from (see Tassoul 1978) are:
|
|
|
|
|
|
$`\left\lbrace \begin{array}{lllll}
|
|
|
\tau _{RR} & & & = & 2\rho \nu ({\partial v_R \over \partial R}-{1\over 3}\nabla \cdot \vec v )\\
|
|
|
\tau _{RR} & & & = & 2\rho \nu ({\partial v_R \over \partial R}-{1\over 3}\nabla \cdot \vec v )\\
|
|
|
\tau _{\varphi \varphi} & & & = & 2\rho\nu ({1\over R}{\partial v_\varphi \over \partial \varphi}+{v_R\over R}-{1\over 3}\nabla \cdot \vec v )\\
|
|
|
\tau_{\theta \theta} & & & = & 2\rho\nu ({1\over R\sin \varphi}{\partial v_\theta \over \partial \theta}+{v_R\over R}+{v_\varphi \cot(\varphi)\over R} -{1\over 3}\nabla \cdot \vec v )\\
|
|
|
\tau _{R\varphi} & = & \tau _{\varphi R} & = & 2\rho\nu[{1\over 2}({1\over R}{\partial v_R \over \partial \varphi}+R{\partial \over \partial R}{v_{\varphi}\over R})] \\
|
... | ... | |