Puchary z mydłem i powierzchnie minimalne: fizyka wykonująca rachunek wariacyjny

Napięcie powierzchniowe jest minimizatorem pola, a nie rysownikiem kształtu

Energia filmu mydłowego jest proporcjonalna bezpośrednio do jego pola powierzchniowego — każdy dodatkowy milimetr kwadratowy filmu kosztuje ustalony ilość energii napięcia powierzchniowego, niezależnie od jego kształtu. Pozwolone na samodzielne działanie, dowolny fizyczny system relaksuje się w kierunku niższej energii, więc film mydłowy rozciągający się na ustalonej granicy (np. krzywą pręta), osiąga kształt o najmniejszym możliwym polu powierzchni wśród wszystkich powierzchni dzielących tę granicę. To dlatego filmy mydlowe wydają się niezwykle eleganckie: są one nie dekoracyjne, a fizyczną reprezentacją rozwiązania zagadnienia zmiennych całkowych, znalezionego przez naturę w czasie rzeczywistym poprzez niczym innym jak tylko cząsteczki przyciągające swoje sasiady.

Zerojaka krzywizna średnia

Matematyczne warunki dla powierzchni bycia minimalnie poleminną (bardziej precyzyjnie, poleminną — lokalny minimum, maksimum lub siodło pola) to, że jej krzywizna średnia H, srednia z dwóch głównych krzywizn w każdym punkcie, musi być dokładnie równa zero wszędzie:

H = (κ₁ + κ₂) / 2 = 0 (warunek powierzchni minimalnej, dla każdego punktu) κ₁, κ₂ = dwie główne krzywizny w punkcie (maksimum i minimum krzywizny, w prostopadłych kierunkach) — dla powierzchni minimalnej są równe i przeciwne: powierzchnia skrzywia się jednym kierunkiem w górę, a drugim w dół o dokładnie tą samą wartość, jak siodło Prosta płaszczyzna prosto spełnia H = 0 wszędzie, dlatego płaskie filmy mydłowe na płaskiej pętli są poprawnym (ale nudnym) przykładem minimalnej powierzchni. Większe wyzwania graniczne zmuszają film do prawdziwych krzywych kształtów, które nadal równoważnie skończą się na zero średniej krzywizny wszędzie — matematyka Jean Baptiste Meusniera, który pierwszy wykazał ten warunek w 1776 roku, łącząc intuicyjne obserwacje fizyczne z geometrią różniczkową ponad sto lat przed experimentami nad filmami mydłowymi.

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

Przykłady klasyczne: catenoid, helicoid, siodło

Najprostsza niepłaskościowa powierzchnia minimalna to catenoid — kształt, jaki formuje film mydlany między dwiema równoległymi osiowymi pęczkami — szyja zciśnięta w kierunku środka, wygenerowana przez obrót krzywej catenary (krzywej łańcuchowej) wokół osi. Jeśli jednak oddzielić te pęczki za daleko, catenoid staje się niestabilny i zanika do dwóch osobnych płaskich dysków — jednego dla każdego pęczka, co jest wyraźnym dowodem na to, że "minimalna" oznacza tylko lokalnie minimalną; po pewnej krytycznej odległości, konfiguracja z dwoma oddzielnymi dyskami ma mniejszą całą powierzchnię. Helicoid, wygenerowany przez spirali prostej wraz ze zmianą jego pozycji wokół osi (kształt schodów spiralnych przybliżony jest krzywą rampy), to drugie klasyczne rozwiązanie, odkryte razem z catenoidem przez Jean Baptiste Meusnier. Widocznie jest on jedynym znanych przykładem powierzchni minimalnej (poza płaszczyzną) która jednocześnie jest powierzchnią rządzoną — wygenerowaną całkowicie przez proste linie. Powierzchnia siodłowa między czterema alternującymi się wysokimi i niskimi barierami to najprostsza metoda do widocznej ilustracji H = 0: krzywo nachyla się w jednym kierunku, a w przeciwnego kierunku, dokładnie taka sama ilość.

Zasady Plateau'a

Joseph Plateau, niezrównany belgijski fizyk, poświęcił kilkadesiąt lat w XIX wieku systematycznym badaniu foli soapowych i przekształcił swoje obserwacje w zasady, które muszą być spełnione przez każdą konfigurację powiązanych folii — te zasady były potem matematycznie dowodzone przez Jean'a Taylor'a w 1976 roku:

1. Folie są gładkimi powierzchniami o stałej średniej krzywizny między połączeniami. 2. TRZY foliowe płaszcze zawsze spotykają się na każdej krawędzi, zawsze pod kątem 120° do siebie (nikt nie ma krawędzi z dwoma lub czterema i więcej foliowymi płaszczami — równoważność sił napinania powierzchni wymaga tego. 3. Krawędzie, gdzie trzy foliowe płaszcze się spotykają, same spotykają się w czwórkę punktów pod kątem tetraedrycznym ≈ 109,47°. Te zasady są przyczyną tego, że grupa bąbelków mydlanych nigdy nie pokazuje przypadkowego skrętu kątów foliowych — każda połączenie, które można znaleźć w rzeczywistej grupie bąbelków, spełnia dokładnie zasady 120°/109,47°, ponieważ dowolna inna konfiguracja ma niezrównaną siłę napinania powierzchni, która natychmiast przekształca połączenie w kąt stabilny.

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°.

Numericzne znalezienie powierzchni minimalnej

Dla małych odstępstw od płaskości równanie powierzchni minimalnej liniarnie redukuje się do typowego równania Laplace’a, ∇²h ≈ 0, gdzie h(x, y) to wysokość folii względem płaszczyzny referencyjnej — taka sama równość, która rządzi spadkiem temperatury w stanie ustalonym i potencjałem elektrycznym. Numeryczne rozwiązanie tej równości na siatkach, przy zadanym warunku brzegowym wysokości, to klasyczna zadanie relaksacyjne: ponownie zamieniaj wysokość każdego punktu wewnętrznego na średnią jego sąsiadów aż cała siatka przestanie się zmieniać.

// 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];
}

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Soap Bubble — Minimal Surface Simulation i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.