Ensuring numerical stability of wave propagation by tuning model parameters using genetic algorithms and response surface methods