Два математика доказали, что при определённых экстремальных условиях уравнения Навье-Стокса выдают бессмыслицу


Уравнения Навье-Стокса при помощи нескольких лаконичных членов описывают одно из самых распространённых явлений физического мира: течение жидкостей. Сегодня эти уравнения, появившиеся ещё в 1820-х, используются для описания всего, от океанских течений и турбулентности, следующей за самолётом до потока крови в сердце.

Хотя физики считают эти уравнения надёжными, как молоток, математики относятся к ним с недоверием. Для математика то, что эти уравнения вроде бы работают, мало что значит. Им нужны доказательства того, что уравнения безошибочны: что для любой жидкости и для долгосрочного прогноза, распространённого сколь угодно далеко в будущее, математика уравнений не подведёт. Такую гарантию оказалось нелегко отыскать. Первый человек или команда, которая сумеет доказать, что уравнения Навье-Стокса будут работать всегда — или представить пример, доказывающий, что они не работают — сможет получить награду за решение одной из "Задач тысячелетия", анонсированных математическим институтом Клэя, и миллионом долларов в придачу [по состоянию на 2017 год только одна из семи задач тысячелетия (гипотеза Пуанкаре) решена Григорием Перельманом / прим. перев.].

Математики разработали множество способов для решения этой задачи. Новая работа, опубликованная в сентябре, ставит серьёзные вопросы по поводу того, сможет ли добиться успеха один из самых популярных подходов к задаче, разрабатываемый в течение многих лет. Работа, которую написали Тристан Бакмастер и Влад Викол из Принстонского университета, представляет собой первый результат, показывающий, как при определённых условиях уравнения Навье-Стокса дают противоречивое описание физического мира.

«Мы пытаемся понять определённые проблемы, присущие этим уравнениям, и то, почему людям, вероятно, придётся их переосмыслить», — говорит Бакмастер.

Работа Бакмастера и Викола показывает, что, если принять при решении уравнений Навье-Стокса очень грубые допущения, они начинают выдавать бессмыслицу: утверждают, что одна и та же жидкость с одними и теми же начальными условиями может прийти в два или более различных состояний. Она может течь одним образом, или же совершенно другим. Если так, то эти уравнения не могут надёжно описывать физический мир, для которого они были разработаны.

Взрывая уравнения


Чтобы понять, как уравнения могут сломаться, представьте себе океанское течение. В его рамках могут существовать локальные течения, в результате чего некоторые его части могут перемещаться в одном направлении и с одной скоростью, а другие — в другом направлении с другой скоростью. Локальные течения взаимодействуют друг с другом в постоянном взаимном действии трения и давления воды, определяющих её поток.

Математики моделируют это взаимодействие при помощи карты, сообщающей вам о направлениях и скорости потока в любой точке жидкости. Эта карта, называемая векторным полем — снимок внутренней динамики жидкости. Уравнения Навье-Стокса берут этот снимок и воспроизводят его, как видео, сообщая, как именно будет выглядеть векторное поле в каждый последующий момент времени.


Карта ветров (windy.com) работает похожим на векторное поле образом. В каждой точке у ветра есть определённое направление и сила

Эти уравнения работают. Они описывают течение жидкости так же надёжно, как уравнения Ньютона предсказывают будущие положения планет; физики постоянно используют их, и они постоянно совпадают с результатами экспериментов. Однако математикам нужно нечто большее, чем эпизодическое подтверждение — им нужно доказательство того, что уравнения не нарушаются, что вне зависимости от того, с какого векторного поля вы начнёте, и от того, как далеко в будущее вы будете его воспроизводить, уравнения всегда дадут вам новое, уникальное векторное поле.

Это и есть тема Задачи тысячелетия, спрашивающей, есть ли у уравнений Навье-Стокса решения (решение, по сути, и есть векторное поле) для всех начальных точек во все моменты времени. Эти решения должны обеспечить точное направление и силу потока в каждой точке жидкости. Решения, дающие информацию с таким бесконечно мелким разрешением, называются «гладкими». У гладкого решения каждая точка поля имеет связанный с ней вектор, позволяющий вам «гладко» путешествовать по полю, не застревая в точках, где вектор отсутствует — в точке, дальнейшее движение из которой вам будет непонятно.

Гладкие решения — полное представление физического мира, но с математической точки зрения они могут существовать не всегда. Математики, работающие над уравнениями, подобными этим, переживают по поводу такой ситуации: вы запускаете уравнения Навье-Стокса и наблюдаете за изменениями векторного поля. По прошествии какого-то конечного времени уравнения говорят вам, что некая частица жидкости двигается с бесконечной скоростью. Тогда у вас будут проблемы. В уравнения входит измерение изменений таких свойств, как давление, трение, скорость жидкости — говоря жаргонным языком, они берут производные этих величин — но производную от бесконечной величины взять не проще, чем поделить на ноль. Так что если уравнения выдают бесконечное значение, можно сказать, что они отказали вам, или «взорвались». Они уже не могут описывать последующие состояния вашей жидкости.

Такой «взрыв» — свидетельство того, что в уравнениях не хватает описания каких-то свойств физического мира, который они должны описывать. «Возможно, уравнения охватывают не все эффекты реальной жидкости, поскольку в реальной жидкости мы не ожидаем» бесконечной скорости движения частиц, как говорит Бакмастер.

Решение Задачи тысячелетия состоит либо в том, чтобы показать, что уравнения Навье-Стокса никогда не взрываются, либо найти условия, при которых это происходит. Одна из стратегий, используемых математиками — смягчить требования к тому, как точно эти уравнения должны описывать требуемые решения.

Нарушение потока


Уравнения Навье-Стокса должны описывать течение любой жидкости, с любыми начальными условиями, и распространять описание бесконечно далеко в будущее. Пытаясь доказать эту их способность, математики иногда «ослабляют», то есть, используют приближённые описания векторных полей, описывающих жидкость. Но с этим возникают трудности.

В идеале, математики хотят доказать, что применение уравнений Навье-Стокса к любой непрерывной, «гладкой» жидкости выдаст один уникальный результат.



Однако проще работать со «слабыми», не такими детализированными векторными полями. И вот математики обнаружили, что некоторые слабые описания выдают неуникальные результаты — позволяют одной и той же жидкости в одних и тех же начальных условиях течь двумя способами.



От слабых к гладким


Когда математики изучают такие уравнения, как эти, они иногда начинают расширять определение того, что считается решением. Гладким решениям требуется максимум информации — в случае с Навье-Стоксом им требуется, чтобы в каждой точке векторного поля, связанного с жидкостью, существовал вектор. Но что, если ослабить требования, и сказать, что вам нужно подсчитывать вектора только для некоторых точек поля, или нужно получить только примерные значения векторов? Такие решения называют «слабыми». Они позволяют математикам почувствовать поведение уравнения без утомительной работы по поиску абсолютно всех решений (что на практике может оказаться и невозможным).


Тристан Бакмастер, математик из Принстонского университета

«С какой-то точки зрения слабые решения ещё легче описать, чем реальные, поскольку знать нужно гораздо меньше», — сказал Камилло Де Леллис, в соавторстве с Лазло Щекелихиди написавший несколько важных работ, заложивших фундамент для работы Бакмастера и Викола.

Слабые решения бывают разной градации. Если представить себе гладкое решение в виде математического изображения жидкости с бесконечным разрешением, то слабые решения будут представлять собой нечто вроде 32-битных, 16-битных или 8-битных версий этого изображения.

В 1934 году французский математик Жан Лере определил важный класс слабых решений. Вместо работы с точными векторами, «решения Лере» берут среднее значение векторов в небольшой окрестности векторного поля. Лере доказал, что всегда можно решить уравнения Навье-Стокса, позволяя вашим решениям принимать форму такого вида. Иначе говоря, решения Лере не взрываются.

Достижение Лере определило новый подход к задаче Навье-Стокса: начать с решений Лере, о существовании которых уже известно, и посмотреть, можно ли превратить их в гладкие решения, существование которых вы хотите доказать. Этот процесс напоминает тот, где вы начинаете с грубой картинки, и смотрите, нельзя ли постепенно подкрутить разрешение, чтобы достичь идеального изображения реальности.

«Одна из возможных стратегий — показать, что эти слабые решения Лере гладкие, и если вы сможете показать, что они гладкие — вы решите Задачу тысячелетия», — сказал Бакмастер.


Влад Вкол представляет собой половину команды, вскрывшей проблемы в подходе к проверке уравнений Навье-Стокса.

Есть и ещё один подвох. Решения уравнений Навье-Стокса соответствуют реальным физическим событиям, а физические события происходят одним возможным образом. Учитывая это, хотелось бы, чтобы у ваших уравнений был только один набор уникальных решений. Если уравнения дают вам множество возможных решений, они не справляются со своей задачей.

Поэтому математики смогут использовать решения Лере для решения Задачи тысячелетия, только если решения Лере уникальны. Неуникальные решения Лере будут означать, что, согласно правилам Навье-Стокса, одна и та же жидкость с одними и теми же начальными условиями может прийти к двум разным физическим состояниям, что не имеет физического смысла, и подразумевает, что уравнения на самом деле не описывают то, что должны.

Новый результат Бакмастера и Викола — первый намёк на то, что для определённых определений слабых решений может происходить именно это.

Множество миров


В своей новой работе Бакмастер и Викол рассматривают ещё более слабые решения, чем решения Лере — решения, в которых используется тот же принцип усреднения, что у и Лере, но ослаблено ещё одно дополнительное требование (известное, как неравенство энергий). Они используют метод «выпуклого интегрирования», берущий начало из работ по геометрии математика Джона Нэша, и позднее привлечённый к изучению жидкостей Де Леллисом и Щекелихиди.

Используя такой подход, Бакмастер и Викол доказывают, что эти очень слабые решения уравнений Навье-Стокса неуникальны. Они, к примеру, демонстрируют, что если начать с полностью спокойной жидкости, к примеру, со стакана с водой рядом с кроватью, возможны два вида развития событий. Первый очевиден: вода начинает со спокойного состояния и остаётся спокойной всегда. Второй фантастичный, но математически возможный: вода начинает со спокойного состояния, взрывается в середине ночи, а затем возвращается в спокойное состояние.

«Это доказывает отсутствие уникальности, поскольку из начальных данных можно сконструировать по меньшей мере два объекта», — говорит Викол.

Бакмастер и Викол доказали существование множества неуникальных слабых решений (не только тех двух, что описаны выше) уравнений Навье-Стокса. Важность этого доказательства ещё предстоит понять. В какой-то момент слабые решения могут стать настолько слабыми, что они перестанут быть связанными с более гладкими решениями, которые должны имитировать. Если так и есть, тогда результат, полученный Бакмастером и Виколом, мало к чему приведёт.

«Такой результат однозначно является предупреждением, но можно спорить о том, что это предупреждение касается самой слабой идеи слабых решений. Существует множество слоёв более сильных решений, на гораздо лучшее поведение которых можно возлагать надежду» в случае уравнений Навье-Стокса, — говорит Де Леллис.

Бакмастер и Викол также мыслят в терминах слоёв, и он нацелились на решения Лере — на доказательство того, что и те допускают множественную физику, в которой одна и та же жидкость из одного и того же состояния может прийти к разным формам в будущем.

«Мы с Тристаном считаем, что решения Лере неуникальны. Мы пока этого не доказали, но наша работа закладывает плацдарм для атаки на эту задачу», — сказал Викол.

Комментарии (27)


  1. vesper-bot
    30.01.2018 15:32

    А что за условие описывает требование неравенства энергий в приложении уравнений Навье-Стокса к физическим объектам? Может, внутренний потенциал вроде радиоактивности, который высвобождается по рандому (т.е. неизвестно, почему), и который вполне в состоянии повлиять на динамику системы нерасчетным образом. Тогда если это требование не выполняется, говорить о валидности решений не приходится, и их работа (может, и важная) оказывается пустой в физическом смысле.


    1. EndUser
      30.01.2018 16:57
      +1

      В физическом смысле вся математика пустая, пока физики своими руками не укажут границы применимости математической формулы. Это ответственность физиков.
      А гладкие уравнения Навье-Стокса в реальном мире порвутся уже на дискретности частиц. Для начала заменить хотя бы кота Шрёдингера на эдакий вентиль — и вот вам нежданчики в гладкости математической модели течения жидкости.
      Но «виноваты» в этом не математики, а физики.
      Потому, что математики ничем никому не обязаны, ничего никому не должны.
      Эти молодые люди просто заточены на приз — показать разрыв в решениях. Математики уже полвека забавляются в специально предусмотренном разделе математики под названием «теория катастроф». И доказывают, что для той или иной функции есть ситуации, где при сколь угодно малом дельта в начальных параметрах мы получаем скачок в её поведении.
      Конкретно эти ищут совсем неустранимый разрыв специально внося дельту. Пытаются доказать, что это не «всюду дифференцируемая» функция.


      1. ilionsd
        31.01.2018 18:24

        Возможно «Теория устойчивости» более распространённое название для этого раздела математики


        1. EndUser
          31.01.2018 19:32

          Это разные разделы.


  1. hurtavy
    30.01.2018 15:35

    А проблема точно в уравнениях, а не в численных методах, которыми получают решения?


    1. Moonrise
      30.01.2018 16:22

      Данная работа теоретическая, в ней нет численных методов.


      1. hurtavy
        30.01.2018 16:27

        В ней векторные поля расчитывали же


        1. Moonrise
          30.01.2018 16:35

          Вот ссылка на работу. Там аналитическое доказательство, нет никаких расчётов.


          1. hurtavy
            30.01.2018 19:11

            Как раз там чётко указано, какое именно приближение они использовали. А в пункте 2.3 указаны параматры метода последовательных приближений, которым и получали решения


            1. Moonrise
              31.01.2018 12:53

              Мы видимо с Вами вкладываем разный смысл в понятие «численные методы». Потому что статья матановая и к численным методам не имеет ни малейшего отношения.


      1. mayorovp
        30.01.2018 16:59

        Сама по себе идея «слабых решений» очень напоминает попытку аналитики численных методов.


        1. Moonrise
          30.01.2018 17:17

          Под «слабыми решениями» традиционно понимают решения в обобщённых функциях.


    1. Victor_koly
      30.01.2018 17:32

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

      Пока Вы не решите уравнение Гамильтона-Якоби для полного гамильтониана системы молекул жидкости, любые «взрывы» будут просто грубым приближением к классической механике (другому грубому приближению к квантовой механике — более точной науке для описания финитного движения частиц).


  1. ivan01
    30.01.2018 15:38
    +2

    Что-то либо перевод некорректный, либо эти математики мало вообще в физике разбираются.

    Неуникальные решения Лере будут означать, что, согласно правилам Навье-Стокса, одна и та же жидкость с одними и теми же начальными условиями может прийти к двум разным физическим состояниям, что не имеет физического смысла, и подразумевает, что уравнения на самом деле не описывают то, что должны.

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


  1. SandroSmith
    30.01.2018 16:31

    Точки бифуркации, емнип, подразумевают наличие неучтённых флуктуаций. В случае же, если математики правы, то 2 исхода возможны без каких либо вмешательств.
    P.S. Пока писал подумал — а ведь эти флуктуации как раз и могут прятаться в неучтённых из-за "ослабления" точках. Так что да, работа, похоже, ни о чём.


    1. ivan01
      01.02.2018 13:18

      Именно, и сами уравнения навье-стокса и все их члены (высосанные не из пальца, а из физического смысла) не подразумевают учета всех флуктуаций всех молекул газа.


  1. Moonrise
    30.01.2018 16:42

    Почитал немного работу. Суть работы, в том, что на начальные условия накладываются слабые значения регулярности, то бишь изначальное векторное поле негладкое, очень ломаное. Но такой эффект есть и в обычных диффурах — если правая часть достаточно гладкая, то решение единственно, а если она ломается, то уже нет. Здесь интереснее нащупать ту границу, за которой теряется единственность.


    1. technic93
      31.01.2018 18:24

      это может иметь смысл для развития численных методов, как думаете?
      Нет ничего удивительного что если мы как то огрубили начальное состояние и не заметили там мльенького торнадо между узлами решетки для численного расчета. Пост слабоват для гиктаймз, без единой формулы о таких вещах не поговоришь. И что такое слабые решения осталось загадкой.


      1. Moonrise
        01.02.2018 08:56

        Эта статья к численным методам не имеет отношения. Одно дело доказать существование решения (или нескольких решений как в данном случае), а совсем другая работа построить толковую разностную схему, которая будет к этому решению сходиться с контролируемой погрешностью.
        Под «слабыми решениями» в работах, подобных этой, понимают решения в интегральном смысле, как правило из какого-нибудь пространства Лебега. Классическое решение дифференциального уравнения должно быть гладким (иначе как брать производную?). Поэтому трюк состоит в том, что от дифференциального уравнения переходят к интегральному, а уж интегрировать мы умеем всё подряд — и разрывные функции, и уходящие на бесконечность и чёрти какие. Вот решение такого уравнения и есть «обобщённое» или «слабое» решение.


  1. sotnikdv
    30.01.2018 23:59

    Уравнения Навье-Стокса это модель. Как и Ньютона, как и Эйнштейна и т.д. Они по определению не включают всех факторов и по определению у модели есть границы применимости. Уверен, ичследователи это знают.

    В общем где-то ученый снова изнасиловал журналиста.


    1. Victor_koly
      31.01.2018 00:42

      Вот как пример, у меня в одной модели выходил бред, если взять шаг по времени 0.05, а при 0.025 — лучше. В совсем другой модели минимальной разумной сеткой расчета я считал 250x250, хотя это не мешало кому-то считать 128x128 — для задачки типа «все интересное происходит в центре плоскости расчета» вполне адекватные результаты выходили. Может не такие «гладкие» (там комплексная величина считается) распределения величины типа амплитуды ЭМ поля, но какие-то общие закономерности увидеть можно.
      Возвращаясь к Навье-Стоксу и этой статье. Не получается ли так, что для любого случая точности сетки «поля скоростей», дающего «неоднозначные решения» можно взять сетку с шагом в 2 раза меньше, адекватно подобрать время поменьше внезапно все проблемы пропадут?


  1. aamonster
    31.01.2018 00:40

    А задачу о верёвке, падающей со стола, они не решали?
    (Кто не в курсе: это классический пример "некорректно поставленной задачи". Верёвка (считаем её бесконечно тонкой и бесконечно гибкой) лежит на столе так, что кончик — на самом краю. За какое время она упадёт со стола? У уравнений, описывающих поведение верёвки, есть два решения).


    1. mayorovp
      31.01.2018 06:38

      Откуда берется второе?


      1. aamonster
        31.01.2018 08:08

        Второе — когда верёвка продолжает лежать на столе.


        1. mayorovp
          31.01.2018 08:28

          Это вообще как?


  1. Ayvan
    31.01.2018 18:24

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


    1. Victor_koly
      31.01.2018 18:48

      Да конечно нет. Конечно при учете квантовых эффектов у Вас любая система выйдет из состояния неустойчивого равновесия, пока не перейдет в любое другое устойчивое.