The low temperature flexural creep properties of Mineral Powder Reinforced Asphalts (MPRAs) are systematically studied using both experimentations and finite element (FE) simulations. First, a new method to obtain the viscoelastic parameters of the pure asphalt is introduced. A digital twin for MPRA is then generated using Matlab, and a reliable simplification of the geometric dimension is adopted and verified. Second, a novel concept of equivalent particle size (EPS) is introduced in the FE models to represent the agglomeration effect of mineral powder in an asphalt. Four different basic EPSs are adopted to simulate the MPRAs with different powder-cement ratios. Third, the exponential function between EPS and simulated creep stiffness is fitted, and a reliable EPS for a specific sample is deduced using the fitted exponential function and experimental result, and it is verified using FE simulations. Finally, the digital twin for the MPRA proposed in this study is further verified by predicting the flexural creep properties at another temperature. The results show that the novel method to obtain the pure asphalt viscoelastic parameters and the digital twin for MPRA are reliable. EPS proposed in this paper could suitably characterize the agglomeration effect of mineral powder in FE simulations. The stress and strain inside the asphalt matrix decrease with smaller EPS. An accurate EPS for a specific MPRA can be deduced using experimental results and the exponential function fitted using different EPSs and simulated creep stiffness. Creep stiffness is larger under bigger powder-cement ratios, whereas the creep rate shows an opposite trend. Mineral powder tends to agglomerate more under bigger powder-cement ratios because the accurate EPS is bigger for MPRA having a bigger powder-cement ratio. The reliable digital twin for MPRA so obtained could be used to consistently predict the flexural creep properties at a different temperature.