most hydrodynamic simulations with dust evolution calculated only the total dust abundance without grain size distribution. [@2010MNRAS.403..620D] post-processed a cosmological hydrodynamic simulation with a dust evolution model which includes stellar dust production and SN destruction, and calculated the expected submillimetre fluxes for high-redshift star-forming galaxies. [@2015MNRAS.451..418Y], using a cosmological zoom-in simulation, calculated the dust distribution in high-redshift galaxies by simply assuming a constant dust-to-metal ratio. [@2015MNRAS.449.1625B] treated dust as a new particle species in addition to the gas, dark matter and star particles in their hydrodynamic simulations. In order to calculate the evolution of dust in galaxies, they included not only the formation and destruction of dust, but also the effect of dust on star formation and stellar feedback. [@2016MNRAS.457.3775M] performed cosmological zoom-in simulations by modelling dust as a component coupled with the gas. They included the processes relevant for the formation and evolution of dust, and revealed the importance of dust growth by the accretion of gas-phase metals. In addition, [@2017MNRAS.468.1505M] ran cosmological simulations to calculate the statistical properties of dust in galaxies. They broadly reproduced the dust abundances in the present-day galaxies, although they tended to underestimate those in high-redshift dusty galaxies. [@2016ApJ...831..147Z] analysed the dust evolution in 