О существовании переходов между стационарными режимами задачи обтекания

Автор: Сазонов Леонид Иванович

Журнал: Владикавказский математический журнал @vmj-ru

Статья в выпуске: 4 т.13, 2011 года.

Бесплатный доступ

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

Система навье -- стокса, оператор озеена, возмущенная полугруппа озеена, стартовая проблема.

Короткий адрес: https://sciup.org/14318366

IDR: 14318366

Текст научной статьи О существовании переходов между стационарными режимами задачи обтекания

Пусть ограниченное тело B С R 3 движется в жидкости, занимающей все пространство вне тела, со скоростью a(t)e i , где e i — единичный орт оси Ox i . Тогда поле скорости жидкости u = u(x, t) удовлетворяет следующей начально-краевой задаче для системы Навье — Стокса

( д? + (u, V )u = A u — V p,

[divu = 0, uldB(t) = —a(t)ei, u^ = 0, где B(t) — область, занимаемая телом в момент времени t. Здесь, не нарушая общности, считаем, что плотность и коэффициент вязкости равны единице. Полагая

u(x, t) = v

t x + j a(s)e1

a(t)e i ,

p(x, t) = r I x

+

t

У a(s)e i

получаем, что v(x,t), r(x,t) удовлетворяют системе

I

^dv + (v, V )v a'(t)e i = Av Vr, divv = 0, v | S Q = 0, v | ^ = a(t)e i ,

в области Q x (0, го ), где Q = R 3 \ B (0).

Пусть v α 1 , v α 2 — стационарные режимы задачи обтекания, т. е. решения стационарной системы Навье — Стокса в области Q

A v Vp = (v, V )v, divv = 0, v | q q = 0,

удовлетворяющие соответственно условиям vU = aiei, v|^ = a2ei.

Задачей о переходе между стационарными режимами vα1 и vα2 будем называть задачу об определении решения нестационарной системы Навье — Стокса (2), удовлетворяющего начальному условию v|t=o = va1 и предельному соотношению limt^^ v = va2. В случае ai = 0 получается задача, впервые рассмотренная Р. Финном [1] и названная им «стартовой проблемой». О работах, посвященных этой задаче и имеющихся здесь трудностях смотрите в [2]. В частности, в [2] при ai = 0 и достаточно малом a2 установлено существование в подходящем банаховом пространстве решения стартовой проблемы v(t), удовлетворяющего оценке llV(v(t) — va2 )kL3 6 ct-i/2, kv(t) - va2 IlLq 6 cq t-(i-3/q^/2, q> 3.

2.    Сведение к интегральному уравнению

Будем предполагать, что стационарная система (3) имеет семейство стационарных решений v a (v a | ^ = ae i ), гладко зависящих от параметра a при a Е [a i , a 2 ]. (В дальнейшем это условие будет уточнено.) При этом предположении решение задачи (2) будем искать в виде v = v a(t) + w, r = p a ( t ) + q, где a(t) — гладкая функция со значениями на отрезке [a i ,a 2 ], такая, что a(0) = a i , a(t) ^ a 2 при t ^ го . Тогда для определения w имеем следующую начально-краевую задачу в области Q при t >  0

f d dw + (v a ( t ) , V )w + (w, V)v a(t) + (w, V )w - a 0 (t)e i = 4 w - V q - a 0 (t) da v a(t) ,

(divw = 0, w|SQ = 0, w|t=o = 0, причем будем требовать выполнения условия w|^ = 0 при всех t > 0.

Введем обозначение va = va — aei. Тогда первое уравнение системы (4) представляется в виде dW + a2diW + (Va2, V)w + (w, V)Va2 + (a(t) — «2) diW + (va(t) — Va2, V)w

-

d                     (5)

+ (w, V )(v a ( t ) v a 2 ) + (w, V )w = 4 w — Vq a' (t) dav a ( t )

Для дальнейшего исследования удобно свести систему (4) к задаче Коши для ОДУ в подходящем банаховом пространстве. Обозначим через S p = S p (^) (1 6 p 6 го ) подпространство в пространстве векторных полей L p (Q), являющееся замыканием множества всех гладких соленоидальных полей с компактным в Q носителем. (Заметим, что здесь и далее мы не различаем обозначения для пространств векторных полей и функций.) Известно (см., например, [3, 4]), что для областей класса C 2 (в дальнейшем это условие предполагается всегда выполненным) существует ограниченный проектор П : L p (Q) ^ S p (^) (1 < p <  го ). Применяя его к системе (4) и учитывая (5), сведем ее к следующей задаче Коши для ОДУ в пространстве S p (Q)

^t. = Aw + B(t)w + Kw + F (t), w|t=o = 0, где операторы A, B(t), K и векторное поле F(t) имеют вид

Aw = П I 4 w a 2 d i w (v a 2 , V )w (w, V )v a 2

B (t)w = - П^(«(<) - a 2 )d i w + (V a(t) - V a 2 , V )w + (w, V )(v a ( t ) - V a 2 )J ,       (8)

Kw = - n(w, V )w, F (t) = - П a (t) dav a(t )j •                  (9)

Рассмотрим действующий в пространстве S p (Q) оператор Озеена

A o = П( 4 - u ^ d i )

с областью определения D(A o ) = S p (Q) П W p2 (Q) П W p (Q).

Здесь Wp(Q), Wp(Q) — соболевские пространства векторных полей, причем элементы из Wp(Q) имеют нулевой след на границе dQ.

Оператор A o порождает в любом пространстве S p (Q) (1 < p <  го ) аналитическую полугруппу T ^ (t), для которой справедливы оценки

/, - Ид' _ 6 Cpq t-|a|/2-3/2(1/p-1/q) , где 1 < p 6 q < го при a = 0; 3/2 6 p 6 q < го при |a| = 1.

При фиксированном значении u ^ = 0 данные оценки установлены в [5] методом гидродинамических потенциалов, в [6] они получены независимо и иными методами, причем доказана их равномерность по параметру u ^ Е (0, r).

Вместе с оператором A o рассмотрим возмущенный оператор Озеена A = A + B , где Bu = - П((u, V )v + (v, V )u).

В [7] для оператора A установлены следующие результаты.

Теорема 1. Если v — соленоидальное поле, причем v,djv Е L^(Q), то возмущенный оператор Озеена A с областью определения D(A) = D(Ao) порождает в любом пространстве Sp(Q), 1 < p < го, аналитическую полугруппу T(t). Если дополнительно v Е Lp(Q), p Е [pi, го], где pi < 3, и при некотором q > 2 возмущенный оператор Озеена A : Sq(Q) ^ Sq(Q) не имеет собственных значений в полуплоскости Re А > 0, то для возмущенной полугруппы Озеена справедлива оценка kT(1)Пд6 ||p_q 6 Cpqt-\6V2-3/2(1/p-1/q>                          (10)

при выполнении следующих условий

|0| 6 1, 1 < p 6 q < го, q > 3/2, причем p > 3/2 при |У| = 1.

Заметим, что из результатов [8] (см. также [9]) следует, что обобщенное решение v a стационарной задачи обтекания (3) удовлетворяет условиям: v a Е W p 2 (Q) для всех p Е (2, го ], V v a Е L p (Q) для всех p Е (3/2, го ]. Далее, как уже отмечалось, мы предполагаем гладкую зависимость решения v a при a Е [a i ,a 2 ] (a i ,a 2 > 0) в норме любого пространства L p (Q) при p Е (2, го ]. Заметим, что локально, т. е. в окрестности фиксированного α данное предположение может быть обосновано методами теории неявных отображений, если для v α выполнены предположения теоремы 1. При малых числах Рейнольдса (в данном случае при малых α) указанное утверждение можно извлечь из результатов [10]. Сошлемся также на работу [11], где исследован вопрос о гладкой зависимости стационарного решения задачи обтекания от числа Рейнольдса в норме некоторого банахова пространства.

В дальнейшем оператор A вида (7) будем рассматривать как возмущение оператора Озеена с той же областью определения D(A) = D(Ao ) и предполагать выполненными предположения теоремы 1. Аналитическую полугруппу, порожденную этим оператором, будем обозначать T(t). В силу наших предположений для нее справедливы оценки (10) и кроме того даva Е Lp(Q) для любого p £ (2, го].

Формально обращая главную часть ОДУ (6), приходим к интегральному уравнению

w(t) = j | t (t s) ( B(s)w(s) + Kw(s) )| ds + F(t),               (11)

где

t

F (t) = 1

T(t s)F(s) ds

а операторы B (t), K и поле F (t) определены в (8), (9).

В дальнейшем решение задачи Коши (6) будем понимать в обобщенном смысле как решение интегрального уравнения (11).

3.    Существование решения задачи о переходе

Предварительно установим ряд оценок для интегральных операторов и поля F(t), входящих в уравнение (11). Всюду в дальнейшем считаем выполненным условие a(t) = a 2 при t > R.

Лемма 1. Справедлива оценка

(1+ t)5kF(t)kq 6 A(p,q,R,5K где

A(p, q, R, 5) = sup

α

∂ vα ∂α

(1 - y)- 1 (1 + R) 5 R 1-Y C p,q p

c p,q — константы из оценки (10), E = max | a 0 (t) | , а параметры 5, y , p, q удовлетворяют соотношениям 0 <5 6 y = (3/2)(1/p 1/q) < 1, 2 < p 6 q < го.

C Используя представления (9), (12) и оценки полугруппы (10), устанавливаем нера- венство

kF (t) k q 6 c pq E suP

α

∂ vα ∂α

p

min( t,R )

У (t s) -Y ds.

Остается получить оценку для интеграла

I (t) = (1 + t) 5

min( t,R ) /  (t

s) Y ds.

В случае t 6 R, осуществляя замену 1 — (1 + s)(1 + t) 1 = u, приводим интеграл I(t) к виду t/(i+t)

I (t) = (1+ t)5+1-Y У u-Y du, из которого ввиду возрастания I(t) на [0, R] получаем

R/ (1+ R )

I (t) 6 (1 + r) 5+1 - y   у .

«7 ... = R 'R

1 Y

.

В случае t > R интеграл I (t) убывает при возрастании t. Поэтому

R

I (t) 6 (1 + R) 5 J (R - s) -Y ds =

R 1-Y (1 + R) 5

1 Y

.

Утверждение леммы следует из установленных оценок (13)–(15). B

Для оценок интегральных операторов введем банахово пространство Xδ,q векторных полей из C([0, го), Sq(Q)) с конечной нормой llu|5,q = SUp(1 + t)5 ku(t)kq.

t

Лемма 2. При 0 <5 6 1/2, q >  3/2 для оператора

t

Bu(t) = У T (t s)B(s)u(s) ds

справедлива оценка

WUkq 6 24C q,q dR 1/2 k u k 5,q , где c q,q — константа из оценки (10),

d = max d(t), d(t) = | a(t) a 2 | + 2 k V a ( t ) V a 2 k ^ -

C Используя оценки (10) с | 0 | = 1, получаем неравенство

t kBu(t)kq 6 3Cq,q У(t — s)-1/2d(s)(1+ S^-ds^uk^.q,   q> 3/2.

Учитывая, что функция d(t) отлична от нуля лишь на отрезке [0, R], приходим к выводу, что достаточно оценить интеграл

min( t,R )

J (t) = (1 + t)5 j (t s) 1/2 (1 + s) 0

δ ds.

Применяя рассуждения, аналогичные доказательству леммы 1, получаем оценку (1 + t) 5 J (t) 6 8R 1/2 .

С учетом (16) приходим к утверждению леммы. B Лемма 3. Для оператора

t

K u(t) = j T(t s)Ku(s) ds, 0

где K определен в (9), при выполнении условий p > 3, q > 3/2, 0 6 2/p - 1/q < 1/3, 0 < в < 1/2, 0 6 5 6 2в + (3/2)(2/p - 1/q) - 1/2

справедлива оценка kKukq 6 C(e,5,p,q)kuke,p, где

C(e,5,p,q) = 3cp/2,qB(1 - a, 1 - 2в), cp/2q — константа из оценки (10), a = (3/2)(2/p - 1/q) + 1/2, B(•, •) — бета-функция Эйлера, определяемая формулой

B(а,в) = j (t - s) a-1 s e-1

ds.

С Вследствие оценок полугруппы T (t) справедливо неравенство

t kK u(t)kq 6 3cp/2,q кикв„ J (t - s)-a(1 + s)"^ ds.                 (17)

Теперь, осуществляя в интеграле

t

I(t) = j(t - s)-a(1 + s)-2e ds замену переменной 1 - (1 + s)/(1 + t) = a, получаем t/(1+t)

I (t) = (1 + t) 1-a-2e j a -a (1 - a) -2e da 6 (1 + t) 1-a- 2 e B (1 - a, 1 - 2в).     (18)

Комбинируя (17) и (18), приходим к заключению леммы. B

Следствие 1. При выполнении условий q > 3 и 1/2 - 3/(2q) 6 5 < 1/2 справедливы неравенства kKukq 6 C(5, 5,9,q)kuk2,q,

||K U 1 - KU 2 k »,q 6 2C(5, 5, q, q) max ku j Ц ^,, k u i - ^k^q .

j

Для исследования вопроса о разрешимости уравнения (11) воспользуемся следующей леммой.

Лемма 4. Пусть X — банахово пространство, B : X ^ X — линейный ограниченный оператор и kBk 6 1/2, Q : X ^ X — нелинейный оператор с оценками kQuk 6 b||u||2,   kQui - Qu2k 6 bmax kukkui - «2k- j

Тогда уравнение u = Bu + Qu + g при выполнении условия kgk 6 1/(16b) имеет единственное решение в шаре {u : |u| 6 1/(4b)}.

C При указанных предположениях шар { u : ||u k 6 1/(4b) } инвариантен относительно отображения u i^ B u + Q u + g, являющегося сжимающим в этом шаре. Поэтому утверждение леммы следует из теоремы Банаха о неподвижной точке. B

Вернемся к уравнению (11). Справедлива следующая

Теорема 2. Пусть существует гладкое семейство стационарных решений v α , α [a i 2 ] задачи обтекания, причем v a 2 таково, что при некотором q >  2 возмущенный оператор Озеена A : S q (Q) ^ S q (Q) вида (7) не имеет собственных значений в полуплоскости Re А > 0. Тогда для любых фиксированных чисел q, p, 6, удовлетворяющих условиям q >  3, p >  2, 1/2 (3/2q) 6 6 < min(1/2, (3/2)(1/p 1/q)), существует такое П = n(q,P, 6), что при | a i a 2 | 6 n уравнение (11) имеет единственное решение w Е X s,q

C Заметим, что уравнение (11) можно представить в виде w = Bw + Kw + F (t), где операторы B, K и поле F определены соответственно в леммах 2 и 3 и в (12). Будем считать, что a(t) = ai + («2 — ai)t при t Е [0, R] и a(t) = «2 при t > R. Тогда при заданных q, p, δ в силу лемм 1, 2 и следствия из леммы 3 выполняются неравенства kF kδ,q 6 A|α1 - α2|, kBukδ,q 6 B|α1 - α2|kukδ,q, kKuks,q 6 Ckuk2q, |Kui — Ku2lkq 6 2Cmax ku ks,q|ui — u^q. j

Вследствие установленных оценок существование решения вытекает из леммы 4.

0 6t6T

Тогда аналогично предыдущим рассуждениям устанавливаем неравенство llwkq,T 6 c(T)kwk^,q,T, в котором c(T) ^ 0 при T ^ 0. Поэтому w(t) = 0 в окрестности точки t = 0.

Пусть £ = max T всех таких T , что w(t) = 0 при t 6 T.

Очевидно, что £ = го , иначе аналогично предыдущему

NW 6 d(T — e)kwk^,q,T при T > (, где d(T — () ^ 0 при T ^ £. Но тогда w(t) = 0 на более широком отрезке чем [0, ^], что противоречит определению ^. Этим завершается доказательство теоремы. B

Замечание. Можно считать a(t) произвольной гладкой функцией и дополнительно требовать малости sup | a 0 (t) | при t Е [0, R].

4.    Оценки сходимости к стационарному режиму

Дальнейшей целью является получение оценок норм разности v(t) — va2, где v(t) — решение системы (2), при t ^ го. Учитывая, что при t > R v(t) = va2 + w(t), приходим к необходимости оценки норм w(t). Для этого будет исследован вопрос о принадлежности решения w уравнения (11) пространствам Xµ,r. В этом случае мы не предполагаем малости решения и, следовательно, малости разности |ai — «2/ а используем лишь его принадлежность одному из пространств Xg,q, в котором существует согласно теореме 2 малое решение.

Ниже существенно используется следующая лемма.

Лемма 5. Оператор I B, где B определен в лемме 2, обратим в любом пространстве X g,q при всех q > 3/2 и всех 6 6 1/2.

C Ввиду соотношения

(Bw)(t) =

/

min( t,R )

T(t s)B(s)w(s) ds

достаточно доказать однозначную разрешимость уравнения

t

w

( t ) /

T(t s)B(s)w(s) ds = f (t)

в пространстве C([0, T],Sq (Q)). Этот факт является следствием того, что оператор B имеет нулевой спектральный радиус в пространстве C ([0,T ],Sq (Q)) при q >  3/2. B

Пусть w — решение уравнения (11) из пространства X g,q , где q — одно из чисел интервала (3, го ), а 6 — одно из чисел, удовлетворяющих неравенству 1/2 3/(2q) <  6 <  min(1/2, (3/2)(1/2 1/q)). Существование малых решений с такими 6, q доказано в теореме 2. Предварительно покажем, что w G X ^,q при заданном q и всех ц < min(1/2, (3/2)(1/2 1/q)). Для элемента Kw имеем в силу леммы 3 включение K w G X pq , где ц — любое число, удовлетворяющее неравенству ц < 6 i = 26 + 3/(2q) 1/2. Применяя лемму 5, получаем, что для w справедливо включение w G X pq при ц < min(6 i , 1/2). Если 6 i > 1/2, то указанное свойство доказано. В противном случае, повторяя приведенное рассуждение, устанавливаем, что Kw и, следовательно, w принадлежат X pq при всех ц <  min(1/2, 6 2 = 26 i + 3/(2q) 1/2). Так как последовательность 6 n = 26 n- i + 3/(2q) 1/2 обладает свойством 6 n 6 n- i = 6 n- i + 3/(2q) 1/2 >  6 + 3/(2q) 1/2 > 0, то при некотором n 6 n становится больше 1/2, что опять влечет выполнение нужного свойства.

Приведенный метод рассуждений обычно называют шнуровкой. Таким образом, шнуровка приводит к следующим выводам: так как при 3 < q 6 6 имеем (3/2)(1/2 1/q) 6 1/2, то w G X ^,q при всех ц < (3/2)(1/2 1/q); в случае q > 6 w G X ^,q при всех ц < 1/2.

Следующие рассуждения показывают, что из принадлежности решения w одному из пространств X g,q следует его принадлежность некоторой совокупности пространств X ^,r . Ввиду предыдущего заключения можно считать, что 6 < min(1/2, (3/2)(1/2 1/q)) и сколь угодно близко к этому числу.

Рассмотрим отдельно два случая 3 < q 6 6 и 6 < q <  го .

В первом случае (3/2)(1/2 1/q) 6 1/2, а F(t) в силу леммы 1 принадлежит всем пространствам X ^,r для всех r > 2 и ц < (3/2)(1/2 1/r). Ввиду леммы 3 Kw принадлежит всем пространствам X ^,r для всех ц и r, удовлетворяющих условиям 0 6 2/q 1/r < 1/3, ц < 26 + (3/2)(2/q 1/r) 1/2. Учитывая выбор 6, получаем, что можно считать, что ц сколь угодно близко к числу 2((3/2)(1/2 1/q)) + (3/2)(2/q 1/r) 1/2 = (3/2)(2/3 1/r). Следовательно, согласно лемме 5 w(t) G X ^,r для всех ц и r, удовлетворяющих неравенствам ц < (3/2)(1/2 1/r), 0 6 2/q 1/r < 1/3. Применение шнуровки приводит к заключению, что в первом случае w(t) G X ^,r при всех r, ц, удовлетворяющих неравенствам 3 < r 6 6, ц < (3/2)(1/2 1/r).

Далее, выбирая q достаточно близким к 3 и повторяя предыдущее рассуждение, устанавливаем, что w G X ^,r для всех r £ (2, 3], ц < (3/2)(1/2 1/r).

Таким образом, установлено, что в первом случае w £ X ^,r для всех r £ (2, 6], ц < (3/2)(1/2 1/r).

Заметим, что этот результат нельзя улучшить.

Обратимся ко второму случаю q £ (6, го ). Так как w £ X pq для всех ц < 1/2, то в силу лемм 3 и 5 w £ X ^,r для всех r >  6, r £ [q/2, го ), ц < 1/2. Следовательно, применение шнуровки, если это необходимо, приводит к включению w £ X ^,r для всех r £ (6, го ), ц < 1/2.

Далее, если w £ X §,q , где q достаточно близко к 6, а 6 >  1/2 3/(2q), то w £ X ^, r для всех ц, r с условиями 0 6 2/q 1r < 1/3, ц < min(1/2, (3/2)(1/2 1/r)). Здесь важно, что r может принимать значения как большие, так и меньшие 6.

Учитывая изложенное, приходим к выводу: из принадлежности w одному из пространств X§,q, где параметры удовлетворяют условиям q £ (3, го), 6 £ (1/2 — 3/(2q), min(1/2, (3/2)(1/2 — 1/q))), следует, что w £ X^,r при всех r > 2, ц < min(1/2, (3/2)(1/2 — 1/r)).

В отличие от первого случая при q £ (6, го ) граница для ц может быть улучшена. Учитывая последнее утверждение, для Bw имеем оценку

k(Bw)(t)k r 6 С , (1 + t) -^                               (19)

для всех r > 2, ц <  1/2 + (3/2)(1/2 1/q).

В то же время согласно лемме 3 справедлива оценка

k (Kw)(t) k r 6 c r,^ (1+ t) -^                              (20)

для всех r >  6, ц < 1/2 + (3/2)(1/3 1/r) = (3/2)(2/3 1/r). Ввиду того, что показатели для ц в (19), (20) больше, чем показатель для F(t) (см. лемму 1) окончательно получаем, что при всех r > 6 w £ X ^,r при всех ц < (3/2)(1/2 1/r).

Таким образом, установлена

Теорема 3. Пусть w — решение уравнения (11), принадлежащее одному из пространств Xg,q при выполнении условий q > 3, 1/2 — 3/(2q) < 6 < min (1/2, (3/2)(1/2 — 1/q)).

Тогда w £ X ^,r для всех r > 2, ц < (3/2)(1/2 1/r) и, следовательно,

||v(t) V a 2 k r 6 C r,^ (1 + t) -^

при r > 2 , ц < (3/2)(1/2 1/r).

Список литературы О существовании переходов между стационарными режимами задачи обтекания

  • Finn R. Stationary solutions of the Navier-Stokes equations//Proc. Symp. Appl. Math.-Providence (R.I.): Amer. Math. Soc., 1965.-Vol. 17.-P. 121-153.
  • Galdi G. P., Heywood J. G., Shibata Y. On the global existence and convergence to steady state of Navier-Stokes flow past an obstacle that is started from rest//Arch. Rational Mech. Anal.-1997.-Vol. 138.-P. 307-318.
  • Galdi G. P. An introduction to the mathematical theory of Navier-Stokes equations. Vol. 1. Linearizad steady problems.-New York: Springer-Verlag, 1994.-448 p.
  • Сазонов Л. И. Гидродинамический проектор во внешней области.-Ростов-на-Дону: РГУ, 2000.-14 с.-Деп. в ВИНИТИ, № 3148.
  • Сазонов Л. И. Обоснование метода линеаризации в задаче обтекания//Изв. РАН. Сер. мат.-1994.-Т. 58, № 5.-С. 85-109.
  • Kobayashi T., Shibata Y. On the Oseen equation in exterior domains//Math. Ann.-1998.-Vol. 310.-P. 1-45.
  • Сазонов Л. И. Оценки возмущенной полугруппы Озеена//Владикавк. мат. журн.-2009.-Т. 11, вып. 3.-С. 51-61.
  • Бабенко К. И. О стационарных решениях задачи обтекания тела вязкой несжимаемой жидкостью//Мат. сб.-1973.-Т. 91, № 1.-С. 3-26.
  • Сазонов Л. И. Об асимптотике решения задачи трехмерного обтекания вдали от обтекаемых тел//Изв. РАН. Сер. мат.-1995.-Т. 59, № 5.-С. 173-196.
  • Сазонов Л. И. Трехмерная стационарная задача обтекания при малых числах Рейнольдса. Построение решения в виде ряда.-Ростов-на-Дону: ЮФУ, 2008.-34 с.-Деп. в ВИНИТИ, № 773.
  • Galdi G. P. Further properties of steady-state solutions to the Navier-Stokes problem past a three-dimensional obstacle//J. of Math. Phys.-2007.-Vol. 48.-P. 1-43.
Еще
Статья научная