Math of Popping Bubbles in a Foam
Berkeley Lab researchers mathematically describe the complex evolution and disappearance of foamy bubbles
May 9, 2013
Bubble baths and soapy dishwater, the refreshing head on a beer and the luscious froth on a cappuccino. All are foams, beautiful yet ephemeral as the bubbles pop one by one. Now, two researchers from the Department of Energy’s (DOE’s) Lawrence Berkeley National Laboratory (Berkeley Lab) and the University of California, Berkeley have described mathematically the successive stages in the complex evolution and disappearance of foamy bubbles, a feat that could help in modeling industrial processes in which liquids mix or in the formation of solid foams such as those used to cushion bicycle helmets.
Applying these equations, they used supercomputers at DOE’s National Energy Research Scientific Computing Center (NERSC) to create mesmerizing computer-generated visualization showing the slow and sedate disappearance of wobbly foams one burst bubble at a time.
The applied mathematicians, James A. Sethian and Robert I. Saye, will report their results in the May 10 issue of Science. Sethian leads the Mathematics Group at Berkeley Lab and is a professor of Mathematics at UC Berkeley. Last month, Sethian was elected to the National Academy of Sciences. Saye will join Berkeley Lab’s Mathematics Group this September as a Luis W. Alvarez Fellow in Computing Sciences, and will graduate from UC Berkeley this May with a PhD in applied mathematics.
“This work has application in the mixing of foams, in industrial processes for making metal and plastic foams, and in modeling growing cell clusters,” said Sethian. “These techniques, which rely on solving a set of linked partial differential equations, can be used to track the motion of a large number of interfaces connected together, where the physics and chemistry determine the surface dynamics.”
The problem with describing foams mathematically has been that the evolution of a bubble cluster a few inches across depends on what’s happening in the extremely thin walls of each bubble, which are thinner than a human hair.
“Modeling the vastly different scales in a foam is a challenge, since it is computationally impractical to consider only the smallest space and time scales,” Saye said. “Instead, we developed a scale-separated approach that identifies the important physics taking place in each of the distinct scales, which are then coupled together in a consistent manner.”
Saye and Sethian discovered a way to treat different aspects of the foam with different sets of equations that worked for clusters of hundreds of bubbles. One set of equations described the gravitational draining of liquid from the bubble walls, which thin out until they rupture. Another set of equations dealt with the flow of liquid inside the junctions between the membranes. A third set handled the wobbly rearrangement of bubbles after one pops. Using a fourth set of equations, the mathematicians solved the physics of a sunset reflected in the bubbles, taking account of thin film interference within the bubble membranes, which can create rainbow hues like an oil slick on wet pavement. Then they used NERSC’s Hopper Systems to solve the full set of equations of motion.
“Solving the full set of equations on a desktop computer would be time-consuming. Instead, we used massively parallel computers at NERSC and computed our results in a matter of days," said Saye.
“Foams were a good test that all the equations coupled together,” Sethian said. “While different problems are going to require different physics, chemistry and models, this sort of approach has applications to a wide range of problems.”
The mathematicians next plan to look at manufacturing processes for small-scale new materials.
“DOE’s longstanding support for core basic applied mathematics has been instrumental in providing the opportunities to develop the mathematics and algorithms behind this work,” said Sethian.
The work is supported by the Department of Energy’s Office of Science, the National Science Foundation, and the National Cancer Institute.
- Multiscale Modeling of Membrane Rearrangement, Drainage, and Rupture in Evolving Foams (May 10, 2013 Science)
- Sethian’s research web site
The Lawrence Berkeley National Laboratory (Berkeley Lab) Computing Sciences organization provides the computing and networking resources and expertise critical to advancing the Department of Energy's research missions: developing new energy sources, improving energy efficiency, developing new materials and increasing our understanding of ourselves, our world and our universe.
ESnet, the Energy Sciences Network, provides the high-bandwidth, reliable connections that link scientists at 40 DOE research sites to each other and to experimental facilities and supercomputing centers around the country. The National Energy Research Scientific Computing Center (NERSC) powers the discoveries of 6,000 scientists at national laboratories and universities, including those at Berkeley Lab's Computational Research Division (CRD). CRD conducts research and development in mathematical modeling and simulation, algorithm design, data storage, management and analysis, computer system architecture and high-performance software implementation. NERSC and ESnet are DOE Office of Science User Facilities.
Lawrence Berkeley National Laboratory addresses the world's most urgent scientific challenges by advancing sustainable energy, protecting human health, creating new materials, and revealing the origin and fate of the universe. Founded in 1931, Berkeley Lab's scientific expertise has been recognized with 13 Nobel prizes. The University of California manages Berkeley Lab for the DOE’s Office of Science.