
В 1934 году Жан Лере доказал, что у уравнений Навье–Стокса всегда есть решения «в обобщённом смысле». Остаётся ли решение гладким или может сорваться в сингулярность, он не узнал, и за девяносто лет не узнал никто. В 2000 году Институт Клэя включил вопрос в семь проблем тысячелетия. Формулировка составителя, Чарльза Феффермана, честная: он просит доказать любое из четырёх утверждений. (A) и (B): без внешней силы гладкость сохраняется всегда, в пространстве и на торе. (C) и (D): существует гладкая сила и гладкие начальные данные, при которых гладкого решения с конечной энергией нет.
Работа OpenAI доказывает (C) и (D). Для любой вязкости \(\nu\gt 0\) построена гладкая сила, компактная в пространстве и во времени, и течение, которое стартует из полного покоя, держит энергию ограниченной и к моменту \(t=1\) выходит на бесконечную скорость. Как это сделано — вот что стоит понять, потому что приём переворачивает саму постановку.
Задачу решают задом наперёд. Уравнение Навье–Стокса можно читать как определение силы. Выпиши любое течение u и давление p, подставь в уравнение, и всё, что не сошлось, объяви внешней силой f. Уравнение выполняется тождественно. Взрывающееся течение придумать легко. Трудность переезжает в одно место: остаток обязан быть гладким. Отдельные слагаемые в нём уходят в бесконечность, и нужно устроить между ними такое сокращение, чтобы сумма и все её производные спокойно прошли через момент взрыва. На это и потрачены 166 страниц.
Вихрь-игла. На вертикальной оси строится вихрь, схлопывающийся к началу координат. Внутри жидкость закручивается по спирали к оси и вытекает вдоль оси вверх и вниз. Разгоняет её сохранение момента импульса: величина \(r\,u_\theta\) у частицы без внешнего момента постоянна, и чем меньше r, тем больше \(u_\theta\). А продолжаться притоку позволяет несжимаемость: втекающую жидкость уносит осевой поток, иначе она забила бы ось.
Ядро сжимается самоподобно: форма профиля не меняется, меняются только масштабы. Пусть \(\tau=1-t\), время до взрыва. Радиус ядра \(\sim\tau^{1/2}\), высота \(\sim\tau^{1/2-h}\) с крошечным \(h\lt 1/100\), скорости \(\sim\tau^{-1/2-h}\). Радиус сжимается быстрее высоты, и ядро вытягивается в тонкую иглу объёмом \(\sim\tau^{3/2-h}\). Отсюда самое изящное место работы: кинетическая энергия ядра \(\sim\tau^{1/2-3h}\) стремится к нулю. Скорость бесконечна, энергия исчезает, объём падает быстрее, чем растёт квадрат скорости.
Вязкость не побеждена, а встроена. Именно вязкость отличает Навье–Стокса от Эйлера, и из-за неё уравнение считали, возможно, безопасным. Здесь считают два числа Рейнольдса. По вращению \(\mathrm{Re}_\theta\sim\tau^{-h}\to\infty\): жидкость успевает сделать всё больше оборотов за время вязкой диффузии, вращение выигрывает. По радиальному притоку \(\mathrm{Re}_r=O(1)\): вязкость остаётся равным партнёром. Скорости всех процессов, осевого и радиального сноса и диффузии поперёк ядра, растут одинаково, как \(\tau^{-1}\), и остаются в главном балансе до самого конца. Вдоль оси диффузия слабее в \(\tau^{2h}\) раз, и этим зазором конструкция пользуется.
Где ломается. Ядро надо пришить к спокойной жидкости снаружи. В кольце между ядром и внешним течением сшивка не сходится, и остаток растёт неограниченно при \(t\to1\). Сила обязана быть гладкой, а тут она должна быть бесконечной. Фоновый вихрь сам себя не держит.
Жидкость производит силу сама. В кольцо подсаживают мелкие колеблющиеся возмущения, каждое в виде полного кольца вокруг оси, локальное по радиусу, высоте и времени, запущенное экспоненциально малой затравкой. Как колебание с нулевым средним даёт ненулевую силу? Когда жидкость идёт наружу, она несёт избыток вращательной скорости; когда внутрь — недостаток. Средняя скорость ноль, но поток момента импульса — произведение двух этих величин — сохраняет знак в обоих полутактах. Переверни оба знака, произведение не изменится. Получается устойчивый средний поток импульса, и его пространственная неоднородность действует на фон как настоящая сила. Физики называют это напряжением Рейнольдса и обычно считают помехой; здесь оно инструмент, и его подгоняют так, чтобы в точности съесть расходящуюся часть остатка.
Импульсы вспыхивают и гаснут сами. Растут они на сдвиге фона, по центробежному механизму: частица с избытком вращения выталкивается наружу, а если угловая скорость падает с радиусом достаточно круто, её избыток относительно соседей растёт. Две вещи подпитывают друг друга, получается экспоненциальный рост, пока он обгоняет вязкое затухание. Гасятся они тоже сдвигом: он сносит разные радиусы с разной скоростью и укорачивает радиальную длину волны, а вязкое затухание растёт как квадрат волнового числа. Начальную длину волны выбирают так, что сперва побеждает усиление, потом неизбежно затухание. Хвост каждого импульса настолько мал, что все его производные обращаются в ноль к моменту взрыва.
Остальная сборка. Импульсов две семьи: заткнуть надо и поток момента импульса, и поток осевого импульса, а у семей разные отношения этих потоков. Ближе к взрыву импульсы ставят на всё более мелкой сетке по месту и времени, как в выпуклом интегрировании: бесконечная лестница поправок, каждая мельче предыдущей. Есть и намеренная асимметрия: при строгой симметрии верх-низ механизм усиления обнуляется в плоскости \(z=0\), поэтому осевой профиль слегка перекошен. Снаружи кольца течение чисто вращательное, не зависит от высоты и точно решает радиальное уравнение теплопроводности: его остаток ноль, силы там не нужно. Это позволяет гладко обрезать течение в пространстве, а компактность переносит результат на тор, отсюда и (D).
Что нового по сравнению с предшественниками. Приём силы и усиления через масштабы принадлежит Кордобе и Мартинесу-Сороа: сперва для трёхмерного Эйлера, потом вместе с Чжэном для Навье–Стокса с ослабленной диссипацией. Тао в 2016 году построил взрыв для усреднённого Навье–Стокса, сохраняющего энергетическое сокращение. Здесь два новых шага: настоящая вязкость и полное уравнение при любом \(\nu\gt 0\), и импульсы не просто усиливаются, а работают как источник нужного напряжения, подставленного точно под дыру фона. Кирпичи все известные: центробежная неустойчивость, снос волновых векторов сдвигом, точные вязкие сдвиговые волны, реализация напряжения колебаниями из конструкций Данери и Секейхиди для Эйлера. Новая только сборка.
Кто и как это сделал. По рассказу OpenAI, доказательство нашла группа примерно из десяти тысяч согласованных агентов на внутренней модели, за 88 часов, отправив 2,7 миллиона сообщений и потратив около 130 миллиардов выходных токенов; ещё 17 часов заняла формализация и проверка в Lean. Работа не подписана людьми: автор указан как OpenAI. За день до этого Тристан Бакмастер и Левент Альпёге выложили три препринта с тем же приёмом гладкой силы для пористой среды, Буссинеска и трёхмерного Эйлера, тоже с Lean. Между сторонами возник спор о первенстве; сами доказательства, по обоюдному признанию, разные, а результат для Эйлера у сторон отличается: у одних с силой, у других без.
Что остаётся открытым. Старт из покоя звучит эффектно, но означает обратное тому, как читается: раз жидкость стояла, всю энергию внесла внешняя сила. Взрыв здесь — не «порядок из ничего», а очень точно рассчитанная накачка. Задача без силы, пункты (A) и (B), не решена, и Институт Клэя числит проблему открытой; OpenAI на приз не претендует. Но приём, при котором старые инструменты доводит до конца машина, а проверяет другая машина, теперь опробован на задаче, которая девяносто лет не поддавалась людям.
Этот раздел адресован автору работы и его коллегам. Соседние работы нашёл смысловой поиск по нашему архиву — не по совпадению слов, а по смыслу. Мы не рецензируем: мы показываем поле вокруг.
Что сделано в работе
Построено течение, у которого скорость уходит в бесконечность за конечное время при гладкой внешней силе, из покоя и с конечной энергией. Сингулярность живёт в игле, сжимающейся быстрее, чем разгоняется, а недостающую силу поставляет сама жидкость через среднее напряжение колебательных импульсов. Для нас это первый строго построенный пример локальной сингулярности сплошной среды при конечной энергии, с которым можно работать дальше.
В чём сила
Механизм явный: вихрь-игла, кольцо, две семьи импульсов, тепловая внешность. Значит, его можно переносить: считать эффективную метрику, спрашивать, что с ним делает исчезающая вязкость, и пробовать воспроизвести в лаборатории.
Куда можно двигаться
Что можно попробовать
Почему это может быть важно
Если сингулярность течения переносится на эффективную метрику, обычная и квантовая жидкости становятся стендом для вопросов, которые в самой гравитации не проверить: как ведёт себя энтропия рядом с сингулярностью и может ли порядок удерживаться там, где непрерывное описание ломается. Работа OpenAI впервые даёт такому стенду явную, строго построенную сингулярность с конечной энергией.
Работы рядом в нашем архиве
Работа не проходила рецензирование. Это наш разбор — мнение, а не приговор.
Сильная сторона
Поворот постановки: течение выписано, сила определена как остаток, и вся трудность честно перенесена в гладкость остатка. Самоограничивающиеся импульсы, которые сами вспыхивают на сдвиге и сами гаснут вязкостью, и подгонка их среднего напряжения под расходимость фона — это законченный, проверяемый механизм, а не оценка «сверху».
Что усилит работу
Вопросы автору
Работу не подписал ни один человек: в графе «автор» стоит OpenAI. Доказательство нашли около десяти тысяч агентов за 88 часов, а проверял его на 17 часов дольше не человек, а Lean.
Работа опубликована не на arXiv. Здесь её страница, PDF и формализация: проверить можно самим, а не только прочитать наш разбор.