Построены математические модели теплообмена при движении несжимаемой вязкой жидкости по цилиндрическим каналам с центральными винтовыми вставками в криволинейных неортогональных координатах. Для ламинарного течения в области тепловой стабилизации трехмерная задача сведена к двухмерной. Исследовано влияние геометрических и физических параметров на распределение температурного поля методом R-функций.