Граничные условия и численный метод решения

Система определяющих уравнений (1)-(4) и замыкающие соотношения (5)-(15) записываются в цилиндрической системе координат.

Краевые условия сводятся: к заданию на входе однородных профилей для осевой скорости осредненного течения, нормальных напряжений (), кинетической энергии и скорости диссипации (,где Tu-интенсивность турбулентности), другие параметры в этом случае равны нулю.

При x= (выход) принято На стенке (r=R) реализуются условия прилипания потока и отсутствия турбулентных пульсаций, а также На оси симметрии (r=0):

, V=W=0,

Численное интегрирование выполняется на неравномерных сетках. Сгущение узлов к твердым поверхностям отвечает замене переменных в исходных уравнениях:

где D-параметр преобразования, обеспечивающий попадание трех- пяти узлов в область . Детали интегрирования преобразованных таким образом уравнений, выполненного с применением экономичных неявных конечно- разностных схем можно найти в [10].