Математическая энциклопедия - навье - стокса уравнения

основные уравнения движения вязкой жидкости, представляющие математическое выражение законов сохранения импульса и массы. Для неустановившегося течения сжимаемой жидкости Н.С. у. в декартовой системе координат могут быть, записаны в виде

где вектор скорости с проекциями на соответствующие оси координат давление,плотность,коэффициент вязкости;проекции вектора массовой силы на координатные оси; субстанциональная производная. При выводе уравнений (1) использован обобщенный закон трения Ньютона, предполагающий, что для движущихся жидкостей и газов напряжения пропорциональны скоростям деформаций. Для исследования сжимаемых течений к уравнениям (1) необходимо добавить уравнение состояния, связывающее между собой давление, плотность и температуру, и уравнение энергии. Уравнения (1), составляющие основу гидродинамики, впервые были получены Л. Навье [1] и С. Пуассоном [2] на основе соображений о действии межмолекулярных сил. Б. Сен-Венан [3] и Дж. Г. Стоке [4] вывели эти уравнения, допуская только, что нормальные и касательные напряжения линейно связаны со скоростями деформаций.

Для течений несжимаемой изотермич. жидкости () уравнения (1) в векторной форме могут быть представлены в виде

При анализе Н.С. у., как правило, рассматриваются в безразмерной форме, к-рая получается путем отнесения всех величин, входящих в уравнения, к соответствующим характерным величинам. Так, в случае стационарных течений несжимаемой жидкости при отсутствии массовых сил в Н.С. у. появляется один определяющий безразмерный параметр, наз. Рейнольдса числом:

где Vи lхарактерные скорости и линейный размер, v кинематическая вязкость.

Для исследования двумерных несжимаемых течений часто используются Н.С. у. в форме Гельмгольца:

где функция тока и вихрь связаны с проек циями скорости и следующим образом:

Основные краевые задачи для стационарных Н.С. у. связаны с исследованием течений в замкнутых полостях, каналах, течений со свободными поверхностями, с обтеканием тел, течений в струях и следах за телами. При этом интегрирование Н.С. у. проводится в областях (конечных или бесконечных), на границе к-рых ставятся условия из соображений физич. характера (условия прилипания или скольжения на поверхности тел, вдува или отсоса на проницаемых поверхностях, условия внешнего потока вдали от обтекаемого тела, условия на свободной границе и др.). Для нестационарных задач помимо граничных условий должны задаваться начальные условия.

Строгий математич. анализ разрешимости краевых задач гидроаэромеханики для Н.С. у. сжимаемого газа отсутствует (1982). Имеются нек-рые результаты в математич. теории динамики вязкой несжимаемой жидкости (см. Гидродинамики математические задачи).

Первоначально усилия исследователей были направлены на отискание точных решений. Напр., для несжимаемой жидкости имеются точные решения для установившихся течений: в плоском канале при заданном постоянном перепаде давления (течение Пуазёйля); между двумя параллельными плоскими стенками, одна из к-рых покоится, а другая движется в своей плоскости с постоянной скоростью (течение Куэтта); в прямолинейной трубе с круглым поперечным сечением при постоянном перепаде давления (течение ХагенаПуазейля). Найдены также нек-рые автомодельные решения, среди них: плоскопараллельное (и осесимметричное) течение вблизи критической точки (течение Хоуарта); течения в суживающемся и расширяющемся каналах (течение Гамеля).

Приближенные решения Н.С. у. основаны на упрощающих предположениях. Здесь следует отметить решения при очень малых числах Рейнольдса , соответствующие так наз. ползущим

движениям, среди к-рых наиболее известно течение Стокса около шара. Предельный случай очень больших чисел Рейнольдса приводит к теории гидродинамического пограничного слоя. Уравнения пограничного слоя позволили решить большой круг практически важных задач на основе широко разработанных приближен ных и численных методов.

Для решения нек-рых классов задач динамики вязких жидкостей и газов разработаны достаточно эффективные алгоритмы, основанные на использовании разностных схем. Напр., для задачи расчета ламинарных течений вязких несжимаемых жидкостей в областях простой формы (или около тел простой формы). Наибольшее распространение здесь получили разностные методы для уравнений в форме (3), хотя для этой системы и имеются трудности, связанные с определением граничных условий для . Первые результаты по решению стационарного варианта системы (3) были получены с помощью простейших явных пятиточечных схем и итерационных методов (см. [6]). Решение стационарных задач динамики вязкой несжимаемой жидкости большей частью основано на использовании метода установления и применении явных и неявных схем для системы (3). Среди явных схем используются двухслойные по времени схемы с симметричной аппроксимацией первых производных центральными разностями и решением второго уравнения из (3) на каждом временном слое с помощью метода Зейделя, а также трехслойная схема, в к-рой конвективные члены аппроксимируются по

схеме "крест", а диффузионные по схеме Дюфорта Франкеля. С помощью уже этих схем были получены нек-рые результаты при решении стационарных задач о ламинарных двумерных течениях в суживающемся и расширяющемся каналах, в прямоугольной выемке с движущейся крышкой, а также нестационарной задаче обтекания в канале плоской пластины, расположенной перпендикулярно к направлению потока.

Неявные схемы, как правило, основаны на применении метода дробных шагов (см. [8]). Общая структура таких схем для уравнения (3) может быть представлена, напр., в виде

где разностные одномерные операторы:

В этих формулах временной шаг,итерационный параметр,итерационные индексы при решении уравнения Пуассона из (3) на временном слое итерационным методом; разностные операторы, аппроксимирующие соответствующие вторые и первые производные. Первое уравнение из (3) используется для нахождения значений , а второе значений на последующем временном слое. Аппроксимация вторых производных, как правило, симметричная, а аппроксимация первых производных в уравнении для проводится или симметричными разностями, или односторонними разностями против потока с учетом знака скорости. Хорошо зарекомендовала себя схема в [9] с примененном монотонной аппроксимации (см. [10]). Эта схема для первого уравнения (3) имеет вид

где hи lшаги сетки по хи у,

Разностные уравнения (4) обычно приводятся к трех-диагональному виду и совместно с соотношениями, аппроксимирующими граничные условия, решаются методом прогонки. При решении стационарных задач методом установления может применяться либо поочередное решение уравнений (4) (без внутренних итераций для определения ), либо одновременное решение соответствующих уравнений из (4) для совместного нахождения и с помощью векторной прогонки. Трудности, связанные с постановкой граничных условий для уравнений (3), заключаются в том, что обычные граничные условия прилипания на твердых стенках для Н.С. у. дают условия только для . Для численного решения уравнения для формально требуют граничные условия для вихря. Эти условия могут быть получены на каждом временном слое либо приближенно на границе области, либо путем интегрирования уравнения для только в области, расположенной внутри основной области интегрирования [9].