Numerical modeling is used to investigate the physical mechanisms of the interaction of ultra-short (sub-picosecond) laser pulses with metallic targets. The laser–target interaction is modeled by using a one-dimensional hydrodynamic code that includes the absorption of laser radiation, the electronic heat conduction, the electron-phonon or electron–ion energy exchange, as well as a realistic equation of state. Laser fluences typical for micromachining are considered. The results of the 1D modeling are then used as the initial conditions for a 2D plasma expansion model. The dynamics of laser plume expansion in femtosecond regime is investigated. Calculations show that the plasma plume is strongly forward directed. In addition, a two-peaked axial density profile is obtained for 400 nm laser wavelength. The calculation results agree with the experimental observations.