These operators are written in different forms in Cartesian, cylindrical and spherical coordinates. For instance, in spherical coordinate system, one has
$$\nabla \cdot \overrightarrow{F}=\frac{1}{r^{2}}\frac{\partial }{\partial r} \left( r^{2}F_{r}\right) +\frac{1}{r\sin \theta }\frac{\partial }{\partial \theta }\left( \sin \theta \cdot F_{\theta }\right) +\frac{1}{r\sin \theta } \frac{\partial F_{\varphi }}{\partial \varphi }.$$
Question: Do identities such as
$$\nabla \cdot \left( \overrightarrow{A}\times \overrightarrow{B}\right) =% \overrightarrow{B}\cdot \nabla \times \overrightarrow{A}-\overrightarrow{A}% \cdot \nabla \times \overrightarrow{B}$$
hold in general when cylindrical and spherical coordinate systems are used or do they have to be adapted?
Added: After having read the comments it occured to me that the invariance of these identities with regard to the coordinate system is a consequence of the definitions of the mentioned operators in terms of integrals, e.g.:
$$\nabla \cdot \overrightarrow{F}=\underset{V\rightarrow 0}{\lim }\frac{1}{V}% \underset{S}{\int \int }\overrightarrow{F}\cdot \overrightarrow{n}\;dA$$
where $V$ is the volume of a bounded closed region $T$, $S$ is the surface of $T$, and $\overrightarrow{n}$ the unit outer normal vector to $S$.