Particle-resolved direct numerical simulations (PR-DNS) play an increasing role in investigating aerosol-cloud-turbulence interactions at the most fundamental level of processes. However, the high computational cost associated with high resolution simulations poses considerable challenges for large domain or long duration simulation using PR-DNS. To address these issues, here we present a digital twin model of the complex physics-based PR-DNS developed by use of the data-driven Fourier Neural Operator (FNO) method. The results demonstrate high accuracy at various resolutions and the digital twin model is two orders of magnitude cheaper in terms of computational demand compared to the physics-based PR-DNS model. Furthermore, the FNO digital-twin model exhibits strong generalization capabilities for different initial conditions and ultra-high-resolution without the need to retrain models. These findings highlight the potential of the FNO method as a promising tool to simulate complex fluid dynamics problems with high accuracy, computational efficiency, and generalization capabilities, enhancing our understanding of the aerosol-cloud-precipitation system.