Research: Surface optimization framework for foam networks

[pjt222]

Summary

Implement minimal surface solver using variational methods for complex foam network topology.

References

• Meng 2026 - Surface optimization governs the local design of physical networks (Nature)
• Witten 1986, Carlip 1988, Saadi & Zwiebach 1989 (mathematical framework)
• Tong lectures (string theory background)

Background

This is a research-level extension exploring the mathematical connection between soap films and string theory worldsheets.

Key Insight from Meng 2026

• Physical networks (including soap film networks) minimize surface area
• Bifurcations naturally occur at 120° angles (Plateau's rules)
• Network topology transitions at critical parameter χ ≈ 0.83
• Riemannian metric tensor describes surface geometry

String Theory Connection

Both soap films and string worldsheets:

• Minimize area (surface tension ↔ string tension)
• Have similar vertex structures (3 films at 120° ↔ cubic string vertex)
• Use quadratic differentials for surface parameterization

Mathematical Framework

Surface Energy Functional

E[X] = ∫∫ √det(γ_αβ) dσ¹dσ²

Where γ_αβ is the induced metric on the surface.

Euler-Lagrange Equations

Minimizing E leads to:

∇²X = 0  (minimal surface equation)

Steiner Network Problem

Find network minimizing total edge length connecting n points.

• Solutions have degree-3 vertices at 120° angles
• Analogous to Plateau borders in foam

Proposed Implementation

Phase 1: Energy Functional

fn surface_energy(mesh: &Mesh) -> f32 {
    mesh.faces.iter().map(|f| f.area()).sum()
}

Phase 2: Gradient Descent Minimization

fn minimize_surface(mesh: &mut Mesh, iterations: usize) {
    for _ in 0..iterations {
        let gradient = compute_area_gradient(mesh);
        for (vertex, grad) in mesh.vertices.iter_mut().zip(gradient) {
            *vertex -= grad * STEP_SIZE;
        }
        // Project to constraints (boundary, volume)
    }
}

Phase 3: Network Topology

struct FoamNetwork {
    bubbles: Vec<Bubble>,
    plateau_borders: Vec<PlateauBorder>,
    vertices: Vec<PlateauVertex>,
}

fn optimize_network_topology(network: &mut FoamNetwork) {
    // T1 transitions: swap neighbor bubbles
    // T2 transitions: bubble disappearance
}

Practical Relevance

• Low for single bubble simulation
• High for complex foam rendering
• Academic interest in physics/mathematics connection

Effort

Very high (research project)

Impact

Advanced foam simulation, theoretical understanding.