I. Метод обобщенных функциональных подстановок в задачах нелинейной диффузии и взаимодействия излучения с веществом

​Аннтотация

Проект направлен на развитие и применение метода обобщенных функциональных подстановок в теории нелинейных волновых и диффузионных уравнений к исследованию процессов в нелинейных средах с диффузией под действием внешнего излучения. В проекте предусматривается исследование процессов образования когерентных структур в нелинейных средах с диффузией под действием внешних излучений и процессов самоорганизации. Одним из основных направлений работ будет теоретический анализ точных решений нелинейных уравнений диффузионного типа с помощью метода обобщенных функциональных подстановок типа Коула-Хопфа и его обобщений

​Цели и задачи

Проект направлен на решение фундаментальной задачи описания процессов нелинейной диффузии в различных конденсированных средах, в том числе кристаллах, под действием внешнего излучения, приводящих к образования когерентных структур за счет нелинейного взаимодействия и самоорганизации отдельных компонент среды. Задачами проекта являются следующие задачи:
1) Разработка новых вариантов метода функциональных подстановок конкретно для задач нелинейной диффузии в средах, находящихся под действием внешнего нейтронного и лазерного облучения.
2) Выяснение условий образования когерентных структур дефектов в материалах под действием внешнего излучения и за счет самоорганизации среды;
3) Построение аналитических моделей процессов образования когерентных структур в материалах под облучением и сопоставление их с экспериментальными данными.
4) Разработка новых методов асимптотического анализа нелинейных диффузионных уравнений с целью выяснения устойчивости когерентных структур в соответствсующих системах под облучением.

Общий план проведения работ по проекту на  2016 год:

Первый год (2016):
1) Развитие метода обобщенных функциональных подстановок Коула-Хопфа для решения задач нелинейной диффузии. Предположительно будут получены специфические типы матричных подстановок для построения уравнений нелинейной диффузии в размерности пространства 1+2 и 1+3.
2) Разработка варианта метода автомодельных решений уравнений нелинейной диффузии для анализа роста когерентных структур на ранних этапах их образования под действием внешнего излучения или самоорганизации в нелинейной среде. Будут проанализированы условия роста когерентных структур за счет самоорганизации в нелинейной среде с диффузией и найдены условия их возникновения при наличии внешних потоков излучения.
3) Развитие методов асимптотического анализа решений уравнений нелинейной диффузии для выяснения устойчивости когерентных структур. Предполагается, что будут получены критерии возникновения устойчивых когерентных структур при наличии внешних потоков излучения и случайных флуктуаций в системе.

Результаты за 2016 год.

1. Метод функциональных подстановок для уравнений диффузионного типа

Этот раздел содерждит результаты в основном представленные в статье: Журавлев В.М. Изв. ВУЗОВ. Поволжский регионю Серия физико-математическая. 2016

​1.1.  Общая идеология метода функциональных подстановок

Опираясь на результаты работ [1, 2, 3, 4] в качестве исходной системы линейных уравнений, рассмотрим уравнения следующего вида $$T_x=A(x,t)T,~~T_t=B(x,t)T\tag{1.1}\label{Eq11}$$ относительно одной вспомогательной комплексной функции $T(x,t)$  двух вещественных переменных $x$  и $t$   и двух вспомогательных комплексных функции $A(x,t)$  и $B(x,t)$. Дифференцируя однократно первое уравнение по $t$ , а втрое - по $x$, получаем вместе с ними замкнутую однородную алгебраическую систему четырех уравнений относительно функции $T$  и трех первых ее производных: $T_x,~T_t$  и $T_{xt}$. Условием совместности этой системы является уравнение: $$A_t=B_x.\tag{1.2}\label{Eq12}$$ 
В силу этого все производные функции $T(x,t)$  можно выразить рекуррентно через саму функцию $T(x,t)$  или одну любую ее производную по формулам:$$T^{[n,k]}=\frac{\partial^{n+k}T}{\partial x^{n}\partial t^{k}}=A^{[n,k]}T,\tag{1.3}\label{Eq13}$$ где функции $A^{[n,k]}(x,t)$ удовлетворяют рекуррентным соотношениям: $$A^{[n+1,k]}=A^{[n,k]}_x+A^{[n,k]}A,~A^{[n,k+1]}=A^{[n,k]}_t+A^{[n,k]}B\tag{1.3a}\label{Eq13a}$$ и $A^{[1,0]}=A,~A^{[0,1]}=B.$
 
К базовой системе  (1) можно добавить произвольное уравнение относительно $T(x,t)$ , которое в итоге с помощью соотношений (1.3) превращается в нелинейное уравнение, относительно функций $A,B$. При этом это уравнение образует замкнутую систему вместе с уравнением  (1.2).  В этом случае базовые соотношения  (1.1) можно рассматривать как обобщенные подстановки Коула-Хопфа. Эти подстановки будем называть подстановками первого уровня.

1.2. Уравнения Лапласа и гидродинамика

1.3. Подстановки второго уровня

Полученная интерпретация уравнений  (1.5) оказывается не единственной. Эти уравнения обладают внутренней скрытой структурой второго уровня, которая выявляется лишь при специальном представлении функций $A$ и $B$. Рассмотрим следующую замену переменных: $$A=u(x,y)\cos(\phi(x,y)),\quad B=u(x,y)\sin(\phi(x,y)).\tag{1.9}\label{Eq19}$$
Подставляя эти соотношения в уравнения  (1.2) и  (1.5) , получаем: $$(u_x+u\phi_y)\sin(\phi)-(u_y-u\phi_x)\cos(\phi)=0$$$$(u_y-u\phi_x)\sin(\phi)-(u_x+u\phi_y)\cos(\phi)+u^2=0.$$

Рассматривая эту систему как систему алгебраических уравнений относительно $p=u_x+u\phi_y$  и $q=u_y-u\phi_x$ , находим их решение: $$\theta_x+\phi_y=-u\cos\phi=-A,\tag{1.10}$$$$\theta_y-\phi_x=-u\sin\phi=-B.\tag{1.11}$$ Здесь $\theta=\ln u$ . Из этой системы перекрестным дифференцированием, используя  (1.5) , получаем следующие уравнения: $$\phi_{xx}+\phi_{yy}=0,\tag{1.12}$$ $$\theta_{xx}+\theta_{yy}=u^2=e^{2\theta}.\tag{1.13}$$ Первое из этих уравнений есть вновь уравнение Лапласа, а второе - уравнение Лиувилля! Эти уравнения и отражают скрытую динамическую симметрию уравнений  (1.5) , которая заключается в наличии прямых обобщенных подстановок Коула-Хопфа для решений уравнений  (1.12) и  (1.13) , связывающих их с решениями уравнения  (1.4) : $$\theta=\frac{1}{2}\ln\left(\left(\frac{\partial\ln T}{\partial x}\right)^2+\left(\frac{\partial\ln T}{\partial y}\right)^2\right),\tag{1.14}$$ $$\phi={\rm arctg}\left(\frac{T_{y}}{T_{x}}\right).\tag{1.15}$$ Эти подстановки будем называть подстановками второго уровня.

3. Подстановки, связанные с уравнением Гельмгольца 

Найденная богатая структура преобразований уравнения Лапласа в уравнение Лиувилля и в себя может распространятся частично и на другие типы линейных уравнений.
    Рассмотрим в качестве такого обобщения замыкающее уравнение Гельмгольца: $$T_{xx}+T_{yy} = \gamma^2 T,\tag{1.24}\label{EqTKG}$$ где $\gamma$ - постоянная. Замыкающее уравнение в виде коэффициентов  $A,B$ теперь будет иметь такой вид: $$A_x+B_x+A^2+B^2 = \gamma^2.$$
По аналогии с предыдущим случаем ищем решение для $A,B$ в том же виде (1.9):
$$A=u\cos(\phi),~~B=u\sin(\phi).$$ В этом случае имеем:$$\Big(u_x+u\phi_y\Big)\sin\phi - \Big(u_y-u\phi_x\Big)\cos\phi =0,$$$$\Big(u_x+u\phi_y\Big)\cos\phi + \Big(u_y-u\phi_x\Big)\sin\phi + u^2-\gamma^2 =0.$$ Полагая $$u_x+u\phi_y = R^2\cos\phi,~~~u_y-u\phi_x = R^2\sin\phi,\tag{1.25}\label{DefR}$$ приходим к следующему соотношению:$$R^2+u^2 = \gamma^2.$$ Отсюда следует, что существует такая функция $\chi$, что:$$R= \gamma\cos\chi, u = \gamma \sin\chi.$$ Это означает, что вещественные решения для $u$ и $R$, соответствующие (\ref{DefR}), являются ограниченными: $|R|<\gamma$, $|u|<\gamma$.  Исключая из (\ref{DefR}) функцию $R$, находим: $$\frac{1}{u}u_x+\phi_y = - \Big(u-\frac{\gamma^2}{u}\Big)\cos\phi,\tag{1.26a}$$$$\frac{1}{u} u_y-\phi_x = - \Big(u-\frac{\gamma^2}{u}\Big)\sin\phi.\tag{1.26b}\label{Eq26}$$ Последние соотношения нельзя привести к (1.10)-(1.11). Однако мы можем их преобразовать к комплексифицированной форме. Именно, умножая  (1.26b) на мнимую единицу и вычитая полученный результат из  (1.26a), приходим к следующему уравнению: $$\left(\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}\right)\Big(\theta+i\phi\Big)=\left(\gamma^2-u^2\right)e^{-\theta-i\phi}.\tag{1.27}\label{EqCu}$$ Введем обозначение: $w=\theta+i\phi$. Тогда уравнение (1.27) можно преобразовать к следующему виду: $$\frac{\partial}{\partial z}w=\left(\gamma^2-u^2\right)e^{-w},$$ где $$\frac{\partial}{\partial z}=\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}$$ - производная по комплексной переменной $z=x+iy$. Введем еще одно обозначенние, полагая: $\Psi=e^{w}$. В этом случае уравнение для $w$ примет такой вид: $$\Psi_{z}=\gamma^2-|\Psi|^2.\tag{1.28}\label{EqPsi}$$