Can Mathematica solve nonlinear, coupled differential equations?

There were syntactic and conceptual problems with your formulation.

Conceptually, NDSolve is a numerical solver, so you need to specify boundary conditions as well as a numerical range of integration for the independent variable, which were missing in your formulation.

Syntactically, NDSolve was complaining about the fact that you had not specified the independent variable for the $u$ and $v$ functions each time. A human reader might be able to infer that from context, but a computer system has to have it all spelled out unequivocally.

Here is an example to show you how you might be able to set this up. Note that I have made no attempt to choose reasonable constraints and boundary conditions; I have just chosen numbers that would give a solution, to give you a jump start on the syntax.

Eq1 = Derivative[2][u][λ] - Cos[u[λ]] Sin[u[λ]] D[v[λ], λ] D[v[λ], λ] == 0
Eq2 = Derivative[2][v][λ] + 2 Cot[u[λ]] D[u[λ], λ] D[v[λ], λ] == 0
NDSolve[
  {Eq1, Eq2, 
   u[0] == 1000, u'[0] == 5, v'[0] == 5, v[0] == 1000}, 
  {u, v}, {λ, 1, 2}
]