The authors would like to acknowledge useful discussions with G. Y. Fu and V. S. Chan. Numerical computations were performed on the ShenMa High Performance Computing Cluster in Institute of Plasma Physics, Chinese Academy of Sciences. This work was supported by the National MCF Energy R&D Program (grant number 2018YFE0309101); the National Natural Science Foundation of China (grant number 12005100); and the Opening Project of Cooperative Innovation Center for Nuclear Fuel Cycle Technology and Equipment, University of South China (grant number 2019KFY15).