Numerical Simulations of Planetesimal Formation