A numerical model for laser ablation with plasma characterization