Thermo-mechanically coupled transport processes of viscoelastic fluids are important components in many applications in mechanical and chemical engineering. The aim of this work is the development of accurate numerical techniques for incompressible, non-isothermal, viscoelastic fluids which take into account the multi scale behavior in space and time, the multiphase character and significant geometrical changes. Based on special CFD techniques including time-space adaptivity and fast hierarchical FEM techniques, the result is a new CFD tool which is validated w.r.t. well-known benchmarks and experimental results. The advantage of such black-box tool is its flexibility in dealing with any rheological models in the coming future.