Доклады Российской академии наук. Математика, информатика, процессы управления, 2020, T. 491, № 1, стр. 29-37

О СТОХАСТИЧЕСКОЙ МОДЕЛИ ВОЛНОВОЙ ТУРБУЛЕНТНОСТИ ЗАХАРОВА–ЛЬВОВА

А. В. Дымов 1*, С. Б. Куксин 234**

1 Математический институт им. В.А. Стеклова Российской академии наук
Москва, Россия

2 Université Paris-Diderot (Paris 7)
Paris, France

3 Shandong University
Jinan, PRC

4 Санкт-Петербургский государственный университет
Санкт-Петербург, Россия

* E-mail: dymov@mi-ras.ru
** E-mail: Sergei.Kuksin@imj-prg.fr

Поступила в редакцию 09.11.2019
После доработки 09.11.2019
Принята к публикации 21.01.2020

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

Аннотация

Авторы обсуждают ряд строгих результатов в стохастической модели волновой турбулентности Захарова–Львова. А именно, рассматривают уравнение Шрёдингера с (модифицированной) кубической нелинейностью и вязкостью на торе большого периода, возмущенное случайной силой, и раскладывают его решение в формальный ряд по амплитуде. Авторы показывают, что в пределе, когда амплитуда стремится к нулю, а период тора – к бесконечности, спектр энергии квадратичной срезки этого разложения сходится к решению волнового кинетического уравнения с вязкостью и внешней силой. Затем обсуждают срезки этого разложения высшего порядка.

Ключевые слова: волновая турбулентность, энергетический спектр, волновое кинетическое уравнение, кинетический предел, нелинейное уравнение Шрёдингера, стохастическое возмущение

1. ВВЕДЕНИЕ

Доказательства утверждений, приведенных ниже без ссылок, можно найти в работах [2, 3].

1.1. Модель

Рассмотрим тор $\mathbb{T}_{L}^{d} = {{\mathbb{R}}^{d}}{\text{/}}(L{{\mathbb{Z}}^{d}})$ периода $L \geqslant 1$ и размерности $d \geqslant 2$. Обозначим через H пространство ${{L}_{2}}(\mathbb{T}_{L}^{d};\mathbb{C})$, снабженное нормированной мерой Лебега:

${{\left\| u \right\|}^{2}} = \left\| u \right\|_{{{{L}_{2}}(\mathbb{T}_{L}^{d})}}^{2} = \left\langle {u,u} \right\rangle ,\quad \left\langle {u,v} \right\rangle = {{L}^{{ - d}}}\Re \int\limits_{\mathbb{T}_{L}^{d}} {u\bar {v}dx} .$

Рассмотрим модифицированное нелинейное уравнение Шрёдингера (НУШ)

$\begin{gathered} \frac{\partial }{{\partial t}}u + i\Delta u - i\nu \left( {{{{\left| u \right|}}^{2}} - {{{\left\| u \right\|}}^{2}}} \right)u = 0, \\ \Delta = {{(2\pi )}^{{ - 2}}}\sum\limits_{j = 1}^d \,\left( {\frac{{{{\partial }^{2}}}}{{\partial x_{j}^{2}}}} \right),\quad x \in \mathbb{T}_{L}^{d}, \\ \end{gathered} $
где $\nu \in \left( {0,\frac{1}{2}} \right]$. Это гамильтоново УРЧП, полученное модификацией стандартного кубического НУШ с помощью другого гамильтонового уравнения $\tfrac{\partial }{{\partial t}}u = - i\nu {{\left\| u \right\|}^{2}}u,$ поток которого коммутирует с потоком кубического НУШ. Модифицированное НУШ может быть получено из стандартного кубического НУШ с помощью подстановки

u =$exp(it\nu {\text{||}}u{\text{'|}}{{{\text{|}}}^{2}})u{\text{'}}.$

Такая модификация достаточно часто используется людьми, работающими с гамильтоновыми УРЧП, так как, сохраняя основные свойства НУШ, она уменьшает некоторые второстепенные технические трудности. Роль этой модификации состоит в том, чтобы убрать из гамильтониана НУШ интегрируемую часть (см. вторую сноску ниже). Если бы вместо кубического уравнения мы рассматривали квадратичное НУШ, описывающее трехволновую систему, то модификация не понадобилась бы.

Далее мы будем записывать решения u как $u(t,x) \in \mathbb{C}$ или $u(t) \in H$. В медленном времени $\tau = \nu t$ рассматриваемое уравнение принимает вид

(1)
$\begin{gathered} \dot {u} + i{{\nu }^{{ - 1}}}\Delta u - i({{\left| u \right|}^{2}} - {{\left\| u \right\|}^{2}})u = 0, \\ \dot {u} = \left( {\frac{\partial }{{\partial \tau }}} \right)u(\tau ,x),\quad x \in \mathbb{T}_{L}^{d}. \\ \end{gathered} $

Задача волновой турбулентности (ВТ) состоит в изучении поведения решений уравнения (1) в пределе

(2)
$\nu \to 0,\quad L \to \infty .$

Далее мы будем писать ряды Фурье функций $u(x)$ в виде

$u(x) = {{L}^{{ - d/2}}}\sum\limits_{s \in \mathbb{Z}_{L}^{d}} {{{{v}}_{s}}{{e}^{{2\pi is \cdot x}}}} ,\quad \mathbb{Z}_{L}^{d} = {{L}^{{ - 1}}}{{\mathbb{Z}}^{d}},$
где ${{v}_{s}}: = \hat {u}(s) = {{L}^{{ - d/2}}}\int\limits_{\mathbb{T}_{L}^{d}} {u(x){{e}^{{ - 2\pi is \cdot x}}}dx} .$ Тогда

(3)
${{\left\| u \right\|}^{2}} = {{L}^{{ - d}}}\sum\limits_s {{{{\left| {{{v}_{s}}} \right|}}^{2}} = :{{{\left\| v \right\|}}^{2}} = \left\| v \right\|_{{{{L}_{2}}(\mathbb{Z}_{L}^{d})}}^{2}} .$

Обозначая через $h$ гильбертово пространство $h = ({{L}_{2}}(\mathbb{Z}_{L}^{d}),\left\| {\, \cdot \,} \right\|)$, мы видим, что разложение в ряд Фурье задает изоморфизм $H \to h$, $u(x)$ $ \mapsto $ $ \mapsto $ $({{v}_{s}} = \hat {u}(s))$.

При исследовании уравнения (1) специалисты, работающие в ВТ, часто говорят о “подкачке энергии в низкие моды и ее диссипации в высоких модах”. Чтобы придать этим словам строгий смысл, Захаров и Львов [12] (см. также в [4], раздел 1.2) предложили рассматривать НУШ с (гипер)вязкостью, возмущенное случайной силой:

(4)
$\begin{gathered} \dot {u} + i{{\nu }^{{ - 1}}}\Delta u - i\rho ({{\left| u \right|}^{2}} - {{\left\| u \right\|}^{2}})u = \\ \, = - {{( - \Delta + 1)}^{{{{r}_{ * }}}}}u + {{{\dot {\eta }}}^{\omega }}(\tau ,x),\quad {{r}_{*}} > 0, \\ \end{gathered} $
где параметр ρ > 0, а случайный процесс η задан своим разложением в ряд Фурье

(5)
${{\eta }^{\omega }}(\tau ,x) = {{L}^{{ - d/2}}}\sum\limits_s {{{b}_{s}}\beta _{s}^{\omega }(\tau ){{e}^{{2\pi is \cdot x}}}} .$

Здесь $\{ {{\beta }_{s}}(\tau ),s \in \mathbb{Z}_{L}^{d}\} $ – стандартные комплексные независимые винеровские процессы11, а постоянные ${{b}_{s}} > 0$ быстро убывают при $\left| s \right| \to \infty $ и получены ограничением на $\mathbb{Z}_{L}^{d}$ положительной шварцевской функции, определенной на ${{\mathbb{R}}^{d}} \supset \mathbb{Z}_{L}^{d}$. Известно, что если ${{r}_{ * }}$ достаточно велико, то задача Коши для уравнения (4) имеет единственное глобальное решение. Применяя к уравнению (4) формулу Ито и обозначая $B = {{L}^{{ - d}}}\sum\limits_s {b_{s}^{2}} $, мы получаем баланс энергии в уравнении (4):

(6)
$\begin{gathered} \mathbb{E}{{\left\| {u(\tau )} \right\|}^{2}} + 2\mathbb{E}\int\limits_0^\tau {{{{\left\| {{{{( - \Delta + 1)}}^{{{{r}_{ * }}}}}u(s)} \right\|}}^{2}}ds} = \\ = \mathbb{E}{{\left\| {u(0)} \right\|}^{2}} + 2B\tau . \\ \end{gathered} $

Мы видим, что $\mathbb{E}{{\left\| {u(\tau )} \right\|}^{2}}$ – средняя энергия решения u, приходящаяся на единицу объема, – оказывается порядка единицы равномерно по L, вне зависимости от размера параметра $\rho $. В дальнейшем $\rho $ будет выбрано зависящим от $\nu $ таким образом, чтобы уравнение для распределения энергии решения по спектру, следующее из уравнения (4), допускало нетривиальный кинетический предел. Из приведенных результатов будет видно, что это требование однозначно определяет размер параметра $\rho $.

Уравнение (4), записанное в Фурье-представлении, имеет вид

$\begin{gathered} {{{\dot {v}}}_{s}} - i{{\nu }^{{ - 1}}}{{\left| s \right|}^{2}}{{{v}}_{s}} + {{\gamma }_{s}}{{{v}}_{s}} = \\ \, = i\rho {{L}^{{ - d}}}\sum\limits_{1,2} {\mathop \delta \nolimits_{3s}^{'12} {{{v}}_{1}}{{{v}}_{2}}{{{\bar {v}}}_{3}}} + {{b}_{s}}{{{\dot {\beta }}}_{s}},\quad s \in \mathbb{Z}_{L}^{d}, \\ \end{gathered} $
где ${{\gamma }_{s}} = {{(1 + {{\left| s \right|}^{2}})}^{{{{r}_{ * }}}}}$. Здесь, следуя традиции ВТ, мы сокращаем ${{v}_{{{{s}_{j}}}}}$ через ${{v}_{j}}$, ${{\gamma }_{{{{s}_{j}}}}}$ через ${{\gamma }_{j}}$ и т.д., а также $\sum\limits_{{{s}_{1}},{{s}_{2}} \in \mathbb{Z}_{L}^{d}} {\kern 1pt} $ через $\sum\limits_{1,2} {\kern 1pt} $, и обозначаем22

$\delta _{{3s}}^{{'12}} = \left\{ \begin{gathered} 1,\quad {\text{если}}\;\;{{s}_{1}} + {{s}_{2}} = {{s}_{3}} + s\;\;{\text{и}}\;\;\{ {{s}_{1}},{{s}_{2}}\} \ne \{ {{s}_{3}},s\} , \hfill \\ 0\quad {\text{в}}\;{\text{противном}}\;{\text{случае}}. \hfill \\ \end{gathered} \right.$

Ввиду множителя $\delta _{{3s}}^{{'12}}$, в двойной сумме выше индекс s3 является функцией индексов ${{s}_{1}},{{s}_{2}},s$, точнее, ${{s}_{3}} = {{s}_{1}} + {{s}_{2}} - s$. В интерактивном представлении

${{v}_{s}}(\tau ) = exp(i{{\nu }^{{ - 1}}}\tau {\text{|}}s{{{\text{|}}}^{2}}){{a}_{s}}(\tau ),\quad s \in \mathbb{Z}_{L}^{d},$
$v$-уравнение принимает вид a-уравнения:

(7)
$\begin{gathered} {{{\dot {a}}}_{s}} + {{\gamma }_{s}}{{a}_{s}} = i\rho {{\mathcal{Y}}_{s}}(a;{{\nu }^{{ - 1}}}\tau ) + {{b}_{s}}{{{\dot {\beta }}}_{s}},\quad s \in \mathbb{Z}_{L}^{d}, \\ {{\mathcal{Y}}_{s}}(a;t) = {{L}^{{ - d}}}\sum\limits_{1,2} {\delta _{{3s}}^{{'12}}{{a}_{1}}{{a}_{2}}{{{\bar {a}}}_{3}}{{e}^{{it\omega _{{3s}}^{{12}}}}}} . \\ \end{gathered} $

Здесь $\{ {{\beta }_{s}}\} $ – новое множество стандартных независимых комплексных винеровских процессов и

$\begin{gathered} \omega _{{3s}}^{{12}} = {{\left| {{{s}_{1}}} \right|}^{2}} + {{\left| {{{s}_{2}}} \right|}^{2}} - {{\left| {{{s}_{3}}} \right|}^{2}} - {{\left| s \right|}^{2}} = \\ = - 2({{s}_{1}} - s) \cdot ({{s}_{2}} - s) \\ \end{gathered} $
(последнее равенство имеет место в силу соотношения ${{s}_{3}} = {{s}_{1}} + {{s}_{2}} - s$, которое выполнено ввиду множителя $\delta _{{3s}}^{{'12}}$). Через ${{\mathcal{Y}}_{s}}({{a}^{1}},{{a}^{2}},{{a}^{3}};t)$ мы будем обозначать естественное полилинейное отображение, соответствующее 3-однородному отображению ${{\mathcal{Y}}_{s}}$; т.е. ${{\mathcal{Y}}_{s}}(a;t) = {{\mathcal{Y}}_{s}}(a,a,a;t)$.

1.2. Энергетический спектр и существующие результаты

Энергетическим спектром решения $u(\tau )$ уравнения (4) называется функция

$\begin{gathered} \mathbb{Z}_{L}^{d} \ni s \mapsto {{n}_{s}}(\tau ) = {{n}_{s}}(\tau ;\nu ,L) = \\ \, = \mathbb{E}{{\left| {{{{v}}_{s}}(\tau )} \right|}^{2}} = \mathbb{E}{{\left| {{{a}_{s}}(\tau )} \right|}^{2}}. \\ \end{gathered} $

Основной интерес для людей, работающих в ВТ, традиционно представляет поведение функции ${{n}_{s}}(\tau )$ и корреляций решений ${{a}_{s}}(\tau )$ в пределе (2). Одним из главных предсказаний ВТ является утверждение о том, что в этом пределе энергетический спектр ${{n}_{s}}(\tau )$ удовлетворяет волновому кинетическому уравнению (ВКУ). Существует множество физических работ, содержащих различные (но согласованные) подходы к анализу поведения энергетического спектра ns в пределе (2) и к выводу соответствующего ВКУ (см., например, [10, 11, 13] и ссылки из этих работ; см. также введение в работу [1]). Несмотря на значительный интерес в математическом сообществе, ни один из этих подходов до сих пор не был строго обоснован.

Точный смысл предела (2) не совсем ясен. Известно (см. в [8]), что для фиксированных $\rho $ и L уравнение (7) допускает предел при $\nu \to 0$. В [8] на физическом уровне строгости было показано, что если выбрать $\rho $ в виде $\tilde {\varepsilon }\sqrt L $, $\tilde {\varepsilon } \ll 1$, то последующий предел $L \to \infty $ ведет к ВКУ. Однако попытки строгого обоснования этого результата до сих пор были неудачными.

Существует лишь несколько строгих работ, посвященных анализу предела (2). В [6] для детерминистского НУШ с d = 2 исследуется предел (2) в случае, когда L стремится к бесконечности много медленнее, чем ${{\nu }^{{ - 1}}}$. Полученное элегантное описание предела далеко от предсказания ВТ и скорее должно рассматриваться в контексте теории усреднения. В недавней работе [1] рассматривается детерминистское НУШ с $d \geqslant 3$ и случайными начальными условиями $u(0,x)$ выбранными так, что аргументы $arg{{v}_{s}}(0)$, $s \in \mathbb{Z} _{L}^{d}$, коэффициентов ряда Фурье для $u(0,x)$ являются равномерно распределенными независимыми случайными величинами. В обозначениях нашей работы, в [1] доказывается, что в пределе (2), когда L стремится к бесконечности гораздо медленнее, чем ${{\nu }^{{ - 1}}}$, но в то же время не слишком медленно, для значений медленного времени $\tau $ порядка ${{L}^{{ - \delta }}}$, δ > 0, энергетический спектр ${{n}_{s}}(\tau )$ приближенно удовлетворяет линеаризации по времени ВКУ в точке τ = 0, где в ВКУ кинетический интеграл следует умножить на малый параметр $\nu $.

Другая подобная задача была рассмотрена в работе [9], также посвященной строгому исследованию детерминистского НУШ со случайными начальными условиями. Для доказательства желаемого результата авторы были вынуждены заменить тор $\mathbb{T}_{L}^{d}$ дискретным тором ${{\mathbb{Z}}^{d}}{\text{/}}(L{{Z}^{d}})$ и вместо дискретного лапласиана на ${{\mathbb{Z}}^{d}}{\text{/}}(L{{\mathbb{Z}}^{d}})$ рассматривать подходящий линейный оператор, диагональный в базисе Фурье. Предполагалось, что начальные условия $u(0,x)$ распределены согласно мере Гиббса системы, так что решение u(t) оказывалось стационарным во времени случайным процессом на H. Другая родственная задача была рассмотрена в [5].

1.3. Результаты

В данной работе мы уточняем предел (2) следующим образом:

(8)
$\begin{gathered} \nu \to 0\;\;{\text{и}}\;\;L \geqslant {{\nu }^{{ - 2 - \epsilon }}}\;\;{\text{для}}\;{\text{некоторого}}\;\;\epsilon > 0, \\ {\text{или}}\;{\text{сперва}}\;\;L \to \infty ,\;\;{\text{а}}\;{\text{затем}}\;\;\nu \to 0. \\ \end{gathered} $

Вторая возможность выше формально соответствует первой с $\epsilon = \infty $. Предположение (8) хорошо согласуется с мнением, широко распространенным в физическом сообществе, что для получения кинетического предела необходимо, чтобы L стремилось к бесконечности очень быстро, а ν–1 – не слишком быстро.

Мы снабжаем уравнение (4)(7) начальным условием

(9)
$u( - T) = 0,$
с каким-нибудь $0 < T \leqslant + \infty $, и в духе ВТ представляем решение (7), (9) в виде формального ряда по ρ:

(10)
$\begin{gathered} a(\tau ) = {{a}^{{(0)}}}(\tau ) + \rho {{a}^{{(1)}}}(\tau ) + \ldots , \\ {{a}^{{(j)}}}(\tau ) = {{a}^{{(j)}}}(\tau ;\nu ,L). \\ \end{gathered} $

Нетрудно видеть, что в случае $T = \infty $ процессы ${{a}^{{(j)}}}$ являются стационарными. Сперва, как и в физических работах (см., например, [10, раздел 6.4]), мы рассматриваем только квадратичную по $\rho $ часть этого разложения. Мы обозначаем ее ${{A}_{s}}(\tau )$, называем квазирешением и исследуем поведение соответствующего ей энергетического спектра Ns(τ) = $E({\text{|}}{{A}_{s}}(\tau ){{{\text{|}}}^{2}})$. В разделах 2–4 мы показываем, что для того чтобы при рассматриваемом пределе (8) энергетический спектр ${{N}_{s}}(\tau )$ обладал нетривиальной асимптотикой, постоянная ρ в уравнении (7) должна вести себя как $\rho \sim {{\nu }^{{ - 1/2}}}$. Соответственно, мы выбираем ее в виде $\rho = {{\nu }^{{ - 1/2}}}{{\varepsilon }^{{1/2}}}$, где параметр 0 < $\varepsilon \leqslant 1$ должен рассматриваться как фиксированная малая постоянная. Тогда энергетический спектр ${{N}_{s}}$ представляется в виде ${{N}_{s}}(\tau ) = N_{s}^{0}(\tau ) + \varepsilon N_{s}^{1}(\tau )$ + O2), где $N_{s}^{0},N_{s}^{1}$ ~ 1, равномерно по $\nu $ и L. Далее в теореме 2 мы доказываем, что при $L \geqslant {{\nu }^{{ - 2 - \varepsilon }}}$ (см. (8)), функция $s \mapsto {{N}_{s}}$ естественно продолжается до функции на ${{\mathbb{R}}^{d}}$, ε2-близкой к решению ${{m}_{s}}(\tau )$ волнового кинетического уравнения с вязкостью и внешней силой, для всех $\tau \geqslant - T$. Оно имеет вид

(11)
$\begin{gathered} {{{\dot {m}}}_{s}}(\tau ) = - 2{{\gamma }_{s}}{{m}_{s}}(\tau ) + 2b_{s}^{2} + \varepsilon {{K}_{s}}(m(\tau )), \\ s \in {{\mathbb{R}}^{d}},\quad m( - T) = 0, \\ \end{gathered} $
где Ks обозначает волновой кинетический интеграл (см. (22)). В разделе 5 мы возвращаемся к формальному ряду (10) для решения ${{a}_{s}}(\tau )$. Мы раскладываем энергетический спектр решения ${{a}_{s}}$ соответственно
(12)
${{n}_{s}}(\tau ) = n_{s}^{0}(\tau ) + \rho n_{s}^{1}(\tau ) + \ldots $
и исследуем это разложение, где, как и прежде, $\rho = {{\nu }^{{ - 1/2}}}{{\varepsilon }^{{1/2}}}$.

Так как характерные временные интервалы в нашей системе имеют порядок τ ~ 1, а медленное время $\tau $ определяется равенством $\tau = {{\nu }^{{ - 1}}}t$, то при ρ ~ ${{\nu }^{{ - 1/2}}}$ характерные быстро-временные интервалы оказываются порядка $t \sim {{\nu }^{{ - 1}}}$ ~ (размер нелинейности)–2. Именно такой временной масштаб обыкновенно рассматривается физиками.

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

Доказательства результатов, представленных в разделах 2–4, содержатся в [2], а результатов из раздела 5 – в [3]. Более подробное обсуждение полученных результатов может быть найдено во введении в [2]. В данной работе все постоянные не зависят от $\nu ,L,\rho ,\varepsilon $ и $\tau ,T$, если иное не оговорено явно.

2. РЕШЕНИЕ В ВИДЕ ФОРМАЛЬНОГО РЯДА ПО ρ

Мы представляем решение ${{a}_{s}}(\tau )$ уравнения (7), (9) в виде формального ряда (10). Тогда

$\dot {a}_{s}^{{(0)}} + {{\gamma }_{s}}a_{s}^{{(0)}} = {{b}_{s}}{{\dot {\beta }}_{s}},\quad a_{s}^{{(0)}}( - T) = 0,$
так что a(0) является гауссовским процессом

$a_{s}^{{(0)}}(\tau )$${{b}_{s}}\int\limits_{ - T}^\tau {{{e}^{{ - {{\gamma }_{s}}(\tau - l)}}}d{{\beta }_{s}}(l)} .$

Процесс a(1) удовлетворяет уравнению

$\begin{gathered} \dot {a}_{s}^{{(1)}}(\tau ) + {{\gamma }_{s}}a_{s}^{{(1)}}(\tau ) = i{{\mathcal{Y}}_{s}}({{a}^{{(0)}}}(\tau );{{\nu }^{{ - 1}}}\tau ), \\ {{a}^{{(1)}}}( - T) = 0, \\ \end{gathered} $
так что
(13)
$a_{s}^{{(1)}}(\tau ) = i\int\limits_{ - T}^\tau {{{e}^{{ - {{\gamma }_{s}}(\tau - l)}}}{{\mathcal{Y}}_{s}}({{a}^{{(0)}}}(l);{{\nu }^{{ - 1}}}l)dl} $
является третьим винеровским хаосом. Аналогично, при $n \geqslant 2$ имеем
$\begin{gathered} \dot {a}_{s}^{{(n)}}(\tau ) + {{\gamma }_{s}}a_{s}^{{(n)}}(\tau ) = \\ \, = i\sum\limits_{{{n}_{1}} + {{n}_{2}} + {{n}_{3}} = n - 1} \,{{\mathcal{Y}}_{s}}({{a}^{{({{n}_{1}})}}}(\tau ),{{a}^{{({{n}_{2}})}}}(\tau ),{{a}^{{({{n}_{3}})}}}(\tau );{{\nu }^{{ - 1}}}\tau ), \\ \end{gathered} $
так что
(14)
$\begin{gathered} a_{s}^{{(n)}}(\tau ) = \\ \, = i\int\limits_{ - T}^\tau {\sum\limits_{{{n}_{1}} + {{n}_{2}} + {{n}_{3}} = n - 1} {{{e}^{{ - {{\gamma }_{s}}(\tau - l)}}}{{\mathcal{Y}}_{s}}({{a}^{{({{n}_{1}})}}}(l),{{a}^{{({{n}_{2}})}}}(l),{{a}^{{({{n}_{3}})}}}(l);{{\nu }^{{ - 1}}}l)dl} } \\ \end{gathered} $
является (2n + 1)-м винеровским хаосом. Итерируя формулу Дюамеля в правой части равенства (14), мы можем выразить ${{a}^{{(n)}}}(l),$ $l \geqslant - T$, через процессы ${{a}^{{(0)}}}(l\,{\text{'}})$, $l\,{\text{'}} \leqslant l$.

Чтобы изучить предельное поведение корреляций решений ${{a}_{s}}(\tau )$ и корреляций энергетического спектра ${{n}_{s}}(\tau )$, записанных в виде формальных рядов (10) и (12), достаточно исследовать предельное поведение корреляций процессов $a_{s}^{{(n)}}(\tau )$. Чтобы продемонстрировать эффекты, которые здесь можно ожидать, предположим сперва, что $T = \infty $, и рассмотрим корреляции процессов $a_{s}^{{(n)}}({{\tau }_{1}})$ и $a_{{s'}}^{{(n')}}({{\tau }_{2}})$ с $n,n{\text{'}} \leqslant 1$. Нетрудно видеть, что

(15)
$\begin{gathered} \mathbb{E}a_{s}^{{(0)}}({{\tau }_{1}})a_{{s{\text{'}}}}^{{(0)}}({{\tau }_{2}}) \equiv 0, \\ \mathbb{E}a_{s}^{{(0)}}({{\tau }_{1}})\bar {a}_{{s'}}^{{(0)}}({{\tau }_{2}}) = \delta _{{s'}}^{s}{{e}^{{ - {{\gamma }_{s}}|{{\tau }_{1}} - {{\tau }_{2}}|}}}\frac{{b_{s}^{2}}}{{{{\gamma }_{s}}}}; \\ \end{gathered} $
кроме того, можно показать, что $\mathbb{E}a_{s}^{{(0)}}({{\tau }_{1}})a_{{s'}}^{{(1)}}({{\tau }_{2}})$ ≡ ≡ $\mathbb{E}a_{s}^{{(0)}}({{\tau }_{1}})\bar {a}_{{s'}}^{{(1)}}({{\tau }_{2}}) \equiv 0$. Обозначим $B(s) = \frac{{b_{s}^{2}}}{{{{\gamma }_{s}}}},$ $s \in {{\mathbb{R}}^{d}}$. Тогда, используя (13), (15) и формулу Вика, мы видим, что корреляции процессов $a_{s}^{{(1)}}(\tau )$ имеют вид
$\begin{gathered} \mathbb{E}a_{s}^{{(1)}}(\tau )a_{{s{\text{'}}}}^{{(1)}}(\tau ) = 0,\quad \mathbb{E}a_{s}^{{(1)}}(\tau )\bar {a}_{{s'}}^{{(1)}}(\tau ) = \delta _{{s{\text{'}}}}^{s}{{J}_{s}}, \hfill \\ {{J}_{s}} = \frac{{2{{\nu }^{2}}}}{{{{\gamma }_{s}}}}{{L}^{{ - 2d}}}\sum\limits_{1,2} \delta _{{3s}}^{{'12}}\frac{{{{\gamma }_{{123s}}}B({{s}_{1}},{{s}_{2}},{{s}_{3}})}}{{{{{(\omega _{{3s}}^{{12}})}}^{2}} + {{{(\nu {{\gamma }_{{123s}}})}}^{2}}}}, \hfill \\ \end{gathered} $
где
$\begin{gathered} {{\gamma }_{{123s}}} = {{\gamma }_{1}} + {{\gamma }_{2}} + {{\gamma }_{3}} + {{\gamma }_{s}}, \\ B({{s}_{1}},{{s}_{2}},{{s}_{3}}) = B({{s}_{1}})B({{s}_{2}})B({{s}_{3}}) \\ \end{gathered} $
(детали этого вычисления см. в [2]). Суммы Js могут быть приближены интегралами

$\begin{gathered} {{I}_{s}} = \frac{{2{{\nu }^{2}}}}{{{{\gamma }_{s}}}}\int\limits_{{{\mathbb{R}}^{d}} \times {{\mathbb{R}}^{d}}} {d{{s}_{1}}d{{s}_{2}}} \times \\ \, \times \frac{{{{\gamma }_{{123s}}}B({{s}_{1}},{{s}_{2}},{{s}_{3}})}}{{4{{{(({{s}_{1}} - s) \cdot ({{s}_{2}} - s))}}^{2}} + {{{(\nu {{\gamma }_{{123s}}})}}^{2}}}}, \\ {{s}_{3}} = {{s}_{1}} + {{s}_{2}} - s. \\ \end{gathered} $

Точнее,

(16)
$\left| {{{J}_{s}} - {{I}_{s}}} \right| \leqslant C_{s}^{\# }{{L}^{{ - 2}}}{{\nu }^{{ - 2}}}\quad \forall s.$

Здесь и далее символом $C_{s}^{\# }$ мы обозначаем различные непрерывные функции от s, убывающие при ${\text{|}}s{\text{|}} \to \infty $ быстрее любой отрицательной степени |s|. Согласно предположению (8), правая часть (16) мала: она не превосходит $C_{s}^{\# }{{\nu }^{{2 + 2\epsilon }}}$.

Асимптотическое поведение интегралов Is известно (см. [7, 2 ]).

Теорема 1. Интеграл Is имеет вид Is = = $\nu I_{s}^{0} + O(C_{s}^{\# }{{\nu }^{2}})$,33 где

$I_{s}^{0} = \frac{\pi }{{{{\gamma }_{s}}}}\int\limits_{{{\Sigma }_{s}}} {\frac{{B({{s}_{1}},{{s}_{2}},{{s}_{1}} + {{s}_{2}} - s)}}{{\sqrt {{{{\left| {{{s}_{1}} - s} \right|}}^{2}} + {{{\left| {{{s}_{2}} - s} \right|}}^{2}}} }}} \,d{{s}_{1}}d{{s}_{2}}{{{\text{|}}}_{{{{\Sigma }_{s}}}}}.$
Здесь ${{\Sigma }_{s}}$квадрика $\{ ({{s}_{1}},{{s}_{2}}){\kern 1pt} :({{s}_{1}} - s) \cdot ({{s}_{2}} - s) = 0\} $, а $d{{s}_{1}}d{{s}_{2}}{{{\text{|}}}_{{{{\Sigma }_{s}}}}}$элемент объема на ней, соответствующий евклидовой структуре на ${{\mathbb{R}}^{{2d}}}$.

Асимптотика, аналогичная найденной в теореме 1, может быть получена для произвольного $ - \infty \leqslant T$ < 0. Подставляя ${{s}_{1}} = s + x,{{s}_{2}} = s + y$ и обозначая $z = (x,y)$, мы переписываем интеграл $I_{s}^{0}$ в виде

$\begin{gathered} \frac{\pi }{{{{\gamma }_{s}}}}\int\limits_\Sigma {\frac{{B(s + x,s + y,s + x + y)}}{{{\text{|}}z{\text{|}}}}dz{{{\text{|}}}_{\Sigma }}} , \\ \Sigma = \{ z = (x,y){\kern 1pt} :x \cdot y = 0\} . \\ \end{gathered} $

Обозначим $F(z): = x \cdot y = - \tfrac{1}{2}\omega _{{3s}}^{{12}}{{{\text{|}}}_{{{{s}_{1}} = s + x,{{s}_{2}} = s + y,{{s}_{3}} = s + x + y}}}$. Тогда $\left| {\nabla F(z)} \right| = \left| z \right|$, так что последний интеграл в точности совпадает с интегралом $\int {B\delta (F)} $ от функции B по дельта-функции от F, см. [13, с. 67]. Так как $F = - \tfrac{1}{2}\omega _{{s + x + ys}}^{{s + xs + y}} = - \tfrac{1}{2}\omega _{{3s}}^{{12}}\delta _{{3s}}^{{12}}$, то пренебрегая минусом, записываем $I_{s}^{0}$ в виде

$\begin{gathered} \frac{\pi }{{{{\gamma }_{s}}}}\int {B\delta (F)dxdy = \frac{{2\pi }}{{{{\gamma }_{s}}}}\int {B\delta (\omega _{{s + x + y\,s}}^{{s + x\,s + y}})dxdy} } = \\ = \frac{{2\pi }}{{{{\gamma }_{s}}}}\int {B\delta (\omega _{{3s}}^{{12}})\delta _{{3s}}^{{12}}d{{s}_{1}}d{{s}_{2}}d{{s}_{3}}.} \\ \end{gathered} $

Рассматривая ${\text{|}}z{{{\text{|}}}^{{ - 1}}}dz{{{\text{|}}}_{\Sigma }} = \delta (F)$ как меру на пространстве ${{\mathbb{R}}^{{2d}}}$, сосредоточенную на квадрике $\Sigma $, в [2] мы показываем, что она дезинтегрируется как ${\text{|}}x{{{\text{|}}}^{{ - 1}}}dx{{d}_{{{{x}^{ \bot }}}}}y,$ где ${{d}_{{{{x}^{ \bot }}}}}$ обозначает меру Лебега на гиперпространстве ${{x}^{ \bot }} = \{ y{\kern 1pt} :y \cdot x = 0\} $. Точнее, для произвольной наблюдаемой f(z) выполнено

(17)
$\int\limits_\Sigma {f(z){\text{|}}z{{{\text{|}}}^{{ - 1}}}dz{{{\text{|}}}_{\Sigma }}} = \int\limits_{{{\mathbb{R}}^{d}}} {{\text{|}}x{{{\text{|}}}^{{ - 1}}}} \left( {\int\limits_{{{x}^{ \bot }}} {f(x,y){{d}_{{{{x}^{ \bot }}}}}y} } \right)dx.$

Это представление оказывается важным инструментом при анализе интегралов вида $\int {f\delta (F)} $. Оно используется ниже.

3. КВАЗИРЕШЕНИЯ

Мы начинаем наше исследование с анализа квадратичной по ρ части разложения (10), которую мы называем квазирешением A(τ) = $({{A}_{s}}(\tau ;\nu ,L),s \in \mathbb{Z}_{L}^{d})$:

${{A}_{s}}(\tau ) = a_{s}^{{(0)}}(\tau ) + \rho a_{s}^{{(1)}}(\tau ) + {{\rho }^{2}}a_{s}^{{(2)}}(\tau ).$

Рассмотрим энергетический спектр квазирешения A, ${{N}_{s}}(\tau ) = \mathbb{E}{{\left| {{{A}_{s}}(\tau )} \right|}^{2}},$ и разложим его по ρ:

(18)
$\begin{gathered} {{N}_{s}}(\tau ;\nu ,L) = n_{s}^{0}(\tau ) + \rho n_{s}^{1}(\tau ) + {{\rho }^{2}}n_{s}^{2}(\tau ) + \\ \, + {{\rho }^{3}}n_{s}^{3}(\tau ) + {{\rho }^{4}}n_{s}^{4}(\tau ), \\ \end{gathered} $
где $n_{s}^{j}(\tau ) = n_{s}^{j}(\tau ;\nu ,L)$. Здесь $n_{s}^{0} = \mathbb{E}{\text{|}}a_{s}^{{(0)}}(\tau ){{{\text{|}}}^{2}} \sim C_{s}^{\# }$44, и нетрудно видеть, что $n_{s}^{1} \equiv 0$, а $n_{s}^{2}$ = $\mathbb{E}{\text{|}}a_{s}^{{(1)}}{{{\text{|}}}^{2}}$ + $2\Re \mathbb{E}a_{s}^{{(0)}}\bar {a}_{s}^{{(2)}}$. Первое слагаемое $\mathbb{E}{\text{|}}a_{s}^{{(1)}}{{{\text{|}}}^{2}}$ в $n_{s}^{2}$ имеет порядок $\nu $ и в случае $T = \infty $ дается теоремой 1; второй член вычисляется аналогично. Таким образом,

(19)
$n_{s}^{0} \sim C_{s}^{\# },\quad n_{s}^{1} \equiv 0,\quad n_{s}^{2} \sim C_{s}^{\# }\nu .$

Мы показываем, что имеет место оценка 55

(20)
${\text{|}}n_{s}^{3}{\text{|}},{\text{|}}n_{s}^{4}{\text{|}} \leqslant C_{s}^{\# }{{\nu }^{2}},$
при $L \geqslant {{\nu }^{{ - 2 - \varepsilon }}}$ (см. (8)).

При каждом $\tau \geqslant - T$, любых ν, L и произвольном $k = 0, \ldots ,4$ функция $s \mapsto n_{s}^{k}(\tau )$ естественным образом продолжается до шварцевской функции на ${{\mathbb{R}}^{d}}$. Предел

$n_{s}^{k}(\tau ;\nu ,\infty ) = \mathop {lim}\limits_{L \to \infty } n_{s}^{k}(\tau ;\nu ,L)$
существует, также является шварцевской функцией аргумента $s \in {{\mathbb{R}}^{d}}$ и удовлетворяет (19), (20). Соответственно, предел энергетического спектра квазирешения ${{N}_{s}}(\tau ;\nu ,\infty )$ = $\mathop {lim}\limits_{L \to \infty } {{N}_{s}}(\tau ;\nu ,L)$ существует и является шварцевской функцией от $s \in {{\mathbb{R}}^{d}}$.

Соотношения (19) и (20) показывают, что правильный скейлинг параметра $\rho $ имеет вид ρ ~ ${{\nu }^{{ - 1/2}}},$ так что мы выбираем ρ в виде

(21)
$\rho = {{\nu }^{{ - 1/2}}}{{\varepsilon }^{{1/2}}},\quad \varepsilon \in (0,1].$

Действительно, при таком скейлинге процесс ${{N}_{s}}(\tau )$ является ε-возмущением линейного процесса $n_{s}^{0}$ и не сходится к $n_{s}^{0}$ в пределе (8): разложение (18) принимает вид ${{N}_{s}} = n_{s}^{0} + \varepsilon N_{s}^{1} + O({{\varepsilon }^{2}}),$ где $N_{s}^{1} = {{\nu }^{{ - 1}}}n_{s}^{2}\sim C_{s}^{\# }$.

4. ВОЛНОВОЕ КИНЕТИЧЕСКОЕ УРАВНЕНИЕ

Для действительной функции ${{\mathbb{R}}^{d}} \ni s \mapsto {{y}_{s}}$ рассмотрим соответствующий ей волновой кинетический интеграл:

(22)
$\begin{gathered} {{K}_{s}}({{y}_{ \cdot }}) = 2\pi \int\limits_{{{\Sigma }_{s}}} {\frac{{d{{s}_{1}}d{{s}_{2}}{{{\text{|}}}_{{{{\Sigma }_{s}}}}}{{y}_{1}}{{y}_{2}}{{y}_{3}}{{y}_{s}}}}{{\sqrt {{{{\left| {{{s}_{1}} - s} \right|}}^{2}} + {{{\left| {{{s}_{2}} - s} \right|}}^{2}}} }}} \times \\ \, \times \left( {\frac{1}{{{{y}_{s}}}} + \frac{1}{{{{y}_{3}}}} - \frac{1}{{{{y}_{1}}}} - \frac{1}{{{{y}_{2}}}}} \right), \\ \end{gathered} $
где для $j = 1,2,3$ мы обозначаем ${{y}_{j}} = {{y}_{{{{s}_{j}}}}}$ и где ${{s}_{3}} = {{s}_{1}} + {{s}_{2}} - s$. В обозначениях, введенных после теоремы 1, интеграл (22) принимает вид

$4\pi \int {{{y}_{1}}{{y}_{2}}{{y}_{3}}{{y}_{s}}\left( {\frac{1}{{{{y}_{s}}}} + \frac{1}{{{{y}_{3}}}} - \frac{1}{{{{y}_{1}}}} - \frac{1}{{{{y}_{2}}}}} \right)} \,\delta (\omega _{{3s}}^{{12}})\delta _{{3s}}^{{12}}d{{s}_{1}}d{{s}_{2}}d{{s}_{3}}.$

Последний интеграл в точности совпадает с кинетическим интегралом, используемым физиками для описания четырехволнового взаимодействия, см. [13, с. 71; 10, с. 91].

Рассмотрим функциональное пространство ${{C}_{r}}({{\mathbb{R}}^{d}}) = \{ x \in C({{\mathbb{R}}^{d}}):$ ${{\left| x \right|}_{r}} = {{\sup }_{s}}{{(1 + \left| s \right|)}^{r}}\left| {x(s)} \right| < \infty \} $. Представление (17) для меры ${{\left| z \right|}^{{ - 1}}}dz\,{\text{|}}{{\,}_{\Sigma }}$ влечет, что волновой кинетический интеграл $K$ определяет непрерывный оператор

$K{\kern 1pt} :{{C}_{r}}({{\mathbb{R}}^{d}}) \to {{C}_{{r + 1}}}({{\mathbb{R}}^{d}}),\quad {{y}_{s}} \mapsto {{K}_{s}}({{y}_{ \cdot }}),$
если $r > d$. Рассмотрим теперь волновое кинетическое уравнение:

(23)
${{\dot {m}}_{s}}(\tau ) = - 2{{\gamma }_{s}}{{m}_{s}}(\tau ) + 2b_{s}^{2} + \varepsilon {{K}_{s}}({{m}_{ \cdot }}(\tau )),\quad s \in {{\mathbb{R}}^{d}}.$

Существует ${{\varepsilon }_{*}} > 0$, такой что при $0 \leqslant \varepsilon \leqslant {{\varepsilon }_{*}}$ уравнение (23) имеет единственное решение $m{\text{*}}(\tau )$, обнуляющееся при $\tau = - T$ и задающее ограниченную непрерывную кривую $m{\text{*}}{\kern 1pt} :[ - T,\infty ) \mapsto {{C}_{r}}({{\mathbb{R}}^{d}})$ в каждом пространстве ${{C}_{r}}({{\mathbb{R}}^{d}})$. Оно может быть записано в виде

$\begin{gathered} m_{s}^{*}(\tau ) = m_{s}^{{*0}}(\tau ) + \varepsilon m_{s}^{{*1}}(\tau ;\varepsilon ), \\ {{m}^{{*0}}}( - T) = {{m}^{{*1}}}( - T) = 0, \\ \end{gathered} $
где ${{m}^{{*0}}},{{m}^{{*1}}} \sim 1$, $m_{s}^{{*0}}(\tau )$ совпадает с $n_{s}^{0}(\tau )$ и удовлетворяет линейному уравнению

(24)
$\dot {m}_{s}^{{*0}}(\tau ) = - 2{{\gamma }_{s}}m_{s}^{{*0}}(\tau ) + 2b_{s}^{2}.$

Далее $\varepsilon $ всегда обозначает фиксированную малую постоянную, не зависящую от $\nu $ и $L$ и удовлетворяющую $\varepsilon \in (0,{{\varepsilon }_{*}}].$ Параметр ε может быть интерпретирован как квадрат амплитуды квазирешения $A(\tau )$, записанного в правильном скейлинге. Следующая теорема – главный результат работы [2].

Теорема 2. Рассмотрим уравнение (4) с $L \geqslant {{\nu }^{{ - 2 - \epsilon }}}$ и $\rho = {{\nu }^{{ - 1/2}}}{{\varepsilon }^{{1/2}}}$. Энергетический спектр ${{N}_{s}} = {{N}_{s}}(\tau ;\nu ,L)$ его квазирешения As приближается с точностью ε2решением m* уравнения (23) в том смысле, что для произвольного $r$

${{\left| {m{\text{*}}(\tau ) - N(\tau )} \right|}_{r}} \leqslant {{C}_{r}}{{\varepsilon }^{2}}\quad \forall \tau \geqslant - T,$
с некоторым ${{C}_{r}} > 0$, если $0 < \nu \leqslant {{\nu }_{{\varepsilon ,r}}}$ для подходящего ${{\nu }_{{\varepsilon ,r}}} > 0$. Предельный энергетический спектр ${{N}_{s}}(\tau ;\nu ,\infty )$ также удовлетворяет этой оценке для каждого r и $0 < \nu \leqslant {{\nu }_{{\varepsilon ,r}}}$.

Уравнение (24) имеет единственное стационарное решение m0, $m_{s}^{0} = \frac{{b_{s}^{2}}}{{{{\gamma }_{s}}}}$, и оно является асимптотически устойчивым. Согласно теореме о неявной функции, для достаточно малого $\varepsilon $ уравнение (23) имеет единственное стационарное состояние mε, близкое к m0, и оно асимптотически устойчиво. Уменьшая при необходимости ${{\varepsilon }_{*}}$, можно считать, что единственное стационарное состояние mε определено для всех $\varepsilon \leqslant {{\varepsilon }_{*}}$. Вместе с теоремой 2 этот результат определяет асимптотическое по времени поведение энергетического спектра ${{N}_{s}}(\tau )$: для каждого r > d,

(25)
${\text{|}}{{m}^{\varepsilon }} - N(\tau ){{{\text{|}}}_{r}} \leqslant {{C}_{r}}({\text{|}}{{m}^{\varepsilon }}{{{\text{|}}}_{r}}{{e}^{{ - \tau - T}}} + {{\varepsilon }^{2}}),\quad \forall \tau \geqslant - T.$

Подходящая модификация теоремы 1 и оценки (16) влечет, что существует и отличен от нуля последовательный предел $\mathop {lim}\limits_{\nu \to 0} \,\mathop {lim}\limits_{L \to \infty } {{\nu }^{{ - 1}}}n_{s}^{2}(\tau ;\nu ,L).$ Трудно усомниться в том, что аналогичный последовательный предел существует и для ${{\nu }^{{ - 2}}}n_{s}^{4}$ (однако это пока не доказано). В этом случае, ввиду оценки (20), при $\rho = {{\nu }^{{ - 1/2}}}{{\varepsilon }^{{1/2}}}$ существует последовательный предел Ns(τ; 0, ∞ ) = $\mathop {lim}\limits_{\nu \to 0} \,\mathop {lim}\limits_{L \to \infty } {{N}_{s}}(\tau ;\nu ,L).$ Тогда он также удовлетворяет утверждению теоремы 2 и асимптотике (25).

В быстром времени t уравнение (4) с $\rho $, выбранным как выше, и $\lambda = {{(\nu \varepsilon )}^{{1/2}}}$ имеет вид

$\begin{gathered} {{u}_{t}} + i\Delta u - i\lambda ({\text{|}}u{{{\text{|}}}^{2}} - {{\left\| u \right\|}^{2}})u = \\ \, = - \nu {{( - \Delta + 1)}^{{{{r}_{ * }}}}}u + \sqrt \nu {{{\dot {\eta }}}^{\omega }}(\tau ,x), \\ \end{gathered} $
где $\left\| {u(t)} \right\| \sim 1$ при $\nu \to 0$, $L \to \infty $, согласно (6). Таким образом, мы получаем, что

1) для существования кинетического предела коэффициент λ при нелинейности должен быть порядка $\sqrt \nu $, при $\nu \to 0,L \to \infty $;

2) время, необходимое, чтобы “сесть” на кинетический режим, имеет порядок $t \sim {{\lambda }^{{ - 2}}}$. Соответствующий временной масштаб (размер нелинейности)–2 совпадает с временным масштабом, обычно рассматриваемым физиками, см. [10, 11, 13].

5. ЭНЕРГЕТИЧЕСКИЙ СПЕКТР РЕШЕНИЙ ${{a}_{s}}(\tau )$, ЗАПИСАННЫХ В ВИДЕ ФОРМАЛЬНОГО РЯДА ПО $\rho $

Вернемся к разложению (10). Как было упомянуто в разделе 2, итерируя формулу (14), можно выразить $a_{s}^{{(n)}}({{l}_{0}})$ через процессы ${{a}^{{(0)}}}(l\,{\text{'}})$ с $l\,{\text{'}} \leqslant {{l}_{0}}$. Тогда $a_{s}^{{(n)}}$ представляется в виде суммы

(26)
$a_{s}^{{(n)}}({{l}_{0}}) = \sum\limits_{\mathcal{T} \in \Gamma (n)} {{I}_{s}}({{l}_{0}};n,\mathcal{T}),$
в которой смысл индекса $\mathcal{T}$, по которому ведется суммирование, будет объяснен ниже, а ${{I}_{s}}(\mathcal{T})$ := := ${{I}_{s}}({{l}_{0}};n,\mathcal{T})$ обозначает интеграл

(27)
${{I}_{s}}(\mathcal{T}) = \int {...\int {{{L}^{{ - nd}}}} } \sum\limits_{{{s}_{1}}, \ldots ,{{s}_{{3n}}}} \,( \ldots )d{{l}_{1}} \ldots d{{l}_{n}}.$

Интегрирование здесь ведется по некоторому выпуклому многограннику в ${{[ - T,{{l}_{0}}]}^{n}}$. Суммирование берется по всем векторам ${{s}_{1}} \ldots ,{{s}_{{3n}}} \in \mathbb{Z}_{L}^{d}$, удовлетворяющим линейным соотношениям, соответствующим множителю $\delta _{{3s}}^{{'12}}$ из определения ${{\mathcal{Y}}_{s}}$. Член, обозначенный через (…), состоит из произведения функций ${{e}^{{ - {{\gamma }_{{s'}}}({{l}_{k}} - {{l}_{j}})}}}$, $exp( \pm i{{\nu }^{{ - 1}}}\omega _{{s_{3}^{'}s_{4}^{'}}}^{{\mathop {s{\text{'}}}\nolimits_1 s_{2}^{'}}})$ и процессов $[a_{{s{\text{''}}}}^{{(0)}}({{l}_{r}})]{\text{*}}$, где $a{\text{*}}$ обозначает либо a, либо $\bar {a}$, с различными индексами $k,j,r$ и различными $s{\text{'}},s_{i}^{'},s{\text{''}}$, принадлежащими множеству $\{ {{s}_{1}}, \ldots ,{{s}_{{3n}}}\} $. Это произведение имеет степень 2n + 1 по отношению к процессу a(0). Каждый интеграл ${{I}_{s}}({{l}_{0}};n,\mathcal{T})$ соответствует ориентированному дереву $\mathcal{T}$ из класса $\Gamma (n)$ деревьев, имеющих корень $a_{s}^{{(n)}}$, случайные величины $[a_{{s'}}^{{(0)}}({{l}_{j}})]{\text{*}}$ в качестве листьев и величины $[a_{{s'}}^{{(n')}}({{l}_{r}})]{\text{*}}$ с $1 \leqslant n{\text{'}} < n$ в качестве вершин, см. рис. 1. В каждую вершину дерева входит одно ребро, а выходит три. Три ребра, выходящие из какой-нибудь вершины $a_{{s{\text{'}}}}^{{(\bar {n})}}(l\,{\text{'}})$, $\bar {n} \geqslant 1$, соответствуют выбору трех членов ${{a}^{{({{n}_{1}})}}},{{a}^{{({{n}_{2}})}}},{{a}^{{({{n}_{3}})}}}$ в разложении (14) для $a_{s}^{{(n)}}(\tau ): = a_{{s{\text{'}}}}^{{(\bar {n})}}(l\,{\text{'}})$.

Рис. 1.

Дерево $\mathcal{T}$ из класса Γ(2).

Запишем энергетический спектр решения a в виде формального ряда (12). Тогда $n_{s}^{0} \sim 1$, члены $n_{s}^{1}$ = 0 и $n_{s}^{2}$ те же, что и в (18), а $n_{s}^{3}$ и $n_{s}^{4}$ немного отличаются от соответствующих членов из (18); эта неаккуратность в обозначениях не должна вызвать трудностей. Новые члены $n_{s}^{3}$ и $n_{s}^{4}$ по-прежнему удовлетворяют оценке (20) (см. ниже). Рассмотрим $n_{s}^{k}(\tau )$ с произвольным $k \geqslant 0$. Имеем $n_{s}^{k}(\tau )$ = = $\mathbb{E}\sum\limits_{{{k}_{1}} + {{k}_{2}} = k} {a_{s}^{{({{k}_{1}})}}(\tau )\bar {a}_{s}^{{({{k}_{2}})}}(\tau )} ,$ где каждый член $a_{s}^{{(k)}}(\tau )$ дается конечной суммой (26), параметризованной деревьями $\mathcal{T} \in \Gamma (k)$. Тогда

$\begin{gathered} \mathbb{E}a_{s}^{{({{k}_{1}})}}(\tau )\bar {a}_{s}^{{({{k}_{2}})}}(\tau ) = \\ \, = \sum\limits_{{{\mathcal{T}}_{1}} \in \Gamma ({{k}_{1}}),{{\mathcal{T}}_{2}} \in \Gamma ({{k}_{2}})} \mathbb{E}{{I}_{s}}(\tau ;{{k}_{1}},{{\mathcal{T}}_{1}})\overline {{{I}_{s}}(\tau ;{{k}_{2}},{{\mathcal{T}}_{2}})} . \\ \end{gathered} $

Здесь

(28)
$\mathbb{E}{{I}_{s}}({{\mathcal{T}}_{1}})\overline {{{I}_{s}}({{\mathcal{T}}_{2}})} = \int {...\int {{{L}^{{ - kd}}}\sum\limits_{{{s}_{1}}, \ldots ,{{s}_{{3k}}}} \,\mathbb{E}( \ldots )d{{l}_{1}} \ldots d{{l}_{k}}} } ,$
а матожидание берется от произведения скобок из формулы (27) для интегралов ${{I}_{s}}({{\mathcal{T}}_{1}})$ и $\overline {{{I}_{s}}({{\mathcal{T}}_{2}})} $. Так как $a_{s}^{{(0)}}(l),\;s \in \mathbb{Z}_{L}^{d}$, – гауссовские случайные величины с корреляциями (15) (в случае $T < \infty $ вторая корреляция немного меняется), то, согласно формуле Вика, каждое матожидание $\mathbb{E}{{I}_{s}}({{\mathcal{T}}_{1}})\overline {{{I}_{s}}({{\mathcal{T}}_{2}})} $ представляется в виде (конечной) суммы по всевозможным виковским спариваниям не сопряженных величин $a_{{{{s}_{j}}}}^{{(0)}}({{l}_{m}})$ с сопряженными величинами $\bar {a}_{{{{s}_{r}}}}^{{(0)}}({{l}_{q}})$. Так как эти случайные величины являются листьями деревьев ${{\mathcal{T}}_{1}}$ или ${{\bar {\mathcal{T}}}_{2}}$, эта сумма может быть параметризована фейнмановскими диаграммами, полученными спариванием деревьев ${{\mathcal{T}}_{1}}$ и ${{\bar {\mathcal{T}}}_{2}}$ по их листьям.66 Так как для $s{\text{'}} \ne s{\text{''}}$ гауссовские величины $a_{{s'}}^{{(0)}}$ и $\bar {a}_{{s{\text{''}}}}^{{(0)}}$ некоррелированы, суммирование $\sum\limits_{{{s}_{1}}, \ldots ,{{s}_{{3k}}}} {\kern 1pt} $ в (28) нужно брать только по таким векторам $({{s}_{1}}, \ldots ,{{s}_{{3k}}})$, для которых у каждой пары спаренных по Вику листьев $a_{{s'}}^{{(0)}}$ и $\bar {a}_{{s{\text{''}}}}^{{(0)}}$ индексы $s{\text{'}}$ и $s{\text{''}}$ равны. Таким образом, во всех фейнмановских диаграммах, соответствующих (28), каждый лист $a_{s}^{{(0)}}({{l}_{m}})$ дерева ${{\mathcal{T}}_{1}}$ спарен либо с листом $\bar {a}_{s}^{{(0)}}({{l}_{q}})$ дерева ${{\mathcal{T}}_{1}}$, либо с листом $\bar {a}_{s}^{{(0)}}({{l}_{q}})$ дерева ${{\bar {\mathcal{T}}}_{2}}$, и т.д. Таким образом, мы показали, что
(29)
$n_{s}^{k}(\tau ) = \sum\limits_{\mathcal{F} \in \mathfrak{F}(k)} {{\mathcal{I}}_{s}}(\tau ;k,\mathcal{F}),$
где суммирование ведется по множеству $\mathfrak{F}(k)$ фейнмановских диаграмм, соответствующих всевозможным спариваниям между листьями деревьев ${{\mathcal{T}}_{1}} \in \Gamma ({{k}_{1}})$ и ${{\bar {\mathcal{T}}}_{2}} \in \overline \Gamma ({{k}_{2}})$, ${{k}_{1}} + {{k}_{2}} = k$.

Разрешая все соотношения, наложенные на индексы sj в (28) с помощью подходящего аффинного преобразования, мы находим, что среди 3k индексов sj имеется в точности k линейно независимых. Обозначая через ${{z}_{1}}, \ldots ,{{z}_{k}} \in \mathbb{Z}_{L}^{d}$ независимые переменные, полученные из индексов sj посредством этого преобразования, мы переписываем сумму в (28) в виде ${{L}^{{ - kd}}}\sum\limits_{{{z}_{1}}, \ldots ,{{z}_{k}} \in \mathbb{Z}_{L}^{d}} {\kern 1pt} $. Приближая последнюю сумму интегралом $\int\limits_{{{\mathbb{R}}^{{kd}}}} { \ldots dz} $, где z = $({{z}_{1}}, \ldots ,{{z}_{k}}),$ мы находим, что интегралы ${{\mathcal{I}}_{s}}(\tau ;k,\mathcal{F})$ принимают сравнительно простой явный вид. А именно, для любого $s \in \mathbb{Z}_{L}^{d}$

(30)
$\begin{gathered} {{\mathcal{I}}_{s}}(\tau ;k,\mathcal{F}) = {{J}_{s}}(\tau ;k,\mathcal{F}) + O({{L}^{{ - 2}}}{{\nu }^{{ - 2}}}C_{s}^{\# }(k)), \\ {{J}_{s}}(\tau ;k,\mathcal{F}) = \\ \, = \int\limits_{{{\mathbb{R}}^{k}}} {dl} \int\limits_{{{\mathbb{R}}^{{kd}}}} {dz{{F}_{\mathcal{F}}}(\tau ,s,l,z)} \exp \left( {i{{\nu }^{{ - 1}}}\sum\limits_{i,j = 1}^k {\alpha _{{ij}}^{\mathcal{F}}({{l}_{i}} - {{l}_{j}}){{z}_{i}} \cdot {{z}_{j}}} } \right), \\ \end{gathered} $
где $l = ({{l}_{1}}, \ldots ,{{l}_{k}})$, функция ${{F}_{\mathcal{F}}}$ – гладкая по (s, z) ∈ ∈ ${{\mathbb{R}}^{{(k + 1)d}}}$ и быстро убывает по s, z и l, а ${{\alpha }^{\mathcal{F}}} = (\alpha _{{ij}}^{\mathcal{F}})$ – кососимметричная (постоянная) матрица без нулевых строк и столбцов. Ее ранг не меньше двух, и для некоторых диаграмм $\mathcal{F}$ он равен двум. Кроме того, все функции $s \mapsto {{\mathcal{I}}_{s}}(\tau ;k,\mathcal{F})$ естественно продолжаются до шварцевских функций на ${{\mathbb{R}}^{d}}$, и после этого продолжения (30) выполняется для каждого $s \in {{\mathbb{R}}^{d}}$. Следовательно, для любого фиксированного $\nu > 0$,

$\begin{gathered} \mathop {lim}\limits_{L \to \infty } {{\mathcal{I}}_{s}}(\tau ;k,\mathcal{F}) = {{J}_{s}}(\tau ;k,\mathcal{F}) \\ \forall \tau \geqslant - T,\quad s \in {{\mathbb{R}}^{d}},\quad k \geqslant 1,\quad \mathcal{F} \in \mathfrak{F}(k). \\ \end{gathered} $

Тогда, ввиду (29),

$n_{s}^{k}(\tau ;\nu ,\infty ): = \mathop {lim}\limits_{L \to \infty } n_{s}^{k}(\tau ;\nu ,L) = \sum\limits_{\mathcal{F} \in \mathfrak{F}(k)} \,{{J}_{s}}(\tau ;k,\mathcal{F}),$
для всех k, s и $\tau $.

Соотношения (19) и (20) наводят на мысль, что

(31)
${\text{|}}n_{s}^{k}{\text{|}} \leqslant C_{s}^{\# }(k){{\nu }^{{k/2}}}\quad \forall \nu $
и всех k, если L достаточно велико в терминах ν–1 и k. В этом направлении мы имеем два нижеследующих результата.

Теорема 3. Для любого k и любого $\mathcal{F} \in \mathfrak{F}(k)$,77

(32)
$\left| {{{J}_{s}}(\tau ;k,\mathcal{F})} \right| \leqslant C_{s}^{\# }(k){{\nu }^{{min(\left\lceil {k{\text{/}}2} \right\rceil ,d}}}\quad \forall \tau \geqslant - T,$
где $\left\lceil {k{\text{/}}2} \right\rceil $ обозначает наименьшее целое $ \geqslant \,\,{\kern 1pt} k{\text{/}}2$.

Если L настолько велико, что

(33)
${{L}^{{ - 2}}}{{\nu }^{{ - 2}}} \leqslant {{\nu }^{{min(\left\lceil {k{\text{/}}2} \right\rceil ,d)}}},$
то, согласно (30), ${{\mathcal{I}}_{s}}$ также удовлетворяет (32), и ввиду (29) член $n_{s}^{k}(\tau )$ ограничен правой частью (32), умноженной на ${\text{|}}\mathfrak{F}(k){\text{|}}$. Таким образом, неравенство (31) выполняется для $k \leqslant 4$, так как $d \geqslant 2$. Однако для больших k (в зависимости от d) оценка (32) существенно слабее, чем желаемая оценка (31). Наш следующий результат показывает, что оценка (32) точна в том смысле, что показатель $min(\left\lceil {k{\text{/}}2} \right\rceil ,d)$ в правой части (32) не может быть заменен на $\left\lceil {k{\text{/}}2} \right\rceil $.

Обозначим через ${{\mathfrak{F}}_{2}}(k) \subset \mathfrak{F}(k)$ множество фейнмановских диаграмм $\mathcal{F}$, для которых матрица ${{\alpha }^{\mathcal{F}}}$ из (30) имеет в точности одну ненулевою строку (и, соответственно, один ненулевой столбец). Это множество не пусто.

Теорема 4. Если $k > 2d$, то для любой диаграммы $\mathcal{F} \in {{\mathfrak{F}}_{2}}(k)$ имеет место оценка ${{J}_{s}}(\tau ;k,\mathcal{F})$ ~ ~ ${{\nu }^{d}}C_{s}^{\# }(k) \gg {{\nu }^{{\left\lceil {k{\text{/}}2} \right\rceil }}}C_{s}^{\# }(k)$. В то же время

(34)
$\left| {\sum\limits_{\mathcal{F} \in {{\mathfrak{F}}_{2}}(k)} \,{{J}_{s}}(\tau ;k,\mathcal{F})} \right| \leqslant {{\nu }^{{k - 1}}}C_{s}^{\# }(k) \ll {{\nu }^{{\left\lceil {k{\text{/}}2} \right\rceil }}}C_{s}^{\# }(k).$

Предположение о том, что сокращения, влекущие неравенство (34), типичны, представляется правдоподобным, и мы ставим следующую задачу.

Задача 1. Доказать, что $\left| {\sum\limits_{\mathcal{F} \in \mathfrak{F}(k)} {{J}_{s}}(\tau ;k,\mathcal{F})} \right|$ ≤ ≤  $C_{s}^{\# }(k){{\nu }^{{k/2}}}$ для всех k и $\nu $. В частности, ${\text{|}}n_{s}^{k}(\tau ;\nu ,\infty ){\text{|}}$ ≤ ≤  $C_{s}^{\# }(k){{\nu }^{{k/2}}}$ и неравенство (31) выполняется, если L достаточно велико в терминах ${{\nu }^{{ - 1}}}$ и $k$.

Если гипотеза, высказанная в задаче, верна, то при подстановке (21) предельное разложение ${{n}_{s}}(\tau ;\nu ,\infty ) = n_{s}^{0}(\tau ;\nu ,\infty ) + \rho n_{s}^{1}(\tau ;\nu ,\infty )$ + ... превращается в формальный ряд по $\sqrt \varepsilon $, равномерно по ν. Тогда его срезка в произвольном порядке $M \geqslant 2$, ${{n}_{{s,M}}}(\tau ;\nu ,\infty )$ = $n_{s}^{0}(\tau ;\nu ,\infty ) + \ldots + {{\rho }^{M}}n_{s}^{M}(\tau ;\nu ,\infty ),$ ε2-близка к ${{N}_{s}}(\tau ;\nu ,\infty )$ и, следовательно, также удовлетворяет утверждению теоремы 2. Мы не знаем, верно ли, что при больших M срезка ${{n}_{{s,M}}}(\tau )$ удовлетворяет уравнению (23) с погрешностью меньшей, чем ε2.

С другой стороны, если гипотеза, высказанная в задаче 1, неверна в том смысле, что

$\mathop {sup}\limits_{\tau \geqslant - T} \left\| {{{n}^{k}}(\tau ;\nu ,L)} \right\| \geqslant C{{\nu }^{{k/2 - \kappa }}},\quad \kappa > 0,$
для некоторого k88 и для всех достаточно малых $\nu $ и больших L, тогда (12) с $\rho = {{\nu }^{{ - 1/2}}}{{\varepsilon }^{{1/2}}}$ не является равномерным по ν формальным рядом по $\sqrt \varepsilon $. Мы не исключаем возможность, что такая ситуация действительно имеет место, так как НУШ появляется в физике при моделировании малых колебаний в разнообразных средах, и при его выводе отбрасываются члены высоких порядков по амплитуде. Поэтому не кажется невозможным, что кинетический предел имеет место для энергетического спектра квазирешений, но не для энергетического спектра точных решений или энергетического спектра срезок ряда (10) высоких порядков.

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

  1. Buckmaster T., Germain P., Hani Z., Shatah J. Onset of the wave turbulence description of the longtime behaviour of the nonlinear Schrodinger equation. 2019. arXiv:1907.03667.

  2. Dymov A., Kuksin S. Formal expansions in stochastic model for wave turbulence. 1: kinetic limit. 2019. arXiv: 1907.04531.

  3. Dymov A., Kuksin S. Formal expansions in stochastic model for wave turbulence. 2: method of diagram decomposition. 2019. arXiv:1907.02279.

  4. Falkovich G. Introduction to turbulence theory / In: J. Cardy, G. Falkovich, K. Gawedzki. Non-equilibrium statistical mechanics and turbulence. Cambridge: Cambridge Univ. Press, 2008. P. 1–43.

  5. Faou E. Linearized wave turbulence convergence results for three-wave systems. 2018. arXiv:1805.11269.

  6. Faou E., Germain P., Hani Z. The Weakly Nonlinear Large-Box Limit of the 2D Cubic Nonlinear Schrodinger Equation // J. Amer. Math. Soc. 2016. V. 29. P. 915–982.

  7. Kuksin S. Asymptotical Expansions for Some Integrals of Quotients with Degenerated Divisors // Russ J. Math. Phys. 2017. V. 24. P. 497–507.

  8. Kuksin S., Maiocchi A. Derivation of a Wave Kinetic Equation from the Resonant-Averaged Stochastic NLS Equation // Physica D. 2015. V. 309. P. 65–70.

  9. Lukkarinen J., Spohn H. Weakly Nonlinear Schrodinger Equation with Random Initial Data // Invent. Math. 2015. V. 183. P. 79–188.

  10. Nazarenko S. Wave Turbulence. Springer, 2011.

  11. Newell A.C., Rumpf B. Wave Turbulence // Annu. Rev. Fluid Mech. 2011. V. 43. P. 59–78.

  12. Захаров В.Е., Львов В.С. Статистическое описание волновых полей // Известия вузов. Радиофизика. 1975. Т. 18. № 10. С. 1470–1487.

  13. Zakharov V., L’vov V., Falkovich G. Kolmogorov Spectra of Turbulence. Springer, 1992.

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

Инструменты

Доклады Российской академии наук. Математика, информатика, процессы управления