Лекция 4. Идеальная несжимаемая жидкость.

advertisement
Лекция 4. Идеальная несжимаемая жидкость.
Жидкость называется идеальной, если коэффициенты вязкости равны нулю. Предположим,
что ρ(t, x) является константой. Тогда уравнения, описывающие движение идеальной
несжимаемой жидкости, имеют вид
µ ¶
f
p
+ , div v = 0.
vt + hv, ∇iv = −∇
(0.1)
ρ
ρ
Будем считать, что жидкость все время находится в области Ω. В частности, для
любых t, t̃ ∈ R и для любой точки x ∈ Ω найдется точка x̃ ∈ Ω такая, что если
частица жидкости в момент времени t расположена в точке x, то в момент t̃ она
перейдет в точку x̃.
Имеет место теорема Лагранжа.
Теорема 1. Предположим, что сила f потенциальна, т.е. f = ∇U . Тогда если
rot v|t=0 = 0, то rot v = 0.
Доказательство. Движение жидкости задается дифференциальным уравнением
ẋ = v(t, x).
(0.2)
Тогда каждому t > 0 можно сопоставить диффеоморфизм g t : Ω → Ω следующим
образом: если y ∈ R3 , то g t (y) — решение уравнения с начальным условием x(0) = y.
При этом g t (y) гладко по совокупности (t, y).
Пусть γ : [0, 1] → R3 — параметризация гладкой замкнутой кривой, γt = g t ◦ γ.
Покажем, что
Z
v(t, x(t)) dl = const.
(0.3)
γt
В самом деле,
d
dt
Z
d
v(t, x(t)) dl =
dt
Z1
γt
v(t, g t ◦ γ(s))
d t
g ◦ γ(s) ds =
ds
0
Z1 µ
=
¶
∂v
d
(t, γt (s)) + hv, ∇iv(t, γt (s))
γt (s) ds+
∂t
ds
0
Z1
+
v(t, γt (s))
d
v(t, γt (s)) ds =
ds
0
1
=
ρ
Z1
0
d
1
(−∇p + ∇U ) γt (s) ds +
ds
2
Z1
1
dv (t, γt (s)) =
ρ
Z
2
0
(−∇p + ∇U ) dl = 0.
γt
1
Если γ является границей поверхности S, то по формуле Стокса из (0.3) получаем
Z
rot v(t, x(t)) dS = const.
gt S
В частности, если rot v(0, x) = 0 для любого x, то для любого t выполнено rot v(t, x(t)) =
0. В самом деле, если rot v(t, x) 6= 0 в некоторой точке x, то возьмем достаточно
малый участок плоскости S, проходящий
через точку x и ортогональный вектору
R
rot v(t, x), и для него получим, что rot v dS 6= 0 — противоречие.
S
Если область, в которой находится жидкость, односвязна, то условие rot v = 0
эквивалентно существованию скалярной функции ϕ такой, что v = ∇ϕ. Из условия
несжимаемости следует, что
∆ϕ = 0.
(0.4)
Если граница области Γ задана и неподвижна, то для функции v, как правило,
пишется условие непроницаемости hv, ni = 0, где n — нормаль к Γ. Значит, ∂ϕ
= 0.
∂n
Учитывая (0.4), получаем задачу Неймана. Если область ограничена, то нетривиальное
(не являющееся константой) решение возможно только при наличии особенностей.
Если область бесконечна, то для корректной постановки задачи еще нужно задать
условие на бесконечности.
Поскольку rot v = 0, то
3
3
∂ |v|2 X ∂v j j X j ∂v i
=
v =
v
= hv, ∇iv i .
i
j
∂xi 2
∂x
∂x
j=1
j=1
Значит, из
следует, что
p
U
vt + hv, ∇iv = −∇ + ∇
ρ
ρ
µ
¶
|v|2 p U
∇ ϕt +
+ −
= 0,
2
ρ
ρ
т.е.
ϕt +
|v|2 p U
+ − = Φ(t).
2
ρ
ρ
Последнее равенство называется интегралом Коши–Лагранжа.
Если движение стационарное, т.е. все функции не зависят от t, то
|v|2 p U
+ − = const.
2
ρ
ρ
Это равенство называется интегралом Бернулли.
0.1
Плоскопараллельное течение.
Движение называется плоскопараллельным, если существует такая декартова система
координат (x, y, z), в которой вектор скорости имеет вид (u(x, y), v(x, y), 0). Условия
несжимаемости и потенциальности записываются соответственно в виде
∂u ∂v
+
= 0,
∂x ∂y
∂u ∂v
−
= 0.
∂y ∂x
2
Тогда существуют такие функции ϕ и ψ, что ∂ϕ
= u, ∂ϕ
= v, ∂ψ
= −v, ∂ψ
= u.
∂x
∂y
∂x
∂y
Отсюда видно, что функция f (x + iy) := ϕ(x, y) + iψ(x, y) удовлетворяет условиям
Коши – Римана, поэтому является аналитической. Заметим также, что линии уровня
функций ϕ и ψ ортогональны, при этом линии уровня функции ψ совпадают с
траекториями движения (поэтому ψ называется функцией тока).
Рассмотрим задачу обтекания. Пусть в R2 задана односвязная ограниченная
область с гладкой границей. Требуется в R2 \Ω найти такую гармоническую функцию
= 0, а на бесконечности —
ϕ, что на ∂Ω выполнено условие непроницаемости ∂ϕ
∂n
∇ϕ(x, y) 2 →
(U∞ , V∞ ).
2
x +y →∞
По теореме Римана и теореме о соответствии границ, найдется конформное отображение
w области R2 \Ω на дополнение единичного круга, непрерывное вплоть до границы,
при этом ∂Ω переходит в единичную окружность. Итак, наша задача распалась на
две:
1. найти это конформное отображение,
2. решить задачу обтекания для круга.
Мы рассмотрим задачу обтекания круга радиуса R. Подберем такую аналитическую
функцию f (z) = ϕ(z) + iψ(z), чтобы ψ||z|=R = const. Тогда окружность переходит
в горизонтальный отрезок. Вспоминаем,
что для R = 1 таким свойством обладает
¡
¢
1
1
функция Жуковского f (z) = 2 z + z . Для круга радиуса R и условия
´
³ на бесконечности
2
lim ∇ϕ(x, y) = (U∞ , 0) искомой является функция f (z) = U∞ z + Rz . Перейдя
снова к вещественным переменным, получаем
¶
¶
µ
µ
R2 x
R2 y
ϕ(x, y) = U∞ x + 2
, ψ(x, y) = U∞ y − 2
.
x + y2
x + y2
x2 +y 2 →∞
Задача о построении конформного отображения круга на односвязную область в
общем случае является сложной. Различные приближенные методы описаны в [1, 2,
3]. Кратко опишем идею некоторых из них.
1. Приблизим область многоугольником. Если ее граница достаточно хорошая, то
соответствующее конформное отображение близко к искомому в равномерной
метрике. Единичный круг на многоугольник отображается с помощью интеграла
Кристоффеля – Шварца
Zz
(ζ − a1 )α1 −1 . . . (ζ − an )αn −1 dζ + B,
w(z) = A
1
где παj — значение угла при j-й вершине, A, B, aj — неопределенные постоянные.
Трудность состоит в нахождении aj . Один из способов их определения состоит
в использовании так называемых обобщенных степенных рядов [2].
2. М.А. Лаврентьевым были найдены приближенные формулы для конформных
отображений областей, близких к тем областям, отображения которых на канонические
известны. При этом использовался вариационный принцип Линделефа и его
обобщения.
3. В [2] описан метод последовательных конформных отображений области, содержащейся
в верхней полуплоскости, на верхнюю полуплоскость. Для широкого класса
3
областей он позволяет получить сколь угодно высокую степень точности приближения.
Сначала строится отображение области
¾
½
x2 y 2
G = x + iy : y > 0, 2 + 2 > 1
a
b
2
2
на верхнюю полуплоскость так, чтобы полуэллипс xa2 + yb2 = 1, y > 0, перешел в
отрезок
|x| 6 a+b, y = 0. Для этого сначала полуэллипс с помощью отображения
√
z+ z 2 +b2 −a2
τ=
переводится в единичную полуокружность, которая с помощью
a+b
функции Жуковского и гомотетии переводится в отрезок. Затем для области
{x+iy : y > ϕ(x)}, где ϕ — положительная гладкая функция, строится последовательность конформных отображений. На каждом шаге берутся три точки на
границе области, через них проводится полуэллипс, и затем строится отображение,
переводящее его в горизонтальный отрезок. Обратное отображение полуплоскости
на область получается как композиция обратных отображений, которые явно
вычисляются по уже найденным точкам.
В [3] рассмотрены методы отображения более сложных областей.
0.2
Осесимметрическое течение.
Пусть r — трехмерный вектор, v — вектор скорости. Течение называется осесимметрическим, если найдется цилиндрическая система координат (x, r, α) со следующим
свойством: все векторы скорости v(r) лежат в полуплоскости, проходящей через ось
x и точку r, и не зависят от координаты α. Если движение потенциально и жидкость
несжимаема, то соответствующая гармоническая функция ϕ зависит только от x и
r.
Рассмотрим одну из таких полуплоскостей и функцию ϕ = ϕ(x, r) в этой полуплоскости. Уравнение Лапласа в цилиндрических координатах имеет вид
µ µ
¶
µ
¶¶
∂ 2 ϕ ∂ 2 ϕ 1 ∂ϕ
1 ∂
∂ϕ
∂
∂ϕ
+ 2 +
≡
r
+
r
= 0.
∂x2
∂r
r ∂r
r ∂x
∂x
∂r
∂r
Это равенство
можно
рассматривать как условие потенциальности для векторного
¢
¡
∂ϕ
поля −r ∂ϕ
,
r
.
Введем
сопряженную функцию ψ, удовлетворяющую равенствам
∂r
∂x
r
∂ϕ
∂ψ
=
,
∂x
∂r
r
∂ϕ
∂ψ
=− .
∂r
∂x
(0.5)
Функция f = ϕ+iψ не является аналитической. Все функции такого вида, являющиеся
решением системы (0.5), задают класс квазиконформных отображений (определение
и свойства таких отображений приведены, например, в [1]; в частности, имеет место
аналог теоремы Римана о конформном отображении).
Рассмотрим некоторые примеры решений уравнения Лапласа. Решения, зависящие
только от r или только от x, имеют соответственно вид ϕ = a ln r + b и ϕ = ax + b.
Подбирая функции ψ, удовлетворяющие (0.5), и выбирая некоторую нормировку,
получаем функции
Z 1 = ln r − ix, Z 2 = 2x + ir2 .
Аналогично можно найти гармонические полиномы ϕ произвольной степени1 и соответствующие функции Z n .
1
Например, через шаровые функции, которые, в свою очередь, находятся методом разделения
переменных — см. [4].
4
Теперь рассмотрим √
решение трехмерного оператора Лапласа с особенностью в
нуле: ϕ = R1 , где R = r2 + x2 . Найдем сопряженную функцию ψ. Учитывая вид
r)
уравнений (0.5) и функции ϕ, удобно искать решение в виде ψ(x, r) = g(x,
. Тогда
R
получается система
∂g
∂g
R2
− rg = −rx, R2
− xg = r2 .
∂r
∂x
Умножим первое уравнение на x, второе на r и вычтем одно из другого. Сократив
∂g
на R2 , получаем r ∂x
− x ∂g
= r. Отсюда видно, что g(x, r) = x является решением.
∂r
Тем самым, ψ = Rx , так что получается комплекснозначная функция Z −1 = 1+ix
.
R
−n
Функции Z определяются по формулам
µ
¶
1
∂n
x
−n
Z =
+i
= ϕ + iψ,
∂xn R
R
2
при этом ∆ϕ = 0 и ψ удовлетворяет (0.5). В частности, Z −2 = −x+ir
.
R3
Посмотрев на эти решения, можно подобрать такую функцию f = ϕ + iψ, что
ψ|R=R0 = const:
f = βZ 2 − µZ −2 .
³ ´1/3
2
Тогда ψ(x, r) = βr2 − µ Rr 3 , поэтому ψ|r=0 = 0, ψ|R=R0 = 0, где R0 = βµ
. Для
функции ϕ(x, r) = 2βx +
µx
R3
∂ϕ
x→∞ ∂x
выполнено lim
∂ϕ
x→∞ ∂r
= 2β, lim
= 0 (т.е. скорость
течения на бесконечности параллельна
оси x и равна 2β).
¡
¢
Если β = µ, то ϕ = βx 2 + R13 . Если ввести полярные координаты x = R cos θ,
.
то касательная к окружности радиуса R составляющая вектора скорости равна R1 ∂ϕ
∂θ
Значит, на единичной окружности v = −3β sin θ.
Тем самым, мы решили задачу обтекания шара и нашли скорость проскальзывания.
Это мы будем использовать при рассмотрении задач о пограничном слое.
СПИСОК ЛИТЕРАТУРЫ
[1] М.А. Лаврентьев, Б.В. Шабат, Проблемы гидромеханики и их математические
модели.
[2] П.Ф. Фильчаков, Приближенные методы конформных отображений. Киев,
Наукова думка, 1964.
[3] В.И. Власов, Краевые задачи в областях с криволинейной границей. М.,
Вычислительный центр АН СССР, 1987.
[4] В.С. Владимиров, В.В. Жаринов, Уравнения математической физики. М.,
Физматлит, 2000.
5
Download