Let's say I have the following nonlinear dynamical system:
\begin{align} \ x' &= -xy \\\\ \ y' &= -xy \\\\ \ z'&=xy \end{align}
For reference, this comes from a chemical system. Namely, it models concentration changes of the three chemical species when reactants X and Y combine to form product Z. Now, I know from conservation of mass that such a dynamical system does have a conserved quantity: X + Y + 2Z. I also know that for an arbitrary dynamical system to be conservative, its divergence has to be zero. In my case, I calculated my system's divergence out to be: \begin{align} \ div(X) &= -x-y \\\\ \end{align} This is clearly not zero except when both X and Y are zero or when X = -Y (but you can't have negative masses so this case is physically impossible). Yet, I know my system is conservative so how can it have this non-zero divergence? Is there something I'm missing?