Мильнуватий піну та мінімальні поверхні: Фізика робить обчислення варіацій

Поверхневий натяг — це мінімізатор площі, а не формувальник

Енергія мильної плівки прямо пропорційна її поверхневій площі — кожну додаткову квадратну міліметр плівки коштує фіксована енергія поверхневого натягу, незалежно від її форми. Залишена сама по собі, будь-яка фізична система прагне до мінімальної енергії, тому мильна плівка, що перекриває фіксовану межу (наприклад, вигнутий дріт), приймає форму, яка має найменшу можливу площу серед усіх поверхонь, які розділяють цю межу. Саме тому мильні плівки виглядають бездоганно: вони не декоративні, а є фізичним втіленням розв’язаної задачі з математичного аналізу (калькулю варіацій), знайденої природою в реальному часі за допомогою нічого складнішого, ніж молекули, що тягнуться один до одного.

Нульове середнє викривлення

Математична умова для поверхні, що мінімізує площу (більш точно — площа-критична — локальний мінімум, максимум або седло) полягає в тому, що її середнє викривлення H, середнє значення двох основних кривин у кожній точці, має бути рівним нулю всюди:

H = (κ₁ + κ₂) / 2 = 0 (умова мінімальної поверхні, в кожній точці) k₁, k₂ = дві основні кривини в точці (максимум і мінімум викривлення, у перпендикулярних напрямках) — для мінімальної поверхні вони рівні та протилежні: поверхня згинається в один бік і в інший на ту ж саму величину, як у випадку седло. Плоский план trivially задовольняє H = 0 всюди, тому плоский плівка мила через плоску петлю є дійсним (якщо нудним) мінімальною поверхнею. Більш цікаві межі змушують плівку набувати справді викривлених форм, які все одно балансуються до нуля середнього викривлення всюди — математика Жана Батьйоля Музньє, який вперше сформулював цю умову в 1776 році, пов'язуючи інтуїтивне фізичне спостереження з диференціальним геометрією більше ніж за століття до того, як плівки мила були вивчені експериментально в серйозні.

H = (κ₁ + κ₂) / 2 = 0        (minimal surface condition, at every point)

κ₁, κ₂ = the two principal curvatures at a point (max and min curvature, in
         perpendicular directions) — for a minimal surface they are equal
         and opposite: the surface curves up one way and down the other by
         exactly the same amount, like a saddle

Класичні приклади: катаноїд, спіраль, ампутація

Найпростіша непланарна мінімальна поверхня – катаноїд, форма, яку набуває мильна плівка між двома паралельними коаксіальними кільцями — тінь, що звужується до середини, яка генерується обертанням катенарії (крива вішаного ланцюга) навколо осі. Якщо від’єднати два кільця занадто сильно, катаноїд стає нестабільним і колапсує на дві окремі плоскі диски, по одному на кожне кільце — вражаюче демонстрація того, що «мінімальний» означає лише локально мінімальний; за критичної відстані між ними від’єднана дводискова конфігурація має меншу загальну площу. Спіраль, яка генерується спіральним рухом прямої лінії вздовж осі (форма сходи з викруткою), є іншим класичним прикладом, відкритим Жан-Батьйолем Мушніє разом із катаноїдом, і вона відома як єдина відома мінімальна поверхня (окрім площини), яка також є рулетовою поверхнею — повністю прокручується прямою лінією. Ампутація – це найпростіший спосіб побачити H = 0 на око: вона вигинається вгору в одному напрямку та вниз у перпендикулярному напрямку, з абсолютно однаковими величинами.

Закони Плоте

Жозеф Платто, сліпий бельгійський фізик, десятиліттями у 19 столітті систематично вивчав мильні плівки та сформулював правила, які повинні дотримуватися будь-яка конфігурація з’єднаних плівок — пізніше математично доведено Жан Тейлор у 1976 році:

1. Плівки є гладкими поверхнями постійної середньої кривизни між точками перетину. 2. Правильно три шари плівок зустрічаються вздовж будь-якого краю, завжди під кутом 120° один до одного (ніколи 2, ніколи 4+) — баланс сили поверхневого натягу вимагає цього. 3. Края (де зустрічаються три шари) збігаються у чотири вершини, завжди під тетраедричним кутом ≈ 109.47°. Ці правила пояснюють, чому кластер мильних бульбашок ніколи не показує випадкового переплетування кутів плівок — кожна точка перетину, яку ви знайдете в реальному кластері бульбашок, точно дотримується правил 120°/109.47°, оскільки будь-яка інша конфігурація має незбалансовану силу натягу, яка миттєво змінює точку перетину, поки вона не стабілізується під цим кутом.

1. Films are smooth surfaces of constant mean curvature between junctions.
2. Exactly THREE film sheets meet along any edge, always at 120° to
   each other (never 2, never 4+) — force balance between equal
   surface tensions demands it.
3. Edges (where three sheets meet) themselves meet in FOURS at a
   point, always at the tetrahedral angle ≈ 109.47°.

Численное вычисление минимальной поверхности

Для малых отклонений от плоской поверхности уравнение минимальной поверхности линеаризуется в уравнение Лапласа, ∇²h ≈ 0, где h(x, y) — высота пленки над опорной плоскостью — то же уравнение, которое определяет установившийся теплопроводность и электростатический потенциал. Его численное решение на сетке, заданном фиксированными высотами границы, представляет собой классическую задачу расслабления: повторно заменяйте высоту каждой внутренней точки средним значением ее соседей до тех пор, пока вся сетка не перестанет изменяться.

// Jacobi расслабление в направлении ∇²h = 0 в внутренней точке for (let iter = 0; iter < maxIterations; iter++) { for each interior point (x, y): hNext[x][y] = 0.25 * (h[x-1][y] + h[x+1][y] + h[x][y-1] + h[x][y+1]); // граничные точки (каркас) остаются FIXED на каждой итерации [h, hNext] = [hNext, h]; }

// Jacobi relaxation toward ∇²h = 0 on an interior grid point
for (let iter = 0; iter < maxIterations; iter++) {
  for each interior point (x, y):
    hNext[x][y] = 0.25 * (h[x-1][y] + h[x+1][y] + h[x][y-1] + h[x][y+1]);
  // boundary points (the wire frame) stay FIXED every iteration
  [h, hNext] = [hNext, h];
}

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Soap Bubble і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.