# Nonlinear finite elements/Balance of angular momentum
## Statement of the balance of angular momentum

The balance of angular momentum in an inertial frame can be expressed as:

## Proof

We assume that there are no surface couples on $\partial {\Omega }$ or body couples in $\Omega$. Recall the general balance equation
    ${\cfrac {d}{dt}}\left[\int _{\Omega }f(\mathbf {x} ,t)~{\text{dV}}\right]=\int _{\partial {\Omega }}f(\mathbf {x} ,t)[u_{n}(\mathbf {x} ,t)-\mathbf {v} (\mathbf {x} ,t)\cdot \mathbf {n} (\mathbf {x} ,t)]~{\text{dA}}+\int _{\partial {\Omega }}g(\mathbf {x} ,t)~{\text{dA}}+\int _{\Omega }h(\mathbf {x} ,t)~{\text{dV}}~.$
In this case, the physical quantity to be conserved the angular momentum density, i.e., $f=\mathbf {x} \times (\rho ~\mathbf {v} )$. The angular momentum source at the surface is then $g=\mathbf {x} \times \mathbf {t}$ and the angular momentum source inside the body is $h=\mathbf {x} \times (\rho ~\mathbf {b} )$. The angular momentum and moments are calculated with respect to a fixed origin. Hence we have
    ${\cfrac {d}{dt}}\left[\int _{\Omega }\mathbf {x} \times (\rho ~\mathbf {v} )~{\text{dV}}\right]=\int _{\partial {\Omega }}[\mathbf {x} \times (\rho ~\mathbf {v} )][u_{n}-\mathbf {v} \cdot \mathbf {n} ]~{\text{dA}}+\int _{\partial {\Omega }}\mathbf {x} \times \mathbf {t} ~{\text{dA}}+\int _{\Omega }\mathbf {x} \times (\rho ~\mathbf {b} )~{\text{dV}}~.$
Assuming that $\Omega$ is a control volume, we have
    $\int _{\Omega }\mathbf {x} \times \left[{\cfrac {\partial }{\partial t}}(\rho ~\mathbf {v} )\right]~{\text{dV}}=-\int _{\partial {\Omega }}[\mathbf {x} \times (\rho ~\mathbf {v} )][\mathbf {v} \cdot \mathbf {n} ]~{\text{dA}}+\int _{\partial {\Omega }}\mathbf {x} \times \mathbf {t} ~{\text{dA}}+\int _{\Omega }\mathbf {x} \times (\rho ~\mathbf {b} )~{\text{dV}}~.$
Using the definition of a tensor product we can write
    $[\mathbf {x} \times (\rho ~\mathbf {v} )][\mathbf {v} \cdot \mathbf {n} ]=[[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]\cdot \mathbf {n} ~.$
Also, $\mathbf {t} ={\boldsymbol {\sigma }}\cdot \mathbf {n}$. Therefore we have
    $\int _{\Omega }\mathbf {x} \times \left[{\cfrac {\partial }{\partial t}}(\rho ~\mathbf {v} )\right]~{\text{dV}}=-\int _{\partial {\Omega }}[[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]\cdot \mathbf {n} ~{\text{dA}}+\int _{\partial {\Omega }}\mathbf {x} \times ({\boldsymbol {\sigma }}\cdot \mathbf {n} )~{\text{dA}}+\int _{\Omega }\mathbf {x} \times (\rho ~\mathbf {b} )~{\text{dV}}~.$
Using the divergence theorem, we get
    $\int _{\Omega }\mathbf {x} \times \left[{\cfrac {\partial }{\partial t}}(\rho ~\mathbf {v} )\right]~{\text{dV}}=-\int _{\Omega }{\boldsymbol {\nabla }}\bullet [[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]~{\text{dV}}+\int _{\partial {\Omega }}\mathbf {x} \times ({\boldsymbol {\sigma }}\cdot \mathbf {n} )~{\text{dA}}+\int _{\Omega }\mathbf {x} \times (\rho ~\mathbf {b} )~{\text{dV}}~.$
To convert the surface integral in the above equation into a volume integral, it is convenient to use index notation. Thus,
    $\left[\int _{\partial {\Omega }}\mathbf {x} \times ({\boldsymbol {\sigma }}\cdot \mathbf {n} )~{\text{dA}}\right]_{i}=\int _{\partial {\Omega }}e_{ijk}~x_{j}~\sigma _{kl}~n_{l}~{\text{dA}}=\int _{\partial {\Omega }}A_{il}~n_{l}~{\text{dA}}=\int _{\partial {\Omega }}{\boldsymbol {A}}\cdot \mathbf {n} ~{\text{dA}}$
where $[~]_{i}$ represents the $i$-th component of the vector. Using the divergence theorem
    $\int _{\partial {\Omega }}{\boldsymbol {A}}\cdot \mathbf {n} ~{\text{dA}}=\int _{\Omega }{\boldsymbol {\nabla }}\bullet {\boldsymbol {A}}~{\text{dV}}=\int _{\Omega }{\frac {\partial A_{il}}{\partial x_{l}}}~{\text{dV}}=\int _{\Omega }{\frac {\partial }{\partial x_{l}}}(e_{ijk}~x_{j}~\sigma _{kl})~{\text{dV}}~.$
Differentiating,
    $\int _{\partial {\Omega }}{\boldsymbol {A}}\cdot \mathbf {n} ~{\text{dA}}=\int _{\Omega }\left[e_{ijk}~\delta _{jl}~\sigma _{kl}+e_{ijk}~x_{j}~{\frac {\partial \sigma _{kl}}{\partial x_{l}}}\right]~{\text{dV}}=\int _{\Omega }\left[e_{ijk}~\sigma _{kj}+e_{ijk}~x_{j}~{\frac {\partial \sigma _{kl}}{\partial x_{l}}}\right]~{\text{dV}}=\int _{\Omega }\left[e_{ijk}~\sigma _{kj}+e_{ijk}~x_{j}~[{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}]_{k}\right]~{\text{dV}}~.$
Expressed in direct tensor notation,
    $\int _{\partial {\Omega }}{\boldsymbol {A}}\cdot \mathbf {n} ~{\text{dA}}=\int _{\Omega }\left[[{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}]_{i}+[\mathbf {x} \times ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})]_{i}\right]~{\text{dV}}$
where ${\mathcal {E}}$ is the third-order permutation tensor. Therefore,
    $\left[\int _{\partial {\Omega }}\mathbf {x} \times ({\boldsymbol {\sigma }}\cdot \mathbf {n} )~{\text{dA}}\right]_{i}==\int _{\Omega }\left[[{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}]_{i}+[\mathbf {x} \times ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})]_{i}\right]~{\text{dV}}$
or,
    $\int _{\partial {\Omega }}\mathbf {x} \times ({\boldsymbol {\sigma }}\cdot \mathbf {n} )~{\text{dA}}==\int _{\Omega }\left[{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}+\mathbf {x} \times ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})\right]~{\text{dV}}~.$
The balance of angular momentum can then be written as
    $\int _{\Omega }\mathbf {x} \times \left[{\cfrac {\partial }{\partial t}}(\rho ~\mathbf {v} )\right]~{\text{dV}}=-\int _{\Omega }{\boldsymbol {\nabla }}\bullet [[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]~{\text{dV}}+\int _{\Omega }\left[{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}+\mathbf {x} \times ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})\right]~{\text{dV}}+\int _{\Omega }\mathbf {x} \times (\rho ~\mathbf {b} )~{\text{dV}}~.$
Since $\Omega$ is an arbitrary volume, we have
    $\mathbf {x} \times \left[{\cfrac {\partial }{\partial t}}(\rho ~\mathbf {v} )\right]=-{\boldsymbol {\nabla }}\bullet [[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]+{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}+\mathbf {x} \times ({\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }})+\mathbf {x} \times (\rho ~\mathbf {b} )$
or,
    ${\mathbf {x} }\times {\left[{\frac {\partial }{\partial t}}(\rho ~\mathbf {v} )-{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}-\rho ~\mathbf {b} \right]}=-{\boldsymbol {\nabla }}\bullet [[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]+{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}~.$
Using the identity,
    ${\boldsymbol {\nabla }}\bullet (\mathbf {u} \otimes \mathbf {v} )=({\boldsymbol {\nabla }}\bullet \mathbf {v} )\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )\cdot \mathbf {v}$
we get
    ${\boldsymbol {\nabla }}\bullet [[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]=({\boldsymbol {\nabla }}\bullet \mathbf {v} )[\mathbf {x} \times (\rho ~\mathbf {v} )]+({\boldsymbol {\nabla }}[\mathbf {x} \times (\rho ~\mathbf {v} )])\cdot \mathbf {v} ~.$
The second term on the right can be further simplified using index notation as follows.
    ${\begin{aligned}\left[({\boldsymbol {\nabla }}[\mathbf {x} \times (\rho ~\mathbf {v} )])\cdot \mathbf {v} \right]_{i}=\left[({\boldsymbol {\nabla }}[\rho ~(\mathbf {x} \times \mathbf {v} )])\cdot \mathbf {v} \right]_{i}&={\frac {\partial }{\partial x_{l}}}(\rho ~e_{ijk}~x_{j}~v_{k})~v_{l}\\&=e_{ijk}\left[{\frac {\partial \rho }{\partial x_{l}}}~x_{j}~v_{k}~v_{l}+\rho ~{\frac {\partial x_{j}}{\partial x_{l}}}~v_{k}~v_{l}+\rho ~x_{j}~{\frac {\partial v_{k}}{\partial x_{l}}}~v_{l}\right]\\&=(e_{ijk}~x_{j}~v_{k})~\left({\frac {\partial \rho }{\partial x_{l}}}~v_{l}\right)+\rho ~(e_{ijk}~\delta _{jl}~v_{k}~v_{l})+e_{ijk}~x_{j}~\left(\rho ~{\frac {\partial v_{k}}{\partial x_{l}}}~v_{l}\right)\\&=[(\mathbf {x} \times \mathbf {v} )({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )+\rho ~\mathbf {v} \times \mathbf {v} +\mathbf {x} \times (\rho ~{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} )]_{i}\\&=[(\mathbf {x} \times \mathbf {v} )({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )+\mathbf {x} \times (\rho ~{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} )]_{i}~.\end{aligned}}$
Therefore we can write
    ${\boldsymbol {\nabla }}\bullet [[\mathbf {x} \times (\rho ~\mathbf {v} )]\otimes \mathbf {v} ]=(\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times ~\mathbf {v} )+({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )(\mathbf {x} \times \mathbf {v} )+\mathbf {x} \times (\rho ~{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} )]~.$
The balance of angular momentum then takes the form
    ${\mathbf {x} }\times {\left[{\frac {\partial }{\partial t}}(\rho ~\mathbf {v} )-{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}-\rho ~\mathbf {b} \right]}=-(\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times ~\mathbf {v} )-({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )(\mathbf {x} \times \mathbf {v} )-\mathbf {x} \times (\rho ~{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} )+{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}$
or,
    ${\mathbf {x} }\times {\left[{\frac {\partial }{\partial t}}(\rho ~\mathbf {v} )+\rho ~{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} -{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}-\rho ~\mathbf {b} \right]}=-(\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times ~\mathbf {v} )-({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )(\mathbf {x} \times \mathbf {v} )+{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}$
or,
    ${\mathbf {x} }\times {\left[\rho {\frac {\partial \mathbf {v} }{\partial t}}+{\frac {\partial \rho }{\partial t}}~\mathbf {v} +\rho ~{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} -{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}-\rho ~\mathbf {b} \right]}=-(\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times ~\mathbf {v} )-({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )(\mathbf {x} \times \mathbf {v} )+{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}$
The material time derivative of $\mathbf {v}$ is defined as
    ${\dot {\mathbf {v} }}={\frac {\partial \mathbf {v} }{\partial t}}+{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {v} ~.$
Therefore,
    ${\mathbf {x} }\times {\left[\rho ~{\dot {\mathbf {v} }}-{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}-\rho ~\mathbf {b} \right]}=-\mathbf {x} \times {\cfrac {\partial \rho }{\partial t}}~\mathbf {v} +-(\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times ~\mathbf {v} )-({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )(\mathbf {x} \times \mathbf {v} )+{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}~.$
Also, from the conservation of linear momentum
    $\rho ~{\dot {\mathbf {v} }}-{\boldsymbol {\nabla }}\bullet {\boldsymbol {\sigma }}-\rho ~\mathbf {b} =0~.$
Hence,
    ${\begin{aligned}0&=\mathbf {x} \times {\cfrac {\partial \rho }{\partial t}}~\mathbf {v} +(\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times ~\mathbf {v} )+({\boldsymbol {\nabla }}\rho \cdot \mathbf {v} )(\mathbf {x} \times \mathbf {v} )-{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}\\&=\left({\frac {\partial \rho }{\partial t}}+\rho {\boldsymbol {\nabla }}\bullet \mathbf {v} +{\boldsymbol {\nabla }}\rho \cdot \mathbf {v} \right)(\mathbf {x} \times \mathbf {v} )-{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}~.\end{aligned}}$
The material time derivative of $\rho$ is defined as
    ${\dot {\rho }}={\frac {\partial \rho }{\partial t}}+{\boldsymbol {\nabla }}\rho \cdot \mathbf {v} ~.$
Hence,
    $({\dot {\rho }}+\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} )(\mathbf {x} \times \mathbf {v} )-{\mathcal {E}}:{\boldsymbol {\sigma }}^{T}=0~.$
From the balance of mass
    ${\dot {\rho }}+\rho ~{\boldsymbol {\nabla }}\bullet \mathbf {v} =0~.$
Therefore,
    ${\mathcal {E}}:{\boldsymbol {\sigma }}^{T}=0~.$
In index notation,
    $e_{ijk}~\sigma _{kj}=0~.$
Expanding out, we get
    $\sigma _{12}-\sigma _{21}=0~;~~\sigma _{23}-\sigma _{32}=0~;~~\sigma _{31}-\sigma _{13}=0~.$
Hence,
    ${{\boldsymbol {\sigma }}={\boldsymbol {\sigma }}^{T}}$