VisheshVerma

Report
Pattern Formation in Physical Systems: From a Snowflake to an Air Bubble
Student: Vishesh Verma, Stevenson High School
Faculty Advisor: Professor Shuwang Li, Applied Math Department, IIT Chicago
Crystal Growth Problem
Crystal growth problem Is a classical example of a phase transformation
from the liquid phase to the solid phase via heat transfer, for example the
formation of snowflakes
Hypothesis
Conclusion
If we change the air pumping rate and the initial
shape, then we will get a variety of patterns.
The patterns found in physical structures depend
on
D
1. The initial shape
Variables
A.
Snowflakes are formed around a nucleus, such as a dust particle
B.
The snowflake forms when water VAPOR freezes around the nucleus, not from frozen
rain
C.
Snowflakes have six sides/branches because the water molecules form hexagonal
crystals when frozen
D.
Snowflakes have different parts (hexagonal pieces and branches) because of the
Mullins-Sekerta Instability
-The Mullins-Sekerta Instability explains how a not-so-smooth part of a crystal
grows into a noticeable bump/branch in the crystal
2. The air pumping rate or other similar driving
forces
Air flow rate, initial shape, number of computational
points, symmetry modes
3. The symmetry of the anisotropy mode
Materials
Implications
Computer cluster (for computation)
1.
Matlab (for post processing data)
Our research project has to understand the
formation of snowflakes and explains why
The shape of the snowflake is affected mainly by the air temperature and the amount of water in
the air (supersaturation)
-All snowflakes are different because no two snowflakes take the exact same path through
the atmosphere, leading them to experience slightly different conditions which make their
appearances different
snowflakes are all different.
Procedure
2.
Our research project also offers insights to
the oil-recovery application in petroleum
All branches grow the same way because “the local conditions are essentially the same for each
arm on a tiny snow crystal” – Kenneth G. Libbrecht
engineering.
INPUT
2 types of forces shape crystal growth
-Micro: surface tension (inversely related to curvature), kinetics
k = f’’(x) / [1+(f’(x))2]3/2
k=1/R
k means curvature, R = radius
-Macro: temperature, humidity
Driving forces for crystal growth
 Air flow rate
 Initial shape
 Number of
computational points
 Physical parameters
References
1. S. Li, J. Lowengrub, P. Leo, A rescaling scheme with
Purpose
application to the long time simulation of viscous
1. To understand the underlying physical mechanism governing the pattern formation
2. To design a method such that a crystal evolves to a desired or predetermined shape
3. To achieve this goal, physicists use a Hele-Shaw cell to study the pattern formation
fingering in a Hele-Shaw cell, Journal of Computational
TOOLBOX
Physics 225, p 554-567, 2007.
2. K. Libbrecht, Morphogenesis on Ice: The Physics of
Snow Crystals, Engineering and Science, p 10-19, 2001.
3. J. Adam, Flowers of Ice-Beauty, Symmetry, and
Complexity: A Review of The Snowflake: Winter’s
A Sister Problem: Hele Shaw bubbles
Secret Beauty, Notices of AMS 52, p 402-416, 2005.
4. E. Ben-Jacob, P. Garik, The formation of patterns in
non-equilibrium growth, Nature 343, p 523-530, 1990.
Results Explanation
Hele Shaw bubble
Set up: two plates with very little gap, viscous liquid (ex. oil) in
between
Drill hole in top plate, put tube in hole, and force air through
Forms a bubble in a shape representing a snow crystal (not
necessarily 6 sides though)
Viscous fingering pattern due to Saffman-Taylor instability
Snow crystals : Mullins-Sekerta Instability :: Hele Shaw
bubble : Saffman-Taylor instability
Side View
OUTPUT
 Different shapes (made
of a collection of points)
 Different fingering
patterns
In the first column, the initial
shape was a four fold star
shape.
We apply a time
decreasing pumping rate. We
observed a tip splitting. The
symmetry mode is four.
In the second column, we
started from a circle with a six
fold symmetry. The pumping
rate is constant
In the third column, we started
from a circle with four fold
symmetry. The pumping rate
was also constant.
Acknowledgements
Vishesh would like to thank the SPARK program
for this opportunity, the computer cluster resource
at IIT, and the Matlab tutorial at the Applied Math
Department of IIT.
S. Li would like to thank the National Science
Foundation for supporting the related research
projects.

similar documents