Numerical simulation of femtosecond laser interaction with silicon