The objectives of this study are mainly two: (1) to validate whether a single grain scratch process can be modeled in two dimensions under the assumption of plane strain, and (2) to select the best discretization approach to model a single grain scratch process. This paper first focuses on the simulation of the orthogonal cutting process (aluminum alloy A2024 T351) using two mesh-based discretization approaches, the pure Lagrangian method (LAG) and the arbitrary Lagrangian–Eulerian method (ALE), and two particle-based approaches, the particle finite element method (PFEM) and smooth particle hydrodynamics (SPH), for both positive and negative rake angles. Benchmarking of the orthogonal cutting models at a rake angle of γ=20∘ is performed with the results of the process forces (cutting and passive forces) of a turning experiment from the literature. It is shown that all models are able to predict the cutting forces, but not the passive force. The orthogonal cutting model is further extended to simulate the cutting mechanism with negative rake tool geometries typically found in grinding and single grit scratching processes. The effects of the negative rake angles on the discretization approaches are studied. The calculated process forces are also compared to the measurements of the single grit scratch process performed at our laboratory. The 2D orthogonal cutting models significantly overestimate the process forces. One of the reasons why the orthogonal 2D cutting model is inadequate is that it cannot describe the complex mechanisms of material removal such as rubbing, plowing, and cutting. To account for these phenomena in LAG, ALE, and SPH discretization approaches, a 3D scratch model is developed. When comparing the process forces of the 3D model with the experimental measurements, all three discretization approaches show good agreement. However, it can be seen that the ALE model most closely matches the process forces with the experimental results. Finally, the ALE 3D scratch model was subjected to sensitivity analysis by changing the cutting speed, the depth of cut and the tool geometry. The results clearly show that the ALE method not only predicts the process forces and form the trends observed in the scratching experiments, but also predicts the scratch topography satisfactorily. Hence, we conclude that a 3D model is necessary to describe a scratch process and that the ALE method is the best discretization method.