Silicone fluids belong to the group of pseudoplastic non-Newtonian fluids with complex rheological characteristics. They are considered in basic and applied researches and in a wide range of industrial applications due to their favorable physical and thermal properties. One of their specific field of applications in the automotive industry is the working fluid of viscous torsional vibration dampers. For numerical studies in the design and development phase of this damping product, it is essential to have thorough rheological knowledge and mathematical description about the silicone oil viscosity. In the present work, adopted rheological measurement results conducted on polydimethylsiloxane manufactured by Wacker Chemie with initial viscosity of 1000 Pas (AK 1 000 000 STAB silicone oil) are processed. As a result of the parameter identification by nonlinear regression, the temperature-dependent parameter curves of the Carreau−Yasuda non-Newtonian viscosity model are generated. By implementing these parameter sets into a Computational Fluid Dynamics (CFD) software, a temperature- and shear-rate-dependent viscosity model of silicone fluid was tested, using transient flow and thermal simulations on elementary tube geometries in the size range of a real viscous torsional vibration damper’s flow channels and filling chambers. The numerical results of the finite volume method provide information about the developed flow processes, with especial care for the resulted flow pattern, shear rate, viscosity and timing.