Rapid transformations through the liquid and vapor phases induced by laser-solid interactions are described by our thermal model with the Clausius-Clapeyron equation to determine the vaporization temperature under different surface pressure condition. Hydrodynamic behavior of the vapor during and after ablation is described by gas dynamic equations. these two models are coupled. Modeling results show that lower background pressure results lower laser energy density threshold for vaporization. the ablation rate and the amount of materials removed are proportional to the laser energy density above its threshold. We also demonstrate a dynamic source effect that accelerates the unsteady expansion of laser-ablated material in the direction perpendicular to the solid. a dynamic partial ionization effect is studied as well. a self-similar theory shows that the maximum expansion velocity is proportional to cs/α, where 1 – α is the slope of the velocity profile. Numerical hydrodynamic modeling is in good agreement with the theory. With these effects, α. is reduced. therefore, the expansion front velocity is significantly higher than that from conventional models. the results are consistent with experiments. We further study the plume propagates in high background gas condition. Under appropriate conditions, the plume is slowed down, separates with the background, is backward moving, and hits the solid surface. then, it splits to be two parts when it rebounds from the surface. the results from the modeling will be compared with experimental observations where possible.