Hi Furiosoreo,
Where has the need for 29 runs come from?
The way that a CCD will work is it will take the number of factors (4 in your case, so k=4) and then produce a number of runs based on the addition of a factorial portion (the 'corners'), the axial values (the extreme values) and the number of centrepoints:
Factorial portion = 2^k = 16
Axial portion = 2 x k = 8
Centrepoints = 2
Sum number of runs = 26
In your case, you can add more centrepoints to make it up to the 29 value. You will also want to consider what type of axial values you are using, in the case above I've mentioned axial values as 'extreme' (usually sticking outside of the cube made by the DoE), but they can on the face, orthogonal or 'rotatable'.
Rotatable is a popular choice because the final model has equal prediction variance, whereas a 'face' choice may be better when you can't go to the extreme values that the design demands.
Here's an old but gold resource from NIST that I used to learn about CCD's back in the day, I'd recommend you give it a read.
Thanks,
Ben
“All models are wrong, but some are useful”