Abstract: An efficient numerical tool for studying shock waves in viscous carbon dioxide flows is proposed. The developed theoretical model is based on the kinetic theory formalism and is free of common limitations such as constant specific heat ratio, approximate analytical expressions for thermodynamic functions and transport coefficients. The thermal conductivity, viscosity and bulk viscosity coefficients are expressed in terms of temperature, collision integrals and internal energy relaxation times. Precomputed in the wide temperature range thermodynamic functions and transport coefficients are implemented to the numerical code which is used for the simulations of the shock wave structure. Including the bulk viscosity to the kinetic model results in the increasing shock width and improves the agreement with experimental data.