Диссертационная работа посвящена разработке численных методов математического моделирования стационарных течений вязкой теплопроводной жидкости в односвязных и многосвязных областях на основе метода R-функций. Рассмотрены линейные (приближение Стокса) и нелинейные задачи. Течения вязкой жидкости описываются нелинейным уравнением четвертого порядка относительно функции тока и системой нелинейных уравнений относительно функции тока и температуры (приближение Буссинеска).