(dm/dt)  =  (Subscript[ρ, air]) A ((Subscript[v, wind])  +  (Subscript[v, out]))/2

(Subscript[P, blade]) ==  
      1/2 (dm/dt) (((Subscript[v, wind])^2)  -  ((Subscript[v, out])^2))

(Subscript[P, turbine]) =  
      1/2 (1/2) A ((Subscript[v, out]) + 
            (Subscript[v, wind])) (Subscript[ρ, air]) (((Subscript[v, wind])^2)  -  
            ((Subscript[v, out])^2))

(Subscript[P, turbine]) = 
      (Subscript[P, turbine])  /.  {(Subscript[v, out])  ->  (Subscript[v, wind])/3}

(* 
      Compute the force on the turbine from decelerating wind *)

(Subscript[f, turbine]) = 
      (1/2) A ((Subscript[v, out]) + (Subscript[v, wind])) (Subscript[ρ, air]) ((Subscript[v, wind])  -  
              (Subscript[v, out]))  /.  {(Subscript[v, out])  ->  (Subscript[v, wind])/3}

(Subscript[f, drag]) =  50 * ((Subscript[v, b])^2)

(Subscript[v, wind])  =  (Subscript[v, b])  +  (Subscript[v, air])

(Subscript[P, prop ]) =  ((Subscript[f, drag])  +  (Subscript[f, turbine])) 
        (Subscript[v, b])

Simplify[(Subscript[P, prop ])]

Simplify[(Subscript[P, turbine])]

Plot[{(Subscript[P, turbine]),  (Subscript[P, prop ])} /.  {A -> 30, 
          (Subscript[v, air])  ->  10,  (Subscript[ρ, air]) ->  1.168},  {(Subscript[v, b]),  
        0,  10}]