The process of laser ablation of carbon in presence of background gas is simulated numerically. The plume dynamics in laser ablation is important to study for many reasons including temperature of plume particles and shielding of target by previously ablated plumes. Shielding leads directly to the change in energy deposition of incident laser pulse at the target surface and in turn influences the ablation dynamics and amount of material removed. Carbon ablation is studied for single and multiple laser hits typical for synthesis of nanotubes. Two models of correction of ablated velocity and pressure resulting from shielding effect are proposed and investigated. Numerical modeling of this plume dynamics and its integral effect of shielding is challenging due to inherent high nonlinearity of the problem. Some of available numerical techniques handles nonlinearity but are dissipative, e.g. Godunov type schemes. Other techniques are less dissipative but fail to account for strong nonlinearity typical for initial stages of ablation, e.g. the ENO-Roe. To effectively model this highly nonlinear plume dynamics a combination of two of above mentioned schemes is developed so as the numerical evaluation of fluxes is close to their physical values and the scheme has minimum dissipation. The non-monotonic behavior of ablated mass as a function of time duration between two laser pulses is studied. (C) 2007 Elsevier B.V. All rights reserved.