Перейти на главную страницу
ОГЛАВЛЕНИЕ
2. Преобразование нестационарных уравнений Навье-Стокса методом последовательных приближений…………………………………………………………………………..…….…………2
3. Общий подход к решению уравнений Навье-Стокса в первом приближении………….…...….5
4. Решение линейного уравнения движения для плоских волн……………………………..……....6
5. Решение линейного стационарного уравнения движения в сферической геометрии…………..8
6. Решение линейного стационарного уравнения движения в цилиндрической геометрии и в других криволинейных геометриях……………………………………………………………….....10
7. Решение линейного нестационарного уравнения движения в трехмерной геометрии…...…..13
8. Решение трехмерного уравнения неразрывности в первом приближении …………………....15
9. Множество решений уравнений Навье-Стокса…………………………………………………..24
10. Общий подход к решению задачи Коши для линейного приближения уравнений Навье-Стокса. Решение с учетом силы гравитации………………………………………………………..28
11. Решение уравнений Навье-Стокса во втором и более высоких приближениях……………...31
12. Подход к решению уравнения неразрывности во втором и более высоких приближениях…………………………………………………………………………………………42
13. Заключение…………….……………………………...…………..………………………………44
Литература…………………………………………………………………………………………….44
Решение получено для сжимаемой и несжимаемой жидкости. Решение нелинейных уравнений ищется методом последовательных приближений. Метод последовательных приближений позволяет свести задачу к решению последовательности систем линейных дифференциальных уравнений для неизвестной функции с учетом нелинейной части, в которую входят известные функции. Затем используется преобразование Лапласа для временной переменной и преобразование Фурье – для пространственных переменных. Таким образом, решение представляется в виде наложения плоских волн с различными частотами и значениями волнового вектора. Движение этих волн рассматривается в трехмерном пространстве. В основу решений задачи положены их представления для плоской волны. Уравнения, полученные в результате использования решений в виде плоских волн, сводятся к уравнениям для амплитуд этих волн.
Показано, что решения уравнения движения и дифференциального уравнения неразрывности для несжимаемой жидкости, построенные в виде плоской волны, становятся идентичными только при учете в дифференциальном уравнении неразрывности вязкости за счет процессов неупругих столкновений. Только при этом условии можно получить решения, описывающие движение турбулентного потока несжимаемой жидкости.
Решение уравнений Навье-Стокса в трехмерных криволинейных системах координат выполняется на основе решения для плоской волны, полученного в трехмерной геометрии. Это осуществляется путем преобразования решения для плоской волны к решению в виде наложения волн другой геометрической формы. В работе подробно рассмотрено получение решения уравнений Навье-Стокса в сферической геометрии, путем перехода от представления решения в виде наложения плоских волн к представлению этого решения в виде наложения сферических волн.
Помимо представления решений в виде плоской волны существуют также и другие представления решений задачи: в виде разложений в ряд по известным аналитическим функциям. Эти функции являются решениями уравнения Гельмгольца и Лапласа в различных ортогональных трехмерных криволинейных системах координат. Путем замены переменных осуществляется переход от представления решений в одних криволинейных координатах к представлениям в других криволинейных координатах. Таким образом, можно получить решения во всех известных криволинейных системах координат. В большинстве случаев переменные легко разделяются. За основу рассмотрения взято представление решений в сферической системе координат.
Подробно рассмотрены решения уравнения движения и уравнения неразрывности с учетом особенностей этих уравнений для различных представлений решений. Рассмотрено решение однородного (без источников жидкости) и неоднородного (со сторонними источниками жидкости) уравнения неразрывности. Раскрыты особенности применения дифференциального и интегрального уравнения неразрывности для несжимаемой и сжимаемой жидкости при построении решений системы. Получено решение с учетом силы гравитации. Полученные аналитические решения позволяют решить классические, корректно поставленные задачи Коши.
Система нестационарных уравнений Навье-Стокса для сжимаемой жидкости в трехмерной геометрии выглядит следующим образом:
где – плотность жидкости в момент времени в точке с координатами , – векторное поле скоростей жидкости, – давление жидкости, – коэффициент динамической вязкости, – коэффициент второй вязкости, – векторное поле массовых сил, – оператор Гамильтона, – оператор Лапласа, .
В дальнейшем примем, что
, (2.3)
где – ускорение свободного падения.
Уравнение (2.1) является уравнением движения, а уравнение (2.2) – уравнением неразрывности.
Краевые и начальные условия выглядят следующим образом:
где – некоторая замкнутая поверхность, – некоторая известная векторная функция.
Уравнения (2.1) и (2.2) будем решать методом последовательных приближений:
, (2.6)
, (2.7)
, (2.8)
где , , , – слагаемые одного порядка малости.
При этом будем полагать, что
, (2.9)
, (2.10)
, (2.11)
где – постоянные величины, независящие от времени и координат . Будем полагать, что эти постоянные величины заданы. Однако в дальнейшем мы допустим также зависимость
Для упрощения записей уравнений введем обозначения
Сначала будем полагать, что жидкость сжимаема. Подставим выражения (2.6)-(2.8) в уравнения (2.1)-(2.2) и (2.4)-(2.5). Учтем обозначения (2.9)-(2.15). Уравнения (2.1) и (2.2) запишем в виде систем уравнений, в которые входят слагаемые одного порядка малости. На первом этапе построения решения не будем использовать функцию . В результате получим первую систему уравнений
Краевые и начальные условия для этой системы выглядят следующим образом:
Вторая система уравнений имеет вид
Краевые и начальные условия для второй системы выглядят следующим образом:
Необходимо отметить, что все системы уравнений следует рассматривать как единую систему уравнений. Поэтому, если мы решаем задачу для скорости в виде приближения
то, нам следует использовать только условия (2.22)-(2.23), не учитывая условия (2.18)-(2.19).
Третья система уравнений записывается как
(2.25)
. (2.26)
Краевые и начальные условия для третьей системы выглядят аналогично предыдущим.
Таким же образом получаются и все остальные системы уравнений. При заданном порядке приближения скорости их число будет конечным.
Для получения уравнений для несжимаемой жидкости во всех системах следует положить
В результате получим первую систему уравнений в виде
Отметим, что в уравнении (2.28) мы не опустили слагаемое .
Вторая система уравнений имеет вид
(2.30)
. (2.31)
Третья система уравнений записывается как
Краевые и начальные условия для всех систем не изменились.
Отметим, что вместо уравнений (2.29), (2.31), (2.33) и т.д., мы можем записать одно уравнение
. (2.34)
Получение уравнений для стационарной задачи тривиально. Особенности этой задачи мы рассмотрим ниже.
Для нестационарной задачи справедлива зависимость (2.12). Для несжимаемой жидкости мы можем записать
, (2.35)
где , – потенциальная энергия поля массовых сил.
массовые силы отсутствуют. Зависимость скорости от времени будет использована при решении задачи о движении сжимаемой и несжимаемой жидкости в поле сил гравитации.
Рассмотрим систему (2.16)-(2.17) с условиями (2.18)-(2.19). Запишем эту систему еще раз.
Эта система справедлива как для сжимаемой, так и для несжимаемой жидкости. Уравнение (3.2) удовлетворяется тождественно. Если жидкость несжимаема, то мы можем использовать вместо него уравнение (2.31)
Но решение уравнения (2.31) требует особого подхода, который будет представлен ниже.
Будем искать решения системы линейных уравнений (3.1)-(3.5) в виде наложений плоских волн с различными частотами и значениями волнового вектора. Отметим специфику. Данные плоские волны задаются в трехмерном пространстве. Временная зависимость включает не только простые периодические функции времени, но и экспоненту с вещественным показателем. Поэтому искомая временная зависимость может быть построена с помощью обратного преобразования Лапласа. Итак,
, (3.6)
, (3.7)
где – мнимая единица, – комплексная переменная, – некоторое вещественное число, – волновой вектор, , , , . Отметим, что вместо функций и можно использовать функции и , где . При этом метод решения уравнений Навье-Стокса не измениться. Но целесообразность использования таких выражений в данной работе рассматриваться не будет.
Положим, что функции , , представляются в виде
где , – вещественные числа, , являются скалярными вещественными функциями векторного аргумента, а , – векторными вещественными функциями векторного аргумента.
Векторный аргумент этих функций выражается набором вещественных переменных . Но решение уравнений Навье-Стокса может быть получено и для комплексной переменной . Более подробно этот вопрос будет рассмотрен ниже. Для упрощения получаемых выражений сначала рассмотрим случай вещественной переменной .
Рассмотрим уравнение движения (3.1). Подставим в него выражения (3.5)-(3.7). В силу того, что переменная интегрирования не зависит от координат, дифференциальные операторы можно внести под знак интеграла. Выполняя элементарные операции дифференцирования экспоненциальной функции, получим векторное линейное алгебраическое уравнение, которое является уравнением для амплитуд плоских волн:
В этом уравнении мы должны сделать подстановки, используя формулы (3.10)-(3.12). В этом случае мы получим систему двух уравнений для функций , , , . Эти функции дополнят друг друга при построении конечного решения, зависящего от времени и координат. Для упрощения изложения решим уравнение (4.1) для комплексных функций и , имея в виду, что все сказанное о них относится к функциям и , а также и .
Определим векторную функцию в виде
, (4.2)
где – любой вектор, перпендикулярный вектору , абсолютное значение которого равно :
Умножим векторное уравнение (4.1) сначала на вектор , а затем на вектор . Получим два скалярных уравнения
Уравнение (4.7) выражает зависимость давления от скорости:
Уравнение (4.8) имеет два решения:
Решение (4.10) означает, что
и величина не является функцией вектора . Тогда
где представляется формулой (4.9).
Решение (4.11) означает, что величина является функцией вектора : . Тогда выражения (3.6) и (3.7) примут вид
, (4.15)
, (4.16)
где . (4.17)
В этих выражениях мы учли, что и являются функциями только вектора .
Итак, мы получили два вида решения нестационарного уравнения движения: (4.9), (4.13), (4.14) и (4.15)-(4.17). Эти два вида решения могут использоваться как отдельно, так и совместно.
Рассмотрим также стационарное уравнение движения:
, (4.18)
Данное уравнение сведется к системе уравнений
Очевидно, что уравнение (4.20) имеет только одно решение:
Тогда решение стационарного уравнения движения имеет вид
При этом зависимость давления от скорости имеет вид
Решение стационарного уравнения соответствует решению нестационарного уравнения вида (4.9), (4.13), (4.14). Решение (4.15)-(4.17) не имеет соответствующего вида решения стационарного уравнения. Однако для случая мнимых значений соответствующее решение стационарного уравнения существует.
Все полученные решения линейного уравнения движения, а именно (4.9), (4.13), (4.14), (4.15)-(4.17), (4.22)-(4.24), строятся в виде наложения плоских волн с различными частотами и значениями волнового вектора. Как известно, такое наложение плоских волн, распространяющихся в пространстве трех измерений, определяет решение волновых уравнений в произвольной трехмерной геометрии. В результате определенных преобразований можно перейти к представлению решения в виде наложения волн другой геометрической формы. Решение уравнения движения в трехмерных системах координат выполняется на основе решения для плоской волны, полученного в трехмерной геометрии. Более подробно вопрос о методе получения решений в различных криволинейных координатах, исходя из решения уравнения для плоской волны, рассмотрен в работах [8,9]. Используем этот метод и в данной работе. Применим основные идеи этого метода к решению уравнений Навье-Стокса.
Сначала получим решения уравнений Навье-Стокса в сферической геометрии, исходя из решений для плоской волны в трехмерной геометрии. Получение решений в сферической геометрии на основе решения для плоской волны очень хорошо разработано в теории специальных функций [5] и используется в различных разделах теоретической физики, в частности, в квантовой механике и теории поля. Решение в других системах координат трехмерного пространства можно получить на основе решений в сферической геометрии путем соответствующей замены координат.
Для определенности мы выбрали знак плюс перед произведением .
Рассмотрим выражение (5.1). Используем известное из теории специальных функций разложение [5]:
, (5.4)
где – сферические функции Бесселя первого рода; .
Воспользуемся теоремой сложения для полиномов Лежандра:
. (5.5)
При этом сферические функции имеют вид
где – произвольный вектор.
Нормировка сферических функций имеет вид
, (5.8)
Вектор представим в виде
где , , – орты прямоугольной декартовой системы координат. Эта систему будем считать главной при решении задачи.
Учтем, что для любого вектора выполнено
, (5.10)
(5.11) . (5.12)
Тогда вектор представляется в виде
В результате получим выражение
Воспользуемся формулами
Полученное выражение приобретает вид
Получим выражение
Введем обозначение
Получим выражение
Получим выражение для давления. Направим ось по направлению вектора :
Согласно формуле (4.24) давление равно
Очевидно, что
Коэффициенты связаны рекуррентными отношениями. Например, имеет место формула
Поэтому
Выражения, полученные для скорости и давления, являются функциями трех пространственных сферических координат. Путем перехода к новым трем пространственным переменным они могут быть представлены в координатах любых известных криволинейных систем, для которых существуют соответствующие формулы перехода. Причем в большинстве случаев все переменные разделяются. Более подробно этот вопрос рассмотрен в работе [9]. В качестве примера получим выражение для вектора скорости в цилиндрических координатах.
Связь между декартовыми и цилиндрическими координатами выражается формулами:
где , , .
Связь между цилиндрическими и сферическими координатами выглядит следующим образом:
, , . (6.2)
Рассмотрим выражение для скорости в упрощенном виде
Мы опустили интегрирование по модулю волнового вектора, поскольку оно не существенно для последующих преобразований. Для удобства обозначения мы положили
Введем в правую часть равенства множитель
Выделим в правой части данного выражения однородный гармонический полином степени
где .
где .
Замечание. Для обоснования данной формулы удобно использовать следующее определение сферических функций:
(6.7)
Используем теорему сложения Гегенбауэра
где - гамма-функция; - полином Гегенбауэра; . ; .
Замечание. В случае разделение переменных не требуется, поскольку величина принимает фиксированное значение.
. (6.9)
, (6.10)
где – полином Якоби.
Первые полиномы Якоби в явном виде выглядят следующим образом:
(6.11)
Введем обозначение:
Также выполнено соотношение:
Учтем также соотношение:
Тогда
Формула справедлива как для , так и для .
Подставим в формулу (6.5) выражения (6.6) и (6.15). Получим
Учтем, что . Получим
Воспользуемся формулами:
где и – постоянные коэффициенты.
Замечание. Данные формулы являются следствиями соотношения
.
В результате получим выражение:
Введем обозначение:
Тогда
Замечание. Функции и можно представить с помощью функций и , где .
Решения в других ортогональных системах криволинейных координат получаются аналогично. Мы можем построить решения в различных ортогональных системах криволинейных координат на основе сферических функций Бесселя. Но возможен и другой подход.
Отметим, что сферические функции Бесселя удовлетворяют частному случаю соответствующего уравнения Ломмеля (см. [9]). Уравнение Ломмеля можно получить из уравнения Бесселя путем замены переменных. Таким образом, решения уравнений Навье-Стокса строятся на основе решений дифференциального уравнения Бесселя, а также уравнения Лапласа в сферической геометрии. Напомним, что частное решение уравнения Лапласа в этой геометрии строится на основе сферических функций. В свою очередь, к решению этих уравнений сводится решение уравнения Гельмгольца в сферической геометрии. Уравнения Лапласа и Гельмгольца в различных трехмерных системах ортогональных криволинейных координат преобразуются к обобщенным уравнениям гипергеометрического типа. Решениями этих уравнений являются специальные функции, характерные для данной криволинейной геометрии. На этих функциях и должно строится решение трехмерных уравнений Навье-Стокса в различных системах криволинейных координат. Решение, представленное в этих функциях, может иметь более простой вид.
Преобразование нестационарных уравнений Навье-Стокса методом последовательных приближений
11 10 2014
4 стр.
Уравнения 1 есть дифференциальные уравнения Навье1 – Стокса2 движения вязкой жидкости, являющиеся математическим описанием полей скоростей и давлений в подвижной среде
17 12 2014
1 стр.
Цель урока: Познакомить учащихся с методом решения уравнений, основанном на применении теоремы Безу. Научить использовать его при решении уравнений
13 10 2014
1 стр.
Единая модель udm удовлетворяет любым аналитическим нуждам, от многомерного анализа и создания отчетов до прогнозирующей аналитики, — все в рамках единственного решения
10 10 2014
1 стр.
Законы или правила Кирхгофа. Делители напряжений и токов. Возможные методы упрощения систем уравнений (метод узловых потенциалов и эквивалентного источника). Машинный метод решения
14 12 2014
1 стр.
Родина геометрии Древний Египет. Первые шаги геометрии в Китае, Индии, Вавилоне, Ассирии были связаны с необходимостью измерять. Основные черты геометрии: в Вавилоне носили арифмет
14 10 2014
1 стр.
В задачах 36-40 исследовать данную систему уравнений на совместимость и решать ее, если она совместна
18 12 2014
1 стр.
Рассмотрена проблема решения плохообусловленных систем линейных алгебраический уравнений (слау). Описан разработанный эволюционного метод решения слау основанный на принципе Парето
24 09 2014
1 стр.