Let's say I want to find the cross product of two gradient fields on a manifold $M$ in any co-ordinate system (not necessarily orthonormal). So I want to define a map $C: T_{p}M$ $\times$ $T_{p}M \to N_pM$ which I guess would be a rank $3$ tensor.
Say f and g are two linear functionals on the manifold, $g_{ij}$ is the metric tensor and $\epsilon$ is the Levi-Civita symbol.
Now, $grad$ $f$ = $g^{ij}$ $\partial_{i}f$ $\partial_{i}$
So, $C(grad$ $f$, $grad$ $c$ $)$ = $C(g^{kl}$ $\partial_{k}f$ $\partial_{l},$ $g^{mn}$ $\partial_{m}f$ $\partial_{n}$ $)$
The cross product is also defined as $*(df \wedge dc)$ = $*($ $\partial_{k}f$ $dx^k$ $\wedge$ $\partial_{m}c$ $dx^m$ $)$
Where $*$ is the hodge dual and $\wedge$ is the skew symmetric tensor product(wedge operator)
Also, according to https://en.wikipedia.org/wiki/Cross_product#Index_notation_for_tensors, it seems like the cross product is:
$(a \times b)^i = C(a, b)^i = g^{ij}\epsilon_{jkl}dx^kdx^l$
Is there any way to derive this general cross product in any coordinate frame? Perhaps according to the two definitions I have above? Or is there anything that I am doing wrong?
The usual cross product is only possible in three dimensions - its not possible in lower or higher dimensions; interestingly it makes the space into a lie algebra.
The generalisation to all (finite-dimensional) spaces is done by (differential) forms with the wedge as product.