Новый метод идентификации свободной поверхности для 3D MPS-метода

Почему важно отслеживать движущиеся поверхности воды

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

Другой подход к разграничению воды и воздуха

В методе Moving Particle Semi-Implicit (MPS) вода моделируется бесчисленными мелкими частицами, которые движутся и взаимодействуют. Традиционные подходы определяют, лежит ли частица на свободной поверхности, простым подсчётом числа её соседей: меньше соседей обычно означает, что частица на грани или рядом с ней. Однако это эмпирическое правило даёт сбои, когда частицы распределены неравномерно, как часто бывает при турбулентных течениях или у острых углов. Авторы в свою очередь сосредотачиваются на том, как именно соседние частицы расположены в пространстве, используя величину, названную дивергенцией относительных позиций, которая измеряет, насколько «разбросаны» соседи вокруг данной частицы. Если соседи сосредоточены в основном с одной стороны, значение дивергенции меняется так, что это сигнализирует о вероятной поверхности.

Добавление направления в тест поверхности

Хотя улучшённая мера дивергенции хорошо согласуется с теорией, авторы показывают, что в трёх измерениях простая пороговая проверка по этой величине не всегда может идеально разделить внутренние и поверхностные частицы, особенно возле острых кромок и внутри полостей. Чтобы исправить это, они вводят дополнительный шаг, который оценивает локальное направление поверхности, или нормаль, исходя из окружающих частиц. Используя это направление, метод задаёт коническую область, обращённую наружу от каждой кандидатной частицы. Если в пределах определённого расстояния внутри этого конуса не находится других частиц, частица считается лежащей на свободной поверхности. Такой комбинированный подход, называемый RPD+NV, учитывает и расположение соседей, и наличие пустого пространства, давая более надёжную картину истинного интерфейса между жидкостью и воздухом.

Тестирование метода на сложных формах и спокойной воде

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

Фиксация волн, ударов и качки

Далее команда испытывает метод в сильно динамичных ситуациях. В трёхмерной симуляции прорыва дамбы столб воды рушится и врезается в стену, создавая быстро изменяющиеся свободные поверхности, всплески и воздушные карманы. По сравнению с традиционными критериями, основанными только на числе соседей или одной лишь дивергенции, метод RPD+NV отслеживает более непрерывные и детализированные поверхности и даёт более плавные, реалистичные поля давления. Существенно, пиковое ударное давление по стене оказывается гораздо ближе к лабораторным измерениям. В случае демпфированной качки, где волны в баке постепенно затухают, новый метод обеспечивает наиболее плавное развитие сил на стенках, указывая на снижение числового шума. Наконец, в тесте, имитирующем качку топлива в кольцевом баке, который раскачивается, формы поверхности и пики давления из симуляции хорошо согласуются с данными высокоскоростной съёмки и датчиков давления, снова без ложных всплесков давления на поверхности.

Что это значит для практических симуляций

Для читателей, не занятых вычислительной гидродинамикой, практический вывод таков: эта работа улучшает способ, которым компьютеры определяют, где заканчивается вода и начинается воздух в моделях на основе частиц. Объединив уточнённую меру окружения соседей каждой частицы с направленным тестом «пустой конической области», новый метод RPD+NV резко сокращает ошибочные классификации поверхностных частиц. В результате получаются более чистые, более устойчивые предсказания давлений для широкого спектра течений — от стоячей воды до сильных ударов и сложной трёхмерной качки. Поскольку метод добавляет небольшие вычислительные затраты, он представляет собой практическое улучшение для инженеров, моделирующих прибрежные сооружения, движение судов, ёмкости для хранения и другие системы, где точная фиксация движущихся поверхностей воды жизненно важна.

Цитирование: Geng, C., Wang, Wh., Heng, My. et al. A new free surface identification method for 3D MPS method. Sci Rep 16, 13829 (2026). https://doi.org/10.1038/s41598-026-44218-9

Ключевые слова: моделирование жидкости на основе частиц, обнаружение свободной поверхности, численное воздействие волн, встряхивание жидкостей в баках, метод движущихся частиц