Simulating the effect of high column density absorbers on the one-dimensional Lyman α forest flux power spectrum