A soap bubble turns surface area into visible geometry. Its film pulls inward, pressure holds it open, curvature connects the two, and gravity, drainage and airflow disturb the ideal shape. Join bubbles and flat-looking films meet at characteristic angles; stretch a film across a wire frame and it searches for a low-area surface. Mathematics explains both the elegance and the limits.
The experiments in this article are educational. Use a teacher-approved mild soap solution, protect eyes, do not ingest it, wash hands after use and wipe slippery spills immediately. Do not use pressurised gas, flames or unknown chemicals. The safest investigation uses small hand-blown films, photographs and calculations.
- Begin with surface energy
- Calculate Laplace pressure
- Explore minimal films
- Understand Plateau borders
- Plan a careful investigation
- Use the FAQ
Why mathematics matters to a soap film
Soap films connect several branches of mathematics. Geometry supplies area, volume and curvature. Calculus describes how area changes when a surface moves. Differential equations express equilibrium. Ratios compare competing forces. Graphs and image analysis turn photographs into measurements. Statistics distinguishes a trend from the scatter of short-lived bubbles.
The first modelling decision is whether the object is a single liquid–gas interface, a soap bubble with two interfaces, or a film on a frame. A liquid droplet has one outer interface. A soap bubble is a thin liquid layer with an inner and an outer gas–liquid interface. A free film between a wire boundary also has two sides, although its pressure conditions can differ.
Confusing these cases creates a factor-of-two error in energy and pressure. Good mathematical reasoning starts with a diagram that labels phases and interfaces before any formula is used.
An observation route
- Look at the boundary: closed bubble, open film or droplet?
- Count interfaces and state the pressure on each side.
- Measure a length scale such as radius or frame width.
- Choose an ideal model: sphere, circular arc, catenoid or numerical surface.
- Calculate area, curvature, pressure or a dimensionless ratio.
- Compare with evidence and list gravity, drainage, airflow and optical effects omitted.
This route keeps a beautiful photograph from becoming an overconfident physical claim.
Surface tension turns area into energy
Surface tension γ has units newtons per metre. It can also be read as joules per square metre because 1 N/m = 1 J/m². In a simple constant-γ model, increasing one interface by area ΔA requires surface energy ΔE = γΔA.
A soap film has two interfaces, so an area A measured across its frame corresponds approximately to surface energy E = 2γA, ignoring the frame edge and thickness variation. A spherical soap bubble of radius R has area 4πR² on each surface and total interfacial area about 8πR². Its surface energy is therefore E = 8πγR² in the thin-film approximation.
This energy interpretation explains why a free film tends to reduce area while respecting its boundary and enclosed volume constraints. “Soap wants the smallest area” is a useful shortcut only after those constraints are stated. A film cannot ignore a wire frame; a closed bubble cannot discard its gas volume instantly.
Area and volume of a sphere
A sphere of radius R has area A = 4πR² and volume V = 4πR³/3. The area-to-volume ratio is A/V = 3/R. Small spheres have more surface area per unit volume than large ones.
If radius doubles, area grows by four and volume by eight. This scaling affects surface energy, pressure and drainage. It also explains why “twice as wide” is not “twice as much bubble”.
Among closed surfaces enclosing a fixed volume, the sphere minimises area. That is an isoperimetric result. It helps explain why a small unconstrained bubble in still conditions is nearly spherical: lower area means lower surface energy for a given volume.
Did You Know? The sphere’s area-minimising property and the pressure–curvature law are related but not identical statements. One is a global optimisation under a volume constraint; the other is a local force balance.
Curvature and the Young–Laplace relation
For one interface with principal radii of curvature R₁ and R₂, the pressure jump is ΔP = γ(1/R₁ + 1/R₂) under a common sign convention. The sum is twice the mean curvature.
For a spherical droplet, R₁ = R₂ = R, so ΔP = 2γ/R. A soap bubble has two interfaces. If both have nearly the same radius and surface tension, the pressure difference between inside and outside is approximately ΔP = 4γ/R.
The factor four is not a new law; it is two spherical interfaces, each contributing about 2γ/R. If a problem says “bubble” without specifying whether it means a gas bubble in liquid, a liquid droplet in gas or a soap bubble film, ask for the physical arrangement.
The relation predicts that smaller radius means larger pressure difference for fixed γ. Curvature, not volume alone, governs the local pressure jump.
Unit check
γ/R has units (N/m)/m = N/m², which is a pascal. If radius is entered in millimetres without conversion, the result will be wrong by a factor of one thousand. Units are especially important because bubble radii are often measured in millimetres while γ is reported in SI units.
Worked example: pressure in small and large bubbles
Take an illustrative surface tension γ = 0.030 N/m. For a soap bubble of radius R = 10 mm = 0.010 m, ΔP = 4×0.030/0.010 = 12 Pa.
For R = 2 mm = 0.002 m, ΔP = 4×0.030/0.002 = 60 Pa. The radius is one fifth, so the pressure difference is five times larger. This inverse relationship is easy to predict before calculating.
For a single-interface liquid droplet of the same radius and γ, the ideal pressure difference would be half those values: 6 Pa at 10 mm and 30 Pa at 2 mm. Always count interfaces.
Surface energy of the 10 mm bubble
Using E = 8πγR², E = 8π×0.030×0.010² ≈ 7.54×10⁻⁵ J. If radius doubles while γ stays constant, surface energy increases fourfold.
This energy is small in everyday terms, yet it dominates the bubble’s shape at small scales. Comparing energy magnitude alone with household energy use would miss the relevant force and length scales.
Two connected bubbles
In the simplest model, a smaller bubble has higher internal pressure than a larger one. If they are connected and surface tension is comparable, gas tends to flow from the smaller to the larger. The small bubble shrinks and its ideal pressure rises further, reinforcing the change.
Real foams can be stabilised by surfactant transport, diffusion limits, geometry and liquid structure. The calculation gives the direction of an ideal tendency, not a complete lifetime prediction.
The NASA educational activity Surface Tension and Bubble Stability uses bubble observations to examine surface tension and notes the higher pressure associated with smaller bubbles. It is a useful primary educational source for connecting measurement with the model.
Minimal surfaces and mean curvature
A minimal surface is a surface that locally minimises area under its boundary conditions. Mathematically, its mean curvature is zero at regular points. Zero mean curvature does not mean the surface is flat. The principal curvatures can be equal in magnitude and opposite in sign, as on a saddle-like patch.
A film spanning a single closed planar loop can be flat because the flat disk has low area. Give the boundary a non-planar shape and the film can become a curved minimal surface. The exact result depends on the wire boundary.
If pressure on both sides of a film is equal, the Young–Laplace pressure jump is zero, so the mean curvature is zero under the ideal constant-tension model. A closed bubble has nonzero pressure difference and positive mean curvature; it is not a minimal surface in the same unconstrained sense. Its sphere minimises area subject to fixed volume.
Local versus global minimum
A surface can be locally stable yet not have the smallest area among every imaginable spanning surface. Boundary geometry may allow several stationary surfaces. As a frame changes, one form can lose stability and jump to another.
This distinction is important in optimisation. A calculus condition such as zero first variation identifies a stationary candidate. A stability or second-variation analysis asks whether nearby changes increase area. A global proof compares all allowed competitors.
NASA’s Minimal Surfaces in a Low-Gravity Environment connects liquid films, low gravity and the classical isoperimetric problem. Reduced gravity is useful because it weakens sagging relative to surface forces, making ideal geometric behaviour easier to study.
The catenoid between two rings
Stretch a soap film between two equal parallel circular rings. One possible ideal surface of revolution is a catenoid. Its radial profile can be written r(z) = a cosh(z/a), where a is the minimum neck radius in an appropriate coordinate system.
The name resembles “catenary” because rotating a catenary about its axis produces a catenoid. Yet a hanging chain and a soap film solve different optimisation problems: the chain balances weight and tension along a curve; the catenoid minimises surface area under ring boundaries.
Not every ring spacing supports a stable catenoid for a given radius. As rings are separated, the neck becomes narrow and a stability limit is reached. The film can collapse into two nearly flat disks. This is a physical way to see that an equation may describe a stationary shape even where that branch is not stable.
Area as an integral
For a surface of revolution r(z), area is A = 2π∫r√[1+(dr/dz)²] dz. Substituting r = a cosh(z/a) gives dr/dz = sinh(z/a), and √[1+sinh²] = cosh because 1+sinh²u = cosh²u.
The integrand becomes 2πa cosh²(z/a). Calculus turns a three-dimensional surface into a one-dimensional integral. Students need not evaluate every expression to appreciate the modelling chain: profile, derivative, surface element, integral, boundary condition.
Why films meet at special angles
In an ideal dry foam with equal surface tension, three films meet along a Plateau border at angles of 120°. The vector sum of three equal tension forces separated by 120° is zero. Draw the vectors head to tail and they form an equilateral triangle.
Four Plateau borders can meet at a vertex at the tetrahedral angle, approximately 109.47°. This geometry distributes directions symmetrically in three dimensions. These are ideal equilibrium rules; gravity, liquid content, unequal tensions and boundaries can modify observed geometry.
Vector proof of 120°
Let three unit vectors have directions 0°, 120° and 240°. Their x components sum to 1 − 1/2 − 1/2 = 0. Their y components sum to 0 + √3/2 − √3/2 = 0. Multiplying each by the same tension magnitude preserves the zero resultant.
If one film has different surface tension, equal 120° spacing no longer balances the vectors. The force triangle must use side lengths proportional to the three tensions. Angles therefore provide evidence about symmetry of the local interfacial forces, not merely a decorative pattern.
Two-dimensional bubble clusters
When bubbles are confined between close transparent plates, their edges can be modelled approximately as circular arcs meeting at vertices. Pressure differences determine arc curvature through a two-dimensional analogue of Laplace balance. Image analysis can fit circles to edges and test whether higher-pressure cells curve toward lower-pressure neighbours.
Plate confinement, wet borders and contact angles complicate the model. A fitted radius is a measurement with uncertainty, not a direct proof of a pressure value unless surface tension and geometry are known.
Gravity, capillary length and the Bond number
Surface effects dominate small bubbles; gravity becomes more noticeable as size grows. The Bond number compares gravity with surface tension: Bo = ΔρgL²/γ, where Δρ is density difference, g is gravitational acceleration and L is a characteristic length.
If Bo is much less than one, surface tension tends to dominate shape. If Bo is much greater than one, gravity is important. Around one, both matter. The choice of L and geometry should be stated.
The capillary length ℓc = √[γ/(Δρg)] is the scale at which Bo is about one. With γ = 0.030 N/m, Δρ ≈ 1000 kg/m³ and g = 9.81 m/s², ℓc ≈ 0.00175 m, or 1.75 mm.
This does not mean every bubble larger than 1.75 mm becomes non-spherical. A thin soap bubble involves gas and a liquid film rather than a bulk droplet, and pressure, drainage and constraints matter. The calculation supplies a comparison scale, not a sharp universal boundary.
Scaling insight
Bond number grows with L². Doubling characteristic length multiplies Bo by four. Small demonstrations can therefore look more surface-dominated than large installations even when they use the same liquid.
This is the same logic used in engineering similarity: preserve relevant dimensionless groups if a model is meant to represent a larger system. Simply enlarging every dimension does not preserve force balance.
Drainage, evaporation and colourful films
A soap film is not a static mathematical sheet. Liquid drains downward under gravity, while capillary suction and viscous resistance redistribute it. Evaporation thins the film. Surfactant concentration changes locally as the surface stretches.
Surface-tension gradients produce Marangoni stresses that can pull liquid along the interface. These stresses help a film resist local thinning in some conditions. A constant-γ model omits that active redistribution.
The changing colours of a thin film arise from optical interference. Light reflected from the front and back surfaces travels different optical paths. Film thickness, refractive index, wavelength and viewing angle determine constructive or destructive interference. Dark regions can indicate very thin film, but colour is not a direct thickness ruler without an optical model.
The mathematics therefore links fluid mechanics with waves. A time-lapse colour pattern can reveal flow qualitatively, yet extracting thickness requires calibration and careful illumination.
Related reading on rainbows, refraction and colour geometry shows another way that wavelength and geometry create visible patterns.
A safe measurement project
Choose one modest question, such as how estimated bubble pressure varies with radius under a fixed assumed surface tension. Blow bubbles gently through a commercial wand or teacher-approved straw procedure without inhaling solution. Photograph them against a scale from a fixed distance.
Use image software to estimate diameter. Repeat at least five bubbles in several size bands. Convert millimetres to metres, calculate ΔP = 4γ/R and plot pressure against 1/R. The ideal model predicts a straight line through the origin with slope 4γ.
This project calculates pressure from radius and assumed γ; it does not independently measure pressure. Label the graph accordingly. To estimate γ experimentally would require an additional calibrated method.
Uncertainty from radius
Because P = 4γ/R, relative sensitivity to radius is approximately δP/P = δR/R in magnitude for small errors. If R = 5.0 ± 0.2 mm, radius uncertainty is about 4%, so calculated pressure has about 4% radius-derived uncertainty before uncertainty in γ is included.
For independent uncertainties, relative uncertainty can be combined approximately by root-sum-square. For conservative bounds, magnitudes may be added. State which interpretation is being used.
Image-analysis checks
- Keep the camera perpendicular to the scale plane to reduce perspective error.
- Place scale and bubble at nearly the same depth.
- Fit a circle to several edge points rather than measuring one noisy diameter.
- Record whether the bubble is obviously distorted or touching a surface.
- Keep original images and note exclusions before seeing the final trend.
A second project photographs three-film junctions in a shallow foam and measures angles. Average repeated junctions and plot a distribution around 120°. Discuss why wet borders, perspective and boundary effects create scatter.
Numerical models and optimisation
Complex wire frames rarely produce a surface described by one elementary equation. A computer can approximate the film using a mesh of triangles. The algorithm moves vertices to reduce total triangle area while holding boundary vertices fixed.
The area of a triangle with edge vectors a and b is ½|a×b|. Sum that quantity over the mesh. A gradient-based method estimates how area changes when each free vertex moves, then updates the vertices in a direction that lowers area.
Mesh quality matters. Large or skewed triangles can distort curvature estimates. A numerical solution should be refined and checked for convergence. Different initial meshes may settle into different local minima.
This is an age-appropriate entrance to computational mathematics. The computer is not “finding the answer by magic”; it is executing an objective function, constraints, representation and stopping rule chosen by people.
Discrete versus continuous claims
A mesh minimises area within its finite representation. The physical film is continuous and also subject to gravity, thickness and material effects. Agreement between simulation and photograph is evidence for the model over that range, not proof that every omitted mechanism is absent.
NASA’s bubble-shape measurements compared with profiles derived from the Laplace formula provide a useful example of testing mathematical profiles against observed geometry rather than treating a curve as self-validating.
Common misconceptions
“Every soap film is a minimal surface.” A film with equal pressure on both sides can approximate a zero-mean-curvature surface. A closed bubble has a pressure difference and minimises area subject to volume, so its mean curvature is not zero.
“Small bubbles have lower pressure because they contain less air.” Under the ideal Laplace model, smaller radius gives a larger pressure difference for fixed surface tension.
“Surface tension is a skin with a fixed strength.” Surface tension is interfacial free energy per area or force per length. Surfactant concentration and dynamics can change it; the film is a flowing liquid layer.
“A sphere has the smallest area, full stop.” A sphere minimises area among closed surfaces with a fixed enclosed volume. Without a volume constraint, shrinking toward zero would reduce area.
“A colourful film has coloured soap.” The colours commonly arise from interference related to film thickness and viewing geometry.
“A simulation proves the physical bubble.” It proves what follows from its equations and numerical representation. Experiments test whether those assumptions are adequate.
Four worked extensions
Extension 1: one large bubble versus eight small bubbles
Suppose one spherical bubble of radius R is replaced by eight bubbles, each enclosing one eighth of the original volume. Because volume scales with radius cubed, each smaller radius is R/2.
The original bubble has one-side area 4πR². Each small bubble has area 4π(R/2)² = πR², and eight have total area 8πR²—twice the original. For soap bubbles, both sides multiply both totals equally, so the factor remains two.
At constant γ, total surface energy doubles. This explains why creating a fine foam requires interfacial energy. Surfactant does not remove the energy requirement; it lowers surface tension and helps stabilise the new interfaces.
The pressure difference in each small bubble is twice that in the original because radius halves. Fragmenting volume therefore increases both total area and ideal Laplace pressure.
Extension 2: work done during slow expansion
For an ideal soap bubble E = 8πγR². Differentiating gives dE/dR = 16πγR. A small radius increase dR requires surface-energy change dE ≈ 16πγR dR when γ is constant.
The volume change is dV/dR = 4πR². The energy change per volume change is (dE/dR)/(dV/dR) = 4γ/R, exactly the soap-bubble pressure difference. This calculus derivation connects energy minimisation with local pressure balance.
The result assumes reversible, quasi-static expansion, constant γ and negligible dissipation. Blowing a real bubble involves airflow, viscosity and changing surfactant concentration, so lung effort is not estimated by this interfacial term alone.
Extension 3: balancing unequal tensions
Suppose three film tensions meeting in a plane have magnitudes 30, 30 and 40 arbitrary units. They cannot meet at equal 120° angles. Treat the tensions as sides of a closed force triangle.
The angle between the two 30-unit force vectors needed so their resultant has magnitude 40 satisfies 40² = 30²+30²+2(30)(30)cos θ. Thus cos θ = (1600−1800)/1800 = −1/9 and θ ≈ 96.38° between those force directions.
The corresponding geometric film angle depends on the force-vector orientation convention, but the central lesson is firm: angle changes encode unequal tensions. A measured non-120° junction is not automatically “bad data”; it may indicate boundaries, dynamics or unequal interfaces.
Extension 4: a lifetime dataset
Imagine bubble lifetimes in seconds: 18, 22, 25, 27, 31, 35, 41, 49 and 110. The mean is raised strongly by the 110 s bubble, while the median is 31 s. A box plot or all individual points communicates the skew better than mean alone.
Lifetime is often variable because film thickness, humidity, contamination and airflow differ. If comparing two solutions, randomise order, use the same wand, control location and collect enough repeated observations. Censor bubbles that remain when observation ends rather than recording them as if they popped at the stopping time.
Survival analysis is designed for time-to-event data with censoring. Even if students do not perform a formal survival model, recognising the data type prevents a biased average.
From two bubbles to a foam network
A foam divides space into cells while trying to manage total interfacial area under volume constraints. In a two-dimensional idealisation, topology offers simple accounting. Let V be vertices, E edges and F faces in a planar network; Euler-type relations connect them once boundaries are defined.
If three edges meet at most interior vertices, counting edge ends gives approximately 3V ≈ 2E away from boundary corrections. Combined with Euler’s relation, this helps explain why six-sided cells are common on average in large planar foams, even though individual cells can have more or fewer sides.
This does not say every bubble is a hexagon. Cell areas vary, boundaries intervene and the network rearranges. When an edge shrinks, a neighbour-switching event can change which cells touch while preserving a force-balanced local geometry.
In three dimensions, the problem is substantially harder. Equal-volume partitions and minimal-area foams lead to rich geometric questions. Famous candidate structures show that an intuitively simple arrangement is not always the best. Mathematics advances by proposing competitors, calculating objective values and testing proofs.
Coarsening over time
Gas can diffuse through films from higher-pressure small bubbles toward lower-pressure large ones. Average bubble size grows while bubble count falls, a process called coarsening. The Laplace pressure relation supplies the driving tendency; diffusion and film properties set the rate.
A time-lapse image sequence can estimate cell areas. Segmenting every bubble is difficult where films merge, glare appears or bubbles leave the frame. State the segmentation rule and manually audit a sample. Automated image analysis is only as credible as its error assessment.
Plot a characteristic size against time on ordinary and logarithmic axes. A straight line on a log–log plot suggests a power law over that range, but it does not prove the same exponent holds forever. Early drainage and late boundary effects can create different regimes.
Surface mathematics in careers and education
The mechanisms in a soap film reappear in emulsion stability, sprays, coatings, porous materials, microfluidics and biological interfaces. Researchers may combine differential geometry with chemistry and fluid mechanics. Process engineers use dimensionless groups and transport models. Imaging specialists quantify curvature and thickness.
At school level, the goal is not to declare a career from one enjoyable experiment. It is to recognise transferable tools: scale analysis, constrained optimisation, vectors, uncertainty and computational meshes. Those tools keep options open across mathematics, physics, chemistry, materials science and engineering.
Singapore students should verify current subject and programme requirements on official institution pages when choosing pathways. A bubble project can demonstrate curiosity and careful investigation, but it does not replace the broader prerequisites for a course or profession.
The strongest portfolio record is an honest one: a clear question, safe method, original data, equations with units, uncertainty, source links and a conclusion limited to the evidence. Beautiful photographs are welcome, but they support rather than substitute for reasoning.
Keeping a dated laboratory notebook also makes unexpected results useful. A failed film or blurred image can reveal a procedural weakness that the next trial corrects.
How students can develop the mathematics
Begin with circles and spheres: radius, area, volume and unit conversion. Use proportional reasoning to predict how area, volume and pressure change when radius doubles. Then distinguish a droplet’s one interface from a soap bubble’s two.
Move to vectors at a three-film junction. Resolve three equal tensions at 120°. Use a spreadsheet to plot ΔP against 1/R and surface energy against R². Add uncertainty bars.
More advanced learners can study mean curvature, the catenoid, calculus of variations, the Bond number and triangular meshes. The progression is valuable because every new idea answers a visible question.
Parents can support curiosity with safe observation: Why is a free bubble nearly round? Why do small bubbles disappear into larger ones in a foam? Why are the colours changing? Encourage the learner to separate observation, calculation and interpretation.
For another measurement pathway, particle counters, sampling volume and uncertainty shows how a small observation becomes a population estimate. The mathematics is different, but the discipline of units and uncertainty is the same.
Frequently asked questions
Why is mathematics important for soap bubbles?
It connects surface area to energy, curvature to pressure, boundary geometry to minimal surfaces, film junctions to vector balance, and photographs to quantitative tests. It also distinguishes ideal rules from real drainage and evaporation.
Why is a bubble spherical?
For a fixed enclosed volume and uniform surface tension, a sphere has the least area and therefore low interfacial energy. Gravity, airflow, contact and motion can distort it.
Why is pressure higher in a smaller bubble?
For an ideal soap bubble, ΔP = 4γ/R. With γ fixed, reducing R increases pressure inversely. This is a pressure difference between inside and outside.
Why does a soap bubble use 4γ/R instead of 2γ/R?
A soap film has inner and outer interfaces. Each approximately spherical interface contributes 2γ/R in the thin-film model, giving a total near 4γ/R.
Is a bubble a minimal surface?
Not in the zero-mean-curvature sense. A closed bubble has a pressure difference and nonzero mean curvature. It minimises area subject to enclosing volume. A frame-spanning film with equal pressure on both sides can approximate a minimal surface.
What is a catenoid?
It is a minimal surface of revolution generated by rotating a catenary. A soap film between suitable circular rings can take this form, within a stability range.
Why do three films meet at 120 degrees?
If their surface tensions are equal, three force vectors balance when separated by 120°. Unequal tensions or external constraints alter the angles.
Can colour tell me exact film thickness?
Only with a calibrated interference model including wavelength, refractive index, incidence and reflection phase. Colour alone is qualitative.
What makes a bubble pop?
Film thinning, evaporation, contamination, contact, airflow and disturbances can contribute. The simple pressure and area models do not predict an exact lifetime.
Is bubble mathematics used beyond play?
Related ideas appear in foams, emulsions, capillarity, microfluidics, materials and optimisation. Career relevance depends on broader science and engineering training; one topic does not guarantee a pathway.
A small film, a large mathematical world
Soap bubbles are optimistic mathematics: a child’s familiar object opens into geometry, pressure, energy, calculus and computation. The formulas are compact, but the reasoning is spacious. Count interfaces, name constraints, check units, measure uncertainty and compare the ideal surface with the changing film.
That combination is the real benefit of learning mathematics. It allows wonder to become a testable question without draining away the wonder. A bubble can remain beautiful while also teaching that every elegant model earns its strength from clearly stated assumptions.