Hello again,
This has been corrected since earlier this morning.
Since we have to consider distance from all the wires to the 'object' to be attracted, we have to think about the shape of the coil as well as how many turns and how much current.
The force varies as:
F=K1*I^2*N^2/L^2 (B varies as 1/L but the force as 1/L^2)
and with a current of 1 ampere we get:
F=K1*N^2/L^2
so the force goes up with the square of the turns N but down with the square of the distance L. This means we need to know the shape of the coil also, or determine it. But if we determine it, it probably wont be the shape we really want since we probably want something with a flat face. We also have to know what size object we intend to attract, or a range of object shapes we might be using this with.
For a simplification into 1 dimension along the axis of the coil, we see that the distance increases with the wire diameter d, but the wire diameter d only changes by:
d=2*sqrt(A1*N/pi)
where A1 is the initial wire cross sectional area, and normalizing for that area, and L=2*d, so we get:
F=K1*N/(16*pi)
so yes the force does increase along the axis when the coil diameter is small compared to the length.
But next we'd have to check the result for when the coil diameter is not small compared to the length (next time). In this case the force probably varies as the square root of the number of turns N cubed (N^(3/2)) (for a square coil with simplification into 1 dimension along the axis).
We could also check other shapes like a pie shape on a point target with the point of the pie facing the target.
Keep in mind that this optimization is for when we increase the number of turns N and at the same time double the cross sectional area of the wire. We are using more turns but bigger wire for each increase in N.