The cryocoolers are being developed for use in space and in terrestrial applications where combinations of long lifetime, high efficiency, compactness, low mass, low vibration, flexible interfacing, load variability, and reliability are essential. Pulse tube cryocoolers are now being used or considered for use in cooling infrared detectors for many space applications. In the development of these systems, as presented in this paper, first the system is analyzed theoretically. Based on the conservation of mass, the equation of motion, the conservation of energy, and the equation of state of a real gas a general model of the pulse tube refrigerator is made. The use of the harmonic approximation simplifies the differential equations of the model, as the time dependency can be solved explicitly and separately from the other dependencies. The model applies only to systems in the steady state. Time dependent effects, such as the cool down, are not described. From the relations the system performance is analyzed. And also we are describing pulse tube refrigeration mathematical models. There are three mathematical order models: first is analyzed enthalpy flow model and heat pumping flow model, second is analyzed adiabatic and isothermal model and third is flow chart of the computer program for numerical simulation. These mathematical reviews describe cryocoolers working and operation.