A review is presented on our recent Vlasov-Fokker-Planck (VFP) simulation code development and applications for high-power laser-plasma interactions. Numerical schemes are described for solving the kinetic VFP equation with both electron-electron and electron-ion collisions in one-spatial and two-velocity (1D2V) coordinates. They are based on the positive and flux conservation method and the finite volume method, and these two methods can insure the particle number conservation. Our simulation code can deal with problems in high-power laser/beam-plasma interactions, where highly non-Maxwellian electron distribution functions usually develop and the widely-used perturbation theories with the weak anisotropy assumption of the electron distribution function are no longer in point. We present some new results on three typical problems: firstly the plasma current generation in strong direct current electric fields beyond Spitzer-Härm’s transport theory, secondly the inverse bremsstrahlung absorption at high laser intensity beyond Langdon’s theory, and thirdly the heat transport with steep temperature and/or density gradients in laser-produced plasma. Finally, numerical parameters, performance, the particle number conservation, and the energy conservation in these simulations are provided.