Why is it that whenever balloons are inflated they converge towards the shape of a sphere regardless of their initial geometry?

On one level this may be a purely geometrical problem due to thermodynamic constraints on a body with finite surface area. This suggests the necessity of solving a global optimisation problem. However, the sequence of deformations undergone may be facilitated by the elastic material the balloons are made of.

In this article I consider the contribution of the latter by analysing the problem in two dimensions and demonstrate that a minimal surface may be entirely due to local mechanical instabilities.

The role of material properties:

Let’s consider an object that is only allowed to extend in one dimension. If you were to elongate such an object it would assume a roughly cylindrical shape.

It follows that we must pay careful attention to the material properties of the balloon.

Reasonable assumptions:

A two-dimensional balloon is essentially an elastic loop that initially has perimeter of length:

\begin{equation} \lvert \partial \mathcal{B}(t=0) \rvert = l_0 \end{equation}

Furthermore, we may make the following reasonable assumptions:

  1. The balloon contains an astronomical number of gas particles that collectively satsify the ideal gas equation.

  2. The balloon is surrounded by a heat bath.

  3. The balloon itself is made of elastic filaments that form a Hamiltonian circuit.

Furthermore, we may assume that the mechanical behaviour of the balloon is largely driven by energy-minimisation processes that I shall detail in the next couple sections.

Isobaric inflation as a consequence of energy minimisation:

If we consider the force required to elongate the elastic boundary of the balloon we may define an associated potential energy function:

\begin{equation} U((\lvert \partial \mathcal{B}(t) \rvert - l_0)^2) \geq 0 \end{equation}

such that:

\begin{equation} U(\cdot) = 0 \iff \lvert \partial \mathcal{B}(t) \rvert = l_0 \end{equation}

Now, if we consider that physical systems tend to minimise potential energy we may infer that the balloon would tend to increase in volume without increasing , the length of its perimeter.

In the case of inflation, after accumulating a pressure difference with respect to its environment the evolution of would be guided by an approximately isobaric process provided that :

\begin{equation} PV = nRT \end{equation}

\begin{equation} \frac{\Delta V}{V} = \frac{\Delta T}{T} \end{equation}

We can go further with this type of reasoning. Not only does the elastic membrane constrain the type of thermodynamic processes that is likely to guide inflation; it also constrains the mechanism for modifying the geometry of the balloon.

Local deformations of elastic filaments lead to minimal surfaces:

If we assume that the balloon constrains an ideal gas that my be modelled as an astronomical number of Newtonian particles, it’s reasonable to suppose that equal pressure is applied to equal areas. Now, if this is the case we may consider pressure-driven deformations of that exploit a local mechanism that is operational everywhere on the boundary. What might such a mechanism look like?

Under a coarse-grained approximation, the boundary consists of a large chain of cylindrical elastic rods. If each individual rod is much larger than the characteristic length where bending occurs any amount of bending will guarantee tensile stress.

It follows that the elastic membrane will try, as much as possible, to increase the enclosed volume while minimising the elongation globally. This global minimisation happens by minimising the bending angle locally. No global coordination is required.

Another way of understanding this process is that deformations of the elastic membrane are mainly driven by local mechanical instabilities that lead to a global minimisation of potential energies.

A polygonal approximation to two-dimensional elastic boundaries:

One approach to modelling the activity of elastic boundaries is to approximate them as polygons with sides of equal length where is large. Given that the sum of the interior angles must add up to we may define the potential energy:

\begin{equation} U = \frac{1}{2} \sum_{i=1}^N (\theta_i - \pi \cdot \big(\frac{N-2}{N}\big))^2 \end{equation}

\begin{equation} \sum_{i=1}^N \theta_i = \pi \cdot (N-2) \end{equation}


\begin{equation} \frac{\partial U}{\partial \theta_i} = \theta_i - \pi \cdot \big(\frac{N-2}{N}\big) \end{equation}

\begin{equation} \Delta \theta_i \propto \frac{\partial U}{\partial \theta_i} \end{equation}

and we find that if we choose the local update with :

\begin{equation} \begin{split} \theta_i^{t+1} & = \theta_i^{t} - \Delta \theta_i \\ & = \theta_i^{t} - \lambda \frac{\partial U}{\partial \theta_i} \\ & = (1-\lambda) \cdot \theta_i^t + \lambda \cdot \pi \cdot \big(\frac{N-2}{N}\big) \end{split} \end{equation}

and we can show that very quickly since:

\begin{equation} x_{n+1} = (1-\lambda) \cdot x_n + \lambda \cdot \alpha \implies x_{n+1} - \alpha = (1-\lambda) \cdot (x_n - \alpha) \end{equation}

\begin{equation} \frac{(x_{n+1}-\alpha)^2}{(x_n - \alpha)^2} = (1-\lambda)^2 \end{equation}

so if we define:

\begin{equation} \epsilon_{n+1}^2 = (x_{n+1}-\alpha)^2 \end{equation}

\begin{equation} \epsilon_{n}^2 = (x_{n}-\alpha)^2 \end{equation}

we find that:

\begin{equation} \lim_{n \to \infty} \epsilon_{n+1}^2 = \epsilon_1^2 \cdot \prod_{n=1}^\infty \frac{\epsilon_{n+1}^2}{\epsilon_{n}^2} = \lim_{n \to \infty} \epsilon_1^2 \cdot (1-\lambda)^{2n} = 0 \end{equation}

so we have exponentially fast convergence to a spherical geometry.


In this article I propose the existence of a local mechanical instability present everywhere in a closed elastic membrane with aspherical geometry. Surprisingly, the net action of this instability leads to exponentially fast convergence to the global minimum. But, this analysis may be refined.

The above analysis is entirely based on phenomenological studies of rubber bands by alternately dropping and manipulating rubber bands on a table. This led me to a useful phenomenological model which may help simulate kinematics but doesn’t approximate the forces involved i.e. dynamics.