The presented thesis deals with a numerical simulation of axially cooled glassmaking parison. Several types of cooling channels have been investigated according to its usage in glassmaking process. Software Star CCM+ was used to describe properties throughout the channels and also the parison itself. The main focus of this thesis lies in comparing of temperature fields, cooling effects and cooling intensity among the chosen channel designs. The boundary condi-tions were kept the same accordingly to glassmaking conditions for every design. Some out-puts were transferred to non-dimensional diagrams in order to recalculation pursuant to other boundary conditions. Properties of diffuser channel designs were compared to an experiment. Transient temperature field of the parison were calculated with Ansys APDL in order to de-termine the accumulation and stationary zone of the parison and to verify the boundary conditions for cooling channels.