The complexity of the phenomena simultaneously occurring, from the very first instants of high-power laser pulse interaction with the target up to the phase explosion, along with the strong changes in chemical-physical properties of matter, makes modeling laser ablation a hard task, especially near the thermodynamic critical regime. In this work, we report a computational model of an aluminum target irradiated in vacuum by a gaussian-shaped pulse of 20ns duration, with a peak intensity of the order of GW/cm□. This continuum model covers laser energy deposition and temperature evolution in the irradiated target, along with the mass removal mechanism involved, and the vaporized material expansion. Aluminum was considered to be a case study due to the vast literature on the temperature dependence of its thermodynamic, optical, and transport properties that were used to estimate time-dependent values of surface-vapor quantities (vapor pressure, vapor density, vapor and surface temperature) and vapor gas-dynamical quantities (density, velocity, pressure) as it expands into vacuum. Very favorable agreement is reported with experimental data regarding: mass removal and crater depth due to vaporization, generated recoil momentum, and vapor flow velocity expansion.