Construction of a self-similar solution to the system of gas dynamics equations describing the outflow of polytropic gas into vacuum from an inclined wall in the inconsistent case

Full Text

Введение

Рассмотрим политропный газ, покоящийся в клиновидной области $BOF$ плоскости $xOy$ (выделено серым цветом и отмечено цифрой 0 на рис. 1, a), образованной двумя непроницаемыми стенками $OB$ и $OF$. Полубесконечная вертикальная стенка $OB$ задается уравнением $x=0$ (при $y⩾0$). Бесконечная косая стенка $OF$ задается уравнением $y=x\mathrm{tg}\alpha$. Линия $BF$ не является непроницаемой стенкой, поэтому на ней никакие условия на газ не накладываются. Часть плоскости $xOy$ слева от вертикальной стенки $OB$ и выше косой стенки $OF$ — вакуум (отмечено цифрой 3 на рис. 1, a).

Газодинамические параметры в покоящемся газе (значения функций $c(t,x,y)$, $u(t,x,y)$ и $v(t,x,y)$) следующие:
$$c=1,u=0,v=0.$$

В момент времени $t=0$ вертикальная стенка $OB$ убирается, после чего начинается истечение газа в вакуум вдоль косой стенки $OF$.

На рис. 1, b приведена конфигурация течения газа в момент времени $t>0$. В области $ABF$ находится покоящийся однородный газ, который отделен звуковой характеристикой $AB$ от области центрированной волны Римана — $ABCD$, помеченной на рис. 1, b цифрой 1. Звуковая характеристика $AB$ является вертикальной прямой, уравнение которой в координатах $t$, $x$, $y$ имеет следующий вид [1, 2]:
$$x=t.$$

Звуковая характеристика $AB$ двигается по покоящемуся газу слева направо со скоростью, равной 1. Центрированная волна примыкает к вакууму через свободную границу, являющуюся вертикальной прямой $CD$. Движение границы $CD$ в вакуум описывается по закону [1, 2]
$$x=-\frac{2t}{\gamma -1}.$$

Скорость движения границы $CD$ равна $-2/(\gamma -1)$. Значения газодинамических параметров течения в области центрированной волны задаются следующими формулами [1, 2]:
$$c=\frac{\gamma -1}{\gamma +1}\frac{x}{t}+\frac{2}{\gamma +1},u=\frac{2}{(\gamma +1)}\frac{x}{t}-\frac{2}{(\gamma +1)},v=0.$$ 

В области $ADE$, помеченной на рис. 1, b цифрой 2, находится двойная волна — искомое двумерное течение. Это течение отделено от центрированной волны звуковой характеристикой $C^{+}$ — линия $AD$. Область двойной волны примыкает к вакууму через свободную границу — линия $DE$. Стенка $AE$ является непроницаемой, поэтому на этой стенке выполняется условие непротекания:
$$v{|}_{AE}=u\mathrm{tg}\alpha {|}_{AE}.$$

Рис. 1. Начальная конфигурация в момент $t=0$ (a) и конфигурация потока в момент $t>0$ (b): 0 — область, в которой находится покоящийся газ; 1 — область течения в виде центрированной волны; 2 — область течения в виде двойной волны; 3 — область вакуума
[Figure 1. (a) Initial configuration$(t=0)$; (b) the flow configuration at $t>0$: 0 — the quiescent gas region; 1 — the flow region in the form of a centered wave;
2 — the flow region in the form of a double wave; 3 — the vacuum region]

Закон движения газа в области двойной волны неизвестен, требуется найти параметры течения газа в области двойной волны как решение следующей начально-краевой задачи:
$$\begin{array}{}\text{(1)}& \{\begin{array}{c}{c}_t+u{c}_x+v{c}_y+\varkappa c(u_x+v_y)=0,\\ u_t+u{u}_x+v{u}_y+\frac{c}{\varkappa }{c}_x=0,\\ v_t+u{v}_x+v{v}_y+\frac{c}{\varkappa }{c}_y=0,\\ c{|}_{C^{+}}=c_0,u{|}_{C^{+}}=u_0,v{|}_{C^{+}}=0,\\ v{|}_{y=x\mathrm{tg}\alpha }=u\mathrm{tg}\alpha {|}_{y=x\mathrm{tg}\alpha },\end{array}\end{array}$$
где $c$ — скорость звука в газе, отн.ед.; $u$ — горизонтальная компонента скорости газа, отн.ед.; $v$ — вертикальная компонента скорости газа, отн.ед; $\varkappa =(\gamma -1)/2$; $\gamma$ — показатель политропы газа, отн.ед. Значения $c_0$ и $u_0$ — параметры газа на характеристике $C^{+}$:
$$c_0=\frac{\varkappa }{\varkappa +1}\frac{x}{t}+\frac{1}{\varkappa +1},u_0=\frac{1}{\varkappa +1}\frac{x}{t}-\frac{1}{\varkappa +1},$$
соответствующие параметрам газа в плоском течении.

Традиционно для этой и аналогичных задач с целью уменьшения объема выкладок искомое двумерное течение строится как решение уравнения для функции $\mathrm{\Phi }(t,x_1,x_2)$ — потенциала скорости газа в двумерном случае [1, 2]. То есть, чтобы описать разлет газа в вакуум, в качестве независимых переменных выбираются
$$\begin{array}{}\text{(2)}& t,u_1={\mathrm{\Phi }}_{x_1},u_2={\mathrm{\Phi }}_{x_2}\end{array}$$
и в результате замены (2) первые три уравнения системы (1) сводятся к одному уравнению:
$$\begin{array}{}\text{(3)}& {\mathrm{\Phi }}_{tt}+2\sum _{i=1}^2[{\mathrm{\Phi }}_{x_i}{\mathrm{\Phi }}_{x_{i}t}+\frac{1}{2}\sum _{k=1}^2(1-{\delta }_{ik}){\mathrm{\Phi }}_{x_i}{\mathrm{\Phi }}_{x_k}{\mathrm{\Phi }}_{x_i{x}_k}-\frac{1}{2}(c^2-{\mathrm{\Phi }}_{x_i}^2){\mathrm{\Phi }}_{x_i{x}_i}]=0,\end{array}$$
где
$$c^2=(\gamma -1)[M-{\mathrm{\Phi }}_t-\frac{1}{2}({\mathrm{\Phi }}_{x_1}^2+{\mathrm{\Phi }}_{x_2}^2)],M=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}.$$

Далее с помощью преобразования Лежандра
$$\begin{array}{}\text{(4)}& \mathrm{\Psi }(t,u_1,u_2)=x_1{u}_1+x_2{u}_2-\mathrm{\Phi }(t,x_1,x_2)+Mt\end{array}$$
делается переход к новой неизвестной функции $\mathrm{\Psi }(t,u_1,u_2)$. Уравнение для $\mathrm{\Psi }(t,u_1,u_2)$ строится из (3) с учетом (4) и введенных обозначений $u=u_1$, $v=u_2$:
$$\begin{array}{}\text{(5)}& {\mathrm{\Psi }}_{tt}({\mathrm{\Psi }}_{uu}{\mathrm{\Psi }}_{vv}-{\mathrm{\Psi }}_{uv}^2)+[c^2-({\mathrm{\Psi }}_{ut}-u{)}^2]{\mathrm{\Psi }}_{vv}+2({\mathrm{\Psi }}_{ut}-u)({\mathrm{\Psi }}_{vt}-v){\mathrm{\Psi }}_{uv}+[c^2-({\mathrm{\Psi }}_{vt}-v{)}^2]{\mathrm{\Psi }}_{uu}=0,\end{array}$$
где $c^2=2\varkappa ({\mathrm{\Psi }}_t-(u^2+v^2)/2)$.

Начальное и граничное условия задачи (1), записанные для функции $\mathrm{\Psi }(t,u,v)$, имеют вид
$$\begin{array}{}\text{(6)}& \mathrm{\Psi }(t,u,0)=t\{\frac{u^2}{2}+\frac{\varkappa }{2}[c_0+\frac{u}{\varkappa }{]}^2\}\end{array}$$
и
$$\begin{array}{}\text{(7)}& {\mathrm{\Psi }}_v(t,u,u\mathrm{tg}\alpha )={\mathrm{\Psi }}_u(t,u,u\mathrm{tg}\alpha )\mathrm{tg}\alpha ,\end{array}$$
где$c_0=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}>0$.

Решение задачи (5)–(7) строится в пространстве годографа $u$, $v$ в виде ряда по степеням $v$. Для нахождения значений газодинамических параметров $c$, $u$, $v$ в пространстве физических переменных $t$, $x$, $y$ необходимо выполнять обратное преобразование Лежандра (4), а затем обратную замену (2).

В работе [3] получено частное точное решение задачи (5)–(7), описывающее двумерное течение газа при выполнении конкретного соотношения между показателем политропы $\gamma$ газа и тангенсом угла $\alpha$ наклона косой стенки:
$$\begin{array}{}\text{(8)}& {\mathrm{tg}}^2\alpha =\frac{\gamma +1}{3-\gamma }.\end{array}$$

Здесь и далее случай выполнения равенства (8) будем называть случаем согласованного течения, и наоборот, когда соотношение (8) не выполняется, тогда рассматриваемый случай двумерного течения будет несогласованным.

Решение (5)–(7) в согласованном случае (8) в координатах $u$, $v$ имеет следующий вид [3]:
$$\begin{array}{}\text{(9)}& c(u,v)=1+\varkappa u+\varkappa \sqrt{\beta }v,\end{array}$$
где $\beta =(\varkappa +1)/(1-\varkappa )$. Видно, что решение (9) линейно относительно независимой переменной $v$, то есть коэффициенты $c_i(u)={\partial }^{i}c/\partial v^i{|}_{v=0}$ ряда
$$c(u,v)=\sum _{i=0}^{\mathrm{\infty }}{c}_i(u)\frac{v^i}{i!}$$
при $i⩾2$ равны $0$.

С. П. Баутиным и С. Л. Дерябиным [2, c. 196–214] в пространстве специальных независимых переменных (решение строилось для функции $\mathrm{\Psi }$ в пространстве годографа) рассмотрена задача об истечении газа в вакуум при произвольном значении угла $\alpha$:$0<\alpha <\pi /2$ — наклона косой стенки, не связанным со значением $\gamma$ (несогласованный случай). Доказано существование и единственность локально-аналитических решений соответствующих начально-краевых задач, построено решение (5)–(7) для функции $\mathrm{\Psi }(t,u,v)$ в виде ряда по степеням $v$. Значения коэффициентов ряда для функции $\mathrm{\Psi }(t,u,v)$ в явном виде не были найдены, значения газодинамических параметров $c$, $u$, $v$ течения газа в пространстве физических переменных $t$, $x$, $y$ не построены.

В цикле работ [4–6] автора настоящей статьи решение исходной системы уравнений газовой динамики (СУГД) рассматриваемой начально-краевой задачи (1) строится для компонент вектора $\mathbf{U}=(c,u,v{)}^{\mathrm{\top }}$ в виде бесконечного ряда по степеням $\vartheta$, где $\vartheta$ — известная функция независимых физических автомодельных переменных $\xi =x/t$, $\eta =y/t$. В работе [4] построено транспортное уравнение для коэффициента $c_1(\xi )$ и найдено решение СУГД в согласованном случае. В работе [5] решение В. А. Сучкова применено к описанию двух способов газодинамического воздействия на специальный призматический объем. В работе [6] исходная начально-краевая задача (1) в результате двух невырожденных замен сводится к характеристической задаче Коши стандартного вида. Для этой новой начально-краевой задачи доказана теорема существования и единственности решения СУГД в виде сходящихся бесконечных рядов. Описан алгоритм построения коэффициентов ряда.

В настоящей работе впервые построено аналитическое решение транспортного уравнения для коэффициента $c_1(\xi )$ решения СУГД, описывающего изэнтропическое истечение политропного газа с косой стенки в общем несогласованном случае в пространстве физических автомодельных переменных $\xi =x/t$, $\eta =y/t$. Для газа с показателем политропы $\gamma =5/3$ (водород) решение транспортного уравнения для коэффициента $c_1(\xi )$ также впервые построено в явном виде. Построенное решение начально-краевой задачи применено к описанию сильного сжатия газа в объеме с непроницаемыми стенками, представляющего собой в поперечном сечении правильный треугольник ($\alpha =\pi /3$).

1. Начально-краевая задача Коши стандартного вида для описания истечения газа в вакуум на косой стенке

В результате невырожденной замены $\xi =x/t$, $\eta =y/t$ исходная начально-краевая задача (1) для СУГД, решение которой при $t>0$ описывает разлет политропного газа в вакуум на косой стенке, записывается в векторно-матричном виде (см. работу [4]) для вектора $\mathbf{U}=(u_1,u_2,u_3{)}^{\mathrm{\top }}=(c,u,v{)}^{\mathrm{\top }}$:
$$\begin{array}{}\text{(10)}& \{\begin{array}{l}A{\mathbf{U}}_{\xi }+B{\mathbf{U}}_{\eta }=\mathbf{0},\\ \mathbf{U}{|}_{C^{+}}={\mathbf{U}}_0,\\ u_3{|}_{\eta =\xi \mathrm{tg}\alpha }=u_2\mathrm{tg}\alpha {|}_{\eta =\xi \mathrm{tg}\alpha }.\end{array}\end{array}$$

Начальное условие задачи (10) задается значением вектора $\mathbf{U}$ на характеристике $C^{+}$ — ${\mathbf{U}}_0=(c_0,u_0,0{)}^{\mathrm{\top }}$. Звуковая характеристика $C^{+}$ задается уравнением $\eta =f(\xi )$, где $f(\xi )$ — неизвестная функция независимой переменной $\xi$. Краевое условие (граничное условие) задачи (10) задается в виде условия непротекания газа на косой стенке, определяемой уравнением $\eta =\xi \mathrm{tg}\alpha$. 

Матрицы $A$ и $B$ задачи (10) следующие:
$$A=\left(\begin{array}{ccc}u-\xi & \varkappa c& 0\\ c/\varkappa & u-\xi & 0\\ 0& 0& u-\xi \end{array}\right),B=\left(\begin{array}{ccc}v-\eta & 0& \varkappa c\\ 0& v-\eta & 0\\ c/\varkappa & 0& v-\eta \end{array}\right).$$

Начально-краевые условия, коэффициенты матриц $A$ и $B$ задачи (10) предполагаются аналитическими функциями.

Для приведения задачи (10) к характеристической задаче Коши стандартного вида сделаем замену переменных
$$\begin{array}{}\text{(11)}& \vartheta =\eta -f(\xi ),{\xi }^{\prime }=\xi ,\end{array}$$
где линия $\vartheta =0$ задает звуковую характеристику $C^{+}$. Якобиан $J$ замены (11) следующий:
$$J=\left|\begin{array}{cc}{\vartheta }_{\eta }& {\vartheta }_{\xi }\\ {\xi }_{\eta }^{\prime }& {\xi }_{\xi }^{\prime }\end{array}\right|=\left|\begin{array}{cc}1& -f^{\prime }\\ 0& 1\end{array}\right|=1-0\cdot f^{\prime }=1,$$
т.е. при выполнении условия $|f^{\prime }|<\mathrm{\infty }$ замена (11) невырожденная.

В результате замены (11) задача (10) преобразуется к виду
$$\begin{array}{}\text{(12)}& \{\begin{array}{l}[B-f^{\prime }(\xi )A]{\mathbf{U}}_{\vartheta }+A{\mathbf{U}}_{\xi }=\mathbf{0},\\ \mathbf{U}{|}_{\vartheta =0}={\mathbf{U}}_0,\\ u_3{|}_{\vartheta =\xi \mathrm{tg}\alpha -f(\xi )}=u_2\mathrm{tg}\alpha {|}_{\vartheta =\xi \mathrm{tg}\alpha -f(\xi )}.\end{array}\end{array}$$

Задача (12) эквивалентна задаче (10). Новая замена переменных $\xi$, $\vartheta$ на переменные $\vartheta$, $\zeta$ по формулам
$$\begin{array}{}\text{(13)}& {\vartheta }^{\prime }=\vartheta ,\zeta =\vartheta +f(\xi )-\xi \mathrm{tg}\alpha \end{array}$$
приводит задачу (12) к следующему виду:
$$\begin{array}{}\text{(14)}& \{\begin{array}{l}[B-f^{\prime }(\xi )A]{\mathbf{U}}_{\vartheta }+[B-\mathrm{tg}\alpha A]{\mathbf{U}}_{\zeta }=\mathbf{0},\\ \mathbf{U}{|}_{\vartheta =0}={\mathbf{U}}_0,\\ u_3{|}_{\zeta =0}=u_2\mathrm{tg}\alpha {|}_{\zeta =0}.\end{array}\end{array}$$

Из соотношения (13), неявно задающего функцию $\xi$ от $(\zeta -\vartheta )$, с учетом неравенства (16) однозначно определяется функция
$$\begin{array}{}\text{(15)}& \xi =\phi (\zeta -\vartheta ).\end{array}$$

Далее для простоты записи будут сохраняться обозначения $f(\xi )$, $f^{\prime }(\xi )$, $c_0(\xi )$, естественно, с учетом наличия связи (15).

Исходное уравнение косой стенки $y=x\mathrm{tg}\alpha$ в координатах $\vartheta$, $\zeta$ имеет вид $\zeta =0$, т.е. косая стенка берется за новую координатную ось. Якобиан замены переменных (13) следующий:
$$J_2=\left|\begin{array}{cc}{\displaystyle \frac{\partial \zeta }{\partial \vartheta }}& {\displaystyle \frac{\partial \zeta }{\partial \xi }}\\ {\displaystyle \frac{\partial {\vartheta }^{\prime }}{\partial \vartheta }}& {\displaystyle \frac{\partial {\vartheta }^{\prime }}{\partial \xi }}\end{array}\right|=\left|\begin{array}{cc}1& f^{\prime }(\xi )-\mathrm{tg}\alpha \\ 1& 0\end{array}\right|=\mathrm{tg}\alpha -f^{\prime }\left(\xi \right).$$

Чтобы замена (13) была невырожденной, необходимо выполнение неравенства
$$\begin{array}{}\text{(16)}& J_2=\mathrm{tg}\alpha -f^{\prime }\left(\xi \right)\ne 0,\end{array}$$
т.е. наклон косой стенки не равен наклону звуковой характеристики, разделяющей области центрированной и двойной волны. Неравенство (16) доказывается в работе [6].

Таким образом, начально-краевая задача (14) для вектора $\mathbf{U}$ записана в координатах $\vartheta$, $\zeta$. Начально-граничные условия (14) определены на прямых, задаваемых уравнениями $\vartheta =0$ и $\zeta =0$. Далее для построения локально-аналитического решения в окрестности точки $(\zeta =0,\vartheta =0)$ задачу (14) необходимо привести к форме характеристической задачи Коши стандартного вида. Для этого СУГД из задачи (14) слева умножается на матрицу $T_1$, а вектор $\mathbf{U}$ заменяется новым вектором:
$$\mathbf{W}=T_2^{-1}\mathbf{U}=(w_1,w_2,w_3{)}^{\mathrm{\top }}=(c+v{\displaystyle \frac{\varkappa }{c_0}}(c_0{f}^{\prime }-f),u+v{f}^{\prime },-{\displaystyle \frac{\varkappa }{c_0}}v{)}^{\mathrm{\top }}.$$

Невырожденные матрицы $T_1$, $T_2$, элементы которых есть аналитические функции независимой переменной $\xi$, имеют вид
$$\begin{array}{l}{T}_1=\left(\begin{array}{ccc}0& 1& 0\\ 1& 0& 0\\ c_0{f}^{\prime }-f& \varkappa c_0{f}^{\prime }& -\varkappa c_0\end{array}\right),T_2=\left(\begin{array}{ccc}1& 0& c_0{f}^{\prime }-f\\ 0& 1& c_0{f}^{\prime }/\varkappa \\ 0& 0& -c_0/\varkappa \end{array}\right).\end{array}$$

Начальное условие для вектора $\mathbf{U}$ заменяется начальным условием для вектора $\mathbf{W}$:
$${\mathbf{W}|}_{\vartheta =0}={\mathbf{W}}_0={\left(\begin{array}{l}{w}_{1,0},w_{2,0},w_{3,0}\end{array}\right)}^{\mathrm{\top }}={\left(\begin{array}{l}{c}_0,u_0,0\end{array}\right)}^{\mathrm{\top }}.$$

Краевое условие задачи (14)
$$u_3{|}_{\zeta =0}=u_2\mathrm{tg}\alpha {|}_{\zeta =0}$$
переписывается через компоненты вектора $\mathbf{W}$ и преобразуется к виду
$$w_3{|}_{\zeta =0}=g(\vartheta )w_2{|}_{\zeta =0},$$
где $g(\vartheta )=-[\frac{\varkappa }{c_0(\phi (\zeta -\vartheta ))}\frac{\mathrm{tg}\alpha }{[1+f_{\xi }^{\prime }(\phi (\zeta -\vartheta ))\mathrm{tg}\alpha ]}{]}_{\zeta =0}$.

Окончательно начально-краевая задача (14) имеет вид
$$\begin{array}{}\text{(17)}& \{\begin{array}{l}{T}_1[B-f^{\prime }A]T_2{\mathbf{W}}_{\vartheta }+T_1[B-\mathrm{tg}\alpha A]T_2{\mathbf{W}}_{\zeta }+\\ +T_1[B-\mathrm{tg}\alpha A]\frac{\partial T_2}{\partial \xi }\frac{\mathbf{W}}{f^{\prime }-\mathrm{tg}\alpha }=\mathbf{0},\\ \mathbf{W}{|}_{\vartheta =0}={\mathbf{W}}_0,\\ w_3{|}_{\zeta =0}=g(\vartheta )w_2{|}_{\zeta =0}.\end{array}\end{array}$$

Задача (17) эквивалентна задаче (12), а та, в свою очередь, задаче (10), так как замены (11) и (13) невырожденные. Для того чтобы задача (17) была характеристической задачей Коши стандартного вида, необходимо выполнение следующего условия при $\vartheta =0$ (см. работу [7]):
$$det(T_1[B-f^{\prime }(\xi )A]T_2){|}_{\genfrac{}{}{0ex}{}{\mathbf{W}={\mathbf{W}}_0}{\vartheta =0}}\equiv 0.$$

Известно, что определитель произведения матриц равен произведению определителей каждой из матриц, отсюда
$$det{T}_1{|}_{\genfrac{}{}{0ex}{}{\mathbf{W}={\mathbf{W}}_0}{\vartheta =0}}\cdot det[B-f^{\prime }(\xi )A]{|}_{\genfrac{}{}{0ex}{}{\mathbf{W}={\mathbf{W}}_0}{\vartheta =0}}\cdot det{T}_2{|}_{\genfrac{}{}{0ex}{}{\mathbf{W}={\mathbf{W}}_0}{\vartheta =0}}\equiv 0,$$ 
т. к. матрицы $T_1$ и $T_2$ невырожденные, их определители отличны от нуля. В результате получаем матричное уравнение
$$\begin{array}{}\text{(18)}& det[B-f^{\prime }(\xi )A]{|}_{\genfrac{}{}{0ex}{}{\mathbf{U}={\mathbf{U}}_0}{\vartheta =0}}\equiv 0,\end{array}$$
решение которого дает выражение для функции $f(\xi )$ в явном виде:
$$f(\xi )=\pm \{\begin{array}{lll}{c}_0\sqrt{{\mathrm{tg}}^2\alpha -2\ln c_0},& \varkappa =1& (\gamma =3),\\ c_0\sqrt{\beta +c_0^{\frac{1-\varkappa }{\varkappa }}({\mathrm{tg}}^2\alpha -\beta )},& \varkappa \ne 1& (\gamma \ne 3).\end{array}$$

В работе [6] для задачи (17) доказывается 

Теорема. Поставленная задача (17) при найденной функции $f(\xi )$ является характеристической задачей Коши стандартного вида и поэтому у нее в некоторой окрестности точки $(\zeta =0,\vartheta =0)$ существует единственное локально-аналитическое решение, представимое в виде сходящегося ряда
$$\mathbf{W}(\zeta ,\vartheta )=\sum _{k=0}^{\mathrm{\infty }}{\mathbf{W}}_k(\zeta )\frac{{\vartheta }^k}{k!},{\mathbf{W}}_k(\zeta )=\frac{{\partial }^k\mathbf{W}}{\partial {\vartheta }^k}{|}_{\vartheta =0}.$$

Алгоритм построения решения задачи (17) в пространстве специальных переменных $\vartheta$, $\zeta$ подробно описан в работе [6]. Для построения решения задачи (17) в пространстве физических автомодельных переменных $\xi$, $\eta$ необходимо осуществлять обратное преобразование согласно заменам (11) и (13), что в общем (несогласованном) случае выполнить в явном виде затруднительно. Поэтому искомое двумерное течение будет строиться как решение задачи (12) в пространстве переменных $\xi$, $\vartheta$ в виде сходящегося ряда по степеням $\vartheta$:
$$\begin{array}{}\text{(19)}& \mathbf{U}(\xi ,\vartheta )=\sum _{k=0}^{\mathrm{\infty }}{\mathbf{U}}_k(\xi )\frac{{\vartheta }^k}{k!},{\mathbf{U}}_k(\xi )=\frac{{\partial }^k\mathbf{U}}{\partial {\vartheta }^k}{|}_{\vartheta =0},\end{array}$$
так как задачи (17) и (12) эквивалентны. Начальное условие для транспортного уравнения задачи (12) будет строиться из граничного условия задачи (17).

2. Построение транспортного уравнения и уравнения начальных условий для коэффициента ${\mathit{c}}_{\mathbf{1}}\mathbf{(}\mathit{\xi }\mathbf{)}$ рассматриваемой начально-краевой задачи

При подстановке в СУГД задачи (12) значения $\vartheta =0$ получаем соотношения между функциями $c_1(\xi )$, $u_1(\xi )$ и $v_1(\xi )$:
$$\begin{array}{}\text{(20)}& u_1={\displaystyle \frac{c_1}{\varkappa }}{\displaystyle \frac{c_0{f}^{\prime }}{c_0{f}^{\prime }-f}},v_1=-{\displaystyle \frac{c_1}{\varkappa }}{\displaystyle \frac{c_0}{c_0{f}^{\prime }-f}}.\end{array}$$

Для построения транспортного уравнения (дифференциальное уравнение для нахождения коэффициента $c_1(\xi )$) необходимо систему (12) продифференцировать по переменной $\vartheta$ и подставить значение $\vartheta =0$. При этом в соответствии с (19) будут справедливы следующие соотношения:
$$\begin{array}{lllll}c{|}_{\vartheta =0}=c_0,& c_{\xi }{|}_{\vartheta =0}=c_0^{\prime },& c_{\vartheta }{|}_{\vartheta =0}=c_1,& c_{\xi \vartheta }{|}_{\vartheta =0}=c_1^{\prime },& c_{\vartheta \vartheta }{|}_{\vartheta =0}=c_2,\\ u{|}_{\vartheta =0}=u_0,& u_{\xi }{|}_{\vartheta =0}=u_0^{\prime },& u_{\vartheta }{|}_{\vartheta =0}=u_1,& u_{\xi \vartheta }{|}_{\vartheta =0}=u_1^{\prime },& u_{\vartheta \vartheta }{|}_{\vartheta =0}=u_2,\\ v{|}_{\vartheta =0}=0,& v_{\xi }{|}_{\vartheta =0}=0,& v_{\vartheta }{|}_{\vartheta =0}=v_1,& v_{\xi \vartheta }{|}_{\vartheta =0}=v_1^{\prime },& v_{\vartheta \vartheta }{|}_{\vartheta =0}=v_2,\end{array}$$
в результате система (12) принимает следующий вид:
$$\begin{array}{}\text{(21)}& \{\begin{array}{l}{u}_1{c}_0^{\prime }-c_0{c}_1^{\prime }+\varkappa c_1{u}_0^{\prime }+\varkappa c_0{u}_1^{\prime }+(v_1-1-u_1{f}^{\prime })c_1+(c_0{f}^{\prime }-f)c_2-\varkappa f^{\prime }{c}_1{u}_1-\varkappa f^{\prime }{c}_0{u}_2+\varkappa c_1{v}_1+\varkappa c_0{v}_2=0,\\ c_1{c}_0^{\prime }+c_0{c}_1^{\prime }+\varkappa u_1{u}_0^{\prime }-\varkappa c_0{u}_1^{\prime }-c_1^2{f}^{\prime }-c_0{f}^{\prime }{c}_2+\varkappa (c_0{f}^{\prime }-f)u_2+\varkappa (v_1-1-u_1{f}^{\prime })u_1=0,\\ -\varkappa c_0{v}_1^{\prime }+c_1^2+c_0{c}_2+\varkappa (v_1-1-u_1{f}^{\prime })v_1+\varkappa (c_0{f}^{\prime }-f)v_2=0.\end{array}\end{array}$$

Из второго и третьего уравнения системы (21) выразим $u_2(\xi )$ и $v_2(\xi )$:
$$\begin{array}{l}{u}_2=-{\displaystyle \frac{c_1{c}_0^{\prime }+c_0{c}_1^{\prime }+\varkappa u_1{u}_0^{\prime }-\varkappa c_0{u}_1^{\prime }-c_1^2{f}^{\prime }-c_0{f}^{\prime }{c}_2+\varkappa (v_1-1-u_1{f}^{\prime })u_1}{\varkappa (c_0{f}^{\prime }-f)}},\\ v_2=-{\displaystyle \frac{c_1^2+c_0{c}_2+\varkappa (v_1-1-u_1{f}^{\prime })v_1-\varkappa c_0{v}_1^{\prime }}{\varkappa (c_0{f}^{\prime }-f)}}.\end{array}$$

Подставим выражения для $u_2(\xi )$ и $v_2(\xi )$ в первое уравнение системы (21), в нем сгруппируем слагаемые, содержащие коэффициент $c_2(\xi )$, в результате получим следующее уравнение:
$$\begin{array}{c}{u}_1{c}_0^{\prime }+\varkappa c_1{u}_0^{\prime }-c_0{c}_1^{\prime }+\varkappa c_0{u}_1^{\prime }+[v_1-1-u_1{f}^{\prime }]c_1-\varkappa f^{\prime }{c}_1{u}_1+\varkappa c_1{v}_1-\varkappa c_0{f}^{\prime }[{\displaystyle \frac{c_1^2}{\varkappa }}{\displaystyle \frac{f^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{(v_1-1-u_1{f}^{\prime })u_1}{c_0{f}^{\prime }-f}}+{\displaystyle \frac{c_0{u}_1^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{c_0}{\varkappa }}{\displaystyle \frac{c_1^{\prime }}{c_0{f}^{\prime }-f}}]-\hfill \\ \hfill -\varkappa c_0{f}^{\prime }[-{\displaystyle \frac{u_1{u}_0^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{c_0^{\prime }}{\varkappa }}{\displaystyle \frac{c_1}{c_0{f}^{\prime }-f}}]+\varkappa c_0[{\displaystyle \frac{c_0{v}_1^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{c_1^2}{\varkappa }}{\displaystyle \frac{1}{c_0{f}^{\prime }-f}}]-\varkappa c_0{v}_1{\displaystyle \frac{v_1-1-u_1{f}^{\prime }}{c_0{f}^{\prime }-f}}+c_2{\displaystyle \frac{(c_0{f}^{\prime }-f{)}^2-(c_0{f}^{\prime }{)}^2-c_0^2}{c_0{f}^{\prime }-f}}=0.\end{array}$$

Множитель при коэффициенте $c_2(\xi )$ тождественно равен нулю, так как справедливо соотношение (18). В результате получаем дифференциальное уравнение относительно коэффициентов $c_1(\xi )$, $u_1(\xi )$ и $v_1(\xi )$:
$$\begin{array}{}\text{(22)}& \begin{array}{c}{u}_1{c}_0^{\prime }+\varkappa c_1{u}_0^{\prime }-c_0{c}_1^{\prime }+\varkappa c_0{u}_1^{\prime }+[v_1-1-u_1{f}^{\prime }]c_1-\varkappa f^{\prime }{c}_1{u}_1+\varkappa c_1{v}_1-\varkappa c_0{f}^{\prime }[{\displaystyle \frac{c_1^2}{\varkappa }}{\displaystyle \frac{f^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{(v_1-1-u_1{f}^{\prime })u_1}{c_0{f}^{\prime }-f}}+{\displaystyle \frac{c_0{u}_1^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{c_0}{\varkappa }}{\displaystyle \frac{c_1^{\prime }}{c_0{f}^{\prime }-f}}]-\hfill \\ \hfill -\varkappa c_0{f}^{\prime }[-{\displaystyle \frac{u_1{u}_0^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{c_0^{\prime }}{\varkappa }}{\displaystyle \frac{c_1}{c_0{f}^{\prime }-f}}]+\varkappa c_0[{\displaystyle \frac{c_0{v}_1^{\prime }}{c_0{f}^{\prime }-f}}-{\displaystyle \frac{c_1^2}{\varkappa }}{\displaystyle \frac{1}{c_0{f}^{\prime }-f}}]-\varkappa c_0{v}_1{\displaystyle \frac{v_1-1-u_1{f}^{\prime }}{c_0{f}^{\prime }-f}}=0.\end{array}\end{array}$$

Используя соотношения (20), можно свести уравнение (22) к дифференциальному уравнению относительно неизвестной функции $c_1(\xi )$ и построить таким образом транспортное уравнение. После упрощений в результате элементарных преобразований исходное уравнение (22) записывается в нормализованном виде:
$$\begin{array}{}\text{(23)}& c_1^{\prime }-[\frac{c_0^2}{f^2}-1-\frac{4}{\varkappa +1}\frac{c_0^2}{c_0^2+f^2}]\frac{c_1}{2c_0}-\frac{\varkappa +1}{\varkappa }\frac{(c_0^2+f^2{)}^2}{4c_0^2{f}^3}{c}_1^2=0.\end{array}$$

Рассмотрим общий несогласованный случай, когда $\varkappa \ne 1$, $\gamma \ne 3$, тогда функцию $f(\xi )$ можно записать в виде произведения
$$\begin{array}{}\text{(24)}& f(\xi )=c_0(\xi )R(\xi ),\end{array}$$
где $R(\xi )=\sqrt{\beta +c_0^{(1-\varkappa )/\varkappa }({\mathrm{tg}}^2\alpha -\beta )}$. Тогда уравнение (23) можно упростить до следующего вида:
$$\begin{array}{}\text{(25)}& c_1^{\prime }-[\frac{1}{R^2}-1-\frac{4}{\varkappa +1}\frac{1}{R^2+1}]\frac{c_1}{2c_0}-\frac{\varkappa +1}{\varkappa }\frac{(1+R^2{)}^2}{4R^3}\frac{c_1^2}{c_0}=0.\end{array}$$

Построим начальное условие для транспортного уравнения (25). Для этого краевое условие из задачи (17) продифференцируем по $\vartheta$ и в получившееся соотношение 
$$w_{3\vartheta }{|}_{\zeta =0}=[g^{\prime }(\vartheta )w_2+g(\vartheta )w_{2\vartheta }{]}_{\zeta =0}$$
подставим значение $\vartheta =0$:
$$w_{3,1}{|}_{\zeta =0}=g_{\vartheta }^{\prime }(0)w_{2,0}{|}_{\zeta =0}+g(0)w_{2,1}{|}_{\zeta =0}.$$
Здесь $g(0)$ и $g^{\prime }(0)$ — значения функции $g(\vartheta )$ и ее производной по $\vartheta$ при $\vartheta =0$. Так как $w_{2,0}{|}_{\zeta =0}=u_0{|}_{\zeta =0}=0$, выполняется
$$\begin{array}{}\text{(26)}& w_{3,1}{|}_{\zeta =0}=g(0)w_{2,1}{|}_{\zeta =0}.\end{array}$$

По определению $w_{3,1}{|}_{\zeta =0}=-\frac{\varkappa }{c_0}{v}_1{|}_{\zeta =0}$. Для определения начального условия для уравнения (25) необходимо найти явный вид для выражения $w_{2,1}{|}_{\zeta =0}$. Для этого в систему из задачи (17) подставим значение $\vartheta =0$. С учетом вида функций $w_1$, $w_2$, $w_3$ и $w_{1\zeta }^{\prime }$, $w_{2\zeta }^{\prime }$, $w_{3\zeta }^{\prime }$ и $w_{1\vartheta }^{\prime }$, $w_{2\vartheta }^{\prime }$, $w_{3\vartheta }^{\prime }$ при $\vartheta =0$ система (17) принимает вид
$$\begin{array}{}\text{(27)}& \{\begin{array}{l}-{\displaystyle \frac{c_0}{\varkappa }}{f}^{\prime }{w}_{1,1}+(c_0{f}^{\prime }-f)w_{2,1}={\displaystyle \frac{u_0^{\prime }f}{f^{\prime }-\mathrm{tg}\alpha }},\\ (c_0{f}^{\prime }-f)w_{1,1}-\varkappa c_0{f}^{\prime }{w}_{2,1}={\displaystyle \frac{c_0^{\prime }f}{f^{\prime }-\mathrm{tg}\alpha }}.\end{array}\end{array}$$

Система (27) представляет собой систему линейных алгебраических уравнений, решив которую, найдем функции $w_{1,1}$, $w_{2,1}$:
$$w_{1,1}=-\frac{c_0^{\prime }}{f^{\prime }-\mathrm{tg}\alpha },w_{2,1}=-\frac{u_0^{\prime }}{f^{\prime }-\mathrm{tg}\alpha }.$$

Подставим в (26) выражения для $w_{3,1}$, $w_{2,1}$, $g(0)$:
$$v_1{|}_{\zeta =0}=-\frac{u_0^{\prime }}{f^{\prime }-\mathrm{tg}\alpha }\frac{\mathrm{tg}\alpha }{1+f^{\prime }\mathrm{tg}\alpha }{|}_{\zeta =0}.$$

С учетом связи (20) между функциями $c_1(\xi )$ и $v_1(\xi )$ и уравнением характеристики (24) получим начальное условие для $c_1(\xi )$ в случае, когда $\varkappa \ne 1$, $\gamma \ne 3$:
$$c_1{|}_{\zeta =0}=-{\displaystyle \frac{2c_0^{\prime }R(1+R^2)}{R^2-2R\mathrm{tg}\alpha -1}}{\displaystyle \frac{\mathrm{tg}\alpha }{2R+\mathrm{tg}\alpha (R^2-1)}}{|}_{\zeta =0}.$$

При $\zeta =0$ и $\vartheta =0$, т. е. когда $\xi =1$, значение $c_0(1)=1$, а значение $R(\xi =1)=\mathrm{tg}\alpha$; отсюда
$$\begin{array}{}\text{(28)}& c_1{|}_{\xi =1}=\frac{2c_0^{\prime }\mathrm{tg}\alpha }{{\mathrm{tg}}^2\alpha +1}.\end{array}$$

Таким образом, с учетом (25) и (28) задача Коши для транспортного уравнения рассматриваемой СУГД в общем несогласованном случае при $\varkappa \ne 1$, $\gamma \ne 3$ будет иметь вид
$$\begin{array}{}\text{(29)}& \{\begin{array}{c}{c}_1^{\prime }-P(\xi )c_1-Q(\xi )c_1^2=0,\\ c_1{|}_{\xi =1}={\displaystyle \frac{2c_0^{\prime }\mathrm{tg}\alpha }{{\mathrm{tg}}^2\alpha +1}},\end{array}\end{array}$$
где 
$$P(\xi )={\displaystyle \frac{1}{2c_0}}[{\displaystyle \frac{1}{R^2}}-1-{\displaystyle \frac{4}{\varkappa +1}}{\displaystyle \frac{1}{R^2+1}}],Q(\xi )={\displaystyle \frac{(1+R^2{)}^2}{4c_0^{\prime }{c}_0{R}^3}}.$$

Первое уравнение задачи (29) — дифференциальное уравнение Бернулли, решение которого можно выписать в квадратурах.

3. Решение транспортного уравнения для коэффициента ${\mathit{c}}_{\mathbf{1}}\mathbf{(}\mathit{\xi }\mathbf{)}$ в несогласованном случае при $\mathit{\varkappa }\mathbf{\ne }\mathbf{1}$, $\mathit{\gamma }\mathbf{\ne }\mathbf{3}$

Запишем функцию $c_1(\xi )$ как произведение двух неизвестных функций $c_1(\xi )=q(\xi )p(\xi )$ и подставим в первое уравнение задачи (29):
$$\begin{array}{}\text{(30)}& p^{\prime }q+(q^{\prime }-P(\xi )q)p-Q(\xi )p^2{q}^2=0.\end{array}$$

Коэффициент при $p(\xi )$ в уравнении (30) зануляется и значение $q(\xi )$ находится из решения дифференциального уравнения
$$\begin{array}{}\text{(31)}& q^{\prime }-P(\xi )q=0.\end{array}$$
Решением (31) будет функция
$$q(\xi )=\exp (\int P(\xi )d\xi ).$$

Найдем значение интеграла, стоящего в показателе экспоненты:
$$\int P(\xi )d\xi ={\displaystyle \frac{1}{2c_0^{\prime }}}\int {\displaystyle \frac{d{c}_0}{c_0(\beta +c_0^{1/\varkappa -1}({\mathrm{tg}}^2\alpha -\beta ))}}-{\displaystyle \frac{1}{2c_0^{\prime }}}\int {\displaystyle \frac{d{c}_0}{c_0}}-{\displaystyle \frac{2}{\varkappa }}\int {\displaystyle \frac{d{c}_0}{c_0(\beta +1+c_0^{1/\varkappa -1}({\mathrm{tg}}^2\alpha -\beta ))}}.$$

Сделаем замену переменных:
$$\begin{array}{}\text{(32)}& z=c_0^{1/\varkappa },c_0=z^{\varkappa },d{c}_0=\varkappa z^{\varkappa -1}dz.\end{array}$$
Здесь и далее предполагается, что $\varkappa \in \mathbb{Q}$. В результате будем иметь
$$\int P(\xi )d\xi ={\displaystyle \frac{\varkappa +1}{2}}\int {\displaystyle \frac{z^{\varkappa -2}dz}{\beta z^{\varkappa -1}+{\mathrm{tg}}^2\alpha -\beta }}-2\int {\displaystyle \frac{z^{\varkappa -2}dz}{(\beta +1)z^{\varkappa -1}+{\mathrm{tg}}^2\alpha -\beta }}-{\displaystyle \frac{1}{2c_0^{\prime }}}\ln c_0=\ln ({\displaystyle \frac{c_0^{(\varkappa -1)/\varkappa }(\beta +1+c_0^{(1-\varkappa )/\varkappa }({\mathrm{tg}}^2\alpha -\beta ))}{c_0\sqrt{\beta +c_0^{(1-\varkappa )/\varkappa }({\mathrm{tg}}^2\alpha -\beta )}}}).$$
В результате
$$\begin{array}{}\text{(33)}& q(\xi )={\displaystyle \frac{c_0^{(\varkappa -1)/\varkappa }(\beta +1+c_0^{(1-\varkappa )/\varkappa }({\mathrm{tg}}^2\alpha -\beta ))}{c_0\sqrt{\beta +c_0^{(1-\varkappa )/\varkappa }({\mathrm{tg}}^2\alpha -\beta )}}}={\displaystyle \frac{c_0^2+f^2}{c_0^{(\varkappa +1)/\varkappa }f}}={\displaystyle \frac{1+R^2}{c_0^{1/\varkappa }R}}.\end{array}$$

Функция $p(\xi )$ является решением уравнения
$${\displaystyle \frac{dp}{p^2}}=\frac{1}{c_0^{\prime }}Q(c_0)q(c_0)d{c}_0,$$
отсюда
$$\begin{array}{}\text{(34)}& \tilde{C}-{\displaystyle \frac{1}{p}}=({\displaystyle \frac{\varkappa +1}{2\varkappa }}{)}^2\int {\displaystyle \frac{(1+R^2{)}^3}{c_0^{(\varkappa +1)/\varkappa }{R}^4}}d{c}_0,\end{array}$$
где $\tilde{C}$ — константа интегрирования. 

Вычислим интеграл, стоящий в правой части выражения (34). После замены (32) интеграл переписывается следующим образом:
$$\begin{array}{}\text{(35)}& \begin{array}{c}\int {\displaystyle \frac{(1+R^2{)}^3}{c_0^{(\varkappa +1)/\varkappa }{R}^4}}d{c}_0=\int {\displaystyle \frac{d{c}_0}{c_0^{(\varkappa +1)/\varkappa }{R}^4}}+3\int \frac{d{c}_0}{c_0^{(\varkappa +1)/\varkappa }{R}^2}+3\int \frac{d{c}_0}{c_0^{(\varkappa +1)/\varkappa }}+\int {\displaystyle \frac{R^{2}d{c}_0}{c_0^{(\varkappa +1)/\varkappa }}}=\hfill \\ \hfill =-\frac{{\mathrm{tg}}^2\alpha -\beta }{c_0}+3\varkappa \int {\displaystyle \frac{z^{\varkappa -3}dz}{(z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta ))}}-\frac{\varkappa (\beta +3)}{c_0^{1/\varkappa }}+\varkappa \int {\displaystyle \frac{dz}{z^{4-2\varkappa }(z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta ){)}^2}}.\end{array}\end{array}$$
Последний интеграл в выражении (35) берется по частям:
$$\begin{array}{}\text{(36)}& \int {\displaystyle \frac{dz}{z^{4-2\varkappa }(z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta ){)}^2}}={\displaystyle \frac{1}{\beta (\varkappa -1)}}\int {\displaystyle \frac{z^{\varkappa -2}d\beta z^{\varkappa -1}}{(z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta ){)}^2}}={\displaystyle \frac{1}{\varkappa +1}}[-{\displaystyle \frac{z^{\varkappa -2}}{z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta )}}+(\varkappa -2)\int {\displaystyle \frac{z^{\varkappa -3}dz}{z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta )}}].\end{array}$$
Выражение (36) подставим в (35):
$$\begin{array}{}\text{(37)}& \int {\displaystyle \frac{(1+R^2{)}^3}{c_0^{(\varkappa +1)/\varkappa }{R}^4}}d{c}_0=-\frac{\varkappa (\beta +3)}{c_0^{1/\varkappa }}-\frac{({\mathrm{tg}}^2\alpha -\beta )}{c_0}+\frac{\varkappa }{\varkappa +1}\frac{z^{\varkappa -2}}{z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta )}+\frac{\varkappa (2\varkappa +5)}{\varkappa +1}\int \frac{z^{\varkappa -3}dz}{z^{\varkappa -1}\beta +({\mathrm{tg}}^2\alpha -\beta )}.\end{array}$$