Слабые решения предельных задач Крокко

-

математическая физика

DOI: 10.18721/JPM.11303 УДК 517

слабые решения предельных задач крокко

М.Р. Петриченко, Д.Д. Заборова, Е.В. Котов, Т.А. Мусорина

Санкт-Петербургский политехнический университет Петра Великого, Санкт-Петербург, Российская Федерация

В работе предлагается прием построения приближенного решения типичной предельной задачи Крокко, который состоит в замене этой исходной задачи интегральным уравнением. Последнее решено прямым вычислением интеграла с использованием теоремы о среднем. Параметр осреднения исключен интегрированием по параметру в промежутке (0, 1). Продемонстрированы расширения способа решения и найдены слабые решения. Для классического случая полученное слабое решение незначительно отличается от точного решения Блазиуса. Приближенное значение постоянной Блазиуса оказывается равным 1/3 и отличается от точного (0,33206) на 0,3 %.

Ключевые слова: задача Коши, интегральное уравнение, теорема о среднем, группа преобразований, уединенная волна

Ссылка при цитировании: Петриченко М.Р., Заборова Д.Д., Котов Е.В., Мусорина Т.А. Слабые решения предельных задач Крокко // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2018. Т. 11. № 3. С. 27-38. DOI: 10.18721/JPM.11303

WEAK SOLuTIONS OF THE CROCCO BOuNDARY PROBLEMS M.R. Petrichenko, D.D. Zaborova, E.V. Kotov, T.A. Musorina

Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russian Federation

A procedure for designing an approximate solution of the Crocco boundary typical problem has been proposed in the paper. The procedure calls for the change of this initial problem by a nonlinear integral equation. The latter was solved by direct calculation of the integral using the mean-value theorem. . The averaging parameter was eliminated by integrating over the parameter in the (0, 1) interval. Widening the scope of the solution procedure was demonstrated and weak solutions were found. For the classical case, the weak solution was not too different from the Blasius exact one. The approximate value of the Blasius constant turned out to be 1/3 and differed from the exact one (0.33206) by 0.3 %.

Key words: Cauchy problem, integral equation, mean value theorem, group of transformations, solitary wave

Citation: M.R. Petrichenko, D.D. Zaborova, E.V. Kotov, T.A. Musorina, Weak solutions of the Crocco boundary problems, St. Petersburg Polytechnical State University Journal. Physics and Mathematics. 11 (3) (2018) 27-38. DOI: 10.18721/JPM.11303

Введение

Известно, что предельные задачи Крок-ко, в основном, связаны с гидродинамическими приложениями, в частности с продольным обтеканием пластины вязким потоком, с нестационарной фильтрацией в однородной и изотропной (скалярной) пористой среде [1 — 3].

Типичная предельная задача Крокко ставится следующим образом [1]:

2 ф —¡~2 + и = 0,

йи

Дф) = (и : 0 < и0 < и < 1), ф е С(2)(Дф)), =Ф(1) = 0,

йи )и=и0

1т(ф) = (0, а), а := ф(и0),

причем в классическом случае Блазиуса

и = 0.

В представленном виде задача (1) широко используется в гидродинамике, где переменная и трактуется как продольная скорость, а распределение ф — как напряжение трения [1].

В задачах теории фильтрации задача (1) возникает при расчете уединенной волны расхода, т. е. при решении предельной задачи для уравнения Буссинеска [2, 3]:

ди д (, ди \ — = — I ки — , д t д51 дs

где и = и(/,5) < 1 — глубина фильтрационного потока > 0, я > 0); к — коэффициент фильтрации.

В частном случае

и(0^) — 1= и(^0) = 0.

В общем случае

к = к(и), с = - к (и )(| ^ ,

где с — скорость фильтрации.

В классическом случае Буссинеска к = 1, с = 1. Тогда уравнение Буссинеска имеет вид

ди д ( ди — = —I и — д t д5 I дs

Наконец, уравнение Крокко возникает в задачах о струйных движениях вязкой

жидкости (свободная конвекция в обогреваемых каналах, свободные затопленные и пристеночные струи и др.) [4, 5].

В отличие от «естественной» постановки прикладных предельных задач, задача (1) удобна тем, что позволяет найти инъектив-ное отображение компакта (и0,1) в компакт

(0, а), ф : (и0,1) ^ (0, а).

Выражаясь точнее, утверждаем, что каждая ветвь решения предельной задачи (1) — это 2-диффеоморфизм ф : (и0,1) (0, а).

Решения предельной задачи (1) приведены в работах [6 — 32]. Эти работы распадаются на два класса:

в первом используются, в основном, аналитические методы и в том числе решения, полученные в виде степенных и расщепляющих (плоских) рядов;

во втором доминируют численные решения.

К аналитическим относятся работы (мы их относим к первому классу), использующие методы теории групп преобразований Ли и разложения в степенные и обвертывающие ряды.

Например, уравнение Буссинеска (1б) допускает линейное преобразование

% = а ± л/а5,

и тогда существует решение уравнения Буссинеска в виде уединенной волны расхода:

% = и + с1 - с21п(и + с2).

В этом решении случай с2 = 0 отвечает центрированной волне расхода, распространяющейся со скоростью ±а1/2 либо вверх, либо вниз по течению.

Особое место занимают работы [26 — 31], содержащие аналитический аппарат для построения решений обыкновенных дифференциальных уравнений (ОДУ) рассматриваемого типа вблизи сингулярной особенности. В указанных публикациях обращается внимание на то, что свойства аналитической функции определяются в основном ее особенностями, которые невозможно изучать, оставаясь на вещественном интервале. «Выход» в комплексную плоскость автоматически означает

выход (отображение) на риманову поверхность решения [31].

«Точные» решения предельной задачи для уравнения Крокко получаются также при использовании степенных рядов по и. Но для этих решений неизвестны тау-беровские теоремы, а именно — ряды для функции ф(и) оказываются «плохо» сходящимися при и ^ 1 — 0. Так, в работе [32] продемонстрирована расходимость ряда для ф(и) на внешней границе пограничного слоя (и ^ 1 — 0) и бифуркация решения во внешней, струйной части пограничного слоя.

Из приведенного обзора опущены так называемые интегральные методы, в которых вместо уравнений на плотности распределения решаются ОДУ для самих распределений (интегральных соотношений). Эти методы идейно ближе к совершенно иным методам — прямым или вариационным.

Таким образом, на сегодняшний день существует единственная альтернатива численным методам решения уравнений типа Блазиуса, Шази, Крокко — плоские ряды для аналитических решений.

Цель настоящего исследования — построить приближенное решение типичной предельной задачи Крокко с использованием процедуры осреднения.

Построение решения задачи

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

Суть этого приема следующая. Пусть

/(и) Е С(1)(0,1), Ли) > 0.

Предельная задача Крокко для промежутка (0, 1) имеет вид

и допускает формальное понижение порядка:

dф U f (v)dv

2— = -Г,

du J

йи 00 ф

Для положительной ветви решения Ф := ф+ функция 1 / ф(и) представляет собой положительное, монотонно убывающее распределение, отображающее промежуток и е (и0, 1) на промежуток фЕ (0, а), где а := ф°0).

Тогда, согласно теореме Бонне, справедливо равенство

^ и

du

j f (v)dv,

2ф ^ф + f (u) = 0, ф '(0) = ф(1) = 0 du

где 0 < 0 < 1 — параметр (правильная дробь).

Пусть в этой формуле u = 1. В этом случае при u ^ 1—0

dфе2 / du = 0(1).

Это означает, что необходимо выполнение равенств

фе (u) = 0(sm), dфе / du = 0(s—m),

m = n,

где m, n — положительные параметры.

Интегрируя равенство (#) еще раз, получим:

1 v

фе (u) = j dv j f (t )dt.

u ev

Это равенство и есть приближенное e -решение уравнения Крокко (отмечено нижним индексом е). Данное решение зависит непрерывным образом от дроби (параметра) е, причем ясно, что

1 u

ф?^) = 0, ф0М = j dv j f (t )dt =

u 0

1u

=j(1 -1 )f (t )dt -j(u -1 )f (t )dt > фе (u) > 0, 0 0

уе e (0,1),

и вообще Эф / Эе < 0.

Исключить параметр е можно, например, осреднением по нему дериватива:

М d е, du { du

что приводит к выражению й ф2

= -1/ и| V/(V)Л.

Окончательно можем записать следующее приближенное решение предельной задачи (1):

1 1 и

ф2 (и) = Г V/(V) 1п1 ¿V - Г V/(V) 1п и dv. (2) 0 V 0 v

Полученный результат не зависит от порядка интегрирования по параметру 0 и по аргументу и. Решение (2) назовем слабым 0 -решением.

Свойства решений предельной задачи

В данном разделе мы перечислим свойства решений предельной задачи (1) (доказательства этих свойств опускаются).

1. Предельные условия в задаче (1) можно заменить одноточечными условиями (Коши):

¿ф | = ф(0) - а = 0, аи

причем в начальных условиях (3) параметр а подбирается так, чтобы ср(1) = 0. Такое наложение условий оправдано в силу непрерывной зависимости ф от параметра а.

2. Существуют две ветви решения предельной задачи (1) и, соответственно, предельной одноточечной предельной задачи (3):

ф+

(рис. 1). Эти ветви связаны следующим образом:

ф+ (и) + ф (и) = 0, 0 < и < 1,

причем

0 <ф+ (и) < а, аф+< 0, ¿ф- < 0; йи йи

а <ф-(и) < 0, аф-> 0, ¿ф- > 0.

йи йи

Предельная задача (1) типична, поскольку к ней сводится, например, однородная предельная задача Крокко. Пусть и0 = 0 и вместо представления (1) рассматривается

Рис. 1. Решение типичной предельной задачи Крокко: положительная (ф+(и)) и отрицательная (ф— (и)) монотонные ветви; вертикальными пунктирами показаны границы интервала

однородная предельная задача Крокко:

т й 2ф п

2ф + и = 0'

ф(0) = ф(1) = 0.

Решение однородной предельной задачи (1в) состоит из двух ветвей:

0 < ф+ (и) и ф- (и) < 0

(положительной и отрицательной), при этом таких, что

ф+ (и) + ф- (и) = 0

для каждого значения и из промежутка 0 < и < 1.

Существует значение и = и* (0 < и*< 1), такое, что йф± / йи = 0 (согласно теореме Ролля). Поэтому каждая из ветвей решения, ф± (и), однородной задачи (1), распадается на два решения типичной предельной задачи Крокко (1):

ф/± (и), Дф; ±) = (0, и *);

фг ± (и), Б(фг ±) = (и *,1).

При этом решения ф(и) = ф±г(и) положительных и отрицательных типичных предельных задач сопрягаются непрерывно и гладко в точке и = и* (рис. 2):

ф± (0) =

й ф±

( ОфГ

йи

йи

= ф± (1) = 0.

и=и +0

Ветви решения однородной предельной задачи распадаются на решения типичных предельных задач:

ф± (и* - 0) -ф± (и* + 0) = 0,

( + Л

й2 ф йи2

< 0,

( Л

йф йи2

> 0.

/и=и V у и=и

3. Решение предельной задачи (1) — (3) удовлетворяет тождеству

1 й"

1 - и2

4. Решение предельной задачи (1) — (3) равносильно задаче на минимум положительного функционала (распределения):

^(ф) = (1/2) Г

1 /v л \2

йф ] , а 1 + и 1п —

йи

Л

ф

йи > 0.

Другими словами, вдоль экстремалей функционала Дф) выполняется условие < , где й¥ — вариация вдоль харак-

Ф

а-

Ф+ г Фн г

и

0 1

Ф" Ф" г

а- н —

Рис. 2 . Зависимости ф±(и), на которых показано, как ветви решения однородной предельной задачи распадаются на решения (##) типичных предельных задач; и" соответствует максимумам функции |ф|

теристики (траектории решения), а 8Е — вариация вдоль допустимой (виртуальной) траектории. Доказательство свойства 4 основано на том, что необходимое условие минимума Е(ф) совпадает с уравнением Крокко, а достаточность условия гарантирована выпуклостью плотности лагранжиана Е (ф).

5. Свойство решения при и0 = 0. В этом случае предельная задача (1) равносильна следующему нелинейному интегральному уравнению:

й ф ¿и

1 и ^¿у

* I

ф(и) = 2

' 2 0 ф(у)'

2I ¿у I

ф

0 ф

} (и - г) tdtл

i

ф

Для решения данного интегрального уравнения можно использовать итерационный процесс.

Пусть нижний индекс обозначает номер итерации, и тогда процесс решения выражается как

¿и

\

vdv

фs-1(у У

ф,(и) = ¿у]- Ш

0с ф5-1

Если 1 / ф Е А(0,1), то ф Е С(1)(0,1). Допустим, что последовательность итераций функции 1 / ф5 образует последовательность Коши. Тогда ф8 ^ ф почти всюду на промежутке 0 < и < 1, ввиду полноты области Х1(0, 1).

6. Решение позволяет сформулировать следствие из теоремы о среднем. Поскольку 1 / ф(и) — монотонно возрастающее распределение, то, согласно теореме Бонне о среднем, справедливы равенства

2ф(и) ¿ф = - -е2),

¿и 2

ф2( и) = 6(1 - и3)(1 -В2),

где 9 — правильная дробь (0 < 0 < 1).

Окончательное выражение для прибли-

женного решения имеет вид

фе (и) = ал/б^а - из)(1 -е2).

Среднее квадратичное значение фе(и), т. е. 9 -аппроксимация решения, определяется равенством

ф2

Отсюда ф(и) = (1/3)л/1 - и3, и это приближение аппроксимирует точное решение Блазиуса, особенно при малых значениях и. Например, значение постоянной Блазиуса составляет а = 1/3. Ее точное значение, недавно анонсированное В.П. Вариным, составляет [29, 30]:

а = 0,33205733621519 629893718006201058 296654709356141267 981810047564019872 417401806440507049 0731855146368... .

Рациональное значение константы отличается от приведенного иррационального менее, чем на 0,3 %.

1

Величина © := ^^¿и в задачах с физиче-

0

ским содержанием представляет диссипацию на отрезке (0, 1). В данном случае ее значение составляет

© = (1/3)1 ^/Г^U3du = ^ „ 0,27.

'I 18 Г(11/6) '

Видно, что за диссипацию «отвечает» распределение ф(и) в окрестности точки и = 0.

7. Сформулируем общее определение нормы 9 -аппроксимации:

(и) = ^=(1 - и3) {|(1 -е2)г/2dе

ч1/г

1 - и

3 (

л/Пг

У/г

^2Г(г/2 + 3/2)) '

где г > 0 — любое положительное вещественное число.

Далее нижний индекс опускается и норма должна быть понятна из контекста. На-

пример, в предыдущем пункте было принято г = 2, и тогда получаем, что

ф^) := ф2(^ =

1 - u3

\f%

S V 2 Г

= 76vl(1 - u3) = (1/3)(1 - u3)-

Далее выполняется равномерно по 0 < u < 1 неравенство Коши — Гельдера

фМ) ^ фг+а (u), Va > 0,

и последовательность норм не убывает с возрастанием индекса r от 0 до да, а точнее,

(1/V6) exp(-с /2 - 0,02) < Щ < 1 / S,

где с — постоянная Маскерони.

8. Первое обобщенное свойство решения. Пусть предельная задача Крокко (1) имеет следующий вид:

2 ф ф + um = 0, du

ф '(0) = ф(0) = 0.

Данное представление уравнения Крок-ко возникает из уравнения Буссинеска, если

к = k(h) = k0(h / H)m-1.

Тогда 0-аппроксимация решения предельной задачи (1а) имеет вид

фе (u) =

и, если применить среднеквадратичное осреднение по е, то слабое решение имеет вид:

фМ =

4Г-

m+2

а =

m + 2 1

m + 2

Тождество (4) записывается следующим образом:

1 idS -

0 V ¿и ) 2

Условие минимума квадратичного функционала практически не меняется и выражается как

F (ф) = (1/2) j

1 (f Т + um ln а ^

du

ф

-u ^ inf > 0.

Величина диссипации в данном случае следует выражению

D =

+2 ^

m + 2

- um+2 du =

1

ri-L- w 3

m + 2 I I 2

(m + 2)2

= O

Г

1

m + 2

3m + 8 2

, m >> 1,

и она уменьшается при увеличении т.

В силу тождества (4а) справедливо выражение

ф(и) := ф(и) ,

а

из которого видно, что с увеличением параметра т степень заполнения профиля

ф

увеличивается.

Вкладывая физическое содержание в решение, предположим, что ф = ф(и) — трение в пограничном слое. Тогда трение на поверхности пластины, обтекаемой продольным вязким потоком, а = ф(0) монотонно снижается и сохраняется от пристеночной к струйной части слоя:

при т ^ да, ф(и) ^ ф(0) = 0, 0 < и < 1.

9. Второе обобщенное свойство решения. Пусть уравнение Буссинеска и предельные условия для него имеют вид

диа

~дТ

д_ ds

(

du

ds

Л

где а, Ь, с — вещественные параметры; Ди) = (г > 0, я > 0), и(0, я) - 1 = и(г, 0) = 0. Тогда соответствующее преобразование Крокко приводит это уравнение к виду

иь =

-О- ф(и) I- й | и^^а ¿ф с + 1 I йи I йи

причем предельные условия ставятся следующим образом:

Ф(1) = Ф'(0) = 0.

В этом случае 9 -решение предельной задачи (7), (8) имеет вид

с

(с + 1) сс+1 Фе (и) = —-—-~ Х

ас+1(Ь + с )с+1(Ь + ас )с+1

{(1 -еЬ/с+1)(1 - и

а+Ь/сс+1

Примечание 1. Если в частном случае а = Ь = 1, то из формулы (9) следует выражение, известное как фильтрационная модель Христиановича для плоского фильтрационного потока:

((

Фе(и) = Тогда

с+1Л

1 -е

2с+1 Л Л

1 - и

/у

с с+1

( с+1 Л с+1

1 -е с

с + 1

и, далее

с

Г| гГс(г + 1) + 1

с + П I с + 1

с(г+1)+1

(с + 1) с

при этом г > 0. Очевидно, что

ф(и) := Ф"

Г

с(г + 2) + 1 с + 1

/ 2с+1Лс+1

1 - и с

и с уменьшением параметра с происходит заполнение профиля безразмерного распределения ф(и).

Например, если с = 1/2, то

ф(и) = (1 - и4)1/3; а если с = 1/3, то

Ф(и) = (1 - и5)1/4.

Примечание 2. Пусть а = 1; Ь, с — свободные параметры. Тогда, в силу е -решения (9), получаем выражения

с

->с+1

фе(и) = (с + 1)^ {(1 - еЬ/с+1)(1 - иЬ/с+2)}с+1, (Ь + с)с+1

с

ф

и при уве__ичении параметра Ь заполнение профиля ф(и) увеличивается.

Например, если с = 1/2, то

<р(и) = (1 - и2(Ь+1))1/3 ^ 0, 0 < и < 1.

Выводы

В результате проведенного исследования установлено следующее.

Типичная предельная задача Крокко допускает положительную и отрицательную ветви решения (ф+ и ф-), такие что

Ф+ (и) + ф- (и) = 0.

Однородная предельная задача Крокко редуцируется на две типичные предельные задачи Крокко, сопрягаемые в критической точке и = и*, такой что

й ф йи

0 < и0 < и* < 1, = ф(и* - 0) - ф(и* + 0) = 0.

Типичная предельная задача Крокко равносильна нелинейному интегральному уравнению. Последнее решается прямым вычислением интеграла с использованием второй теоремы о среднем. Параметр осреднения исключается интегрированием по параметру в промежутке (0, 1).

В данной статье демонстрируются расширения предложенного способа решения. Для классического случая а = Ь = с = 1 слабое е -решение незначительно отличается от точного. Приближенное значение постоянной Блазиуса а = ф(0) оказывается равным 1/3. При этом точное значение ф(0) = 0,33206.

с

ае =

СПИСОК ЛИТЕРАТУРЫ

1. Crocco L. Sulla strato limite laminare nei gas lungo una lamina plana // Rend. Math. Appl. Ser. 5. 1941. Vol. 21. No. 2. Pp. 138-152.

2. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. 2-е изд. М.: Наука, 1977. 664 с.

3. Петриченко М.Р. Предельная задача Бус-синеска и ее гидравлическое приложение // Сб. трудов заочной конф. «Научное обозрение физико-математических и технических наук в XXI веке». М.: «Prospero», 2015. С. 3-7.

4. Петриченко М.Р. Предельные задачи Крокко в теории струй вязкой жидкости // Сб. трудов заочной конф. «Научное обозрение физико-математических и технических наук в XXI веке». М.: «Prospero», 2016. С. 10-14.

5. Nemova D., Reich E., Subbotina S., Khayrutdinova E, Petrichenko M.

Heat and mass transfer in a vertical channel under heat-gravitational convection conditions // Experimental Fluid Mechanics Proceedings of the International inference. Prague, 2015. Pp. 604-616.

6. Akdi M., Sedra M.B. Numerical solution of the Blasius problem //The African Review of Physics. 2014. Vol. 9. No. 0022. Pp. 165-168;

7. Allan F.M., Syam M.I. On the analytic solutions of the nonhomogeneous Blasius problem // Journal of Computational and Applied Mathematics. 2005. Vol. 182. No. 2. Pp. 362-371.

8. Aminikhah H. Analytical approximation to the solution of nonlinear Blasius viscous flow equation by LTNHPM // International Scholarly Research Network. ISRN Mathematical Analysis. 2012. Vol. 2012. Article ID 957473 (10 p.).

9. Aminikhah H., Kazemi S. Numerical solution of the Blasius viscous flow problem by quartic B-spline method // Hindawi Publishing Corporation. International Journal of Engineering Mathematics. 2016. Vol. 2016. Article ID 9014354 (6 p.). http:// dx.doi.org/10.1155/2016/9014354.

10. Andrzejczak G., Nockowska-Rosiak M., Przeradzk B. A note on Blasius type boundary value problems // Opuscula Math. 2013. Vol. 33. No. 1. Pp. 5-17.

11. Bataller R.C. Numerical comparisons of Blasius and Sakiadis flows // Matematika. 2010. Vol. 26. No. 2. Pp. 187-196.

12. Cortell R. Numerical solutions of the classical Blasius flat-plate problem // Applied Mathematics and Computation. 2005. Vol. 170. No. 1. Pp. 706-710.

13. Chang C.-W., Liu C.-S. The Lie-group shooting method for boundary layer equations in fluid mechanics // Proceedings of the Conference

of Global Chinese Scholars on Hydrodynamics. Beijing, 2006 -07. Pp. 103-108.

14. Fang T., Lee C.F. A moving-wall boundary layer flow of as slightly rarefied gas free stream over a moving flat plate // Applied Mathematics Letters. 2005. Vol. 18. No. 5. Pp. 487-495.

15. Tiegang Fang, Wei Liang, Chia-fon F. Lee. A new solution branch for the Blasius - A shrinking equation sheet problem // Computers and Mathematics with Applications. 2008. Vol. 56. No. 12. Pp. 3088-3095.

16. Fazio R. Blasius problem and Falkner — Skan model: Tnpfer's algorithm and its extension // Computers & Fluids. 2013. Vol. 73. 15 March. Pp. 202-209.

17. Goh J., Majid A.A., Ismail A.I.M. A quatic B-spline for second-order singular boundary value problems // Computers and Mathematics with Applications. 2012. Vol. 64. No. 12. Pp. 115-120.

18. He J.-H. A simple perturbation approach to Blasius equation // Applied Mathematics and Computation. 2003. Vol. 140. No. 2-3. Pp. 217-222.

19. Chein-Shan Liu, Jiang-Ren Chang. The Lie-groups shooting method for multiple-solutions of Falkner - Skan equation under suction-injection conditions // International Journal of Non-Linear Mechanics. 2008. Vol. 43. No. 9. Pp. 844-851.

20. Liu C.-S. Cone of non-linear dynamical system and group preserving schemes // International Journal of Non-Linear Mechanics. 2001. Vol. 36. No. 7. Pp. 1047-1068.

21. Liu C.-S. The Lie-group shooting method for nonlinear two-point boundary value problems exhibiting multiple solutions // Comput. Model. Eng. Sci. 2006. Vol. 13. No. 2. Pp. 149-163.

22. Ding Xu, Xin Guo. Fixed point analytical method for nonlinear differential equations // Journal of Computational and Nonlinear Dynamics. 2013. Vol. 8. No. 1. Р. 011005

23. Варин В.П. Плоские разложения и их приложения. М.: Препринты ИПМ РАН им. М.В. Келдыша, 2014. № 23. 25 с. URL: http:// keldysh.ru/papers/2014/prep2014_23. pdf.

24. Варин В.П. Плоские разложения и их приложения // Журнал вычислительной математики и математической физики. 2015. Т. 55. № 5. С. 807-821.

25. Варин В.П. Плоский асимптотический ряд Блазиуса // Математический форум (сборник). Сер. «Итоги науки. Юг России». Владикавказ, 2015. С. 34-47.

26. Варин В.П. Специальные решения уравнения Чейзи. Препринты ИПМ РАН им. М.В. Келдыша. 2015. № 43. 39 с. http://keldysh.ru/ papers/2015/prep2015_43_eng.pdf.

27. Варин В.П. Специальные решения уравнения Шази // Журнал вычислительной математики и математической физики. 2017. Т. 57. № 2. С. 210-236.

28. Варин В.П. Интегрирование ОДУ на Ри-мановых поверхностях как вычислительный инструмент // Теоретические основы конструирования численных алгоритмов и решение задач

математической физики. Тез. докл. XXI Всерос. конф. и молод. школы-конф., посв. памяти К.И. Бабенко. Институт прикладной математики им. М.В. Келдыша РАН. 2016. С. 73-74.

29. Фаиз Ахмад. Применение уравнение Крокко — Ванга к решению задачи Бласиуса // Электронный журнал «Техническая акустика». http://ejta.org, 2007, 2.

Статья поступила в редакцию 23.05.2018, принята к публикации 04.06.2018.

СВЕДЕНИЯ ОБ АВТОРАХ

ПЕтРИЧЕНКО Михаил Романович — доктор технических наук, заведующий кафедрой гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.

195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 fonpetrich@mail.ru

ЗАБОРОВА Дарья Дмитриевна — ассистент кафедры гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация. 195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 zaborova-dasha@mail.ru

КОтОВ евгений Владимирович — аспирант кафедры гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация. 195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 ekotov.cfd@gmail.com

МуСОРИНА татьяна Александровна — аспирантка кафедры гидравлики и прочности Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.

195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 flamingo-93@mail.ru.

REFERENCES

[1] L. Crocco, Sulla strato limite laminare nei gas lungo una lamina plana, Rend. Math. Appl., Ser. 5. 21 (2) (1941) 138-152.

[2] P.Ya. Polubarinova-Kochina, Teoriya dvizheniya gruntovykh vod [Groundwater movement theory], 2nd ed., Nauka, Moscow, 1977.

[3] M.R. Petrichenko, Predelnaya zadacha Bussineska i yeye gidravlicheskoye prilozheniye [The Boussinesq limit problem and its hydraulic application], In the collection of scientific papers "The Scientific Review of Physics, Mathematics and Engineering Sciences for the 21th Century", "Prospero" publ. house, Moscow (2015) 3-7.

[4] M.R. Petrichenko, Predelnyye zadachi Krokko v teorii struy vyazkoy zhidkosti [The Crocco limit problems in the viscous-fluid jet theory], In the collection of scientific papers "The Scientific Review of Physics, Mathematics and Engineering Sciences for the 21th Century", "Prospero" publ. house, Moscow, (2016) 10-14.

[5] D. Nemova, E. Reich, S. Subbotina, et al., Heat and mass transfer in a vertical channel under heat-gravitational convection conditions, Experimental Fluid Mechanics Proceedings of the International Conference, Prague (2015) 604-616.

[6] M. Akdi, M.B. Sedra, Numerical solution of the Blasius problem, The African Review of Physics 9 (0022) (2014) 165-168.

[7] F.M. Allan, M.I. Syam, On the analytic solutions of the nonhomogeneous Blasius problem, Journal of Computational and Applied Mathematics. 182 (2) (2005) 362-371.

[8] H. Aminikhah, Analytical approximation to the solution of nonlinear Blasius viscous flow equation by LTNHPM, International Scholarly Research Network, ISRN Mathematical Analysis. 2012 (2012) ID 957473 (10 p).

[9] H. Aminikhah, S. Kazemi, Numerical solution of the Blasius viscous flow problem by

quartic B-spline method, Hindawi Publishing Corporation, International Journal of Engineering Mathematics. 4 (2016) ID 9014354 (6 p), http:// dx.doi.org/10.1155/2016/9014354..

[10] G. Andrzejczak, M. Nockowska-Rosiak, B. Przeradzk, A note on Blasius type boundary value problems, Opuscula Math. 33 (1) (2013) 5-17.

[11] R.C. Bataller, Numerical comparisons of Blasius and Sakiadis flows, Matematika. 26 (2) (2010) 187-196.

[12] R. Cortell, Numerical solutions ofthe classical Blasius flat-plate problem, Applied Mathematics and Computation. 170 (1) (2005) 706-710.

[13] C.-W. Chang, C.-S. Liu, The Lie-group shooting method for boundary layer equations in fluid mechanics, Proceedings of the Conference of Global Chinese Scholars on Hydrodynamics, Beijing, 2006 -07, 103-108.

[14] T. Fang, C.F. Lee, A moving-wall boundary layer flow of as slightly rarefied gas free stream over a moving flat plate, Applied Mathematics Letters. 18 (5) (2005) 487-495.

[15] Tiegang Fang, Wei Liang, Chia-fon F. Lee, A new solution branch for the Blasius - A shrinking equation sheet problem, Computers and Mathematics with Applications. 56 (12) (2008) 3088-3095.

[16] R. Fazio, Blasius problem and Falkner — Skan model: T^fer's algorithm and its extension, Computers & Fluids. 73 (15 March) (2013) 202-209.

[17] J. goh, A.A. Majid, A.I.M. Ismail, A quatic B-spline for second-order singular boundary value problems, Computers and Mathematics with Applications. 64 (12) (2012) 115-120.

[18] J.-H. He, A simple perturbation approach to Blasius equation, Applied Mathematics and Computation. 140 (2003) 217-222.

[19] Chein-Shan Liu, Jiang-Ren Chang, The Lie-groups shooting method for multiple-solutions of Falkner - Skan equation under suction-injection conditions, International Journal of Non-Linear Mechanics. 43 (2008) 844-851.

[20] C.-S. Liu, Cone of non-linear dynamical system and group preserving schemes, International

Received 23.05.2018, accepted 04.06.2018.

Journal of Non-Linear Mechanics. 36 (7) (2001) 1047-1068.

[21] C.-S. Liu, The Lie-group shooting method for non-linear two-point boundary value problems exhibiting multiple solutions, Comput. Model. Eng. Sci. 13 (2) (2006) 149-163.

[22] Ding Xu, Xin Guo. Fixed point analytical method for nonlinear differential equations, Journal of Computational and Nonlinear Dynamics. 8 (1) (2013) 011005 (9 p).

[23] V.P. Varin, Flat expansions and their applications, Moscow, KIAM Preprint (23) (2014). URL: http://keldysh.ru/papers/2014/prep2014_23. eng.pdf.

[24] V.P. Varin, Flat expansions and their applications, Computational Mathematics and Mathematical Physics. 55 (5) (2015) 797-810.

[25] V.P. Varin, Ploskiy asimptoticheskiy ryad Blaziusa [The Blasius's flat asymptotic expansion], In: Mathematics Forum Ser. "Sciences summary, The South of Russia" (2015) 34-47.

[26] V.P. Varin, Special solutions to Chazy equation, Moscow, KIAM Preprint (43) (2015). http://keldysh.ru/papers/2015/prep2015_43_eng.pdf.

[27] V.P. Varin, Special solutions to Chazy equation, Computational Mathematics and Mathematical Physics. 57 (2) (2017) 211-235.

[28] V.P. Varin, Integrirovaniye ODU na Rimanovykh poverkhnostyakh kak vychislitelnyy instrument, [ODE integration on the Riemann surfaces as numerical instrument], In: Teoreticheskiye osnovy konstruirovaniya chislennykh algoritmov i resheniye zadach matematicheskoy fiziki [Theoretical foundation of synthesizing the numerical algorithms and the solution of math-physical problems], Tezisy dokladov XXI Vserossiyskoy konferentsii i Molodezhnoy shkoly-konferentsii, posvyashchennoy pamyati K.I. Babenko [Sci. abstracts of the 21th All-Russian Conf. and the Youth School-Conf. Dedicated to the Memory of K.I. Babenko], KIAM (2016) 73-74.

[29] Faiz Ahmad, Application of Crocco - Wang equation to the Blasius problem, Electronic Journal "Technical Acoustics", http://www.ejta.org, 2007, 2.

THE AuTHORS

PETRICHENKO Mikhail R.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

fonpetrich@mail.ru

ZABOROVA Dariya D.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

zaborova-dasha@mail.ru

KOTOV Eugeniy V.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

ekotov.cfd@gmail.com

MUSORINA Tatiana A.

Peter the Great St. Petersburg Polytechnic University

29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation

flamingo-93@mail.ru

© Санкт-Петербургский политехнический университет Петра Великого, 2018