Based on iterated Crank–Nicolson (CN) procedure, an alternative algorithm with perfectly matched layer (PML) formulation is proposed in the body‐of‐revolution (BOR) finite‐difference time‐domain (FDTD) lattice for the simulation of rotational symmetric geometrics. For the nonuniform domain simulation, an alternative subgridding method is employed to during the simulation. The iterated CN procedure improves the efficiency through preventing the calculation of tri‐diagonal matrices. The alternative subgridding method enhances the accuracy in nonuniform domains by the calculation of subregions. Numerical example is carried out for the demonstration of effectiveness including efficiency, accuracy and absorption. Through the results, the proposed scheme shows considerable absorption and accuracy improvement in nonuniform domains. Compared with the other CN schemes, the iterated CN procedure can significantly increase the efficiency with small time steps. In conclusion, the advantages and novelty of the proposed algorithm can be described as follows: (1) The iterated CN procedure is proposed for rotational symmetric geometrics. (2) Absorption boundary condition for iterated CN is proposed in BOR‐FDTD. (3) An alternative subgridding method for iterated CN procedure is proposed in BOR‐FDTD lattice. Thus, the proposed algorithm shows potential in nonuniform rotational symmetric geometrics open region simulation.