They shouldn't commute, that they seem like they should is mostly a superficial result of the notation.
In the following we rewrite these derivatives in a way such that it's manifest that they're not the same thing. The idea is simply to make sense of these in terms of operations that are a little more global than shifting around terms in expressions on a tangent space.
Claim
Let $(M, g)$ be a Riemannian manifold, $f \in C^\infty(M)$. On some local coordinate $x^\mu$ defined on open $U \in M$:
$$
\partial_v \partial^\mu f = d\left(~dx^\mu[\nabla f]~\right)[\partial_\nu]
$$
$$
\partial^\mu \partial_v f = dx^\mu[~\nabla(~df[\partial_v]~)~]
$$
where
- $d: C^\infty(M) \to \Omega^1(M) \equiv \Gamma(T^*M)$ ($\Gamma$ denote space of smooth sections in a bundle) is the De Rham differential.
- $\nabla: C^\infty(f) \to \Gamma(TM)$ is the usual gradient, defined by the property
$$
df[v] = \langle\nabla f, v\rangle
$$
Proof of the claim is just a direct computation writing the RHS of the expressions above out in components. A useful identity is the local expression for the gradient: $dx^\mu[\nabla f] = g^{\mu\nu}\partial_v f$ (which can prove as an exercise or look it up in any standard Riemannian geometry text).
On the sum If you're trying to compute a Laplacian, the correct version is $\partial_\mu \partial^\mu$, see, for example, the brief derivation of local formula for the Laplace-Beltrami operator on wikipedia.