Математические модели и методы механики сплошных сред. Учебное пособие
.pdf
4.2. НЕЛИНЕЙНАЯ ЗАДАЧА |
81 |
является линейным непрерывным функционалом, определенным на пространстве H(Ω). Действительно, линейность функционала F очевидна. Непрерывность следует из неравенства Коши — Буняковского и неравенства Пуанкаре — Фридрихса (лемма 2.2). Более точно, справедлива оценка
|(f, Φ)| ≤ f 2,Ω Φ 2,Ω ≤ cp(Ω) f 2,Ω Φ H ,
где cp(Ω) — постоянная из неравенства Пуанкаре — Фридрихса. Из этой оценки по определению нормы линейного функционала (см. определение 2.2) следует, что
F H(Ω) ≤ cp(Ω) f 2,Ω, |
(4.25) |
где через H(Ω) обозначено пространство линейных непрерывных функционалов, сопряженное к H(Ω). Значит, функционал F ограничен и, следовательно, непрерывен.
В силу этого по теореме Рисса (теорема 2.1) существует единственный элемент F H(Ω) со свойствами
[F, Φ] = F, Φ H(Ω) ≡ (f, Φ)2 |
Φ H(Ω), |
(4.26a) |
||
F H = F H(Ω) . |
|
|
|
(4.26b) |
Далее, оценим модуль второго слагаемого в левой части (4.22) |
||||
n |
|
≤ |v|2|Φ|, затем по нера- |
||
сначала по неравенству k=1 vkv · Φxk |
||||
|
|
|
|
|
венству Коши — Буняковского, |
после этого |
по первому или второму |
||
(в зависимости от размерности n) неравенству Лерэ и, наконец, по
неравенству Пуанкаре — Фридрихса, и получим |
||||||
|
n |
|
|
Ω |
|
|
|
k=1(vkv, Φxk )2 |
|
≤ |
|
|v|2|Φ| dx ≤ v 42 Φ 2 ≤ c1 v H2 Φ H , |
|
|
|
|||||
|
|
|
(4.27) |
|||
где постоянная c1 зависит только от постоянных в неравенствах Лерэ и Пуанкаре — Фридрихса.
В силу произвольности пробной функции Φ на основании этой
оценки |
3заключаем, что при каждом фиксированном v H(Ω) ли- |
|
нейный |
функционал |
|
|
n |
|
|
G: Φ → (vkv, Φxk )2 |
(4.28) |
|
k=1 |
|
3 Линейность очевидна.
82 4. СТАЦИОНАРНЫЕ ЗАДАЧИ ДИНАМИКИ ВЯЗКОЙ ЖИДКОСТИ
является непрерывным, т. е. принадлежит пространству H(Ω) . Снова применяем теорему Рисса и заключаем, что существует единственный элемент G H(Ω) со свойствами
|
n |
|
[G, Φ] = G, Φ H(Ω) ≡ |
(vkv, Φxk )2 Φ H(Ω), |
(4.29a) |
|
k=1 |
|
(4.27) |
|
|
G H = G H(Ω) ≤ |
c1 v H1/2. |
(4.29b) |
Очевидно, функционал G определяется по функции v однозначно, и значит элемент G определяется по v однозначно в силу утверждения единственности в теореме Рисса. Таким образом по формулам (4.28) и (4.29a) конструируем оператор A: H(Ω) → H(Ω), который переводит элемент v в элемент G, т. е.
G = Av v H(Ω). |
(4.30) |
Легко видеть, что оператор A является нелинейным.
С помощью представлений (4.26a), (4.29a) и (4.30) записываем интегральное равенство (4.22) в виде
[νv − Av − F, Φ] = 0 Φ H(Ω). |
(4.31) |
Введем в рассмотрение оператор AF: H(Ω) → H(Ω) по формуле
AFv := |
1 |
(Av + F) v H(Ω). |
(4.32) |
ν |
С учетом (4.32) и в силу произвольности пробной функции Φ интегральное равенство (4.31) сводится к операторному уравнению
v = AFv. |
(4.33) |
Таким образом, мы свели обобщенную формулировку задачи SNE в виде интегрального (вариационного, другими словами) уравнения (4.22) к эквивалентному ей операторному уравнению (4.33). Как очевидно из структуры этого уравнения, его решение — это неподвижная точка оператора AF. Для доказательства существования неподвижной точки прим´еним теорему Лерэ — Шаудера.
Этап 2. Сначала заметим, что H(Ω) является сепарабельным банаховым пространством: сепарабельность следует, например, из существования базиса в H(Ω), см. теорему 4.4. Значит, H(Ω) удовлетворяет требованиям теоремы Лерэ — Шаудера. Проверим, что оператор AF также удовлетворяет всем требованиям этой теоремы.
4.2. НЕЛИНЕЙНАЯ ЗАДАЧА |
83 |
Во-первых, покажем, что оператор A компактен, откуда сразу вытекает компактность оператора AF. Произвольно возьмем слабо сходящуюся в H(Ω) последовательность {vm}m=1,2,...:
vm v в H(Ω). |
(4.34) |
m→+∞ |
|
Вспомним из курса функционального анализа, что если последовательность сходится слабо в некотором пространстве X, то она ограничена в этом пространстве. В силу этого свойства, в частности, имеем, что существует постоянная c2, не зависящая от номера члена последовательности m, такая что
vm 2 ≤ c2, vm 2 ≤ c2, |
(4.35a) |
vm − v 2 ≤ c2, vm − v 2 ≤ c2. |
(4.35b) |
Еще вспомним теорему Реллиха (теорема 2.5), в силу которой из (4.34) следует сильное предельное соотношение
|
vm |
− |
v |
2 |
m−→ |
0. |
(4.36) |
||
|
|
|
|
|
|||||
|
|
|
|
|
→ |
+ |
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
Теперь, используя (4.35), неравенство Коши – Буняковского и оценку (2.18) из замечания 2.8, конструируем следующую цепочку соотношений:
|
|
|
|
|
(4.29a) |
k=1 vkmvm |
− vkv, Φxk 2 ≤ |
|
|
|
||||||
0 ≤ [Avm − Av, Φ] = |
|
|
|
|||||||||||||
|
|
n |
|
|
(4.30) |
|
n |
|
n |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
||||
≤ k=1 (vkm − vk) vm, Φxk |
+ |
k=1 (vm − v) vk, Φxk |
≤ |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Коши-Б. |
|
m |
− v 4 v |
m |
|
|
|
|
|
m |
|
|
|
|
||
≤ v |
|
|
4 Φ 2 + v |
|
− v 4 v 4 Φ 2 ≤ |
|
||||||||||
(2.18) |
|
|
|
|
|
|
|
(2.18) |
|
|
|
|
|
|
||
(4.35a) |
|
|
|
|
|
|
(4.35b) |
|
|
|
|
|
|
|||
≤ c3 Φ H vm − v 4 |
|
≤ c4 Φ H vm − v 2, |
|
|
|
|||||||||||
в которой Φ H(Ω) — это произвольная пробная функция, а постоянные c3 и c4 зависят от постоянной c2 и постоянной cr(Ω) из оценки (2.18).
Полагая Φ := Avm −Av, сокращая на Avm −Av H и учитывая (4.36), из этой цепочки неравенств выводим, что
0 |
≤ |
Avm |
− |
Av |
H |
≤ |
c4 |
|
vm |
− |
v |
2 |
m−→ |
0, |
(4.37) |
||
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
→ |
+ |
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
84 4. СТАЦИОНАРНЫЕ ЗАДАЧИ ДИНАМИКИ ВЯЗКОЙ ЖИДКОСТИ
то есть последовательность {Avm}m=1,2,... сильно сходится к Av при m → ∞. Итак, оператор A переводит слабо сходящуюся последовательность в сильно сходящуюся, а значит, является компактным (см. замечание 2.2). Отсюда следует, что AF также является компактным.
Во-вторых, построим априорную4 оценку для решения уравнения v = λAFv в H(Ω), где λ [0, 1].
Умножая скалярно в H(Ω) обе части этого уравнения на v, с учетом (4.32), (4.29a) и (4.30), получаем интегральное равенство
v H2 |
= [v, v] = |
λ |
n (vkv, vxk )2 + |
λ |
[F, v]. |
(4.38) |
||||
ν |
ν |
|||||||||
|
|
|
|
|
|
k=1 |
|
|||
Здесь для слагаемых в правой части имеем |
|
|||||||||
и |
|
|[F, v]| ≤ F H v H |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
n |
1 n |
|
∂v2 |
|
||||||
k=1(vkv, vxk )2 = |
|
k=1 Ω |
vk |
|
dx = |
|
||||
2 |
∂xk |
|
||||||||
ф-ла |
|
|
Ω |
n ∂vk |
|
|
|
Ω |
|
|
Грина |
|
1 |
|
1 |
|
|||||
= |
− |
|
v2 k=1 |
|
dx = − |
|
|
|v|2 div v dx = 0. |
||
2 |
∂xk |
2 |
||||||||
В силу этого из (4.38) выводим неравенство v 2H ≤ λν F H v H , откуда
|
λ |
(4.26b) |
λ |
(4.25) |
λ cp(Ω) |
|
(λ≤1) |
cp(Ω) |
|
||
v H ≤ |
|
F H |
= |
|
F H(Ω) |
≤ |
|
f 2 ≤ |
|
f 2. |
|
ν |
ν |
ν |
ν |
||||||||
|
|
|
|
|
|
|
|
|
|
|
(4.39) |
Таким образом получили, что любое возможное решение урав- |
|||||||||||
нения v = λAFv (λ [0,1]) лежит в шаре |
|
H(Ω). |
|
||||||||
B = v H(Ω): v H ≤ (cp(Ω)/ν) f 2 |
|
||||||||||
Итак, оператор AF: H(Ω) → H(Ω) удовлетворяет всем требованиям теоремы Лерэ — Шаудера и, следовательно, имеет неподвижную точку v H(Ω). Согласно построениям, проведенным на первом этапе доказательства, v является с.о.р. задачи SNE в смысле
4 Обратим внимание на то, что в формулировке теоремы Лерэ – Шаудера не требуется гарантий, чтобы решение уравнения v = λAFv существовало, и речь таким образом идет об условной, то есть априорной, оценке.
4.2. НЕЛИНЕЙНАЯ ЗАДАЧА |
85 |
определения 4.2, причем оценка Лерэ (4.23) следует из цепочки соотношений (4.39).
Остается заметить, что оценка Лерэ также несложно выводится непосредственно из интегрального равенства (4.22) при подстановке в это равенство решения v на место пробной функции Φ. Этот вывод аналогичен выводу оценки Лерэ из интегрального равенства (4.38). Значит, в принципе оценка Лерэ выполняется для любого возможного с.о.р. задачи SNE, а не только для того, существование которого установили выше с помощью теоремы Лерэ — Шаудера5.
Теорема доказана.
Замечание 4.5. Обратим внимание на то, что мы доказали на самом деле несколько более общий результат, чем теорема 4.7. А именно доказательство теоремы 4.7 служит обоснованием существования с.о.р. задачи SNE и в том случае, когда заданной функцией в уравнении количества движения (4.21a) является функционал F H(Ω) , не обязательно связанный с какой-либо функцией f L2(Ω)n формулой (4.24). При этом, конечно, интегральное равенство (4.22) в определении 4.2 и оценка Лерэ (4.23) принимают соответственно вид
n |
|
||
ν[v, Φ] − (vkv, Φxk )2 − F, Φ H(Ω) = 0 Φ H(Ω) |
(4.40) |
||
k=1 |
|
||
и |
|
||
v H ≤ |
1 |
F H(Ω) . |
(4.41) |
ν |
|||
Перейдем к исследованию вопроса о единственности с.о.р. задачи SNE. Докажем следующую теорему.
Теорема 4.8. (О единственности медленных течений.)
Пусть заданная распределенная массовая сила f L2(Ω)n удовлетворяет дополнительному условию
f 2 < |
ν2 |
(4.42) |
cr(Ω)2cp(Ω) , |
в котором cr(Ω) — это постоянная из оценки (2.18) в замечании 2.8 и cp(Ω) — постоянная из неравенства Пуанкаре – Фридрихса. Тогда существует не более одного с.о.р. задачи SNE в смысле определения 4.2.
5 В литературе известны и другие способы доказательства теоремы 4.7, не связанные с применением теоремы Лерэ — Шаудера, см., например, [24, 34].
86 4. СТАЦИОНАРНЫЕ ЗАДАЧИ ДИНАМИКИ ВЯЗКОЙ ЖИДКОСТИ
Доказательство. Возьмем комбинацию оценки Лерэ (4.23) и условия (4.42). Получим, что каждое возможное с.о.р. v подчиняется оценке
v H ≤ |
cp(Ω) |
f 2 < |
ν |
(4.43) |
|
|
|
. |
|||
ν |
cr(Ω)2 |
||||
Для разности u двух предполагаемых с.о.р. v1 и v2 из (4.22) выводим интегральное равенство
n n
ν[u, Φ] − (ukv2, Φxk )2 − (v1ku, Φxk )2 = 0,
k=1 k=1
в котором Φ H(Ω) — это произвольная пробная функция.
Предположим, что u = 0. Подставляя Φ := u и учитывая тож-
n
дество (v1ku, uxk )2 = 0, с помощью неравенства Коши – Буня-
k=1
ковского, оценки (2.18) в замечании 2.8 и оценки (4.43), находим
n |
n |
ν u H2 = |
(ukv2, uxk )2 + (v1ku, uxk )2 ≤ |
k=1 |
k=1 |
Коши-Б. |
(2.18) |
≤ |
u 4 v2 4 u H ≤ cr(Ω)2 u H2 v2 H |
откуда вытекает, что ν u 2H < ν u 2H .
Это неравенство абсурдно. Значит, u = 0, т. е. v1 рема доказана.
(4.43)
< ν u 2H ,
= v2, и тео-
Замечание 4.6. По аналогии с классическим понятием в гидродинамике (см., например, [12, гл. II, § 19; 25, § 9.4]), величинуv H cr(Ω)2/ν можно назвать обобщенным числом Рейнольдса. Соответственно, оценка (4.43) означает, что обобщенное число Рейнольдса меньше единицы и это можно стандартно трактовать как случай медленного течения .
Завершим изложение настоящей главы кратким обсуждением вопроса о внутренней регулярности с.о.р. задачи SNE.
Имеет место следующая теорема, по своей сути аналогичная теоремам 3.2 и 4.2.
Теорема 4.9. Пусть выполнены условия теоремы 4.7. Пусть Ω C2 — строго внутренняя подобласть области Ω. Тогда с.о.р. v задачи SNE принадлежит пространству Соболева H2(Ω )n.
4.2. НЕЛИНЕЙНАЯ ЗАДАЧА |
87 |
Доказательство теоремы 4.9 можно найти в [9, гл. V, § 5, теорема 6]. Так же как и обоснование теорем 3.2 и 4.2, оно основано на специальном выборе пробной функции Φ в вариационной формулировке задачи, т. е. в интегральном равенстве в определении с.о.р. Однако стоит отметить, что доказательство теоремы 4.9 не сводится к модификации доказательств теорем 3.2 и 4.2: оно технически гораздо более сложное и существенно базируется на теории гидродинамических потенциалов, которая лежит за пределами нашего пособия.
Наконец, заметим, что в силу теорем 4.7 и 4.9 результат теоремы 4.5 перестает быть условным, и в условиях теоремы 4.7 можем восстановить давление p H1(Ω ) в любой строго внутренней подобласти Ω Ω. Соответственно, система уравнений Навье – Стокса (4.21a), (4.21b) выполняется почти всюду в Ω .
ГЛАВА 5
Нестационарные задачи динамики вязкой несжимаемой жидкости
В этой главе исследуем начально-краевые задачи для уравнений Стокса и Навье – Стокса вязкой несжимаемой жидкости. В этих задачах из числа независимых переменных существенно выделяется временн´ая переменная t, и свойства регулярности решений и заданных входных данных различаются для t и пространственных переменных xi. В связи с этим начнем изложение с двух параграфов, имеющих техническое значение: в них мы представим широко применяющиеся в теории нестационарных задач (часто говорят: эволюционных задач) леммы Гронуолла и понятие о пространствах Бохнера — пространствах с выделенной переменной времени t.
5.1. Леммы Гронуолла
Приведем две простые леммы, которые часто используются при получении априорных оценок решений нестационарных задач и для обоснования единственности решений таких задач.
Лемма 5.1. (Дифференциальная форма леммы Гронуолла.) Пусть неотрицательная абсолютно непрерывная функция y: [0, T ] → R удовлетворяет при почти всех t [0, T ] дифференциальному неравенству
|
|
dy(t) |
≤ C1(t)y(t) + C2(t), |
|
||
|
|
|
|
|
||
|
|
|
dt |
|
||
где C1 и |
C2 — неотрицательные функции |
из |
||||
L1(0, T ). Тогда справедлива оценка |
|
|
||||
y(t) ≤ e |
0t C1(τ)dτ |
y(0) + t C2(τ) e− 0τ C1 |
(ξ)dξ dτ |
, |
||
|
|
|
0 |
|
|
|
(5.1)
пространства
t [0, T ]. (5.2)
89
90 5. НЕСТАЦИОНАРНЫЕ ЗАДАЧИ ДЛЯ ВЯЗКОЙ ЖИДКОСТИ
Доказательство. Действительно, (5.1) можно записать в ви-
де
dtd e− 0t C1(τ)dτ y(t) ≤ C2(t) e− 0t C1(ξ)dξ ,
а затем проинтегрировать в пределах от нуля до t, что и дает (5.2).
Лемма 5.2. (Интегральная форма леммы Гронуолла.)
Пусть y: [0, T ] → R — непрерывная неотрицательная функция, C1 и C2 — неотрицательные функции из пространства L1(0, T ) и C — неотрицательная постоянная. Тогда из неравенства
y(t) ≤ C + t |
C1(τ)y(τ) + C2(τ) dτ |
(5.3) |
|||||
следует неравенство |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||
y(t) ≤ e 0t C1(τ)dτ C + t C2(τ) e− 0τ C1 |
(ξ)dξ dτ , t [0, T ]. |
(5.4) |
|||||
|
0 |
|
|
|
|
||
Доказательство. Обозначим правую (б´ольшую) часть нера- |
|||||||
венства (5.3) через x(t). Тогда (5.3) принимает вид |
|
||||||
dx(t) |
≤ C1(t)x(t) + C2(t). |
|
|||||
|
|
|
|
|
|||
|
dt |
|
|
|
|||
Остается применить лемму 5.1. |
|
|
|||||
5.2. Пространства Бохнера
Введем определения и изложим основные свойства пространств функций с выделенной переменной t, которые далее систематически используются в этой главе.
5.2.1.Определение пространств Бохнера. Пусть Ω Rn
иT > 0. Через QT обозначим пространственно-временной цилиндр
Ω × (0, T ). Для u: QT → Rm (m N) отображение u(t): x → u(x, t) при фиксированном значении t (0, T ) — это элемент некоторо-
го функционального пространства, например, пространства Лебега или Соболева, т. е. функция t → u(t) отображает интервал (0, T ) в некоторое пространство функций. Пространства функций, отображающих интервал времени (0, T ) в некоторое банахово пространство, называются пространствами Бохнера.
