Jacobian evolution of a smooth flow (${\partial }_{t}J[{\phi }_t]=(\nabla \cdot F)({\phi }_t)J[{\phi }_t]$)

For a smooth flow map ${\phi }_t$ generated by a vector field $F$, the chain rule and $\det (I+tDF)=1+t\mathrm{tr}(DF)+o(t)$ give ${\partial }_{t}J[{\phi }_t](x)=(\nabla \cdot F)({\phi }_t(x))J[{\phi }_t](x)$. With initial value one, its solution is the exponential of the time integral of this divergence. Thus a divergence-free generator preserves Euclidean volume and orientation. This connects the Jacobian determinant with a volume-preserving vector field.