A computational model for two-stage 4K-pulse tube cooler: Part I. Theoretical model and numerical method