In the present paper, we formulate a new mathematical model for the dynamics of COVID-19 with quarantine and isolation. Initially, we provide a brief discussion on the model formulation and provide relevant mathematical results. Then, we consider the fractal-fractional derivative in Atangana–Baleanu sense, and we also generalize the model. The generalized model is used to obtain its stability results. We show that the model is locally asymptotically stable if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathcal{R}_{0}<1$\end{document} . Further, we consider the real cases reported in China since January 11 till April 9, 2020. The reported cases have been used for obtaining the real parameters and the basic reproduction number for the given period, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathcal{R}_{0}\approx 6.6361$\end{document} . The data of reported cases versus model for classical and fractal-factional order are presented. We show that the fractal-fractional order model provides the best fitting to the reported cases. The fractional mathematical model is solved by a novel numerical technique based on Newton approach, which is useful and reliable. A brief discussion on the graphical results using the novel numerical procedures are shown. Some key parameters that show significance in the disease elimination from the society are explored.