Wild guess: Are you after this?
Normal[Series[{
{-((U^2 w0 (CD w0 + CL\[Alpha] U ArcSin[Abs[w0]/U]))/(U^2 + w0^2)^(3/2))},
{0},
{(U^2 w0 (CD U - CL\[Alpha] w0 ArcSin[Abs[w0]/U]))/(U^2 + w0^2)^(3/2)}},
w0 -> 0, Assumptions -> w0 > 0 && U > 0]]
(* {{-(((CD + CL\[Alpha]) w0^2)/U)}, {0}, {CD w0}} *)
It's slightly different than the result your substitutions imply, so my guess may be wrong. However, asymptotic expansions are common operations, and you seem to be treating w0
as if it were small and approaching zero. I thought it worth suggesting.
In my experience, algebraic manipulation via substitution rules is tricky and bug-prone. Rules of the form variable -> expr
that replace all instances of the variable
by the same expr
are generally safe. But I find it's more often better to ask Mathematica to calculate the mathematical operation I want. There are exceptions, though.