{\alpha} \angg{ \frac{\partial\eps}{\partial
k_{\alpha}}\varphi}   +m{\bf a}_{\alpha} \ang{
\frac{\partial\eps}{\partial k_{\alpha}}\dot{\varphi} } +
m{\bf a}_{\alpha} \angg{ \frac{\partial\dot{\eps}}{\partial
k_{\alpha}}\varphi }  + \dot{\bf g} .$$ Multiplying (62) by $-\dot{\bf u}\,$ and adding the result to the right hand side of (56) yield $$\begin{aligned}
\dot{E}  =  \dot{\bf u}\dot{\bf J} - \frac{1}{2}\dot{\rho}\dot{\bf u}^2 +
\sigma_{\alpha\beta}a_i^{\alpha}a_k^{\beta} \frac{\partial\dot{u}_i}{\partial
x_k}  - m\ddot{\bf  u} {\bf a}_{\alpha}
\angg{
\varphi\frac{\partial\eps}{\partial k_{\alpha}}}
  -  \dot{\bf  u}\nabla E_0
   +
   \dot{W}  -
\dot{\bf    u}\dot{\bf g} + \frac{\partial}{\partial t}
\langle \langle \eps\varphi\rangle \rangle .\end{aligned}$$ The last term is considered in Appendix 2. Substituting the time derivatives ${\dot{\bf J}}$ and ${\dot \rho}$ by means of (53) and (51) one obtains after cumbersome calculations: $$\begin{aligned}
\dot{E} & + & \nabla_k\left\{\frac{1}{2}\rho \dot{u}^2\dot{u}_k     +
\dot{u}_i\left(\Pi_{ik} - \rho\dot{u}_i\dot{u}_k  +  E_0\delta_{ik}  +
\langle \langle \eps f\rangle \rangle \delta_{ik}\right) \right.
\nonumber\\ & - &\left.\frac{1}{2}m\dot{u}^2
\ang{\frac{\partial\eps}{\partial p_k}f}  + \ang{\eps
\frac{\partial\epsilon}{\partial p_k}f} \right\}  \nonumber\\ & = &
\frac{\partial\dot{u}_i}{\partial  x_k}\left\{\Pi_{ik}        -
\rho\dot{u}_i\dot{u}_k + \sigma_{i k}
- \langle \lambda_{i k}  f\rangle    +
E_0\delta_{ik}   \right.\nonumber\\
&  -  &\left.  \dot{u}_i   j_{0i}   -\dot{u}_kj_{0i}
 \frac{}{}\right\} + \langle \eps\hat{I}f\rangle  +
\ang{{\bf F}\frac{\partial\eps}{\partial {\bf p}}f}  +
\dot{W} - \dot{\bf u}\dot{\bf g}\end{aligned}$$ where $$\lambda_{ik} = 2 \frac{\p \eps}{\p g^{\alpha\beta}}a^\alpha_i a^\beta_k
\qquad
\sigma_{ik} = -2 \frac{\p E_0}{\p g^{\alpha\beta}}a^\alpha_i a^\beta_k .
\ee             %65
     The last three terms 