Доклады Российской академии наук. Физика, технические науки, 2021, T. 496, № 1, стр. 24-26

САМОИНДУЦИРОВАННЫЕ ПРОЦЕССЫ ВЯЗКО-НЕВЯЗКОГО ВЗАИМОДЕЙСТВИЯ В ЛАМИНАРНОМ ПОГРАНИЧНОМ СЛОЕ НАД ПОРИСТОЙ ПОВЕРХНОСТЬЮ

Член-корреспондент РАН И. И. Липатов 12*

1 Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского
Жуковский, Московская обл., Россия

2 Московский физико-технический институт (национальный исследовательский университет)
Жуковский, Московская обл., Россия

* E-mail: igor_lipatov@mail.ru

Поступила в редакцию 13.10.2020
После доработки 13.12.2020
Принята к публикации 15.12.2020

Полный текст (PDF)

Аннотация

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

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

При обтекании пористых поверхностей может реализоваться вдув газа, в условиях, если давление под пористой поверхностью превосходит давление над поверхностью.

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

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

Существенно, что граничное условие на стенке, описывающее связь между вертикальной скоростью и перепадом давления, является линейной и не меняет однородность линейной задачи. Граничное условие этого вида позволяет пересмотреть ранее полученные результаты [14], описывающие самоиндуцированные процессы взаимодействия течения в пограничном слое с внешним сверхзвуковым потоком.

Используя полученные в [1] результаты, можно записать задачу следующим образом:

(1)
$\begin{gathered} u\frac{{\partial u}}{{\partial x}} + {v}\frac{{\partial u}}{{\partial y}} + \frac{{\partial p}}{{\partial x}} = \frac{{{{\partial }^{2}}u}}{{\partial {{y}^{2}}}}, \\ \frac{{\partial u}}{{\partial x}} + \frac{{\partial {v}}}{{\partial y}} = {\text{0,}}\quad p = - \frac{{\partial A}}{{\partial x}}, \\ \end{gathered} $
$\begin{gathered} u = {\text{0,}}\quad {v} = - \beta p\quad {\text{при}}\quad y = {\text{0;}} \\ u = y + A(x)\quad {\text{при}}\quad y \to \infty ; \\ u = y\quad {\text{при}}\quad x \to - \infty , \\ \end{gathered} $
где для безразмерных параметров, соответствующих декартовым координатам, отсчитываемым вдоль поверхности и по нормали к ней, для компонентов вектора скорости, давления, толщины вытеснения и коэффициента пористости приняты следующие обозначения: $x$, $y$, $u$, ${v}$, $p$, A, $\beta $.

Эта задача отличается от задачи, описывающей возмущенное течение около непроницаемой стенки вследствие появления условия для вертикальной скорости на стенке [1, 2].

Для малых возмущений решение представимо в виде $u = y + {{u}_{1}}$, $p = {{p}_{1}}$. Такое решение приводит к следующему виду уравнений в первом приближении:

(2)
$\begin{gathered} y\frac{{\partial {{u}_{1}}}}{{\partial x}} + {{{v}}_{1}} + \frac{{\partial {{p}_{1}}}}{{\partial x}} = \frac{{{{\partial }^{2}}{{u}_{1}}}}{{\partial {{y}^{2}}}}, \\ \frac{{\partial {{u}_{1}}}}{{\partial x}} + \frac{{\partial {{{v}}_{1}}}}{{\partial y}} = 0. \\ \end{gathered} $

Решение задачи представимо в виде $({{u}_{1}},{{{v}}_{1}},{{p}_{1}},A)$ = = ${{e}^{{\alpha x}}}(U,V,P,B)$:

(3)
$\begin{gathered} y\alpha U + V + \alpha P = U\,{\text{''}}, \\ \alpha U + V\,{\text{'}} = 0,\quad U(\infty ) = B. \\ \end{gathered} $

После ряда преобразований можно получить уравнение Эйри,

(4)
$y = {{\alpha }^{{ - 1/3}}}Y,\quad U\,{\text{'}} = F,\quad YF = F\,{\text{''}},$
имеющему следующее решение $F = CAi(Y)$:

(5)
$U(\infty ) = \int\limits_0^\infty {FdY} = C{{\alpha }^{{ - 1/3}}}\int\limits_0^\infty {Ai(Y)dY} = \frac{C}{{3{{\alpha }^{{1/3}}}}},$
$\begin{gathered} P = - \alpha B = - \frac{{C{{\alpha }^{{2/3}}}}}{3}, \\ (\alpha - \beta )P = {{\alpha }^{{2/3}}}U_{w}^{{''}} = C{{\alpha }^{{1/3}}}Ai'(0), \\ \end{gathered} $
$ - \frac{{C{{\alpha }^{{2/3}}}}}{3}(\alpha - \beta ) = C{{\alpha }^{{1/3}}}Ai'(0).$

Вводя новые переменные ${{\alpha }^{{1/3}}}(\alpha - \beta ) = - 3Ai{\text{'}}(0)$, α = ${{\delta }^{3}}\beta $, окончательно приходим к соотношению, связывающему инкремент роста возмущений и параметр проницаемости стенки:

$\delta ({{\delta }^{3}} - 1) = - \frac{{3Ai{\text{'}}(0)}}{{{{\beta }^{{4/3}}}}}.$

Можно рассмотреть два предельных случая. Первый соответствует малому коэффициенту пористости $\beta \ll 1$

(6)
$\delta \approx {{\left[ { - \frac{{3Ai'(0)}}{{{{\beta }^{{4/3}}}}}} \right]}^{{1/4}}} \approx \frac{1}{{{{\beta }^{{1/3}}}}}{{\left[ { - 3Ai{\text{'}}(0)} \right]}^{{1/4}}}.$

Тогда поскольку $\alpha = {{\delta }^{3}}\beta $, можно получить

(7)
$\alpha \approx [ - 3Ai{\text{'}}0){{]}^{{3/4}}}.$

Аналогичный результат получен для непроницаемой стенки [1, 2].

В другом предельном случае при больших значениях коэффициента пористости $\beta \to \infty $

(8)
$\delta \approx 1,\quad \alpha \approx \beta ,\quad \beta \gg 1.$

Тогда большие величины инкремента роста α соответствуют малым величинам протяженности области возмущенного течения.

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

В настоящей работе изучены режимы взаимодействия, связанные с разрежением. Такого рода течения могут реализовываться, если, например, в донной области задан отрицательный перепад давления или на пограничный слой падает волна разрежения. В этом случае при наличии пористого участка поверхности реализуется вдув газа.

В общем случае задача (1) должна быть решена численно. Численный метод описан в [3, 4]. В результате расчетов получены распределения индуцированного градиента давления и поверхностного трения как функции индуцированного возмущения давления. Отметим, что и распределение толщины вытеснения A(x) неизвестно заранее и определяется в результате решения. Представленные на рис. 1, 2 результаты соответствуют двум значениям коэффициента пористости $\beta $ = 0.25, 0.75.

Рис. 1.

Зависимость индуцированного возмущения давления p от продольной координаты x при различных значениях коэффициента пористости (1$\beta = 0.25$, 2$\beta = 0.75$).

Рис. 2.

Зависимость поверхностного трения $\tau = \frac{{\partial u}}{{\partial y}}$ от продольной координаты x при различных значениях коэффициента пористости (1$\beta = 0.25$, 2$\beta = 0.75$).

Следует отметить, что исследуемый режим взаимодействия течения в пограничном слое с внешним сверхзвуковым потоком характеризуется существенным влиянием, по крайней мере, двух факторов. Один связан с влиянием донного перепада давлений и разгоном течения в результате передачи возмущений в пограничном слое вверх по потоку от донного среза. Второй связан с наличием пористого участка поверхности, вдув через который приводит к росту толщины вытеснения и ослабляет эффект разгона. При этом увеличение отрицательного значения перепада давления приводит к росту скорости поверхностного вдува. Для неразгонного течения это рано или поздно привело бы к появлению отсоединения пограничного слоя [6] или отрыва пограничного слоя, инициируемого вдувом. Но в данном случае в условиях разгона процесс оказывается более сложным, и отрыв пограничного слоя не происходит. Хотя в общем случае можно предположить, что разгонное течение будет плавно переходить в течение сжатия с появлением или обычного отрыва, или отсоединения пограничного слоя.

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

Список литературы

  1. Нейланд В.Я. К теории отрыва ламинарного пограничного слоя в сверхзвуковом потоке // Изв. АН СССР. МЖГ. 1969. № 4. С. 53–57.

  2. Stewartson K., Williams P.G. // Self-induced separation. Proc. Roy. Soc. A. 1969. V. 312. № 1509. P. 181–206.

  3. Нейланд В.Я., Боголепов В.В., Дудин Г.Н. Асимптотическая теория сверхзвуковых течений вязкого газа. М.: Наука, 2004. 536 с.

  4. Neyland V.Ya., Bogolepov V.V., Dudin G.N., Lipatov I.I. Asymptotic Theory of Supersonic Viscous gas Flows. Elsevier Ltd, 2008. 536 p.

  5. Данилов Д.С., Толкачев Г.Ю., Липатов И.И. Самоиндуцированный отрыв ламинарного пограничного слоя и процессы вязко-невязкого взаимодействия в течении около пористой поверхности // Письма в ЖТФ. 2010. Т. 36. № 19. С. 68–75.

  6. Липатов И.И. Отсоединение ламинарного пограничного слоя при вдуве в сверхзвуковой поток // Изв. АН СССР. МЖГ. 1987. № 6. С. 16–20.

  7. Lipatov I.I. Some models of viscous-inviscid interaction // Phil.Trans. of the Royal Soc. Ser. A. 2005. V. 363. P. 1830.

Дополнительные материалы отсутствуют.