Научная статья на тему 'Метод построения решений уравнений Навье – Стокса'

Метод построения решений уравнений Навье – Стокса Текст научной статьи по специальности «Математика»

CC BY
1382
225
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЧАСТНАЯ ПРОИЗВОДНАЯ / ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ / ДВИЖЕНИЕ / ЖИДКОСТЬ / ИНТЕГРАЛ / ВЯЗКОСТЬ / ВИХРЬ / PARTIAL DERIVATIVE / DIFFERENTIAL EQUATION / MOTION / FLUID / INTEGRAL / VISCOSITY / VORTEX

Аннотация научной статьи по математике, автор научной работы — Коптев Александр Владимирович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Method of Solution Construction for Navier – Stokes Equations

Navier — Stokes equations for viscose incompressible fluid flow are under the consideration. Analytical solution method is suggested. The method is based on the first integral of NavierStokes equations and Riccati partial differential equations. New solutions of Navier — Stokes eq uations have been constructed.

Текст научной работы на тему «Метод построения решений уравнений Навье – Стокса»

1. Bolotov A. A., Gashkov S. B., Frolov A. B. Elementamoe vvedenie v jellipticheskuju kriptografiju: Protokoly kriptografii na jellipticheskih krivyh. M.: KomKniga, 2006. S. 76-81.

2. Glusko K. L., Titov S. S. Normal'nye bazisy i derevo kvadratichnyh rasshirenij binarnyh polej // Ne-kotorye aktual'nye problemy sovremennoj matematiki i matematicheskogo obrazovanija: Materialy nauchnoj konferencii «Gercenovskie chtenija — 2012». SPb.: BAN, 2012. S. 221-226.

3. Glusko K. L., Titov S. S. Reshenie kvadratnyh uravnenij v konechnyh poljah harakteristiki dva // Problemy teoreticheskoi i prikladnoi matematiki: Trudy 43-i Vserossiiskoi molodezhnoi konferencii. Ekaterinburg: UrO RAN, 2012. S. 23-25.

4. Glusko K. L., Titov S. S. Specifika problem svjazi i upravlenija na transporte // Innovacionnyj transport. Ekaterinburg: Izd-vo UrGUPS. 2012. № 2 (3). S. 44-50.

5. Demkina O. E., Titov S. S., Torgashova A. V. Rekurrentnoe vychislenie neprivodimyh mnogochlenov v zadachah dvoichnogo kodirovanija // Molodye uchenye — transportu: Trudy IV nauchno-tehnicheskoj konferencii. Ekaterinburg: Izd-vo UrGUPS, 2003. S. 391-404.

6. Logachev O. A., Sal'nikov A. A., Jawenko V. V. Bulevy funkcii v teorii kodirovanija i kriptologii. M.: MCNMO, 2004. S. 41-53.

7. Parshin A. V. Klassicheskij protokol paketnoj kommutacii: Monografija. Ekaterinburg: Izd-vo UrGUPS, 2007. 242 s.

8. Rozhnev A. Ju. Teorija zapretov mnogobitovyh funkcij i ee primenimost' v sistemah svjazi na zhelez-nodorozhnom transporte // Sovremennye problemy nauki i obrazovanija. 2012. № 2; URL: www.science-education.ru/102-5 849.

9. Sergeev I. S. O realizacii nekotoryh operacij v konechnyh poljah shemami logarifmicheskoj glubiny: Avtoref. dis. ... kand. fiz.-mat. nauk. M., 2007. 96 s.

10. Titov S. S., Torgashova A. V. Generacija neprivodimyh mnogochlenov, svjazannyh stepennoj zavi-simost'ju kornej // Doklady Tomskogo gosudarstvennogo universiteta sistem upravlenija i radiojelektroniki. 2010. N 2 (22). Ch. 1. S. 310-318.

11. Biryukov A., Shamir A., Wagner A. Real Time Cryptanalysis of A5/1 on a PC. Fast Software Encryption Workshop. 2000. April 10-12. P. 1-18.

12. Koblitz N. Hyperelliptic cryptosystems // Journal of cryptology. 1989. N 1. P. 139-150.

13. Koblitz N. Algebraic aspects of Cryptography: Springer. 2004. 109 p.

14. Menezes A. J., Vanstone S. Elliptic Curve Cryptosystems and their implementation // Journal of Cryptology. 1993. N 6. P. 209-224.

15. Rosing M. Implementing Elliptic Curve Cryptography. Greenwich: Manning Publication, 1998.

338 p.

А. В. Коптев

МЕТОД ПОСТРОЕНИЯ РЕШЕНИЙ УРАВНЕНИЙ НАВЬЕ — СТОКСА

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

Ключевые слова: частная производная, дифференциальное уравнение, движение, жидкость, интеграл, вязкость, вихрь.

METHOD OF SOLUTION CONSTRUCTION FOR NAVIER — STOKES EQUATIONS

Navier — Stokes equations for viscose incompressible fluid flow are under the consideration. Analytical solution method is suggested. The method is based on the first integral of Navier-Stokes equations and Riccati partial differential equations. New solutions of Navier — Stokes equations have been constructed.

Keywords: partial derivative, differential equation, motion, fluid, integral, viscosity, vortex.

Введение

Уравнения Навье — Стокса — это известный тип нелинейных дифференциальных уравнений в частных производных, которые описывают движение жидких и газообразных сред при наличии вязкости. Они представляют и чисто математический интерес, и имеют многочисленные приложения к практическим задачам.

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

ди ди д( p + Ф) 1

U ~дх + V ду _ дх + Re

dv dv _ д(p + Ф) 1

U дх +V ду ду + Re

+

ду2 )

^ô2v ô2vл дх2 ду2 )

(1)

(2)

_ ». (3)

дх ду

Здесь u, v, p — основные неизвестные, скорости и давление; Ф — потенциал внешних сил, заданная функция; Re — число Рейнольдса, положительный параметр.

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

Сегодня многие вопросы изучены недостаточно и требуют дополнительного исследования. Нет общего метода решения и нет общего подхода к решению как самих уравнений Навье — Стокса, так и граничных и начальных задач для них. Нет окончательного решения проблемы существования [6, 7]. Недостаточно изучены свойства гладкости решения, нет четкого представления об асимптотике решения при Re ^ +œ .

Нерешенных математических проблем, связанных с уравнениями Навье — Стокса, множество. Их сложность оценена по достоинству. Математический институт Клэя (Clay Mathematical Institute, USA) определил изучение уравнений Навье — Стокса как одну из семи основных проблем третьего тысячелетия [8].

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

В частности, в работах [3, 4] показано, что для случая несжимаемой жидкости уравнения Навье — Стокса допускают интеграл. Этот результат имеет место в общем случае, без каких-либо дополнительных допущений о граничных и начальных условиях, характере области, значении числа Рейнольдса и других предположениях частного характера.

Дальнейшее интегрирование

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

Попытаемся проанализировать пути дальнейшего интегрирования и построения решений для этого случая.

Первый интеграл уравнений (1)-(2), соответствующих указанному случаю, представляет три соотношения. Они связывают основные неизвестные и, V, р, ассоциированное неизвестное ^ 2 и произвольные аддитивные функции одного переменного [3].

В безразмерных переменных эти соотношения можно представить в виде

р = -ф -

1

+ — 2

2

д2¥2 д 2¥

Л

дх2

+

ду 2

+ а( у ) + в(х ) =

2 2 и - V

( ди

Re

дх

дv ^ ду

д2¥2 д2¥2

ду2

Яе

ди дv

---------1-------

ду дх

Л

дх2

д2^ 2 дхду

■ + .

2 (а( У )-в( х ))

(4)

(5)

(6)

где ¥ 2 — ассоциированное неизвестное; а(у), Р(х) — произвольные функции одного переменного, у или х соответственно.

К этим соотношениям следует добавить уравнение неразрывности (3) и произвести интегрирование, рассматривая (3)-(6) как систему дифференциальных уравнений относительно и, V, р, ¥ 2.

Анализ уравнений (3)-(6) приводит к выводу, что определяющую систему фактически образуют только три уравнения, а именно уравнения (3) и (5)-(6). Эти три уравнения определяют функции и, V, ¥2. Неизвестное р фигурирует в соотношении (4) аддитивно, и оно легко находится из (4), если найдены и, V, ¥ 2.

Рассмотрим определяющую систему (3), (5)-(6). Обращает на себя внимание следующее обстоятельство. В уравнениях (5)-(6) присутствуют квадратически нелинейные чле-

2 2

ны и , V , и • V вместе с первыми производными скоростей по координатам и некоторыми другими величинами, не зависящими от и и V.

Можно констатировать, что относительно и и V уравнения (5)-(6) представляют аналоги уравнения Риккати в частных производных.

Уравнения Риккати — это хорошо изученный тип обыкновенных дифференциальных уравнений [1]:

/' + г1( х)■ / + г2 (х)-/2 + г3 (х) = 0,

(7)

где / (х) — неизвестное, г (х) — заданные коэффициенты.

Для (5)-(6) также имеем уравнения типа (7), но с частными производными, и неизвестные представляют уже функции двух переменных.

Уравнение Риккати (7) обладает многими примечательными свойствами. В частности, эти уравнения допускают решение в квадратурах при некоторых определенных значениях коэффициентов [2].

2

2

1

Логично предположить, что аналогичное свойство, пусть и в модернизированном варианте, также имеет место и для аналога уравнения Риккати в частных производных. Примем эту гипотезу за основу и попытаемся исходя из этого проинтегрировать систему (3),

(5)-(6) и построить некоторые новые решения уравнений Навье — Стокса.

Уравнения (5)-(6) удобно рассматривать совместно, если ввести комплексную скорость и = и - ¡V, как функцию комплексного переменного 2 = х + ¡у [5].

При таком подходе для неизвестных и и V должны выполняться условия Коши — Ри-

мана:

ди дv . ди & ~

дХ + ду = °- дУ-5 = °' (8)

Первое из этих условий совпадает с уравнением неразрывности (3), и это уравнение в дальнейшем будет выполнено автоматически.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Преобразуем уравнения (5)-(6) в комплексных переменных. Для удобства дальнейшего исследования умножим уравнение (6) на -2/. Это уравнение примет вид

2/

-2/ • ^ +----------

Re

ди дv

- + -

д2 ¥

= 2/д—2. (9)

ду дх ) дхду

Положим также в (4)-(5) для простоты а (у) = р° и р(х) = 42°°, где р° — аддитивная постоянная давления.

Обратим внимание на следующую закономерность. Все члены уравнений (5) и (9), рассмотренные попарно, по одному из каждого уравнения, представляют вещественную и мнимую части некоторой функции комплексного переменного.

Рассмотрим нелинейные члены уравнений (5) и (9). Они представляют соответственно вещественную и мнимую части функции и2 (г). Это справедливо на основании равенства

и2 = и2 - V2 - 2/ • т. (1°)

Рассмотрим далее члены уравнений (5), (9), содержащие первые производные и и V.

Предварительно найдем как функции комплексного переменного. Справедливы ра-

аг

венства

аи = 1 (ди дv ^ / (дм ди Л = ди / (ди дv ^ (11)

ёг 2 ду) 2 + ду) дх 2 ^5у + дх ) . ( )

Отсюда ясно, что члены с производными скоростей в уравнениях (5), (9) есть веще-

4 аи

ственная и мнимая части величины ----------—.

ке аг

Перейдем теперь к рассмотрению правых частей уравнений (5), (9).

Введем функцию

тт ( ) ■ д^2 2 (12)

и*( г) = и*- ы*=^т-А - ¡^т-2. (12)

дх ду

Пусть ¥2 (2) есть некоторая функция комплексного переменного. Тогда справедливо

д2У 2 д2У 2 дх2 ду2

и.(г). Производная и*(г) может быть преобразована как

2 2 2

равенство ^ 22 +-------2^ = °. Выполнены также условия Коши — Римана для функции

ёи* 1 ( ди* дуЛ і (ду* ди* і 1

• +

д2х¥2 д2х¥2 і .д2х¥2

- і . (13)

dz 2 ^ дх ду) 2 ^дх ду) 2 ^ дх2 дУ2 ) дхду

Из уравнения (13) следует, что правые части уравнений (5) и (9) есть соответственно

« ёи *

вещественная и мнимая части величины -2—-—.

ёг

В результате два уравнения, (5) и (9), равносильны одному комплексному уравнению. Оно имеет вид

4 аи + и2 = -2^. (14)

Яе ёг ёг

Имеем обыкновенное дифференциальное уравнение известного типа. Уравнение (14) есть специальное уравнение Риккати [1] с правой частью F(г) = -2 (~и •

Оно допускает аналитическое решение при многих известных вариантах правой части F(г) . Некоторые из них приведены в работе [2]. Создается благоприятная ситуация для построения новых решений.

Построение решений

Применим описанный выше подход и построим некоторые решения уравнений Навье — Стокса, представляющие практический интерес. Для этого достаточно разрешить уравнение (14).

Самое простое решение получается, если в правой части уравнения (14) положить нуль. Тогда приходим к уравнению

4 ёи- + и2 = °. (15)

Re ёг

Его решением является

41

и = Я- •—!-, (16)

Яе ¿1 - г

где г1 = А + ¡В — некоторая фиксированная точка комплексной плоскости. В этой точке имеем простой полюс.

Данное решение допускает простую гидромеханическую трактовку. Точку г1 можно рассматривать как точку расположения стока жидкости с обильностью

2

т1=пке. (17)

Данное решение может представлять интерес для гидрологии и для теории нефтедобычи.

Найдем скорости, соответствующие этому решению. Достаточно отделить вещественную и мнимую части.

В результате получаем выражения

и = А_______А- х______ у = _4_______В1 - У_____ (18)

“ке‘ (А - х)2 +(В1 - у)2, _ке' (А - X)2 +(В1 - У)2' ( }

С учетом уравнения (4) находим давление

р - Ро = -ф-•( а - X)2 + (В, - У )2 • (19)

где ро — аддитивная постоянная давления, Ф — потенциал внешних сил.

С практической точки зрения может быть интересно решение в виде суммы нескольких простых полюсов. В простейшем случае — двух. Пусть, например, в каждой из точек г1 = А1 + ¡В1 и г2 = А2 + В имеем простой полюс. Следуя гидромеханической трактовке, можно считать, что в этих точках расположены стоки с одинаковой обильностью (17). Решение уравнения (14) ищем в виде

^ке-Ьг-г+ггг). (20)

Подстановка уравнения (20) в (14) приводит к соотношению

16 1 ( 1 1 і ёи*

-------- —

Яе2 г2 - г1 ^ г1 - г г2 - г) ёг

Это уравнение допускает интегрирование. В результате получаем

(21)

и* = -^--------- • 1п . (22)

Яе (г2 - г1) г2 - г

На основании уравнения (12) выражение (22) определяет частные производные функции ¥ 2. Возможно и последующее интегрирование (22). В результате может быть найдена функция ¥ 2. Однако эта функция не представляет особого интереса. В выражении (4) для давления фигурирует лишь лапласиан ¥2, а последнее заведомо равно нулю.

Таким образом, для рассматриваемого случая определяющая система разрешена. Отделяем в уравнении (20) вещественную и мнимую части и находим решение уравнений Навье — Стокса, соответствующее двум изолированным стокам:

4 ( А\ — х А2 — х

и = —• --------^----------;т + - 2

Яе ^(4 - х)2+(В - у)2 (А - х)2+(В - у)2)'

у=ЯЧ------------------2 +-------Вг-у-----------------------зі (23)

Яе I (А - х)2+(В - у)2 (А - х)2+(В - у)2)

Используя выражение (4), находим давление:

- = _________• (__________1_________ __________1_________

Р Ро Ке2 ((А - х)2+(В - у)2)2 + ((А2 - х)2 + (В2 - у)2)2 +

+2 (А1 - х) •(А2 - х) + (В1 - У) • (В2 - у) ) ( )

((А1 - х) + (В1 - у) ) •((А2 - х)2 + (В2 - у)2 )

Представляется возможным путем некоторой вариации коэффициентов получить решение и для другого важного случая.

Это случай двух изолированных и различно ориентированных вихрей с одинаковыми по модулю интенсивностями. Такая схема может быть интересна для метеорологии, так как дает простейшую модель движения двух различно ориентированных атмосферных вихрей — циклона и антициклона.

Для этого случая комплексная скорость получается из уравнения (20) умножением правой части на i. Для этой величины имеем выражение:

-г+ф;) • (25)

Отделяя вещественную и мнимую части, находим скорости

и = 4 /---------В1-У----, +--------В2-У

Яе ^ (А1 - х)2 + (В - У)2 (А2 - х)2 + (В - У)2 )’

= 4 Г А1 - х А2 - х ^ (26)

I (А - х)2 + (В - У)2 (А2 - х)2 + (В2 - У)2)' ( ’

Выражение для давления совпадает с уравнением (24).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

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

Выводы

Предложен новый метод построения решений 2D уравнений Навье — Стокса. Метод основан на первом интеграле уравнений Навье — Стокса и на свойствах (возможностях) уравнения Риккати в частных производных.

Каждый известный случай интегрируемого в квадратурах уравнения Риккати генерирует новые решения уравнений Навье — Стокса.

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

СПИСОК ЛИТЕРАТУРЫ

1. Егоров А. И. Обыкновенные дифференциальные уравнения. М.: Изд-во физ.-мат. лит-ры, 2007. 448 с.

2. Зайцев В. Ф., Полянин А. Д. Обыкновенные дифференциальные уравнения: Справочник. М.: Физматлит, 2001. 576 с.

3. Коптев А. В. Интегралы уравнений Навье — Стокса. Саранск: Труды Средне-волжского математического общества. 2004. № 1. Т. 6. С. 215-225.

4. Коптев А. В. Общий интеграл уравнений Навье — Стокса // Математическая физика и ее приложения: Материалы второй международной конференции. Самара: МИАН им. В. А. Стеклова, СамГУ, 2010. С.179-180.

5. Лаврентьев М. А., Шабат Б. В. Методы теории функции комплексного переменного. М.: Наука, 1987. 688 с.

6. Ладыженская О. А. Математические вопросы динамики вязкой несжимаемой жидкости. М.: Наука, 1970. 288 с.

7. Темам Р. Уравнения Навье — Стокса. Теория и численный анализ. М.: Мир, 1981. 408 с.

8. Charles L. Fefferman. Existence and Smoothness of the Navier — Stokes equation. Preprint, Princeton Univ., Math. Dept. Princeton, NJ, 2000. Р. 1-5.

REFERENCES

1. Egorov A. I. Obyknovennye differencial'nye uravnenija. M.: Izd-vo fiz.-mat. lit-ry, 2007. 448 s.

2. Zajcev V. F., Poljanin A. D. Obyknovennye differencial'nye uravnenija. Spravochnik. M.: Fizmatlit, 2001. 576 s.

3. Koptev A. V. Integraly uravnenij Nav'e — Stoksa. Saransk: Trudy Sredne-volzhskogo matemati-cheskogo obwestva. 2004. № 1. T. 6. S. 215-225.

4. Koptev A. V. Obwij integral uravnenij Nav'e — Stoksa // Matematicheskaja fizika i ee prilozhenija: Materialy vtoroj mezhdunarodnoj konferencii. Samara: MIAN im. V. A. Steklova, SamGU, 2010. S. 179-180.

5. Lavrent'evM. A., Shabat B. V. Metody teorii funkcii kompleksnogo peremennogo. M.: Nauka, 1987.

688 s.

6. Ladyzhenskaja O. A. Matematicheskie voprosy dinamiki vjazkoj neszhimaemoj zhidkosti. - M.: Nauka, 1970. 288 s.

7. Temam R. Uravnenija Nav'e — Stoksa. Teorija i chislennyj analiz. M.: Mir, 1981. 408 s.

8. Charles L. Fefferman. Existence and Smoothness of the Navier — Stokes equation. Preprint, Princeton Univ., Math. Dept. Princeton, NJ, 2000. Р. 1-5.

Ю. Н. Ловягин

О ПРОБЛЕМЕ НОРМИРУЕМОСТИ БУЛЕВЫХ АЛГЕБР

Рассматриваются вопросы, связанные с существованием меры на полной булевой алгебре. Доказано, что на регулярные булевы алгебры допускают полумеру. Приведен пример, показывающий, что не всякая регулярная булева алгебра является нормируемой.

Ключевые слова: булева алгебра, регулярная булева алгебра, мера, нормируемость, полунормируемость.

Yu. Lovyagin

ON THE PROBLEM OF MEASURABLIZABILITY OF BOOLEAN ALGEBRAS

This paper discusses the issues related to the existence of measures on a complete Boolean algebra. It is argued that for regular Boolean algebras half-measures are possible. An example is given showing that not every regular Boolean algebra is measurable.

Keywords: Boolean algebra, regular Boolean algebra, the measure, measurable Boolean algebra, seminormativity.

i Надоели баннеры? Вы всегда можете отключить рекламу.