“…Although there has been extensive research on convergence and superconvergence of FEMs for various parabolic OCPs, mostly focused on linear or semilinear parabolic cases (see, e.g., [6,10,16,26,30]), the results on convergence and superconvergence are O(h + k) and O(h 3 2 + k), respectively. Recent years, VD are used to deal with different OCPs in [7,13,14].…”