close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
УДК 517.711.3
ЗАДАЧА ТЕРМОКОНВЕКЦИИ
ДЛЯ ЛИНЕАРИЗОВАННОЙ МОДЕЛИ НЕСЖИМАЕМОЙ
ВЯЗКОУПРУГОЙ ЖИДКОСТИ НЕНУЛЕВОГО ПОРЯДКА
Т.Г. Сукачева
THE THERMOCONVECTION PROBLEM
FOR THE LINEARIZIED MODEL OF THE
INCOMPRESSIBLE VISCOELASTIC FLUID
OF THE NONZERO ORDER
T.G. Sukacheva
Рассматривается линеаризованная модель термоконвекции несжимаемой вязкоупругой жидкости ненулевого порядка. На основе теории относительно p-секториальных операторов и вырожденных полугрупп операторов доказана теорема существования единственного решения задачи
Коши–Дирихле для соответствующей системы дифференциальных уравнений в частных производных, и получено описание расширенного фазового пространства указанной задачи.
Ключевые слова: уравнение соболевского типа, несжимаемая вязкоупругая жидкость, относительно p-секториальный оператор, расширенное фазовое пространство.
The Cauchy – Dirichlet problem for the linearized system modeling
thermoconvection of the incompressible viscoelastic fluid of the nonzero order
is considered. This problem is investigated on the base of the theory of
relatively p-sectorial operators and degenerative semi-groups of operators. The
theorem of the existence of the unique solution of this problem is proved and
the description of its extended phase space is received.
Keywords: Sobolev type equation, an incompressible viscoelastic fluuid,
relatively p-sectorial operator, extended phase space.
Введение
Система уравнений
40

(1 − λ∇2 )vt = ν∇2 v − (ṽ · ∇)v − (v · ∇)ṽ+



k

X



βl ∇2 wl − gγθ − ∇p + f ,




 l=1
0 = ∇ · v,




∂wl



= v + αl wl , αl ∈ R− , l = 1, k,


∂t



θt = κ∇2 θ − v · ∇θ + v · γ
(1)
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
получена на основе линеаризованной системы Осколкова ненулевого порядка [1, 2] и приближенного уравнения теплопроводности, моделирующего термоконвекцию несжимаемой
вязкоупругой жидкости. Она моделирует эволюцию скорости v = (v1 , v2 , . . . , vn ), vi =
vi (x, t), градиента давления ∇p = (p1 , p2 , . . . , pn ), pi = pi (x, t) и температуры θ = θ(x, t)
неньютоновской жидкости — несжимаемой вязкоупругой жидкости Кельвина – Фойгта ненулевого порядка. Параметры λ ∈ R, ν ∈ R+ и κ ∈ R+ характеризуют упругость, вязкость
и теплопроводность жидкости соответственно; g ∈ R+ — ускорение свободного падения;
вектор γ = (0, . . . , 0, 1) — орт в Rn . Параметры βl ∈ R+ , l = 1, k определяют время
ретардации (запаздывания) давления. Свободный член f = (f1 , . . . , fn ), fi = fi (x, t), отвечает внешнему воздействию на жидкость, а вектор-функция ṽ = (ṽ1 , . . . , ṽn ), ṽk = ṽk (x)
соответствует стационарному решению исходной системы [1, 2]. В [3] содержится обоснование линеаризованной модели несжимаемой вязкоупругой жидкости Кельвина – Фойгта
нулевого порядка.
Рассмотрим первую начально-краевую задачу
v(x, 0) = v0 (x),
wl (x, 0) = wl0 (x),
θ(x, 0) = θ0 (x),
∀x ∈ Ω;
v(x, t) = 0,
wl (x, t) = 0,
θ(x, t) = 0,
∀(x, t) ∈ ∂Ω × R+
(2)
для системы (1). Здесь Ω ⊂ Rn , n = 2, 3, 4 — ограниченная область с границей ∂Ω класса
C∞.
Впервые задачу термоконвекции для несжимаемых вязкоупругих жидкостей Кельвина
– Фойгта поставил А.П.Осколков [4]. Им же была исследована разрешимость задачи (2) для
соответствующей нелинейной модели нулевого порядка в случае λ−1 > −λ1 ( λ1 – наименьшее собственное значение задачи Дирихле для оператора Лапласа в области Ω)[5]. Первая
начально-краевая задача для этой модели рассматривалась в [3, 6], а для ее модификации
на случай плоско-параллельного течения в [7]. В этих работах изучалась ситуация, когда
свободный член f не зависит от времени, а в [8] — указанная неавтономная задача. Нестационарная линеаризованная модель нулевого порядка изучалась в [9], а ее обобщение в [10].
Нашей целью является изучение разрешимости задачи (1), (2) при нестационарном свободном члене f = f (x, t). Эту задачу мы исследуем в рамках теории полулинейных уравнений соболевского типа. Основным инстументом исследования является понятие относительно p-секториального оператора и порожденной им разрешающей вырожденной полугруппы
операторов. Доказана теорема существования единственного решения указанной задачи и
получено описание ее расширенного фазового пространства. Статья состоит из трех частей.
В первой части приводятся известные результаты из теории полулинейных уравнений соболевского типа, необходимые нам в дальнейшем [3, 8]. Во второй части проводится редукция
задачи (1), (2) к полулинейному уравнению соболевского типа. В третьей части устанавливается существование квазистационарных полутраекторий указанной задачи и описывается
ее расширенное фазовое пространство. Отметим, что результаты этой статьи обобщают результаты [11] .
1. Полулинейные уравнения соболевского типа
Пусть U и F — банаховы пространства, оператор L ∈ L(U; F), т.е. линеен и непрерывен,
ker L 6= {0}; оператор M : dom M → F линеен, замкнут и плотно определен в U, т.е.
M ∈ Cl(U; F). Через UM обозначим линеал dom M, снабженный нормой графика k| · k| =
kM · kF + k · kU . Пусть оператор F ∈ C ∞ (UM ; F), функция f ∈ C ∞ (R̄+ ; F).
Серия
≪
Математическое моделирование и программирование≫, вып. 10
41
Т.Г. Сукачева
Рассмотрим задачу Коши
u(0) = u0
(3)
для полулинейного нестационарного уравнения соболевского типа
L u̇ = M u + F (u) + f (t).
(4)
Определение 1. Локальным решением (далее просто — решением) задачи (3), (4) назовем
вектор-функцию u ∈ C ∞ ((0, T ); UM ), удовлетворяющую уравнению (4) и такую, что u(t) →
u0 при t → 0 + .
Введем в рассмотрение L-резольвентное множество ρL (M ) = {µ ∈ C : (µL − M )−1 ∈
L(F; U)} и L-спектр σ L (M ) = C \ ρL (M ) оператора M.
Определение 2. Оператор M называется (L, p)-секториальным, если существуют конπ
станты a ∈ R, k ∈ R+ , Θ ∈ ( , π) такие, что
2
L (M ) = {µ ∈ C : | arg(µ − a)| < Θ , µ 6= a } ⊂ ρL (M ) ;
(i) SΘ,a
k
L
Qp
(ii) max{ k R(µ,p)
(M ) kL(U ) ,
k LL
(µ,p) (M ) kL(F ) } ≤
q=0 |µq − a|
L
при любых µ, µ0 , µ1 , . . . , µp ∈ SΘ,a (M ).
Здесь
L
R(µ,p)
(M ) = Πpq=0 RµLq (M ),
p
L
LL
(µ,p) (M ) = Πq=0 Lµq (M )
соответственно правая и левая (L, p)-резольвенты оператора
−1 [12].
M )−1 L, LL
µ (M ) = L(µL − M )
M,
RµL (M ) = (µL −
Определение 3. Оператор M называется сильно (L, p)-секториальным, если он (L, p)L (M )
секториален и при всех µ, µ0 , . . . , µp ∈ SΘ,a
L
(i) k M R(µ,p)
(M ) ( µL − M )−1 f kF ≤
при всех f из некоторого плотного в F линеала;
(ii) k ( µL − M )−1 LL
(µ,p) (M ) kL(F ; U ) ≤
const(f )
Q
|µ − a| pq=0 |µq − a|
|µ − a|
const
Qp
.
q=0 |µq − a|
Замечание 1. Если p = 0, то (L, p)- и сильно (L, p)-секториальный оператор M называется
соответственно L- и сильно L-секториальным [3].
Будем рассматривать задачу (3), (4) в предположении, что оператор M сильно (L, p)секториален. При условии сильной (L, p)-секториальности оператора M решение задачи (3),
(4) может быть неединственным, что показывает пример, приведенный в [13]. Поэтому сузим
понятие решения уравнения (4). Также известно [14] – [16], что решения задачи (3), (4)
существуют не для всех u0 ∈ UM . Поэтому введем два определения.
Определение 4. Множество B t ⊂ UM × R̄+ назовем расширенным фазовым пространством уравнения (4), если для любой точки u0 ∈ UM такой, что (u0 , 0) ∈ B0 существует
единственное решение задачи (3), (4), причем (u(t), t) ∈ Bt .
Замечание 2. Если B t = B × R̄+ , где B ⊂ UM , то множество B называется фазовым пространством уравнения (4). Ранее вместо термина ≪расширенное фазовое пространство≫
42
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
использовался термин ≪конфигурационное пространство≫ [8], что вносило некоторую путаницу в терминологию [17].
Определение 5. Пусть пространство U расщепляется в прямую сумму U = U0 ⊕ U1
так, чтобы ker L ⊂ U0 . Решение u = v + w, где v(t) ∈ U0 , а w(t) ∈ U1 при всех t ∈ (0, T ),
уравнения (4) назовем квазистационарной полутраекторией, если Lv̇ ≡ 0.
Замечание 3. Понятие квазистационарной полутраектории обобщает понятие квазистационарной траектории, введенное для динамического случая [13, 15, 16].
В силу того, что оператор M сильно (L, p)-секториален, пространства U и F расщепляются в прямые суммы U = U 0 ⊕ U 1 , F = F 0 ⊕ F 1 [12], где
U 0 = {ϕ ∈ U : U t ϕ = 0 ∃t ∈ R+ },
F 0 = {ψ ∈ F : F t ψ = 0 ∃t ∈ R+ } —
ядра, а
U 1 = {u ∈ U : lim U t u = u},
t→0+
F 1 = {f ∈ F : lim F t f = f } —
t→0+
образы аналитических разрешающих полугрупп
Z
Z
1
1
L
µt
t
t
R (M )e dµ,
F =
LL (M )eµt dµ,
U =
2πi Γ µ
2πi Γ µ
L (M ) — контур такой, что arg µ → ±Θ при |µ| → +∞) линейного однородного
(Γ ⊂ SΘ,a
уравнения Lu̇ = M u.
Обозначим через Lk (Mk ) сужение оператора L(M ) на U k (U k ∩dom M ), k = 0, 1. Тогда
Lk : U k → F k , Mk : U k ∩ dom M → F k , k = 0, 1, причем сужения M0 и L1 операторов M
и L на пространства U 0 ∩ dom M и U 1 соответственно являются линейными непрерывными
операторами и имеют ограниченные обратные операторы. Эти утверждения следуют из
соответствующих результатов [12]. Поэтому приведем уравнение (4) к эквивалентной форме
Ru̇0 = u0 + G(u) + g(t)
u0 (0) = u00 ,
u̇1 = Su1 + H(u) + h(t)
u1 (0) = u10 ,
(5)
−1
где uk ∈ U k , k = 0, 1, u = u0 + u1 , операторы R = M0−1 L0 , S = L−1
1 M1 , G = M0 (I −
−1
−1
−1
Q)F, H = L1 QF, g = M0 (I − Q)f, h = L1 Qf. Здесь Q ∈ L(F )(≡ L(F ; F )) — проектор,
расщепляющий пространство F требуемым образом.
Определение 6. Систему уравнений (5) назовем нормальной формой уравнения (4).
Замечание 4. В случае, когда оператор M сильно L-секториален, нормальная форма
уравнения (4) (в случае f (t) ≡ 0) имеет вид (5.1) в [7].
В дальнейшем ограничимся изучением таких квазистационарных полутраекторий уравнения (4), для которых Ru̇0 ≡ 0. Для этого предположим, что оператор R — бирасщепляющий [18], т.е. его ядро ker R и образ im R дополняемы в пространстве U. Положим
U 00 = ker R, а через U 01 = U 0 ⊖ U 00 обозначим некоторое дополнение к подпространству
U 00 . Тогда первое уравнение нормальной формы (5) редуцируется к виду
Ru̇01 = u00 + u01 + G(u) + g(t),
(6)
Математическое моделирование и программирование≫, вып. 10
43
где u = u00 + u01 + u1 .
Серия
≪
Т.Г. Сукачева
Теорема 1. Пусть оператор M сильно (L, p)-секториален, а оператор R – бирасщепляющий. Пусть существует квазистационарная полутраектория u = u(t) уравнения (4).
Тогда она удовлетворяет соотношениям
0 = u00 + u01 + G(u) + g(t),
u01 = const.
(7)
Доказательство. Первое соотношение вытекает из (6) в силу требования квазистационарности Ru̇0 = Ru̇01 ≡ 0. Второе соотношение вытекает из тождества Ru̇01 ≡ 0, так как по
теореме Банаха об обратном операторе сужение оператора QR R(I − PR ) на U 01 есть непрерывно обратимый оператор. Здесь QR и PR — проекторы на im R и ker R соответственно,
ker PR = U 01 .
Замечание 5. Второе соотношение в (7) поясняет смысл термина ≪квазистацинарные
полутраектории≫, т.е. это такие полутраектории, которые ≪стационарны по некоторым переменным≫. Другими словами, квазистационарная полутраектория обязательно лежит в
некоторой плоскости (I − PR )u0 = const.
Теорема 1 устанавливает необходимые условия существования квазистационарной полутраектории уравнения (4). Перейдем к рассмотрению достаточных условий. Известно,
что при условии сильной (L, p)-секториальности оператора M оператор S секториален [12].
Значит, он порождает на U 1 аналитическую полугруппу, которую мы обозначим через
{U1t : t ≥ 0}, так как в действительности оператор U1t есть сужение оператора U t на U 1 .
Из того, что U = U 0 ⊕ U 1 следует, что существует проектор P ∈ L(U), соответствующий
данному расщеплению. Можно показать, что P ∈ L(UM ) [12]. Тогда пространство UM рас0 ⊕ U 1 так, что вложение U k ⊂ U k ,
щепляется в прямую сумму UM = UM
k = 0, 1, плотно
M
M
′
и непрерывно. Символом Av обозначена производная Фреше в точке v ∈ V оператора A,
определенного на некотором банаховом пространстве V.
Теорема 2. Пусть оператор M сильно (L, p)-секториален, оператор R — бирасщепляющий, оператор F ∈ C ∞ (UM ; F), а вектор-функция f ∈ C ∞ (R+ ; F). Пусть
(i) в некоторой окрестности Ou0 ⊂ UM точки u0 выполнено соотношение
00
01
1
0 = u01
0 + (I − PR )(G(u + u0 + u ) + g(t));
(8)
0 ), и оператор I + P G′ : U 00 → U 00 — топлинейный изомор(ii) проектор PR ∈ L(UM
R u0
M
M
0
00 = U ∩ U 00 );
физм (UM
M
(iii) для аналитической полугруппы {U1t : t ≥ 0} выполнено соотношение
Z
τ
0
kU1t kL(U 1 ; U 1 ) dt < ∞
M
∀τ ∈ R+ .
(9)
Тогда существует единственное решение задачи (3), (4), являющееся квазистационарной полутраекторией уравнения (4).
Доказательство. Рассмотрим окрестность Ou0 точки u0 . В этой окрестности первое уравнение (5) приобретет вид
1
0 = u00 + PR (G(u00 + u01
0 + u ) + g(t))
44
(10)
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
в силу условия (i). Далее, из (i) в силу теоремы о неявной функции существуют окрестности
00 , (U 00 = U 00 ∩ U ) O
1
1
1
00
1
Ou00
⊂ UM
M
u10 ⊂ UM (UM = U ∩ UM ) точек u0 = PR (I − P )u0 , u0
M
0
∞
соответственно и отображение δ : Ou10 → Ou00
класса C такое, что уравнение
0
u00 = δ(u1 , t)
(11)
эквивалентно уравнению (10).
Теперь в силу (11) второе уравнение (5) в окрестности Ou10 приобретет вид
1
u̇1 = Su1 + H(δ(u1 ) + u01
0 + u ) + h(t),
(12)
1
∞ по построению.
где оператор H((I + δ)(·) + u01
0 ) : Ou10 → U принадлежит классу C
Для доказательства однозначной разрешимости задачи u1 (0) = u10 для уравнения (12)
воспользуемся методом Соболевского–Танабэ, изложенным в [19, глава 9]. В силу (iii), гладкости оператора H и вектор-функции h все условия теорем 9.4, 9.6 и 9.7 в [19] выполне1 , то при некотором T ∈ R существует единственное решение
ны. Поэтому если u10 ∈ UM
+
1
1
1 .
u = u (t), t ∈ [0, T ) уравнения (12) такое, что u1 (t) → u10 при t → 0+ в топологии UM
Итак, решение задачи (3), (4) в данном случае будет иметь вид u = u1 + δ(u1 ) + u01
0 , и
это решение будет квазистационарной полутраекторией по построению.
Замечание 6. Для любой квазистационарной полутраектории уравнения (4) соотношение
(8) непосредственно вытекает из первого уравнения (7).
Замечание 7. Условие (9) для обычных аналитических полугрупп, имеющих оценку
kU1t kL(U 1 ; U 1 ) < t−1 const, не выполняется. В дальнейшем мы собираемся использовать теоM
рему 2 именно в такой ситуации, и потому необходимо сделать некоторые пояснения. Пусть
1 ] , α ∈ [0, 1] — некоторое интерполяционное пространство, построенное по опеUα1 = [U 1 ; UM
α
1 ; F), вектор-функция f ∈ C ∞ (R̄ ; F)
ратору S. В теореме 2 условие ≪оператор F ∈ C ∞ (UM
+
∞ (U 1 ; U 1 ), h ∈ C ∞ (R ; U 1 ),≫ а соотношение (9)
≪дополним условием≫ оператор H ∈ C
+
α
α
M
заменим соотношением
Z τ
∀τ ∈ R+ .
(13)
kU1t kL(U 1 ; Uα1 ) dt < ∞
0
Тогда утверждение теоремы 2 не изменится. Обсуждение этого круга вопросов см. в [19,
глава 9 ].
Замечание 8. Пусть выполнены условия теоремы 2 (возможно, с учетом замечания 7).
Построим плоскость A = {u ∈ UM : (I − PR )(I − P )u = u01
0 } и множество M = {u ∈ UM :
PR ((I − P )u + G(u) + g(t)) = 0}. По условию теоремы их пересечение A ∩ M 6= ∅, так
как содержит по крайней мере точку u0 . Более того, существует C ∞ -диффеоморфизм I + δ,
отображающий окрестность Ou10 на некоторую окрестность Ou0 ⊂ A ∩ M. Следовательно, в
качестве начального значения можно брать не только точку u0 , но и любую из некоторой ее
окрестности Ou0 . Это значит, что Ou0 является частью расширенного фазового пространства
B t уравнения (4).
Теперь пусть Uk и Fk — банаховы пространства, операторы Ak ∈ L(Uk , Fk ), а операторы
Bk : dom Bk → F линейны и замкнуты с областями определений dom Bk плотными в
Uk , k = 1, 2. Построим пространства U = U1 × U2 , F = F1 × F2 и операторы L = A1 ⊗
A2 , M = B1 ⊗ B2 . По построению оператор L ∈ L(U; F), а оператор M : dom M → F
линеен, замкнут и плотно определен, dom M = dom B1 × dom B2 .
Теорема 3. [20] Пусть операторы Bk сильно (Ak , pk )-секториальны, k = 1, 2; причем
p1 ≥ p2 . Тогда оператор M сильно (L, p1 )-секториален.
Серия
≪
Математическое моделирование и программирование≫, вып. 10
45
Т.Г. Сукачева
2. Редукция к полулинейному уравнению соболевского типа
Рассмотрим задачу (2) для системы (1), представленной в виде


(1 − λ∇2 )vt = ν∇2 v − (ṽ · ∇)v − (v · ∇)ṽ+



k

X



βl ∇2 wl − gγθ − p + f ,




 l=1
0 = ∇(∇ · v),





∂wl


= v + αl wl , αl ∈ R− , l = 1, k,


∂t



 θ = κ∇2 θ − v · ∇θ + v · γ.
(14)
t
Здесь p = ∇p, т.к. во многих гидродинамических задачах знание градиента давления
предпочтительнее, чем знание давления [21]. Впервые такая замена уравнения неразрывности сделана в [22]. Нас будет интересовать локальная однозначная разрешимость задачи
(14), (2), эквивалентной исходной задаче (1), (2). Эту задачу удобно рассматривать в рамках
теории полулинейных уравнений соболевского типа, изложенной вкратце в п.1.
Для того, чтобы редуцировать задачу (14), (2) к задаче (3), (4) введем, следуя [22], пространства H2σ , H2π , Hσ и Hπ . Здесь H2σ и Hσ — подпростанства соленоидальных функций в
n
◦
пространствах (W22 (Ω)) ∩(W21 (Ω))n и (L2 (Ω))n соответственно, а H2π и Hπ — их ортогональные (в смысле (L2 (Ω))n ) дополнения. Через Σ обозначим ортопроектор на Hσ , причем его
n
◦
сужение на пространство (W22 (Ω)) ∩(W21 (Ω))n будем обозначать тем же символом. Положим
Π = I − Σ.
Формулой A = ∇2 En : H2σ ⊕ H2π → Hσ ⊕ Hπ , где En — единичная матрица порядка n,
зададим линейный непрерывный оператор с дискретным конечнократным спектром σ(A) ⊂
R, сгущающимся лишь на −∞. Формулой B : v → ∇(∇ · v) зададим линейный непрерывный
сюръективный оператор B : H2σ ⊕ H2π → Hπ с ядром ker B = H2σ .
◦
Положим U10 = H2σ × H2π × Hp , F10 = Hσ × Hπ × Hp , где Hp = Hπ ; U1i = H2 ∩ H1 =
H2σ ×H2π , и F1i = L2 = Hσ ×Hπ , i = 1, k. Тогда пространства U1 = ⊕kl=0 U1l , F1 = ⊕kl=0 F1l .
Операторы A1 и B1 : U1 → F1 определим формулами A1 = diag [Â1 , Ek ] , где
Â1 =
Ǎ1 O
O O
,
Ǎ1 =
Σ(I − λA)Σ ΣA(I − λA)Π
Π(I − λA)Σ ΠA(I − λA)Π
;
2
B1 = (B1ij )i,j=1 , где
B111


νΣA νΣA O
=  νΠA νΠA −I  ,
O
B
O
B112


β1 ΣA . . . βk ΣA
=  β1 ΠA . . . βk ΠA  ,
O
...
O
B121 содержит k строк вида (I, I, O), B122 = diag [α1 , . . . , αk ].
Замечание 9. Пространство U1 (F1 ) определяется точно так же, как пространство U (F)
в модели [23], а оператор A1 (B1 ) совпадает с оператором L (M1 ) в [23].
Замечание 10.
Обозначим через Aσ сужение оператора ΣA на H2σ . По теореме
Солонникова-Воровича-Юдовича спектр σ(Aσ ) вещественен, дискретен, конечнократен и
сгущается лишь на −∞.
46
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
Теорема 4. (i) Операторы A1 , B1 ∈ L(U1 ; F1 ), и, если λ−1 6∈ σ(A), то оператор A1 —
бирасщепляющий, ker A1 = {0} × {0} × Hp × {0} × . . . × {0}, im A1 = Hσ × Hπ × {0} × F1 ×
|
{z
}
k
F2 × . . . × Fk .
(ii) Если λ−1 6∈ σ(A) ∪ σ(Aσ ), то оператор B1 (A1 , 1)-ограничен.
Замечание 11. Доказательство теоремы 4 приведено в [8], только в другой терминологии.
Впервые понятие относительно ограниченного оператора введенно в [24]. Случай относительно секториального оператора рассматривался в [7, 25, 26].
Далее положим U2 = F2 = L2 (Ω) и формулой B2 = κ∇2 : dom B2 → F2 определим
◦
линейный замкнутый и плотно определенный оператор B2 , dom B2 = W22 (Ω) ∩ W21 (Ω). Если оператор A2 положить равным I, то в силу секториальности оператора B2 [27, гл. 1]
справедлива
Теорема 5. Оператор B2 сильно A2 -секториален.
Положим U = U1 × U2 , F = F1 × F2 . Вектор u пространства U имеет вид u =
col (uσ , uπ , up , w1 , . . . , wk , uθ ), где col (uσ , uπ , up , w1 , . . . , wk ) ∈ U1 , а uθ ∈ U2 . Здесь uσ =
Σv , uπ = (I − Σ)v = Πv , up = p̄. Аналогичный вид имеет вектор f ∈ F. Операторы L и
M определим формулами L = A1 ⊗ A2 и M = B1 ⊗ B2 . Оператор L ∈ L(U; F), а оператор
M : dom M → F линеен, замкнут и плотно определен, dom M = U1 × dom B2 .
Из теоремы 4 и замечания 2.1.1 [3] следует, что оператор B1 сильно (A1 , 1)-секториален.
В силу этого и теорем 3, 5 справедлива
Теорема 6. Пусть λ−1 6∈ σ(A), тогда оператор M сильно (L, 1)-секториален.
Перейдем к построению нелинейного оператора F. В данном случае его удобно представить в виде F = F1 ⊗ F2 , где F1 = F1 (uσ , uπ , uθ ) = col (−Σ(((ũσ + ũπ ) · ∇)(uσ + uπ ) − ((uσ +
uπ ) · ∇)(ũσ + ũπ ) + gγuθ ), −Π(((ũσ + ũπ ) · ∇)(uσ + uπ ) − ((uσ + uπ ) · ∇)(ũσ + ũπ ) + gγuθ ), 0, . . . , 0),
| {z }
k+1
а F2 = F2 (uσ , uπ , uθ ) = (uσ + uπ ) · (γ − ∇uθ ).
Формально найдем производную Фреше Fu′ оператора F в точке u,

Σa(uσ , uπ )
Σa(uσ , uπ )

Πa(u
,
u
)
Πa(u
σ
π
σ , uπ )
Fu′ = 

O
O
(γ − ∇uθ ) · (∗) (γ − ∇uθ ) · (∗)

O
−gΣγ

O
−gΠγ
,

O
O
O −(uσ + uπ ) · (∗)
где a(uσ , uπ ) = −((∗) · ∇)(ũσ + ũπ ) − ((ũσ + ũπ ) · ∇)(∗) , а на место символа * следует ставить
соответствующую координату вектора v в случае, когда мы хотим найти вектор Fu′ v.
Далее, в нашем случае пространство UM = U1 ×dom B2 (в силу непрерывности оператора
B1 ). Используя стандартную технику (см., например, [15, 16]), нетрудно показать, что
при любых u ∈ UM оператор Fu′ ∈ L(UM ; F). Аналогично устанавливается, что вторая
производная Фреше Fu′′ оператора F — непрерывный билинейный оператор из UM × UM в
F, а Fu′′′ ≡ O. Таким образом, справедлива
Теорема 7. Оператор F ∈ C ∞ (UM ; F).
Вектор-функцию f представим в виде f = f1 ⊗ f2 , где f1 = col(Σf , Πf , 0, . . . , 0), f2 = 0.
| {z }
k+1
Будем предполагать, что f ∈ C ∞ (R̄+ ; F). Итак, редукция задачи (14), (2) к задаче (3), (4)
Серия
≪
Математическое моделирование и программирование≫, вып. 10
47
Т.Г. Сукачева
закончена. В дальнейшем всюду отождествляем задачи (14), (2) и (3), (4).
3. Расширенное фазовое пространство
и квазистационарные полутраектории
Теперь перейдем к проверке выполнения условий теорем 1 и 2.
В силу теоремы 6 и соответствующих результатов [12] существует аналитическая полугруппа {U t : t ∈ R+ } разрешающих операторов уравнения (4), которую в данном случае
естественно представить в виде U t = V t ⊗W t , где V t (W t ) — сужение оператора U t на U1 (U2 ).
Поскольку оператор B2 секториален, то W t = exp(tB2 ), что влечет за собой W ◦ = {0} и
W 1 = U2 .
Рассмотрим полугруппу {V t : t ∈ R+ }. В силу теорем 4 и 6 и цитируемой монографии
[12] данная полугруппа продолжима до группы {V t : t ∈ R}. Ее ядро V 0 = U100 ⊕ U101 , где
−1
2
2
U100 = {0} × {0} × Hp × {0} × . . . × {0}(= ker A1 по теореме 5), а U101 = ΣA−1
λ Aλπ [Hπ ] × Hπ ×
−1
{0} × . . . × {0} . Здесь Aλ = I − λA, Aλπ — сужение оператора ΠAλ на Hπ . Известно , что
{z
}
|
k+1
если λ−1 6∈ σ(A) ∪ σ(Aσ ), то оператор Aλπ : Hπ → H2π — топлинейный изоморфизм (см.,
например, [8]). Обозначим через U11 образ V 1 . Тогда пространство U1 разлагается в прямую
сумму подпространств: U1 = U100 ⊕ U101 ⊕ U11 .
−1
Построим оператор R = B10
A10 ∈ L(U100 ⊕ U101 ), где A10 (B10 ) — сужение оператора
−1
00
01
A1 (B1 ) на U1 ⊕ U1 . (Оператор B10
существует в силу теоремы 6 и соответствующих результатов [12]). По построению ker R = U100 , а в [22] показано, что im R = U100 . Значит,
оператор R — бирасщепляющий. Обозначим через PR проектор пространства U100 ⊕ U101
0 ), где
на U100 вдоль U101 . В силу конструкции пространства UM проектор PR ∈ L(UM
0 = U ∩ (U 00 ⊕ U 01 )(≡ U 00 ⊕ U 01 ). Зафиксируем это в следующем утверждении.
UM
M
1
1
1
1
Лемма 1. Пусть λ−1 6∈ σ(A) ∪ σ(Aσ ). Тогда оператор R — бирасщепляющий, причем
0 ).
PR ∈ L(UM
Введем в рассмотрение проекторы
Pk = diag [Pˆk , 0],
Qk = diag [Q̂k , 0],
k = 0, 1.
(Подробное описание этих проекторов см. в [23]. Из результатов [23] и в силу того, что ядро
W 0 = {0}, следует, что I − P = (P0 + P1 ) × O, Q = (I − Q0 − Q1 ) × I, P : U → U 1 , Q :
F → F 1 . Применяя проектор I − P к уравнению (4) в данной транскрипции, получаем
Π(νA(uσ + uπ ) − ((ũσ + ũπ ) · ∇)(uσ + uπ ) − ((uσ + uπ ) · ∇)(ũσ + ũπ )+
k
X
βl ∇2 wl − up − gγuθ + f (t)) = 0,
(15)
l=1
Buπ = 0.
Отсюда в силу теоремы 1 и свойств оператора B получаем необходимое условие квазистационарности полутраектории uπ ≡ 0. Другими словами, все решения нашей задачи (если они
существуют) с необходимостью должны лежать в плоскости B = {u ∈ UM : uπ = 0}. А так
как Πup = up , то из первого уравнения (15) получаем соотношение (7) в нашей транскрипции
up = Π(νAuσ − ((ũσ + ũπ ) · ∇)uσ − (uσ · ∇)(ũσ + ũπ ) +
48
k
X
l=1
βl ∇2 wl − gγuθ + f (t)).
(16)
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
Очевидно, P0 ≡ PR , поэтому второе уравнение (15) есть соотношение (8) применительно к
нашей ситуации. Итак, справедлива
Лемма 2. В условиях леммы 1 любое решение задачи (3), (4) лежит во множестве
At = {u ∈ UM : uπ = 0, up = Π(νAuσ − ((ũσ + ũπ ) · ∇)uσ − (uσ · ∇)(ũσ + ũπ )+
k
X
l=1
βl ∇2 wl − gγuθ ) + fπ (t)}.
00 (≡
Замечание 12. Из (16) сразу следует условие (iii) теоремы 2 для любой точки u00 ∈ UM
U100 × {0}). Поэтому ввиду замечания 8 множество At — простое банахово многообразие C ∞ диффеоморфное подпространству U11 × U2 — является кандидатом на роль расширенного
фазового пространства задачи (14), (2).
◦
Приступим к проверке условий (9) и (13). Построим пространство Uα = U1 × W21 (Ω). Данное пространство, очевидно, будет интерполяционным пространством для пары [U, UM ]α ,
причем α = 1/2. Как отмечено выше, полугруппа {U t : t ∈ R+ } продолжается до группы
1 = U
1
{V1t : t ∈ R} на U11 , где V1t — сужение оператора V t на U11 . Поскольку UM
M ∩ U1 (по
построению) и оператор B1 непрерывен (теорема 4 ), то в силу равномерной ограниченности
полугруппы {U t : t ∈ R+ } имеем
Z τ
kV1t kL(U 1 ; U 1 ) dt ≤
1
0
constkB1 kL(U1 ; F1 )
Z
τ
0
M
kV1t kL(U 1 ) dt < ∞
1
∀τ ∈ R+ .
(17)
−
Далее, в силу неравенства Соболева [19, гл.9] полугруппа {W t : t ∈ R+ } удовлетворяет
оценке
Z τ
kW t k
dt < ∞.
(18)
◦
0
L(dom B2 ; W21 (Ω))
Положим Uα1 = Uα ∩ U 1 , где U 1 = U11 × U2 . Тогда из (17) и (18) вытекает
Лемма 3. В условиях леммы 1 выполняется соотношение (9).
И наконец, выполняя требование (13), найдем оператор H и вектор-функцию h. Оператор H естественно представить в виде H = H1 × H2 , где H1 = A−1
11 (I − Q0 − Q1 )F1 , а
1 ; U 1 ), показываетH2 ≡ F2 (A11 — сужение оператора A1 на U11 ). Включение H ∈ C ∞ (UM
α
∞
ся аналогично тому, как было показано включение F ∈ C (UM ; F). Вектор-функцию h(t)
определим как h = A−1
11 (I − Q0 − Q1 )f. Определение операторов Q0 и Q1 см. в [23]. В силу
бесконечной гладкости f h ∈ C ∞ (R̄+ ; Uα1 ).
Таким образом, все условия теоремы 2 выполнены. Поэтому справедлива
Теорема 8. Пусть λ−1 6∈ σ(A) ∪ σ(Aσ ). Тогда при любом u0 ∈ A0 и некотором T ∈
R+ существует единственное решение u(t) ∈ C ∞ ( (0, T ); UM ) задачи (1), (2), являющееся
квазистационарной полутраекторией.
Работа выполнена при поддержке АВЦП Развитие научного потенциала высшей школы (2009 – 2011 годы), проект № 2.1.1/2301.
Автор выражает искреннюю признательность профессору Г.А.Свиридюку за внимание
и интерес к данным исследованиям.
Серия
≪
Математическое моделирование и программирование≫, вып. 10
49
Т.Г. Сукачева
Литература
1. Осколков, А.П. Начально-краевые задачи для уравнений движения жидкостей
Кельвина-Фойгта и жидкостей Олдройта / А.П. Осколков // Труды мат. ин-та АН
СССР. – 1988. – №179. – С. 126 – 164.
2. Осколков, А.П. Нелокальные проблемы для одного класса нелинейных операторных
уравнений, возникающих в теории уравнений типа С.Л.Соболева / А.П. Осколков //
Зап. науч. семин. ЛОМИ. – 1991.– Т. 198.– С. 31 – 48.
3. Свиридюк, Г.А. К общей теории полугрупп операторов / Г.А. Свиридюк // Успехи мат.
наук. – 1994. – Т. 49, №4. – С. 47 – 74.
4. Осколков, А.П. О некоторых нестационарных линейных и квазилинейных системах,
встречающихся при изучении движения вязких жидкостей / А.П. Осколков // Зап.
науч. семин. ЛОМИ АН СССР. – 1976.– Т. 59. – С. 133 – 177.
5. Осколков, А.П. К теории жидкостей Фойгта / А.П. Осколков // Зап. научн. сем. ЛОМИ.
– 1980. – Т. 96. – С. 233 – 236.
6. Свиридюк, Г.А. Разрешимость задачи термоконвекции вязкоупругой несжимаемой жидкости / Г.А. Свиридюк // Изв. вузов. Матем. – 1990. – №12. – С. 65 — 70.
7. Свиридюк, Г.А. Фазовые пространства полулинейных уравнений типа Соболева с относительно сильно секториальным оператором / Г.А. Свиридюк // Алгебра и анализ. –
1994. – Т. 6, №5. – С. 216 – 237.
8. Сукачева, Т.Г. Исследование математических моделей несжимаемых вязкоупругих жидкостей: дис. ... д-ра физ.-мат. наук / Т.Г. Сукачева; Новгород. гос. ун-т. – Великий
Новгород, 2004. – 249 с.
9. Сукачева, Т.Г. Нестационарная линеаризованная модель движения несжимаемой вязкоупругой жидкости / Т.Г. Сукачева // Вестн. Челяб. гос. ун-та. Сер. Математика.
Механика. Информатика. – 2009. – №20 (158), вып. 11. – С. 77 – 83.
10. Сукачева, Т.Г. Нестационарная линеаризованная модель движения несжимаемой вязкоупругой жидкости высокого порядка / Т.Г. Сукачева // Вестн. Юж.-Урал. гос. ун-та,
серия ≪Математическое моделирование и программирование≫. – 2009. – № 17 (150),
вып. 3. – С. 86 – 93.
11. Сукачева, Т.Г. Задача термоконвекции для линеаризованной модели движения несжимаемой вязкоупругой жидкости / Т.Г. Сукачева // Вестн. Юж.-Урал. гос. ун-та, серия
≪Математическое моделирование и программирование≫. – 2010. – № 16 (192), вып. 5. –
С. 83 — 93.
12. Sviridyuk, G.A. Linear Sobolev type equations and degenerate semigroups of operators /
G.A. Sviridyuk, V.E. Fedorov. – Utrecht; Boston; Köln; Tokyo: VSP, 2003.
13. Свиридюк, Г.А. Квазистационарные траектории полулинейных динамических уравнений типа Соболева / Г.А. Свиридюк // Изв. РАН. Сер. Математика. – 1993. – Т. 57, №3.
– С. 192 – 207.
14. Levine, H.A. Some nonexistance and instability theorems for solutions of formally parabolic
equations of the form Dut = −Au + F (u) / H.A. Levine // Arch. Rat. Mech. Anal. – 1973.
– V. 51, № 5. – P. 371 – 386.
15. Свиридюк, Г.А. Задача Коши для одного класса полулинейных уравнений типа Соболева / Г.А. Свиридюк, Т.Г. Сукачева // Сиб. мат. журн. – 1990. – Т. 31, №5. – С. 109 – 119.
50
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
16. Свиридюк, Г.А. Фазовые пространства одного класса операторных уравнений /
Г.А. Свиридюк, Т.Г. Сукачева // Дифференц. уравнения. – 1990. – Т. 26, №2. –
С. 250 – 258.
17. Свиридюк, Г.А. Некоторые математические задачи динамики вязкоупругих несжимаемых сред / Г.А. Свиридюк, Т.Г. Сукачева //Вестник МаГУ. Математика. – 2005. –
Вып. 8. – С. 5 – 33.
18. Борисович, Ю.Г. Нелинейные фредгольмовы отображения и теория Лере-Шаудера /
Ю.Г. Борисович, В.Г. Звягин, Ю.И. Сапронов // Успехи матем. наук. – 1977. – Т. 32,
№4. – С. 3 – 54.
19. Марсден, Дж. Бифуркация рождения цикла и ее приложения / Дж. Марсден , М. МакКракен. – М.: Мир, 1980. – 368 с.
20. Бокарева, Т.А. Исследование фазовых пространств уравнений типа Соболева с относительно секториальными операторами: дис. ... канд. физ.-мат. наук / Т.А. Бокарева. –
Санкт-Петербург, 1993. – 107 с.
21. Ладыженская, О.А. Математические вопросы динамики вязкой несжимаемой жидкости
/ О.А. Ладыженская. – Изд. 2. – М.: Наука, 1970. – 288 с.
22. Свиридюк, Г.А. Об одной модели слабосжимаемой вязкоупругой жидкости / Г.А. Свиридюк // Изв. вузов. Математика. – 1994. – № 1. – С. 62 – 70.
23. Сукачева, Т.Г. Об одной модели движения несжимаемой вязкоупругой жидкости
Кельвина-Фойгта ненулевого порядка / Т.Г. Сукачева // Дифференц. уравн. – 1997.
– Т. 33, №4. – С. 552 – 557.
24. Свиридюк Г.А. Полулинейные уравнения типа Соболева с относительно ограниченным
оператором / Г.А. Свиридюк // ДАН СССР. – 1991. – Т.318, № 4. – С. 828 – 831.
25. Свиридюк, Г.А. Полулинейные уравнения типа Соболева с относительно секториальными операторами / Г.А. Свиридюк // Докл. РАН. – 1993. – Т. 329, №3. – С. 274 – 277.
26. Свиридюк, Г.А. Аналитические полугруппы с ядрами и линейные уравнения типа Соболева / Г.А. Свиридюк, В.Е. Федоров // Сиб. мат. журн. – 1995. – Т. 36, №5. –
С. 1130 – 1145.
27. Хенри, Д. Геометрическая теория полулинейных параболических уравнений / Д. Хенри. — М.: Мир, 1985. – 376 c.
References
1. Oskolkov A.P. Initial-value problems for equations of motion Kelvin-Voight and Oldroyd
fluids [Nachal’no-kraevye zadachi dlya uravneniy dvizheniya zhidkostey Kel’vina-Foygta i
zhidkostey Oldroyta] Trudy mat. in-ta AN SSSR, 1988, no. 179, pp. 126 – 164.
2. Oskolkov A.P. Nonlocal problems for a class of nonlinear operator equations arising in the
theory of Sobolev type equations [Nelokal’nye problemy dlya odnogo klassa nelineynykh
operatornykh uravneniy, voznikayushchikh v teorii uravneniy tipa S.L.Soboleva] Zap. nauch.
semin. LOMI, 1991, vol. 198, pp. 31 – 48.
3. Sviridyuk G.A. On the general operator semigroups theory [K obshchey teorii polugrupp
operatorov] Uspekhi mat. nauk., 1994, vol. 49, no. 4, pp. 47 – 74.
4. Oskolkov A.P. Some nonstationary linear and quasilinear systems occurring in the study
of movement viscous fluids [O nekotorykh nestatsionarnykh lineynykh i kvazilineynykh
sistemakh, vstrechayushchikhsya pri izuchenii dvizheniya vyazkikh zhidkostey] Zap. nauch.
semin. LOMI AN SSSR, 1976, vol. 59, pp. 133 – 177.
Серия
≪
Математическое моделирование и программирование≫, вып. 10
51
Т.Г. Сукачева
5. Oskolkov A.P. On the theory of Voigt liquids [K teorii zhidkostey Foygta] Zap. nauchn. sem.
LOMI, 1980, vol. 96, pp. 233 – 236.
6. Sviridyuk G.A. Solubility of the thermal convection of viscoelastic incompressible fluid
[Razreshimost’ zadachi termokonvektsii vyazkouprugoy neszhimaemoy zhidkosti] Izv. vuzov.
Matem., 1990, no. 12, pp. 65 — 70.
7. Sviridyuk G.A. Phase spaces of semilinear Sobolev type equations with relatively strong
sectorial operator [Fazovye prostranstva polulineynykh uravneniy tipa Soboleva s otnositel’no
sil’no sektorial’nym operatorom] Algebra i analiz., 1994, vol. 6, no. 5, pp. 216 – 237.
8. Sukacheva T.G. The study of mathematical models of incompressible viscoelastic fluids: dis.
... Dr. Sci. Science [Issledovanie matematicheskikh modeley neszhimaemykh vyazkouprugikh
zhidkostey: dis. ... d-ra fiz.-mat. nauk]. Velikiy Novgorod, 2004. 249 p.
9. Sukacheva T.G. Unsteady linearized model of the motion of an incompressible viscoelastic
fluid [Nestatsionarnaya linearizovannaya model’ dvizheniya neszhimaemoy vyazkouprugoy
zhidkosti] Vestn. Chelyab. gos. un-ta. Ser. Matematika. Mekhanika. Informatika, Vyp. 11,
2009, no. 20 (158), pp. 77 – 83.
10. Sukacheva T.G. Unsteady linearized model of the motion of an incompressible viscoelastic
fluid of the high order [Nestatsionarnaya linearizovannaya model’ dvizheniya neszhimaemoy
vyazkouprugoy zhidkosti vysokogo poryadka] Vestnik Yuzhno-Ural’skogo gosudarstvennogo
universiteta. Seriya ≪Matematicheskoe modelirovanie i programmirovanie≫, 2009, no. 17
(150), vyp. 3, pp. 86 – 93.
11. Sukacheva T.G. The problem of thermal convection for a linearized model of
the motion of an incompressible viscoelastic fluid [Zadacha termokonvektsii dlya
linearizovannoy modeli dvizheniya neszhimaemoy vyazkouprugoy zhidkosti] Vestnik
Yuzhno-Ural’skogo gosudarstvennogo universiteta. Seriya ≪Matematicheskoe modelirovanie
i programmirovanie≫, 2010, no. 16 (192), vyp. 5, pp. 83 — 93.
12. Sviridyuk G.A., Fedorov V.E. Sobolev type equations and degenerate semigroups of operators,
Utrecht-Boston: VSP, 2003. 179 p.
13. Sviridyuk G.A. Quasi-stationary trajectories of semilinear dynamical Sobolev type equations
[Kvazistatsionarnye traektorii polulineynykh dinamicheskikh uravneniy tipa Soboleva] Izv.
RAN. Ser. Matematika, 1993, vol. 57, no. 3, pp. 192 – 207.
14. Levine H.A. Some nonexistance and instability theorems for solutions of formally parabolic
equations of the form Dut = −Au + F (u) [Some nonexistance and instability theorems for
solutions of formally parabolic equations of the form Dut = −Au + F (u)] Arch. Rat. Mech.
Anal., 1973, vol. 51, no. 5, pp. 371 – 386.
15. Sviridyuk G.A., Sukacheva T.G. Cauchy problem for a class of semilinear Sobolev type
equations [Zadacha Koshi dlya odnogo klassa polulineynykh uravneniy tipa Soboleva] Sib.
mat. zhurn., 1990, vol. 31, no. 5, pp. 109 – 119.
16. Sviridyuk G.A., Sukacheva T.G. Phase spaces of class of operator equations [Fazovye
prostranstva odnogo klassa operatornykh uravneniy] Differents. uravneniya, 1990, vol. 26,
no. 2, pp. 250 – 258.
17. Sviridyuk G.A., SukachevaT.G. Some mathematical problems of the dynamics of viscoelastic
incompressible media [Nekotorye matematicheskie zadachi dinamiki vyazkouprugikh
neszhimaemykh sred], Vestnik MaGU. Matematika, 2005, Vyp. 8, pp. 5 – 33.
18. Borisovich Yu.G., Zvyagin V.G., Sapronov Yu.I. Nonlinear Fredholm maps and LeraySchauder theory [Nelineynye fredgol’movy otobrazheniya i teoriya Lere-Shaudera] Uspekhi
matem. nauk., 1977, vol. 32, no. 4, pp. 3 – 54.
52
Вестник ЮУрГУ, №37 (254), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
19. Marsden Dzh., Mak-Kraken M. Hopf bifurcation and its applications [Bifurkatsiya rozhdeniya
tsikla i ee prilozheniya], Moscow: Mir, 1980. 368 p.
20. Bokareva T.A. Investigation of phase space of Sobolev type equations with relatively sectorial
operators: Dis. ... cand. Sci. Science [Issledovanie fazovykh prostranstv uravneniy tipa
Soboleva s otnositel’no sektorial’nymi operatorami: dis. ... kand. fiz.-mat. nauk], SanktPeterburg, 1993. 107 p.
21. Ladyzhenskaya O.A. The mathematical theory of dinamic of viscous incompressible fluid
[Matematicheskie voprosy dinamiki vyazkoy neszhimaemoy zhidkosti, izd. 2.], Moscow:
Nauka, 1970. 288 p.
22. Sviridyuk G.A. A model of weakly viscoelastic fluid [Ob odnoy modeli slaboszhimaemoy
vyazkouprugoy zhidkosti] Izv. vuzov. Matematika, 1994, no. 1, pp. 62 – 70.
23. Sukacheva T.G. A model of motion of an incompressible viscoelastic Kelvin-Voigt fluid of
nonzero order [Ob odnoy modeli dvizheniya neszhimaemoy vyazkouprugoy zhidkosti Kel’vinaFoygta nenulevogo poryadka] Differents. uravn., 1997, vol. 33, no. 4, pp. 552 – 557.
24. Sviridyuk G.A. Semilinear Sobolev type equation with relatively bounded operator
[Polulineynye uravneniya tipa Soboleva s otnositel’no ogranichennym operatorom] DAN
SSSR, 1991, vol. 318, no. 4, pp. 828 – 831.
25. Sviridyuk G.A. Semilinear Sobolev type equation with relatively sectorial operators
[Polulineynye uravneniya tipa Soboleva s otnositel’no sektorial’nymi operatorami] Dokl. RAN,
1993, vol. 329, no. 3, pp. 274 – 277.
26. Sviridyuk G.A., Fedorov V.E. Analytic semigroup with kernels and linear Sobolev type
equations [Analiticheskie polugruppy s yadrami i lineynye uravneniya tipa Soboleva] Sib.
mat. zhurn., 1995, vol. 36, no. 5, pp. 1130 – 1145.
27. Khenri D. Geometric theory of semilinear parabolic equations [Geometricheskaya teoriya
polulineynykh parabolicheskikh uravneniy], Moscow.: Mir, 1985. 376 p.
Тамара Геннадьевна Сукачева, доктор физико-математических наук, доцент, кафедра
алгебры и геометрии, Новгородский государственный университет им. Ярослава Мудрого
(Россия, Великий Новгород), tamara.sukacheva@novsu.ru.
Tamara G. Sukacheva, Doctor of Physical and Mathematical Sciences, Associate Professor,
Department of Algebra and Geometry, Novgorod State University (Russia, Velikiy Novgorod),
tamara.sukacheva@novsu.ru.
Поступила в редакцию 16 июня 2011 г.
Серия
≪
Математическое моделирование и программирование≫, вып. 10
53
Документ
Категория
Без категории
Просмотров
3
Размер файла
300 Кб
Теги
линеаризованной, ненулевом, термоконвекции, задачи, модель, порядке, жидкости, вязкоупругих, несжимаемой
1/--страниц
Пожаловаться на содержимое документа