The components of the exterior derivative of a p-form are
\left(d\tilde{\omega}\right)_{ki\cdots j}=\partial_{[k}\omega_{i\cdots j]}
Where the square brackets denote anti-symmetrisation. The second exterior derivative is then
\left[d\left(d\tilde{\omega}\right)\right]_{lki\cdots j}=\partial_{[l}\partial_{[k}\omega_{i\cdots j]]}
Here there is a double anti-symmetrisation. Now
\left[d\left(d\tilde{\omega}\right)\right]_{lki\cdots j}=-\partial_{[k}\partial_{[l}\omega_{i\cdots j]]}
By antisymmetry. But partials commute so
\left[d\left(d\tilde{\omega}\right)\right]_{lki\cdots j}=\partial_{[l}\partial_{[k}\omega_{i\cdots j]]}=-\partial_{[l}\partial_{[k}\omega_{i\cdots j]]}=0
For example, take the exterior derivative of a function f
(df)_i=\partial_if
Which is a one-form. Now take the second exterior derivative
\left[d(df)\right]_{ji}=\partial_j\partial_if-\partial_i\partial_jf=0
By commutativity of partials.
[I've missed the (p+1) etc. multiples off the exterior derivatives because I couldn't be arsed to keep putting them, as they're irrelevant. The factorials are accounted for in the anti-symmetrisation brackets if you're not familar with that notation.]
Oxymoron said:
d^2\omega = \frac{1}{r!}\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{\lambda}\partial x^{\nu}}\right)dx^{\lambda}\wedge dx^{\nu} \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}
I'll show you what I've just said but including covariant bases as you're doing.
d^2\omega = \frac{1}{r!}\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{\lambda}\partial x^{\nu}}\right)dx^{\lambda}\wedge dx^{\nu} \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}
= \frac{1}{r!}\left[\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{1}\partial x^{2}}\right)dx^{1}\wedge dx^{2}+\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{2}\partial x^{1}}\right)dx^{2}\wedge dx^{1}+\cdots \right] \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}
But dx^{2}\wedge dx^{1}=-dx^{1}\wedge dx^{2}, so
\frac{1}{r!}\left[\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{1}\partial x^{2}}\right)dx^{1}\wedge dx^{2} \wedge dx^{\mu_1}+\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{2}\partial x^{1}}\right)dx^{2}\wedge dx^{1}+\cdots \right] \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}
= \frac{1}{r!}\left[\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{1}\partial x^{2}}\right)dx^{1}\wedge dx^{2}-\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{2}\partial x^{1}}\right)dx^{1}\wedge dx^{2}+\cdots \right] \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}
Partials commute so
\frac{1}{r!}\left[\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{1}\partial x^{2}}\right)dx^{1}\wedge dx^{2}-\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{2}\partial x^{1}}\right)dx^{1}\wedge dx^{2}+\cdots \right] \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}
= \frac{1}{r!}\left[\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{1}\partial x^{2}}\right)dx^{1}\wedge dx^{2}-\left(\frac{\partial^2\omega_{\mu_1\dots\mu_r}}{\partial x^{1}\partial x^{2}}\right)dx^{1}\wedge dx^{2}+\cdots \right] \wedge dx^{\mu_1} \wedge \dots \wedge dx^{\mu_r}=0
You had all that you needed, you just didn't put it together. Unless what you're actually asking is how one proves the operator identity d^2=0. This follows from the above because d is defined by its operation on differential forms, so if it holds for a general p-form then one is allowed to just right the operator identity d^2=0. Or it can be taken with the other properties of the exterior derivative that you gave as axioms which together uniquely determine the form of the exterior derivative operator, however I've not seen this proof, but I of know.