Well, I think there won't be any closed-form solution for Δv or thrust involving the prop shape and pitch angle. It's just too messy - real aerodynamics. Maybe if those get condensed into a coefficient of lift, and you have RPM, but still messy. Is this an assigned problem or just something you're thinking about?
In the real world you either do complicated CFD simulations, or use rules of thumb (perhaps from software - try
http://www.drivecalc.de/PropCalc/PCHelp/Help.html ), or you use something off-the-shelf that has test data and specifications, or you just try things out with prototypes and measure the thrust. Usually you'll be limited in diameter that will fit on the plane, the engine power and RPM, and you'll know about what speed you will be going (e.g. static thrust). These, together with typical prop. design performance factors give a rule of thumb estimate of thrust.
If you have the ambient pressure and the pressure behind the prop, of course it's easy, just pressure differential times swept area.
***
F force, N
m mass, kg
a acceleration, m/s^2
ρ rho, density, kg/m^3
V volume, m^3
Δt time, a short interval, should really be dt
v velocity, m/s
A propeller swept area, m^2
π pi
D prop. diameter, m
Δv change in velocity, should really be dv
v_avg average v over the period that the prop. is accelerating the air
v_0 initial v of the incoming air
T thrust, N
v_f in my last post meaning the v of the air after it has been accelerated, not the final velocity of the aircraft.
Prove:
T= π/4*D^2*(v+Δv/2)*ρ*Δv
(substitutions indented)
F=ma (Newton)
m= ρ*V
V/Δt=v*A (volume per unit time equals the velocity of the incoming stream times the area)
| V=v*A*Δt
| m=ρ*v*A*Δt
A=π/4*D^2 (circular area expressed as diameter rather than radius)
| m=ρ*v*π/4*D^2*Δt
| F=ρ*v*π/4*D^2*Δt*a
a=Δv/Δt (definition of acceleration)
| F=ρ*v*π/4*D^2*Δt*Δv/Δt
| F=ρ*v*π/4*D^2*Δv
v_avg = v_0 +1/2*a*Δt (constant acceleration formula from every standard textbook, using Δt because of the short acceleration period and to get the Δt canceled out of the final equation while keeping the (v+Δv/2) term. Yes, it's sort of cheating a little.)
| = v_0 +1/2*Δv/Δt*Δt
| = v_0 + Δv/2
v= v_0 (more cheating to get the answer I want)
| v_avg = v + Δv/2
| F=ρ*(v + Δv/2)*π/4*D^2*Δv
F=T
So
| T=π/4*D^2*(v+Δv/2)*ρ*Δv
as stated in the original equation.
Edit: Here's a reference
http://www.scribd.com/doc/68114590/Journal-04-MDOpropulsion to show you just how complicated designing real propellers is.