Numerical modeling for CW multimode laser beams transformed by time-varying random phase plated is presented. The method is based on diffraction integral formula. The smoothing of multimode laser beams and the energy distribution dependence on the element size and time-varying speed are discussed in detial. The numerical results account for the experimental results satisfactorily.