Parts’ processing on circular grinding machines is carried out according to the specified stepwise control cycles of radial and axial feeds, which are laid down in the control program (CP). Despite the precise positioning systems, modern CNC circular grinding machines have a variable compliance of the technological system along the processing length. It causes different cutting depth and occurrence of various errors in dimensions of the diametric surface. Also, the dimensions error increases due to the fact that the processing of parts batch is carried out in unstable conditions and under the influence of different variable factors, which can include allowance fluctuation, blunting of the wheel grains, changes of the wheel diameter and the contact area of the wheel with the workpiece, etc. However, there are still no means of control developed CP, which include cycles of cutting modes, to ensure the accuracy of parts batch` processing. This article describes the digital twin (DT) of the circular plunge grinding with CNC, which allows simulating layer-by-layer metal removal during the whole grinding cycle by calculating the cutting depth on each workpiece revolution under various combinations of unstable processing conditions and variable compliance of TS. Modeling of layer-by-layer allowance removal is performed in several sections of the treated surface; these sections have different compliance. At the end of the grinding cycle in the DT the processing accuracy is estimated on the basis of the obtained values of radii dimensions in all sections of the processed surface. The introduction of DT allows speeding up the CP development, reduce the time of launching CP in production, and increase the productivity of circular grinding operations with CNC by optimizing the cycles of the cutting modes while ensuring the specified processing accuracy.