Слабые пределы решений задач о движении неоднородной жидкости тема диссертации и автореферата по ВАК РФ 01.01.02, кандидат физико-математических наук Саженков, Сергей Александрович

• Саженков, Сергей Александрович
• кандидат физико-математических наук
• 1998

Оглавление диссертации кандидат физико-математических наук Саженков, Сергей Александрович

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Введение с.З

Глава 1. Задача о движении твердых тел в неньютоновской несжимаемой жидкости.

1. Постановки задач. с

2. Разрешимость задачи А о движении неньютоновской жидкости. с

3. Разрешимость задачи Б о движении твердых тел в неньютоновской жидкости. с

4. Приложение. с

Глава 2. Преобразование Лагранжа. Теоремы существования и единственности решений кинетических уравнений.

1. Постановки задач и результаты главы 2. с

2. Оператор Лагранжа. с

3. Доказательство теоремы 1.3 о лагранжевом преобразовании уравнения Л.Тартара. с

4. Доказательство теоремы 1.4 о лагранжевом преобразовании транспортного уравнения. с

5. Доказательство теоремы 1.1 о разрешимости

уравнения Л.Тартара. с

6. Доказательство теоремы 1.2 о единственности

решения транспортного уравнения. с

7. Приложение. Задача о движении неоднородной

вязкой несжимаемой жидкости с быстроосциллирующими

начальными данными. с

ЛИТЕРАТУРА с

Введение диссертации (часть автореферата) на тему «Слабые пределы решений задач о движении неоднородной жидкости»

ВВЕДЕНИЕ.

1.1 Общие положения и обзор известных результатов.

Движение сплошной среды, занимающей объем О С -К" (ограниченный или неограниченный), описывается системой динамических уравнений [1], состоящей из уравнения неразрывности (баланса массы)

Ар + а^рг?) = 0, (ж,*) еПх [о,т], (1.1)

уравнения количества движения

рдгг+ргг-= + (я,*) € О х [0,Т], (1.2)

уравнения ¿момента количества движения

=/од + Уф + (я ® ж) : Т, (ж, ¿) е О х [0,Т]. (1.3)

Здесь х - радиус-вектор точки пространства Кп, [0 ,Т] - промежуток времени, в течение которого происходит движение, р = р(х^) - плотность, й = й(х^) - вектор скорости, Т = ||гу||^=1 ~ тензор напряжений, / = /(#,£) - вектор массовых сил, Л4 - внутренний момент количества движения, д, = С^п - распределенные массовые и поверхностные пары сил (п - вектор нормали к поверхности).

При изучении движения определенной сплошной среды уравнения (1.1)-(1.3) конкретизируются заданием специфических свойств среды и вектора массовых сил, постулированием реологического уравнения

Т = Т (*,аг,гГ,(1.4)

то есть заданием зависимости тензора напряжений от кинематических величин t, х, й и некоторого набора скалярных констант {рг-}.

Математические исследования уравнений (1.1)-(1.4) мотивируются потребностями развития технологий промышленного производства, стимулируются совершенствованием численных методов решения задач математической физики и постоянным улучшением применяемой для расчетов вычислительной техники. Результаты и методы решений, получаемые при изучении проблем механики сплошных сред находят свое место в теории дифференциальных уравнений, а потому представляют и самостоятельный научный интерес.

В первой главе настоящей диссертации рассматриваются задачи, относящиеся к проблемам движения неньютоновской несжимаемой жидкости при отсутствии внутреннего момента количества движения и распределенных пар сил, то есть М, $ и тождественно равны нулю. Условие несжимаемости и уравнение неразрывности (1.1) эквивалентны уравнению переноса

Dtp + й■Чp = 0, (я,*) 6 V х [0,Т] (1.5)

и условию соленоидальности поля скоростей

аьт = о, (ж,*)' е п х [о,г]. (1.6)

Вследствие М = 0, д = 0 и ф = О уравнение (1.3) специализируется в равенство вида

Тц — Т^гч

которое выражает свойство изотропности (невыделенности направлений) среды.

Реологическое уравнение вязкой несжимаемой жидкости имеет вид

т = -р,1 + ту, (ж,*)еПх[0 ,Т], (1.8)

1, « -Ь

где I = \\8ij||"-=1, = < ' . , р* - давление в жидкости, тензор

( и, г 7=

ту - так называемая "вязкая" часть тензора напряжений. Если ту определяется законом Стокса [1],

ту = 2рВ{й), (1.9)

где £>(«) = \\liDiUj + _ тензор скоростей деформаций, ¡1 > О

- коэффициент вязкости (ц не зависит от и в явном виде), то уравнения (1.2), (1.5), (1-6), (1.8), (1.9) (заметим, что равенство (1.7) выполняется автоматически) описывают движение изотропной линейной вязкой несжимаемой жидкости и называются уравнениями Навье-Стокса. Теория этих уравнений очень обширна и включает в себя, в частности, ответы на вопросы о корректности задач в ряде функциональных классов. Отметим, что основные результаты о существовании и единственности решений в пространствах Соболева содержатся в [2], [3].

Вызывают интерес также задачи о движении твердых тел, погруженных в жидкость с реологическим уравнением в виде (1.8), (1.9). В [4], [5] рассмотрен случай, когда жидкость заполняет все пространство R2 и в ней плавает одно тело. В [6], [7] предполагается, что жидкость заполняет ограниченную область Q с гладкой границей, вектор массовых сил равен нулю и допускается, что погруженных в жидкость тел может быть несколько. Во всех этих работах доказано существование слабого обобщенного решения и в [б], [7] изучены процессы столкновения тел друг с другом и с границей В [8] проведен численный анализ экспериментов по взаимодействию твердого тела, плавающего в жидкости, и стены.

Уравнения Навье-Стокса получили широкое применение при исследовании проблем гидромеханики, но тем не менее они не являются универсальными по отношению ко всему многообразию изотропных вязких несжимаемых сред. Это приводит к появлению моделей, в которых зависимость тензора ту от D{u) нелинейна. Жидкости с реологическим уравнением (1.8) и нелинейной зависимостью rv от D(u) называются ненъютоновскими. Физическим аспектам неньютоновских жидкостей посвящены публикации [9] - [11]. Первые математические исследования по неньютоновской жидкости принадлежат О.А.Ладыженской (1966-1970 гг).

В [3], [12], [13] в предположении, что оператор D(u) —> tv(D(u)) является монотонным и удовлетворяет условиям возрастания

\tv(D{u))\< Ci(l + \D(u)\)p-\ Ci > 0, 1 < p < ос, (1.10)

rv(D(Ü)) : D(u) > C2\D(u) C2 > 0, p > Pl > p - 1, Pl > 0, (1.11)

доказано существование слабых решений при р > 1 + и единственность при г?(ж,0) 6 #(fi), р > В [14] (1970 г) доказаны существование и единственность слабых решений при таких же числах р, п в случае, когда

r^ = IVul^V«. (1.12)

Заметим, что при таком ту тензор напряжений не удовлетворяет принципу изотропности (1.7), однако, доказательства из [14] технически несложно перекладываются для изотропного тензора (см., например,

tv = \D(u)\P~2D(u)

(1.13)

который, кстати, удовлетворяет условиям (1.10), (1-11).

В [16] (1970 г.) предложена разностная схема решения задачи о движении неньютоновской жидкости.

В публикациях последних лет [17], [18] (1993-1994 гг.) рассмотрены задачи с тензором напряжений, определяемым формулами (1.8), (1.10), (1.11), и доказаны при определенных р и п разрешимость и единственность в смысле слабого обобщенного решения в случае, когда жидкость занимает все пространство Rn и начальные данные й(х, 0) и р(х,0) есть периодические функции. В этих же работах предложено понятие мерозначного решения уравнений неньютоновской жидкости. В [15] приведен весьма полный обзор этих и других результатов о наличии мерозначных, слабых и сильных обобщенных решений.

Глава 1 настоящей работы посвящена продолжению исследований математических вопросов о движении неньютоновской жидкости. Доказаны теоремы существования слабых обобщенных решений задач о движении дилатантной (dilatory (англ.) - медленный, запоздалый) жидкости и о движении твердых тел в этой жидкости. Точный вид формулировок и результатов первой главы приведен в п.1.2 этого введения.

Во второй главе настоящей диссертации исследуются уравнения вида (1.1). Они имеют общее название "кинетические уравнения" [19] и входят составной частью не только в уравнения механики сплошной среды (1.1)-(1.3), но и во многие другие модели математической физики. В силу этого, важное значение имеет изучение задачи Коши для уравнения (1.1), которая формулируется следующим образом: Пусть Qt = U х Rn - пространство переменных u(x,t), Pq(x) -

заданные функции. Требуется найти обобщенную функцию р = р(х, t), удовлетворяющую уравнению

Dtp + div(pu) = 0, (х, t) Е П х [0, Т]

(1.1)

и начальному условию

р(х, 0) = ро(х).

(1.14)

Теоремы существования и единственности классического решения задачи (1.1), (1-14) опубликованы во многих учебниках по теории обыкновенных дифференциальных уравнений и уравнений в частных производных (например, [20]). Исследованию вопросов о существовании и единственности решений в пространствах Лебега Ьр(0,), р 6 [1,оо] ■ посвящены работы. [2], [21], [22]. В [23] изучены вопросы корректности задачи (1.1), (1-14) в пространствах Орлича. В [21] предложено понятие ренормализ о ванного решения, доказаны теоремы существования и единственности ренормализованных решений и, в дальнейшем, эти результаты нашли применение при исследовании уравнений Больцмана [24], [25] и Власова-Максвелла [26].

Полученные во второй главе диссертации результаты являются продолжением исследования вопросов корректности задач Коши для кинетических уравнений вида (1.1). Подробные формулировки и постановки содержатся в п. 1.3 настоящего введения.

Структура диссертации.

Диссертация состоит из введения, главы 1, которая содержит 4 раз. дела, главы 2 (7 разделов) и списка литературы (47 названий).

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

В п. 1.4 введения приведен список основных обозначений, используемых в работе.

1.2 Содержание главы 1.

Рассматривается задача о движении неньютоновской несжимаемой жидкости в следующей постановке.

Задача А.

Жидкость занимает ограниченную область О пространства В?. Требуется найти поле скоростей и : давление р* : Я? 4 Л и вязкость жидкости ¡л : Ят —> Я, Ят = О х [0,Т], удовлетворяющие следующим уравнениям и начальным и граничным условиям:

Dtu + £ Mi A« - di vfaW) = f- (x, t) e QT (1.15)

i= 1

3

Ajt + I>iA/x = 0, (x,t)eQT (1.16)

i= 1

divw = 0, (x,t)£QT (1.17)

«(s,i)|<=o = щ(х), x e Q (1.18)

=0, ju(x,t)lt=0 = Mx) (1-^)

0 < 771 < juo(x) < M < oo, m, M = const , X e 0. (1.20)

Относительно W предполагается выполнение следующего условия: /iW £ дФ(Ю(и)), дФ(Б(й)) - субдифференциал функционала Ф(х) = \hQ(x,x(x))dx, где

0(х у(х)) = { ПРИ 1x1 - М

v 'AV \ +оо, при 1x1 > м, М = const < +оо

в точке х = D{u).

Определение 1.1. Обобщенным решением задачи А называется пара функций {u(x,t), fi(x,t)} таких, что

и £ Lp(0,T; V(Q)) П ^(О.Т;^)), \D(u(t))\ < M при п.в. t £ [0,T], ¡1 e Loo(Qt), rh < n(x,t) < M при п.в. (x,t) 6 Qt и выполняются вариационное неравенство

JQt Dt<p(<p - u)dxdt + J D{u)\p~2D(u) : D(0 - u)dxdt

- JQtU®u:V(0-ujdxdt > J f(0 - u)dxdt (1.21)

и интегральное тождество

f fi(Dtip + uV,ip)dxdt + Jn {¿(x, 0>)ф{х, 0)dx = 0. (1.22)

Здесь $ e Lp(0,T;V(tt)) - произвольное поле, удовлетворяющее условиям Dt<p G £Р(0,Г;7(П)), \D($)\ < M, £|t=0 = Щ, Ф € C\QT) -произвольная функция такая, что ip\t=r = 0.

Доказывается теорема существования обобщенного решения задачи

А:

Теорема 1.1. Пусть

/ е 1/(0,Т;У'(0)), < М 20>е = /,

гГ0€Я(П), р(гГ0)М|<М, р"1 + (У)-1 = 1,

Тогда существует обобщенное решение задачи А.

Далее, эта задача рассматривается с начальными данными для вязкости ¡л в виде

„Г* 01-1 1. *€П\Уо,

Здесь Уо есть объединение непересекающихся подобластей У0^, I = 1,..., ТУ с гладкими (липшицевыми) границами Е^.

Доказывается, что при е —> 0 последовательность обобщенных решений задачи А сходится к обобщенному решению задачи Б о движении твердых тел в неньютоновской жидкости под действием гидродинамических реакций, которая состоит в следующем. Задача Б.

Требуется найти области I = 1,..., N, моделирующие твер-

дые тела, поле скоростей й : Ят \ Уг -» Я3 и давление внутри жидкости р, : где^г = {(У(*),*),*е[0,Т]}, У{1) = ¿¡У««, удовлетворяющие системе, включающей в себя уравнения движения жидкости

= М) £ От \ Ут, (1-23)

г=1

<н™ = 0, (м) (1-24)

и уравнения Эйлера движения твердых тел под действием гидродинамических реакций в неподвижной декартовой системе координат [27]:

Ж

с

тп

йг

= + С1"25)

(«7^ + 2т^х^рг^РI — т^-и^ ® х^

-т<4° ® + (J{1) + - ®

tit

= jvV)p®{xxj)dx + J^xx(Tn)da: / = 1,..., iV. (1.26)

Замечание 1.1. Вязкость жидкости считается постоянной и равной 1 во всем объеме жидкости.

К уравнениям (1.23)-(1.26) добавляются начальные и граничные условия

И')(о) = V0{1\ Z = 1,..., iV, (1.27)

u(z,i)|i=o = u0(s), x£Q\Vq, (1.28)

u(s,i)U2 = 0, (1.29)

«(ж,*)^ x ix~xc(t))) ldvT, (1-30)

= u\t=Q = = я?, ar G F0. (1.31)

В формулах (1.25)-(1.31) обозначено E^ - поверхность твердого тела ft - нормаль к Е^ , внешняя по отношению к р^ > 0 -плотность твердого тела, в рассматриваемой задаче положим р^ = 1 для любого /, 1 < I < N, т® > 0 - масса тела, J^ - тензор инерции тела в главных осях инерции, jf-> = Sij fyw p^(xf_1 + xf+1)dx, i,j = 1,2,3, v^ - скорость центра инерции тела, х^ - радиус-вектор центра инерции тела, и^ - угловая скорость тела, Т = ||rij||ij=i _ тензор напряжений вязкой жидкости, гг-;- = — + W, р* - внутреннее давление жидкости, / - вектор внешних массовых сил.

Введем в рассмотрение функцию, характеризующую расположение твердых тел

Mxty=ll> хеУ {х* } \ о, х eci\v{t).

В терминах этой функции начальное условие (1.27) можно записать в виде

л(м)=лом={;; (1-з2)

Следуя [6], введем в рассмотрение некоторые специальные функциональные классы и сформулируем понятие обобщенного решения задачи Б :

Char(E') — класс характеристических функций подмножеств множества Е,

К(х) = [Ф в НЦП)\В(ф)(х) = О, X £ S(x)h ^ * G Char(fi), S(x) = {xett\x(x) = i}.

Определение 1.2. Обобщенным решением задачи Б называется пара функций {и,Л} таких, что й G 2^(0, Т; Н(С1)) П LP(Q,T;V(Q)), и G К (А), |£>(ЗД)| < М при п.в. t G [О, Г]; Л G Char(QT); Л G С(0,Т; L#(Q)), i9 < оо и для которых имеют место интегральное неравенство

[ Dt<p(<p - ujdxdt + [ \D{u)\p~2D{u) : D(<p - ujdxdt

J QT J QT

— I и ® и : V(<p — ujdxdt > f f(<p — ujdxdt (1.33)

J Qf J Qx

и интегральное тождество

/ A(D^ + u4ii)dxdt+( = (1.34)

jqt jq

Здесь ф G Lp(0,T;V(Q)) - произвольное векторное поле, удовлетворяющее условиям Dt<p G Lp(0,T; V(O)) ф G K(Л), < M, (p\t=Q = u0, ф G C1{Qt) ~ произвольная функция такая, что ф^=т = 0.

Теорема существования обобщенного решения задачи Б формулируется следующим образом:

Теорема 1.2. Пусть щ G Я(О) П К(А0), \D(uQ)| < М п.в. в QT,

fe Lp/(0,T;F'(Q)), \D(F)| < M, г<?е DtF = /,

f1 + (/r1 = i) p>j.

Тогда существует обобщенное решение задачи Б.

Замечание 1.2. Из того факта в определении обобщенного решения, что и G К(Л), следует, что эволюция во времени I = 1,..., N является твердотельным движением, поскольку решениями уравнения D(uj(x) = 0 являются функции и(х) = vc + и х (х — хс), где vc, хс и и не зависят от х [3, с.70].

Замечание 1.3. Класс пробных функций для неравенства (1.33) зависит от решения задачи. Как будет видно, такой выбор пробных

функций не вызывает недоразумений при доказательстве существования обобщенного решения. Заметим, что пробные функции, зависящие от решения, уже рассматривались в формулировках обобщенных решений в работах [4,6].

О характере взаимодействия тел и границы между собой. При доказательстве теоремы 1.2 показывается, что поле скоростей й(х, t) почти липшицево, откуда по теореме Осгуда следует, что траектории движения точек жидкости и тела определяются единственным образом в Qt. Таким образом, несложно видеть, что если два тела (или тело и граница) соприкасались в начальный момент времени, то они будут соприкасаться при любом t £ [О ,Т] и иметь друг относительно друга нулевую скорость, а если два тела (или тело и граница) не имели общих точек, то они не "столкнуться" ни при каком t £ (О, Т].

1.3 Содержание главы 2.

Одной из задач, связанных с уравнения переноса вида

Dtp(x, t) + v(x, t) ■ Vxp(x, t) = 0, (x, t) £ QT, (1.35)

Qt £ fi x [0, T], Q С -R2, является задача об эволюции особенностей бы-строосциллирующих решений уравнения вида (1.35), которая может быть сформулирована следующим образом. Имеется некоторая сла-босходящаяся последовательность ps решений уравнения (1.35), р* = w-\imp£. Требуется описать множество слабых пределов f*(x,t) = ии- lim f(p£(x,t)) для всевозможных непрерывных функций /. В работе

[28] Л.Тартара была предложена конструкция Я-меры, которая содержит информацию о слабых пределах последовательностей реА\р^ где А - произвольный псевдодифференциальный оператор нулевого порядка. В плоском случае ü-мера определяется при п.в. t £ [0, Т] посредством равенства

< Ни а<Р№ >= ^mfQ(pi(pe ~ Р*)АМРе ~ P*)]dx, (1.36)

в котором <pi, ip2 £ Co(Q), А : L2(jR2) -»• L2(R2) - псевдодифференциальный оператор нулевого порядка с символом а £ С1^1), где S1 -единичная окружность.

Так как множество функций {^1^2«} плотно в Со(0х5:), то, как показано в [28], эта формула определяет при п.в. t £ [0,Т] неотрицатель-

ную борелевскую меру на Q х 51. Семейство ii-мер {fit}, зависящих от t как от параметра, служит решением кинетического уравнения JI. Тартара

DtiH+ diYx(fit v) + Dy(fit Y :V^) = 0, (t,x,y) G [0,Г] x Q x S\

(1.37)

где Y(y) = \\Yij{y)\\i,j=i,i ~ матрица, Уц G С°°(Д).

Возникает вопрос об отыскании минимальных условий на поле скоростей обеспечивающих корректность задачи Коши для уравнения (1.37) в классе борелевских мер. В настоящей работе дается ответ для соленоидальных векторных полей v(x,t):

Теорема 1.1. Пусть v G 1^(0,Т; Vo(х>у) есть неотрица-

тельная борелевская мера на Q х S1, принадлежащая пространству распределений

(1/2(0, С(5Х)))* и абсолютно непрерывная относительно меры Лебега на Q,, тогда существует единственное решение pit(x,y) задачи Коши для уравнения (1.37), имеющее следующие свойства:

1). щ{х,у)\г=ъ = цо(х,у), (х,у) G Q х S1;

2). fjneHo^iHQMS1)))*);

3). fit при п.в. t G [0,Т] есть неотрицательная мера на Q х 51, абсолютно непрерывная относительно меры Лебега на Q.

Главная проблема при доказательстве теоремы 1.1 заключается в следующем. Представляя уравнение (1.37) в эквивалентном виде

DtiH + divx,y(lHV) = 0, где V - {иъ v2, Y : Wxv},

несложно видеть, что норма ||div^i2/V^jjx2(o,T';i00(П)) не ограничена, а значит невозможно применить известную технику, основанную на получении оценок для решений с использованием леммы Гронуолла (см., напр., [21]).

Аналогичная трудность возникает при изучении вопроса о единственности решения транспортного уравнения

Dtu+ divx(vu) + си = 0, (x,t) G Q х [О,Т], (1.38)

О, С Rn - ограниченная область, если норма ||c||£1(o!T;Lco(fi)) не ограничена.

Это уравнение в настоящей работе рассмотрено наряду с уравнением (1.37) и доказана теорема единственности:

Теорема 1.2. Пусть

уе 17(0,Т;<а(0)ПЯ(П)), 1 < а < оо, 1 < 7 < оо; (*)

ийеЬр(П), 1<р<оо; (**)

с€17(0,Г;£,(П)), 1 < ? < оо; (* * *)

и выполнены условия согласования д > — а?)-1, р-1 < 1 — (]У —

а)(аУУ)-1 (если а) или д-1 < 1 — р~1, р > 1 - любое (если N < а).

Тогда, если существует некоторое решение и(х^) £ £¿(0, Г; Ьр(0,))}

+ 7-1 < 1 задачи Коши для уравнения (1.38) с начальной функцией и0 £ Ьр(0,), то это решение единственное.

Идея доказательств теорем 1.1 и 1.2 состоит в использовании ла-гранжевого представления для уравнений (1.37) и (1.38) и основана на том, что более простой вид уравнений по сравнению к исходному позволяет преодолеть отмеченные выше сложности. Возникает, однако, другая проблема. Как оказалось, минимальным условием гладкости на поле скоростей у(х, ¿), при котором изучена и доказана возможность перехода от эйлеровых координат к лагранжевым, является принадлежность г; пространству 1а(0, Т; Ж22(^)) (П £ Л3) [29], а в случае менее гладкого поля скорости этот вопрос остается открытым.

В настоящей работе предложено понятие лагранжевого преобразования - обобщение понятия лагранжевого представления на случай соленоидальных полей скоростей, принадлежащих пространству £7(0,Т; ТФ^(О)), 1 < а, 7 < оо. Лагранжево преобразование реализуется посредством оператора Лагранжа, определенного следующим образом.

Рассматривается задача Коши для уравнения переноса

В^) + сПу», = 0, (*,*)€ <3г, п FV(x,s)\s=t = f(x,t), хеп, te%т}

с данными Коши / £ 1^(0,Т; £р(0)).

Согласно [21, с.516] при почти каждом £ £ [0,Т] существует единственное решение принадлежащее как функция переменных

X, 5 пространству С([0,Т]; ЬР(Щ) (если р < ос) или Ьоо(Ят) ПС([0,Т7];Ьр/(0)), р1 < со - любое (если р = со).

Ассоциированный с векторным полем V оператор Лагранжа С : Ьр(0,) —>• определяется для почти каждого £ € [О,Т] по закону

£[/]0М)=^М), (1.40)

где - решение задачи Коши (1.39).

Образ функции / при отображении С называется лагранжевым преобразованием этой функции.

Устанавливается, что

Предложение 1.1

1). £[/](ж,£) есть измеримая в Ят функция;

2). С[ДеЬ,(0,Т;Ьр(П)) ;

3). имеют место оценки:

если р < оо, то ||£[/]|Ь(0,т;£р(П)) = IIЯксод^М«))'

если р = со, то \\С[/]\\ь,{0гТ.Мп)) < Шыо.т-мп))',

4). если V в СЩО^СЦП) П Я(О)), отображение / : фт Л, / 6 С^фг), задано в эйлеровых координатах, то £[/](£,£) = [/]$(£,£)> г<9е есть представление функции / б лагранжевых координатах;

5). существует оператор С~1, обратный к оператору Лагранжа С, то есть С о Сг1 и С"1 о С есть тождественные преобразования, действующие из Ьр(0) в Ьр{0) при п.в. I Е [0,Т\, при любых р £ [1,оо]. Для оператора С~1 справедливы утверждения пп. 1).-3). настоящего предложения (в формулировках которых на месте С фигурирует С~1);

6). если у £ С1([0,Г];С01(^) П Н(П)), отображение / : дТ Я, / 6 С1(Ят), задано в лагранжевых координатах, то £_1[/](ж,= [$]х(х,{), где [/]х есть представление функции / в эйлеровых координатах.

На основании предложения 1.1 доказываются теоремы о лагранже-вом преобразовании кинетических уравнений.

Теорема 1.3. Пусть Щ(х^) = ¿), г',; = 1,2 (будем обо-

значать и(х^) = || Ом) 11?, ¿=1/

Уравнение (1.37) имеет своим решением неотрицательную абсолютно непрерывную относительно меры Лебега на П борелевскую меру щ(х,у), определенную на О х 51 при п.в. t Е [О,Т] такую, что ¡¿г Е ¿2(0, Г; (¿2(0; С(51)))*); ¿^{х^у) = <1хс1ц>х(у) тогда и только тогда, когда уравнение

+2^(17 : Ущ) = 0, (£, ж, у) Е [О, Г] х О х 51 (1.41)

имеет своим решением неотрицательную абсолютно непрерывную относительно меры Лебега на О, борелевскую меру щ Е £2(0,Т; (Ь2№такую, что й^(х,у) = <1х<1\х(у), ¡31 ф(у)с1А^х(у) = £ [/51 ^(у)^.,-(у)] ОМ) для любой функции ф Е С(Я1), при п.б. (ж, ¿) Е фт-

((1,х) —У и есть слабоизмеримые относительно

меры Лебега на С^т отображения С^т в пространство мер на 51.)

В формулировке теоремы 1.3 использован тот факт, что если /и(х,у) есть некоторая борелевская мера, определенная на произведении множеств X х У, абсолютно непрерывная относительно меры Лебега на Х1 то при любой измеримой по мере ц функции / по теореме Лебега -Никодима [30, с. 301-303] справедливо представление

< / >= ¡х<1х ¡у /(х,у)(1их(у), (т.е. <1ц(х,у) = йхйрх(у)),

где х —> рх есть слабоизмеримое относительно меры Лебега на X отображение X в пространство мер на У.

Теорема 1.4. Пусть С(х,Ь) = £[с](х,1), показатели а, -у, 8, р, д удовлетворяют условиям теоремы 1.2.

Уравнение (1.38) имеет решение и{х,Ь) Е Ьр(0,)), тогда и

только тогда, когда уравнение

Б& + Си = 0, (я,*)€фг (1-42)

имеет решение V Е Ьз(0,Т] Ьр(0,)) такое, что = С[и](х^).

Уравнения (1.41) и (1.42) будем называть лагранжевыми преобразованиями уравнений (1.37) и (1.38) соответственно.

Справедливость теорем 1.3 и 1.4 позволяет на основании результатов о существовании или единственности решений уравнений (1.41) и

(1.42) делать заключения о справедливости точно таких же результатов для уравнений (1.37) и (1.38).

В приложении (раздел 7) рассматривается начально-краевая задача для системы классических уравнений Навье-Стокса, описывающих движение неоднородной вязкой несжимаемой жидкости с быстроосцил-лирующими начальными данными в ограниченной области Q G R2. При этом считается, что значения плотности р(х, t) переносятся вдоль траекторий движения частиц жидкости со скоростью v(x,t).

Быстроосциллирующие начальные данные моделируются как пределы при £ —у 0 последовательностей начальных распределений р0£(х) = Ps{x, 0) G Lco(Q), vQe(x) = ve(x,0) G Vq(Q):

т < Рое < m, Pqs —>■ Ро * -слабо в Loo(n),

Vq£ —> Vq сильно в Vq^Q).

Известно [2], что для каждого г > 0 существует обобщенное решение {v£(x,t), p£(x,t), Vp£(x,t)}, с начальными данными {щ£(х), ро£(х)}. Доказывается, что р£ —» р *-слабо

х [0,Г]),

v£ V слабо в L2(0,T]V2(Q)), Vp£ Vp слабо в L2(Q х [0,Т]), где v(x,t), p(x,t) и Vp(x,t) составляют сильное обобщенное решение исходной задачи с начальными распределениями щ(х) и ро(х). Устана-• вливается, что Н-мера Тартара ассоциированная с последовательностью p£(x,t), служит решением кинетического уравнения (1.37), в

,, / ¿sin 2?/ —cos2 у \ котором Y — 2 • 9 1 • « • у sin- у —^sm2у J

1.4 Основные обозначения. Геометрические объекты.

Область Q - открытое множество в йп, в диссертации рассматриваются ограниченные области; х = (xi,... ,хп) - точка пространства Rn; дО, - граница области Q;

Qt = О х [0, Т] - пространственно-временной цилиндр, t G [0, T] - время;

..S1 - сфера радиуса 1 в R2: 51 = {у в R21 у2 + у\ — 1}.

Функциональные пространства.

В диссертации используются обозначения, общепринятые в литературе (см., например, [14], [19]):

Р(П), Т>((^}т) ~ пространства основных функций в О,, От; ЬР(П), Ьр(С}т), Р > 1 - пространства Лебега; ]¥р(П), / = 1,2,...,р>1,- пространства Соболева; ]¥10'р(О,) - замыкание £>(0) в норме пространства И^П);

- сопряженное пространство к И1 < р < оо; Нт(П) = ]¥?(П), #0™(П) = Ж0т'2(П), т = 1,2,...; Н-т(П) = ]¥Гт(П);

С^(Г2), Ск([0,Т]),... - пространства к раз непрерывно дифференцируемых функций в О, [О,Т],...;

Ьр(0,Т]Х) - пространство измеримых отображений / отрезка [О, Т] в нормированное пространство X, ||/||£р(0дуг) = (¡о \\1\\х^)1/р'^ ЬР11Р2(Е) - пространство функций /(#,£), £ Е, суммируемых с

разными показателями по I и £ (соответственно, с показателями р2 и Р1);

С* ([О, Т]; X) - пространство к раз непрерывно дифференцируемых отображений [О ,Т] в X;

V{Q) = {0\ipieV(Q), 1 = 1 =

= {<Д I ^ € Я*(П), г = 1,..., п, = 0};

- замыкания У(П) в нормах (И^(0))п, (Нк(П))п, £ = 1,2,..., соответственно; Т/'(0), - сопряженные к У(О) и

Я(О) - замыкание У(П) в норме (¿г^))"-

Рассматриваются также некоторые специфические пространства: СЬаг(Е), К$(х)-, Яо(х)> К(х) ~ определения этих пространств приведены в п. 1.2 введения и разделе 1 главы 1.

Обозначения норм в нормированных пространствах соответствуют общепринятым: норма элемента пространства X обозначается символом || ■ \\х (в символом || • \\р&). В диссертации также рассматриваются три специальные нормы:

|| • || - в главе 1: специальная норма в в главе 2: норма неотри-

цательной меры;

|| • Ц^) - в главе 1: специальная норма в Ьр, р > 2. Обозначения некоторых математических операторов и операций.

Г) — й. п. - А.

~ т ' и% ~ дХг '

А : В = вуЬу - свертка двух тензоров А = и В =

а ® 6 = ||аг-6у - диада двух векторов а и 6;

V, - дифференциальный оператор градиент: УгТ = ||-градиент вектора V, V?; = (И^,..., Д^) - градиент скалярной функции v]

сИу, сНу^ - дифференциальный оператор дивергенция: сНу^ = Е"=х А^ - дивергенция вектора V, (НуА = || Е"=х Аау ||"=1 - дивергенция тензора

А = \\а^=1.

1.5 Апробация работы.

Материалы диссертации докладывались на Сибирской школе-семинаре "Математические проблемы механики сплошных сред", г.Новосибирск, 1997 г.; на III Сибирском конгрессе ИНПРИМ, г.Новосибирск, 1998 г.; на I Краевой конференции по математике, г.Барнаул, 1998 г.; на семинаре кафедры дифференциальных уравнений НГУ "Теоретические и вычислительные проблемы задач математической физики", руководитель семинара - д.ф.-м.н., проф. А.М.Блохин; на семинаре кафедры прикладной математики НГУ, руководитель семинара - д.ф.-м.н., проф. Т.И.Зеленяк.

Основные результаты диссертации опубликованы в работах [43]-[47].

Автор диссертации выражает глубокую признательность члену -корреспонденту РАН профессору П.И.Плотникову за научное руководство и доценту кафедры теоретической механики НГУ В.Н.Старовойтову за полезные обсуждения при подготовке диссертации.

Список литературы диссертационного исследования кандидат физико-математических наук Саженков, Сергей Александрович, 1998 год

ЛИТЕРАТУРА.

1. Седов Л.И. Механика сплошной среды. М.: Наука, 1970.

2. Антонцев С.Н., Кажихов A.B., Монахов В.Н. Краевые задачи механики неоднородных жидкостей. Новосибирск: Наука, 1983.

3. Ладыженская O.A. Математические вопросы динамики вязкой несжимаемой жидкости. М.: Наука, 1970.

4. Юдаков Н.В. О разрешимости задачи движения твердого тела в вязкой несжимаемой жидкости. // Динамика сплошной среды. Новосибирск, 1974. Вып.18.

5. D.Serre. Chute Libre d'iin Solide dans un Fluide Visqueux Incompressible. Existence. // Japan Journal of Applied Mathematics, V.4, N1, pp.99-110, 1987.

6. K.-H. Hoffmann, V.N.Starovoitov. On a motion of a solid body in a viscous fluid. Two-dimensional case. // Advances in Mathematical Sciences and Applications. V.10, N1, 1998. или Preprint M9680, November 1996, Technische Univ. München.

7. K.-H. Hoffmann, V.N.Starovoitov. Zur Bewegung einer Kugel in einer zähen Flüssigkeit. Preprint M9681, November 1996, Technische Univ. München.

8. R. Hsu, R. Ganatos. The motion of a rigid body in viscous fluid bounded by a plane wall. // J. of Fluid Mechanics, V.207, pp.29-72, 1989.

9. W.R.Schowalter. Mechanics of Non-Newtonian Fluids. Pergamon Press, Oxford, 1978.

10. R.R.Huilgol. Continuum Mechanics of Viscoelastic Liquids. Hindusthan Publishing Corp., Dehli, 1975.

11. K.R.Rajagopal. Mechanics of Non-Newtonian fluids. В сборнике Recent Developments in Theoretical Fluid Mechanics под ред. G.P. Galdi, J.Necas. Pitman Research Notes in Mathematics, Series 291, Longman Scientific and Technical, Essex, pp.129-162, 1993.

12. O.A.Ladizhenskaya. New equations for description of the viscous incompressible fluids and solvability in the large of the boundary value problems for them. В сборнике Boundary Value Problems of Mathematical Physics V, Amer. Math. Soc., Providence, Rhode Island, 1970.

13. O.A.Ladizhenskaya. Modification of Navier-Stokes equations for the large velocity gradients. В сборнике Boundary Value Problems of Mathematical Physics and Related Aspects of Function Theory II, Consultants Bureau, New York, pp.57-59, 1970.

14. Лионе Ж.-Л. Некоторые методы решения нелинейных краевых задач. М.: Мир, 1972.

15. J.Malek, J.Necas, M.Rokyta, M.Ruzicka. Weak and Measure-valued Solutions to Evolutionary PDEs. Chapman and Hall, London, etc., 1996.

16. S.Kaniel. On the initial value problem for an incompressible fluid with nonlinear viscosity. // J. Math. Tech. Vol.19, N 8, pp.681-707, 1970.

17. H.Bellout, F.Bloom, J.Necas. Young measure-valued solutions for Non-Newtonian incompressible fluids. // Comm. in Partial Diff. Eq., N 11-12, pp.1763-1803, 1994.

18. J.Malek, J.Necas, M.Ruzicka. On the Non-Newtonian incompressible fluids. Math. Models and Methods in Appl. Sci. Vol.3, N 1, pp.35-63, 1993.

19. Владимиров B.C. Уравнения математической физики. M.: Наука, 1976.

20. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М.: Наука, 1964.

21. R.J.DiPerna, P.L.Lions. Ordinary differential equations, transport theory and Sobolev spaces. // Inventiones mathematicae 98, p.511-547, 1989.

22. B.Desjardins. Linear transport equations with values in Sobolev spaces and application to the Navier-Stokes equations. // Arch. Rational Mech. Anal., 1995.

23. Кажихов А.В., Мамонтов A.E. Об одном классе выпуклых функций и точных классах корректности задачи Коши для уравнения переноса в пространствах Орлича. // Сиб. мат. журнал, Т.39, 1998, N 4, с.831-850.

24. R.J.DiPerna, P.L.Lions. On the Cauchy problem for Boltzmann equations: global existence and weak stability. // C.R. Acad. Sci. Paris 306, pp.343-346, 1988.

25. R.J.DiPerna, P.L.Lions. On Fokker-Planck-Boltzmann equations.

Commun. Math. Phys., 1989.

26. R.J.DiPerna, P.L.Lions. Global weak solutions of Vlasov-Maxwell system. // Commun. Pure Appl. Math.

27. 1. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Механика. М.: Наука, 1988. Т.1.

28. L.Tartar. H-measures, a new approach for studying homogenisation oscillations and concentration effects in partial differential equations. // Proceedings of the Royal Society of Edinburgh, 115A, p.193-230, 1990.

29. R. Temam. Lagrange representation of Navier-Stokes Equations. // J. Diff. Eq., 1985.

30. Бурбаки H. Интегрирование. Меры, интегрирование мер. М.: Наука, 1967.

31. Решетняк Ю.Г. Теоремы устойчивости в геометрии и анализе. Новосибирск: Наука, 1982.

32. Т. Kato. On classical solution of the two-dimensional nonstationary Euler equations. // Arch. Rat. Mech. and Analysis, V.25, N.3, p.188-200, 1967.

33. Иосида К. Функциональный анализ. М.: Мир, 1967.

34. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Наука, 1976.

35. Лойцянский Л.Г. Механика жидкости и газа. М.: Физматгиз, 1970.

36. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1981.

37. Стейн И. Сингулярные интегралы и дифференциальные свойства функций. М.: Мир, 1973.

38. R.Temam. Navier-Stokes Equations. Theory and Numerical Analysis. 1977, North-Holland Pub. Co., Amsterdam - New York - Tokyo.

39. Экланд И., Темам P. Выпуклый анализ и вариационные проблемы. М.: Мир, 1979.

40. Канторович Л.В., Акилов Г.П. Функциональный анализ. М.: Наука, 1984.

41. Шилов Г.Е., Гуревич Б.Л. Интеграл, мера и производная (общая теория). М.: Наука, 1967.

42. Сансоне Дж. Обыкновенные дифференциальные уравнения. Т.1. М.: Издательство иностранной литературы, 1953.

43. Саженков С.А. Задача о движении твердых тел в неньютоновской несжимаемой жидкости. // Сибирский математический журнал, 1998 г., Т.39, N 1, стр.146 - 160.

44. Саженков С.А. Решения задачи о движении вязкой несжимаемой жидкости с быстро осциллирующими начальными данными. // Динамика сплошной среды, 1998 г., Вып. 113, 12 с. (принята к печати).

45. Саженков С.А. Решения задачи о движении вязкой несжимаемой жидкости с быстро осциллирующими начальными данными. // Тезисы докладов Сибирской школы-семинара "Математические проблемы механики сплошных сред", г.Новосибирск, 1997 г., с. 122.

46. Саженков С.А. Решения задачи о движении вязкой несжимаемой жидкости с быстро осциллирующими начальными данными. // Тезисы докладов III Сибирского конгресса ИНПРИМ, г.Новосибирск, 1998 г., с.118 - 119.

47. Саженков С.А. Представление в лагранжевых координатах в случае негладкого поля скоростей. О единственности решения транспортного уравнения. // Тезисы докладов I Краевой конференции по математике, г.Барнаул, 1998 г., с. 27 - 28.