Second integral generalization of the Crocco invariant for 3D flows behind detached bow shock wave

Full Text

Введение

Полная энергия (полная энтальпия) газа не меняется при переходе через скачок уплотнения. Поэтому если набегающий сверхзвуковой поток однороден, то течение за скачком остается изоэнергетическим. Отошедший головной скачок имеет искривленную форму, и течение за ним становится вихревым (кроме лидирующей линии тока, на которой завихренность равна нулю [1]). В осесимметричном (без закрутки) случае для изоэнергетических течений Л. Крокко показал [2], что вдоль линий тока сохраняется отношение завихренности к давлению, умноженному на расстояние до оси симметрии \((I_K=\mathrm{\Omega}/pr)\). Для закрученных осесимметричных течений инвариант Л. Крокко обобщен в [3]. Оказалось, что окружная составляющая завихренности \(\mathrm{\Omega}_{\mathrm{\varphi}}\) представима в виде \(\mathrm{\Omega}_{\mathrm{\varphi}}=r^{-1}\mathrm{\rho}C_1 + rp C_2\), где \(C_1\) и \(C_2\) — функции линий тока (\(\mathrm{\rho}\) — плотность), причем в отсутствие закрутки \(C_1 \equiv 0\).

В общем пространственном случае неинтегральный инвариант Л. Крокко был обобщен в [4] интегральным инвариантом. В [4] показано, что величина криволинейного интеграла по замкнутой вихревой линии
\[
I_1=\displaystyle \int_{\mathrm{\gamma}} (p/\mathrm{\Omega})\,\mathrm{d}l
\]
одинакова для всех вихревых линий \(\mathrm{\gamma}\), лежащих на одной изоэнтропийной поверхности. Это был первый инвариант, пригодный для верификации расчетов течений за скачком в общем пространственном случае.

Появлению инварианта \(I_1\) предшествовало обнаружение факта замкнутости вихревых линий за отошедшим скачком [1]. Если рассматривать векторные поля, связанные с течением (то есть такие, которые выражаются через скорость, плотность и давление газа и через производные этих параметров), то следует отметить, что недавно в [5] найдено еще одно (кроме завихренности) векторное поле с замкнутыми векторными линиями. Речь идет о векторе \(\mathbf{a}=\mathbf{V}\times\boldsymbol{\nabla}\mathrm{\sigma}\), где \(\mathrm{\sigma}=p\mathrm{\rho}^{-k}\) (\(\mathbf{V}\), \(\mathrm{\rho}\) и \(p\) — скорость, плотность и давление газа, \(k\) — показатель адиабаты).

Данная работа посвящена поиску интегрального инварианта, связанного с замкнутыми линиями вектора \(\mathbf{a}\).

1. Уравнения движения

Рассмотрим течение идеального (вязкость и теплопроводность отсутствуют) газа, в частицах которого выполняется соотношение \(p=R\mathrm{\rho}T\), где \(T\) — температура, \(R\) — отношение универсальной газовой постоянной к молярной массе, за отошедшим головным скачком, возникшим при обтекании тела равномерным сверхзвуковым потоком. Энтропийная функция \(\mathrm{\sigma}=p\mathrm{\rho}^{-k}\) постоянна вдоль линий тока:
\((\mathbf{V}\mathbf{\cdot}\boldsymbol{\nabla})\mathrm{\sigma}=0\), но может иметь разные значения на различных линиях тока за скачком. Течение за скачком вследствие однородности набегающего потока является изоэнергетическим, и уравнение Эйлера, записанное в форме Крокко [6], имеет вид
\begin{equation}\label{eq:a1}
\mathbf{\Omega}\times\mathbf{V}=(k-1)^{-1}RT\boldsymbol{\nabla}\mathbf{\ln \mathrm{\sigma}}, \quad \mathbf{\Omega}=\operatorname{\mathbf{rot}}\mathbf{V}.
\end{equation}
Замыкает систему уравнение неразрывности
\begin{equation}\label{eq:a2}
\operatorname{\mathrm{div}}(\mathrm{\rho}\mathbf{V})=0.
\end{equation}

2. Воображаемые частицы

Критерий Гельмгольца—Зоравского [7, 8] применительно к стационарному и соленоидальному (см. формулу \eqref{eq:a2}) полю вектора \(\mathrm{\rho}\mathbf{V}\) можно сформулировать следующим образом.

Если в области \(G\) существует такое векторное поле \(\mathbf{q}\)  (имеющее размерность скорости), что во всей области выполнено равенство
\begin{equation}\label{eq:a3}
\mathrm{\rho}\mathbf{V}\times\operatorname{\mathbf{rot}}(\mathbf{q}\times(\mathrm{\rho}\mathbf{V}))=0,
\end{equation}
то воображаемые частицы, составляющие в некоторый момент времени сегмент векторной линии \(\mathrm{\rho}\mathbf{V}\), лежащий в области \(G\), двигаясь со скоростью \(\mathbf{q}\), будут составлять сегмент одной из векторных линий \(\mathrm{\rho}\mathbf{V}\) в каждый последующий момент времени (до тех пор, пока эти частицы находятся в области \(G\)).

Такие воображаемые частицы, следуя [8], будем называть q-частицами. Очевидно, что векторные линии \(\mathrm{\rho}\mathbf{V}\) совпадают с линиями тока. Поэтому если скорость \(\mathbf{q}\) существует, можно считать, что сегменты линий тока переносятся q-частицами. Докажем существование в течении за скачком поля \(\mathbf{q}\) и найдем его выражение через параметры течения и их производные.

В течении газа за скачком будем искать \(\mathbf{q}\) в виде \(\mathbf{q}=\mathrm{\lambda}\mathbf{a}\), где \(\mathrm{\lambda}\) — некоторое скалярное поле. Тогда выражение \(\mathbf{q}\times(\mathrm{\rho}\mathbf{V})\), стоящее под знаком ротора в \eqref{eq:a3}, примет вид \(\mathrm{\lambda}\mathbf{a}\times(\mathrm{\rho}\mathbf{V})=\mathrm{\lambda}\mathrm{\rho}(\mathbf{V}\times\boldsymbol{\nabla}\mathrm{\sigma})\times\mathbf{V}\). Раскроем двойное векторное произведение и учтем ортогональность скорости \(\mathbf{V}\) и градиента \(\boldsymbol{\nabla}\mathrm{\sigma}\), вытекающую из \eqref{eq:a1}. Получим \(\mathbf{q}\times(\mathrm{\rho}\mathbf{V})=\mathrm{\lambda}\mathrm{\rho}\mathbf{V}^2 \boldsymbol{\nabla}\mathrm{\sigma}\).

Следовательно, если положить \(\mathrm{\lambda}=\mathrm{\lambda}_0 \mathrm{\rho}^{-1}\mathbf{V}^{-2}\), где \(\mathrm{\lambda}_0\) — произвольная ненулевая константа, обеспечивающая для вектора \(\mathbf{q}=\mathrm{\lambda}\mathbf{a}=\mathrm{\lambda}_0 \mathrm{\rho}^{-1}\mathbf{V}^{-2}\mathbf{a}\) размерность скорости, то уравнение \eqref{eq:a3} окажется выполненным (поскольку в этом случае под знаком ротора окажется градиент \(\boldsymbol{\nabla}\mathrm{\sigma}\), умноженный на константу \(\mathrm{\lambda}_0\)). Таким образом, существование скорости \(\mathbf{q}\) доказано, и можно считать, что линии тока переносятся q-частицами, движущимися со скоростью
\begin{equation}\label{eq:a4}
\mathbf{q}=\mathrm{\lambda}_0 \mathrm{\rho}^{-1}\mathbf{V}^{-2}\mathbf{a}=\mathrm{\lambda}_0 \mathrm{\rho}^{-1}\mathbf{V}^{-2}(\mathbf{V}\times\boldsymbol{\nabla}\mathrm{\sigma}).
\end{equation}

3. О ненулевом значении \(\mathbf{q}\)

В следующем разделе будет интегрироваться величина, обратная к \(|\mathbf{q}|\). Поэтому докажем, что всюду, за исключением лидирующей линии тока (линии торможения), величина \(\mathbf{q}\) отлична от нуля. Будем исходить из общепринятого [4, 5, 9–11] допущения о том, что по крайней мере в некоторой окрестности выпуклой носовой части скорость \(\mathbf{V}\) обращается в нуль только в передней точке торможения. Поэтому исходя из \eqref{eq:a4} и из ортогональности \(\mathbf{V}\) и \(\boldsymbol{\nabla}\mathrm{\sigma}\) приходим к выводу, что величина \(\mathbf{q}\) может обращаться в нуль только вместе с \(\boldsymbol{\nabla}\mathrm{\sigma}\).

Докажем, что всюду, за исключением лидирующей линии тока, величина \(\boldsymbol{\nabla}\mathrm{\sigma}\), а вместе с ней и величина \(\mathbf{q}\), отличны от нуля. Поскольку \((\mathbf{V}\mathbf{\cdot}\boldsymbol{\nabla}\mathrm{\sigma})=0\), градиент этого скалярного произведения также будет равен нулю, то есть \(\boldsymbol{\nabla}(\mathbf{V}\mathbf{\cdot}\boldsymbol{\nabla}\mathrm{\sigma})=0\). Используя известное векторное тождество для градиента скалярного произведения и учитывая, что ротор градиента равен нулю, получим
\begin{equation}\label{eq:a5}
(\mathbf{V}\mathbf{\cdot}\boldsymbol{\nabla})(\boldsymbol{\nabla}\mathrm{\sigma})+((\boldsymbol{\nabla}\mathrm{\sigma})\mathbf{\cdot}\boldsymbol{\nabla})\mathbf{V}+(\boldsymbol{\nabla}\mathrm{\sigma})\times\mathbf{rot}\,\mathbf{V}=0.
\end{equation}

Это выражение замечательно тем, что компоненты \(\boldsymbol{\nabla}\mathrm{\sigma}\) дифференцируются только в первом слагаемом \((\mathbf{V}\mathbf{\cdot}\boldsymbol{\nabla})(\boldsymbol{\nabla}\mathrm{\sigma})\), а в другие слагаемые \eqref{eq:a5} компоненты \(\boldsymbol{\nabla}\mathrm{\sigma}\) входят линейно как коэффициенты при различных производных компонент скорости. Поэтому можно применить известный [1, 10, 11] способ — представить \eqref{eq:a5} в виде
\[
\frac{d}{d \, l}(\boldsymbol{\nabla}\mathrm{\sigma})=\frac{1}{|\mathbf{V}|}A(\boldsymbol{\nabla}\mathrm{\sigma}),
\]
где \(l\) — переменная длина дуги вдоль линии тока, \(A\) — матрица размером 3\(\times\)3 с коэффициентами, непрерывно зависящими от производных компонент скорости \(\mathbf{V}\). При заданном поле скорости эту систему уравнений можно считать автономной относительно компонент \(\boldsymbol{\nabla}\mathrm{\sigma}\). Поэтому, согласно известным свойствам автономных систем обыкновенных дифференциальных уравнений [12], получаем следующее

Утверждение 1. В течении за отошедшим головным скачком на любой линии тока на всем ее участке, на котором \(|\mathbf{V}|>0\), либо \(|\boldsymbol{\nabla}\mathrm{\sigma}|\equiv 0\), либо \(|\boldsymbol{\nabla}\mathrm{\sigma}|>0\).

Сразу за скачком искривленной формы завихренность равна нулю только в начале лидирующей линии [13, 14]. При этом завихренность на скачке с дозвуковой стороны лежит в касательной к скачку плоскости [1], а скорость на скачке с дозвуковой стороны имеет ненулевую нормальную к скачку составляющую. Поэтому всюду, кроме начала лидирующей линии (где \(\mathbf{\Omega}=0\)), на дозвуковой стороне скачка вектор \(\mathbf{\Omega}\times\mathbf{V}\) отличен от нуля. С учетом \eqref{eq:a1} это значит, что в начале всех линий тока на скачке, кроме начала лидирующей линии, \(|\boldsymbol{\nabla}\mathrm{\sigma}|>0\). Согласно утверждению 1, это значит, что во всем течении за скачком, кроме лидирующей линии тока, градиент энтропийной функции отличен от нуля \((|\boldsymbol{\nabla}\mathrm{\sigma}|>0)\).

Как замечено в начале данного раздела, величина \(\mathbf{q}\) может обращаться в нуль только вместе с \(\boldsymbol{\nabla}\mathrm{\sigma}\). Поэтому во всем течении за скачком, кроме лидирующей линии тока, скорость \(\mathbf{q}\), определяемая формулой \eqref{eq:a4}, отлична от нуля.

4. Интегральный инвариант

В [5] показано, что по крайней мере вблизи передней точки торможения линии вектора \(\mathbf{a}=\mathbf{V}\times\boldsymbol{\nabla}\mathrm{\sigma}\) замкнуты и один раз охватывают лидирующую линию тока (линию тока, которая пересекает отошедший скачок по нормали и которая, как доказано в [9], совпадает с линией торможения). Поскольку траектории q-частиц лежат на замкнутых линиях \(\mathbf{a}\), эти траектории (как и линии \(\mathbf{a}\)) не только лежат на изоэнтропийных поверхностях, но и один раз опоясывают их.

Применим теперь идею доказательства, предложенную в [4], для замкнутых траекторий q-частиц, лежащих на изоэнтропийных поверхностях. Если \(\mathrm{\gamma}\) — замкнутая линия вектора \(\mathbf{a}\), \(l\) — переменная длина дуги на \(\mathrm{\gamma}\), то q-частица проходит по кривой \(\mathrm{\gamma}\) расстояние \(\mathrm{d}l\) за время \(\mathrm{d}l/|\mathbf{q}|\). Поэтому время, за которое q-частица сделает полный оборот по замкнутой линии \(\mathrm{\gamma}\), равно криволинейному интегралу \(\displaystyle \int_{\mathrm{\gamma}}|\mathbf{q}|^{-1}\,\mathrm{d}l\). Если q-частицы составляли в начальный момент времени одну из линий тока, то, сделав один оборот, они должны составлять ту же самую линию тока. Поэтому время полного оборота у всех q-частиц, лежащих на одной изоэнтропийной поверхности, одинаково, и величина \(I=\displaystyle \int_{\mathrm{\gamma}}|\mathbf{q}|^{-1}\mathrm{d}l\) есть новый (второй) интегральный инвариант изоэнтропийных поверхностей. Используя \eqref{eq:a4} и \eqref{eq:a1}, получим
\begin{equation}\label{eq:a6}
\mathbf{q}=\mathrm{\lambda}_0 \mathrm{\rho}^{-1}\mathbf{V}^{-2}\mathbf{V}\times\boldsymbol{\nabla}\mathrm{\sigma}=\mathrm{\lambda}_0 (k-1)\mathrm{\sigma}\frac{1}{p}(\mathbf{\Omega}-\mathbf{V}(\mathbf{\Omega}\mathbf{\cdot}\mathbf{V})/\mathbf{V}^2).
\end{equation}

Легко проверить, что вектор \(\mathbf{\Omega}_{\mathbf{a}}=\mathbf{\Omega}-\mathbf{V}(\mathbf{\Omega}\mathbf{\cdot}\mathbf{V})/\mathbf{V}^2\), входящий в выражение \eqref{eq:a6}, есть проекция завихренности на вектор \(\mathbf{a}\), то есть, что \(\mathbf{\Omega}_{\mathbf{a}} =\mathbf{a}(\mathbf{\Omega}\mathbf{\cdot}\mathbf{a})/\mathbf{a}^2\). Любая замкнутая линия \(\mathrm{\gamma}\), по которой производится интегрирование в новом инварианте, лежит на изоэнтропийной поверхности. Следовательно, множитель \(\mathrm{\lambda}_0 (k-1)\mathrm{\sigma}\) в \eqref{eq:a6} остается постоянным вдоль \(\mathrm{\gamma}\) и одинаков для всех \(\mathrm{\gamma}\), лежащих на одной изоэнтропийной поверхности. Поэтому выражение для нового инварианта можно упростить:
\begin{equation}\label{eq:a7}
I_2=\int_{\mathrm{\gamma}}\frac{p}{|\mathbf{\Omega}_{\mathbf{a}}|}\mathrm{d}l.
\end{equation}

В незакрученном осесимметричном течении любая линия \(\mathrm{\gamma}\) совпадает с вихревой линией, то есть \(\mathbf{\Omega}_{\mathbf{a}}=\mathbf{\Omega}\), а отношение \({p}/{|\mathbf{\Omega}|}\) не меняется вдоль \(\mathrm{\gamma}\), и интеграл \eqref{eq:a7} равен \({2\mathrm{\pi}rp}/{|\mathbf{\Omega}|}\). Таким образом, в осесимметричном случае инвариант \eqref{eq:a7} совпадает с неинтегральным инвариантом Крокко \(I_K=\mathrm{\Omega}/(pr)\), а в общем пространственном случае обобщает его.

Заключение

Течение за отошедшим головным скачком исследовано с использованием полных (без каких-либо упрощений) уравнений Эйлера. Ранее было известно, что векторные линии \(\mathbf{a}=\mathbf{V}\times\boldsymbol{\nabla}\mathrm{\sigma}\) замкнуты, лежат на изоэнтропийных поверхностях и один раз опоясывают их. (Такие линии обозначены символом \(\mathrm{\gamma}\).) В данной работе обнаружено, что криволинейный интеграл \eqref{eq:a7} сохраняет свое значение для всех линий \(\mathrm{\gamma}\), лежащих на одной изоэнтропийной поверхности. Это второй из известных интегральных инвариантов, пригодных для проверки численных расчетов несимметричных течений за скачком.

About the authors

Grigory Sizykh

Author for correspondence.
Email: o1o2o3@yandex.ru
ORCID iD: 0000-0001-5821-8596
SPIN-code: 5348-6492
Scopus Author ID: 6508163390
ResearcherId: ABI-3162-2020
http://www.mathnet.ru/person112378

Cand. Phys. & Math. Sci; Associate Professor; Dept. of Applied Mathematics

References

Supplementary files