Симуляция 30 000 частиц в потенциале Леннарда-Джонса.
На левую половину действует внешняя сила, направленная вниз, справа - направленная вверх.

Допилил другой вариант рендера и управление с клавиатуры. Внезапно оказалось очень залипательно разгоняться и врезаться. Цвет частицы зависит от ее скорости.

Как такового размера у частиц в этой модели нет.
Вместо этого есть пара параметров потенциала. Первый - позиция потенциальной ямы, по сути - расстояние между центрами частиц, на котором им энергетически максимально выгодно находиться.
Все частицы в этой симуляции изначально расположены на расстоянии 1.0 друг от друга. Там, где этот параметр меньше единицы (верхняя ячейка, r=0.7), частицы слабо притягиваются друг к другу до тех пор, пока не расположатся на расстоянии, задаваемым этим параметром, собираясь в более плотную субстанцию. С продвижением по ячейкам сверху вниз и слева направо, r увеличивается до значения 2.0 (нижняя правая ячейка). В тех ячейках, где r > 1.0 частицы обладают положительной потенциальной энергией, отчего и разлетаются. Потенциальной энергии частиц в нижних ячейках хватает на то, чтобы чтобы часть вещества испарилась, разорвав связи частиц с соседями.

Второй параметр задает глубину потенциальной ямы - то, насколько крепко частицы оказываются связаны друг с другом, когда расположены на максимально энергетически выгодном расстоянии r.
Сверху вниз, слева направо увеличиваем этот параметр, зафиксировав r=1 и получаем пару ячеек с газом, который моментально конденсировался в жидкость, потому что был очень холодным (частицы изначально были в потенциальных ямах), несколько ячеек, в которых присутствует и жидкость и ее пары. В последней ячейке жидкость начала застывать, потенциальной энергии гравитации не хватает на то, чтобы ее растопить.

Этим такие симуляции и притягательны - простая модель покрывает широкий спектр явлений

В смысле, этим же кодом можно промоделить испарение и конденсацию, например. Как вот тут. Сначала греем установку виртуальным термостатом, пока леннард-джонзиум не начнет испаряться, потом справа начинаем охлаждать. И он конденсируется)

Забыл поменять параметры после вчерашних экспериментов перед тем, как запустил интерактивный сим. Удивился, когда группа частиц захватила меня и не отпускала, пока не поднял температуру всей симуляции)

А что-то в этом есть. Довольно прикольно собирать эти "капли" и пытаться разогнать. Или самому разгоняться и врезаться в них. Приходится следить за температурой в симуляции, при слишком низкой не разогнать "каплю", при слишком высокой она испаряется - частицы легко преодолевают притяжение и разлетаются, не поспевая за мной.

Начал реализовывать динамический размер шага интегрирования.
С динамическим шагом можно разгоняться до больших скоростей и красиво врезаться.