Pattern Formation and Dynamics in Nonequilibrium Systems 1st Edition Michael Cross download
Pattern Formation and Dynamics in Nonequilibrium Systems 1st Edition Michael Cross download
https://ebookname.com/product/pattern-formation-and-dynamics-in-
nonequilibrium-systems-1st-edition-michael-cross/
https://ebookname.com/product/competing-interactions-and-pattern-
formation-in-nanoworld-1st-edition-elena-vedmedenko/
https://ebookname.com/product/universality-in-nonequilibrium-
lattice-systems-geza-odor/
https://ebookname.com/product/nonequilibrium-molecular-dynamics-
theory-algorithms-and-applications-1st-edition-billy-d-todd/
https://ebookname.com/product/information-systems-the-state-of-
the-field-1st-edition-john-leslie-king/
Dealing with Feeling Lucky Duck Books 2nd Edition Tina
Rae
https://ebookname.com/product/dealing-with-feeling-lucky-duck-
books-2nd-edition-tina-rae/
https://ebookname.com/product/bad-connection-1st-edition-d-s-
weissman/
https://ebookname.com/product/speaking-of-god-in-thomas-aquinas-
and-meister-eckhart-beyond-analogy-1st-edition-anastasia-
wendlinder/
https://ebookname.com/product/labour-conditions-for-construction-
decent-work-building-cities-and-the-role-of-local-
authorities-1st-edition-roderick-lawrence/
https://ebookname.com/product/holy-delight-typology-numerology-
and-autobiography-in-donne-s-devotions-upon-emergent-occasions-
kate-gartner-frost/
The Adrenal medulla 1986 1988 1st Edition Carmichael
https://ebookname.com/product/the-adrenal-medulla-1986-1988-1st-
edition-carmichael/
PAT T E R N F O R M AT I O N A N D D Y N A M I C S
IN NONEQUILIBRIUM SYSTEMS
Many exciting frontiers of science and engineering require understanding of the spa-
tiotemporal properties of sustained nonequilibrium systems such as fluids, plasmas,
reacting and diffusing chemicals, crystals solidifying from a melt, heart muscle, and
networks of excitable neurons in brains.
This introductory textbook for graduate students in biology, chemistry, engineer-
ing, mathematics, and physics provides a systematic account of the basic science
common to these diverse areas. This book provides a careful pedagogical motivation
of key concepts, discusses why diverse nonequilibrium systems often show similar
patterns and dynamics, and gives a balanced discussion of the role of experiments,
simulation, and analytics. It contains numerous illustrative worked examples, and
over 150 exercises.
This book will also interest scientists who want to learn about the experi-
ments, simulations, and theory that explain how complex patterns form in sustained
nonequilibrium systems.
M I C HAE L C R O S S
California Institute of Technology
HENRY GREENSIDE
Duke University
Published in the United States of America by Cambridge University Press, New York
www.cambridge.org
Information on this title: www.cambridge.org/9780521770507
A catalogue record for this publication is available from the British Library
vii
5 Models 172
5.1 Swift–Hohenberg model 174
5.1.1 Heuristic derivation 175
5.1.2 Properties 178
5.1.3 Numerical simulations 182
5.1.4 Comparison with experimental systems 184
5.2 Generalized Swift–Hohenberg models 186
5.2.1 Non-symmetric model 186
5.2.2 Non-potential models 187
5.2.3 Models with mean flow 187
5.2.4 Model for rotating convection 189
5.2.5 Model for quasicrystalline patterns 191
5.3 Order-parameter equations 191
5.4 Complex Ginzburg–Landau equation 195
5.5 Kuramoto–Sivashinsky equation 196
5.6 Reaction–diffusion models 198
5.7 Models that are discrete in space, time, or value 200
5.8 Conclusions 200
5.9 Further reading 201
6 One-dimensional amplitude equation 206
6.1 Origin and meaning of the amplitude 209
6.2 Derivation of the amplitude equation 212
6.2.1 Phenomenological derivation 212
6.2.2 Deduction of the amplitude-equation parameters 215
6.2.3 Method of multiple scales 216
6.2.4 Boundary conditions for the amplitude equation 217
6.3 Properties of the amplitude equation 219
6.3.1 Universality and scales 219
6.3.2 Potential dynamics 222
6.4 Applications of the amplitude equation 224
6.4.1 Lateral boundaries 224
6.4.2 Eckhaus instability 228
6.4.3 Phase dynamics 232
6.5 Limitations of the amplitude equation formalism 235
6.6 Conclusions 236
6.7 Further reading 237
7 Amplitude equations for two-dimensional patterns 242
7.1 Stripes in rotationally invariant systems 244
7.1.1 Amplitude equation 244
7.1.2 Boundary conditions 246
xiii
electrical control techniques to prevent epilepsy or a heart attack. In these and other
cases such as forecasting the weather or predicting earthquakes, improvements in
the design, control, and prediction of nonequilibrium systems are often limited by
our incomplete understanding of sustained nonequilibrium dynamics.
Finally, a third reason for learning the material in this book is to develop specific
conceptual, mathematical, and numerical skills for understanding complex phenom-
ena. Many nonequilibrium systems involve continuous media whose quantitative
description is given in terms of nonlinear partial differential equations. The solu-
tions of such equations can be difficult to understand (e.g. because they may evolve
nonperiodically in time and be simultaneously disordered in space), and questions
such as “Is the output from this computer simulation correct?’’, “Is this simulation
producing the same results as my experimental data?’’ or “Is experimental noise
relevant here?’’ may not be easily answered. As an example, one broadly useful
mathematical technique that we discuss and use several times throughout the book
is multiscale perturbation theory, which leads to so-called “amplitude equations’’
that provide a quantitatively useful reduction of complex dynamics. We also discuss
the role of numerical simulation, which has some advantages and disadvantages
compared to analytical theory and experimental investigation.
To help the reader master the various conceptual, mathematical, and numerical
skills, the book has numerous worked examples that we call etudes. By analogy
to a musical etude, which is a composition that helps a music student master a
particular technique while also learning a piece of artistic value, our etudes are
one- to two-page long worked examples that illustrate a particular idea and that
also try to provide a non-trivial application of the idea.
Although this book is intended for an interdisciplinary audience, it is really a
physics book in the following sense. Many of the nonequilibrium systems in the
Universe, for example a germ or a star, are simply too complex to analyze directly
and so are ill-suited for discovering fundamental properties upon which a general
quantitative understanding can be developed. In much of this book, we follow
a physics tradition of trying to identify and study simple idealized experimental
systems that also have some of the interesting properties observed in more complex
systems.
Thus instead of studying the exceedingly complex dynamics of the Earth’s
weather, which would require in turn understanding the effects of clouds, the solar
wind, the coupling to oceans and ice caps, the topography of mountains and forest,
and the effects of human industry, we instead focus our experimental and theoretical
attention on enormously simplified laboratory systems. One example is Rayleigh–
Bénard convection, which is a fluid experiment consisting of a thin horizontal layer
of a pure fluid that is driven out of equilibrium by a vertical temperature difference
that is constant in time and uniform in space. Another is a mixture of reacting and
In this opening chapter, we give an informal and qualitative overview – a pep talk –
to help you appreciate why sustained nonequilibrium systems are so interesting and
worthy of study.
We begin in Section 1.1 by discussing the big picture of how the Universe is filled
with nonequilibrium systems of many different kinds, a consequence of the fact
that the Universe had a beginning and has not yet stopped evolving. A profound and
important question is then to understand how the observed richness of structure in
the Universe arises from the property of not being in thermodynamic equilibrium.
In Section 1.2, a particularly well studied nonequilibrium system, Rayleigh–Bénard
convection, is introduced to establish some vocabulary and insight regarding what
is a nonequilibrium system. Next, in Section 1.3, we extend our discussion to
representative examples of nonequilibrium patterns in nature and in the labora-
tory, to illustrate the great diversity of such patterns and to provide some concrete
examples to think about. These examples serve to motivate some of the central ques-
tions that are discussed throughout the book, e.g. spatially dependent instabilities,
wave number selection, pattern formation, and spatiotemporal chaos. The humble
desktop-sized experiments discussed in this section, together with theory and sim-
ulations relating to them, can also be regarded as the real current battleground for
understanding nonequilibrium systems since there is a chance to compare theory
with experiment quantitatively.
Next, Section 1.4 discusses some of the ways that pattern-forming nonequilib-
rium systems differ from the low-dimensional dynamical systems that you may
have seen in an introductory nonlinear dynamics course. Some guidelines are also
given to determine qualitatively when low-dimensional nonlinear dynamics may
not suffice to analyze a particular nonequilibrium system. In Section 1.5, a strategy
is given and explained for exploring nonequilibrium systems. We explain why fluid
dynamics experiments have some advantages over other possible experimental sys-
tems and why certain fluid experiments such as Rayleigh–Bénard convection are
especially attractive. Finally, Section 1.6 mentions some of the topics that we will
not address in this book for lack of time or expertise.
While this book will explain some of what is known about these questions, espe-
cially at the laboratory level which allows controlled reproducible experiments, we
can say at a hand-waving level why the Universe is interesting rather than boring:
the Universe was born in a cosmological Big Bang and is still young when measured
in units of the lifetime of a star. Thus the universe has not yet lasted long enough to
come to thermodynamic equilibrium: the Universe as a whole is a nonequilibrium
system. Because stars are young and have not yet reached thermodynamic equilib-
rium, the nuclear fuel in their core has not yet been consumed. The flux of energy
from this core through the surface of the star and out into space drives the complex
dynamics of the star’s plasma and magnetic field. Similarly, because the Earth is
still geologically young, its interior has not yet cooled down and the flux of heat
from its hot core out through its surface, together with heat received from the Sun,
drives the dynamics of the atmosphere, ocean, and mantle. And it is this same flux
of energy from the Earth and Sun that sustains Earth’s intricate biosphere.
This hand-waving explanation of the origin of nonequilibrium structure is
unsatisfactory since it does not lead to the quantitative testing of predictions by
experiment. To make progress, scientists have found it useful to turn to desktop
experimental systems that can be readily manipulated and studied, and that are also
easier to analyze mathematically and to simulate with a computer. The experiments
and theory described in this book summarize some of the systematic experimental
and theoretical efforts of the last thirty years to understand how to predict and to
analyze such desktop nonequilibrium phenomena. However, you should appreciate
that much interesting research remains to be carried out if our desktop insights are
to be related to the more complex systems found in the world around us. We hope
that this book will encourage you to become an active participant in this challenging
endeavor.
Cold T2
Warm T1
the upper plate is cool. As an example to visualize (but a bit impractical for actual
experimentation as you will discover in Exercise 1.5), consider a square room whose
lateral width L is larger than its height d , and in which all furniture, doors, win-
dows, and fixtures have been removed so that there is only a smooth flat horizontal
floor, a smooth flat horizontal ceiling, and smooth flat vertical walls (see Fig. 1.1).
The floor and ceiling are then coated with a layer of copper, and just beneath the
floor and just above the ceiling some water-carrying pipes and electronic circuits
connected to water heaters are arranged so that the floor is maintained at a constant
temperature T1 and the ceiling is maintained at a constant temperature T2 .1 Because
copper conducts heat so well, any temperature variations within the floor or within
the ceiling quickly become negligible so that the floor and ceiling can be considered
as time-independent constant-temperature surfaces. The supporting sidewalls are
made of some material that conducts heat poorly such as wood or Plexiglas.
A typical nonequilibrium experiment for the room in Fig. 1.1 would then be
simply to fix the temperature difference T = T1 − T2 at some value and then
to observe what happens to the air. “Observe what happens’’ can mean several
1 Uniformly warming the floor and cooling the ceiling is not the usual way that a room is heated. Instead, a
convector – a localized heat source with a large surface area – is placed somewhere in the room, and heat is lost
through the windows instead of through the ceiling. (What we call a convector everyone else calls a radiator
but this is poorly named since the air is heated mainly by convection, not by radiation.) But this nonuniform
geometry is more complicated, and so less well suited, than our idealized room for experiment and analysis.
things depending on the questions of interest. By introducing some smoke into the
room, the pattern of air currents could be visualized. A more quantitative obser-
vation might involve recording as a function of time t some local quantity such
as the temperature T (x0 , t) or the x-component of the air’s velocity vx (x0 , t) at a
particular fixed position x0 = (x0 , y0 , z0 ) inside the room. Alternatively, an experi-
mentalist might choose to record some global quantity such as the total heat H (t)
transported from the floor to the ceiling, a quantity of possible interest to mechan-
ical engineers and architects. These measurements of some quantity at successive
moments of time constitute a time series that can be stored, plotted, and analyzed.
A more ambitious and difficult observation might consist of measuring multivariate
time series, e.g. measuring the temperature field T (x, t) and the components of the
velocity field v(x, t) simultaneously at many different spatial points, at successive
instants of time. These data could then be made into movies or analyzed statisti-
cally. All of these observations are carried out for a particular fixed choice of the
temperature difference T and over some long time interval (long enough that
any transient behavior will decay sufficiently). Other experiments might involve
repeating the same measurements but for several successive values of T , with each
value again held constant during a given experiment. In this way, the spatiotemporal
dynamical properties of the air in the room can be mapped out as a function of the
parameter T , and various dynamical states and transitions between them can be
identified.
The temperature difference T is a particularly important parameter in a convec-
tion experiment because it determines whether or not the fluid is in thermodynamic
equilibrium. (It is precisely the fact that the nonequilibrium properties of the entire
room can be described by a single parameter T that constitutes the idealization
of this experiment, and that motivated the extra experimental work of coating the
floor and ceiling with copper.) If T = 0 so that the ceiling and floor have the
same common temperature T = T1 = T2 , then after some transient time, the air will
be in thermodynamic equilibrium with zero velocity and the same uniform temper-
ature T throughout. There is typically a transient time associated with approaching
thermodynamic equilibrium because the air itself is rarely in such equilibrium with-
out taking special precautions. For example, there might be a small breeze in the
air when the door to the experimental room is closed or some part of the air may
be a bit warmer than some other part because someone walked through the room.
But as long as the room is sealed and the floor and ceiling have the same tempera-
ture, all macroscopic motion in the air will die out and the air will attain the same
temperature everywhere.
As soon as the temperature difference T becomes nonzero (with either sign),
the air can no longer be in thermodynamic equilibrium since the temperature is
spatially nonuniform. One says that the air is driven out of equilibrium by the
2 More precisely, the fluid velocity at a wall is zero in a frame of reference moving with the surface. Exercise 1.9
suggests a simple experiment using an electric fan to explore this point.
T1
Fig. 1.2 Illustration of the driving and dissipative forces acting on small parcels
of air near the floor and ceiling of the experimental room in Fig. 1.1 whose floor
is warmer than its ceiling. The parcels are assumed to be small enough that their
temperatures are approximately constant over their interiors. The acceleration of
the parcels by buoyancy forces is opposed by a friction arising from the fluid
viscosity and also by the diffusion of heat between warmer and cooler regions of
the fluid. Only when the temperature difference T = T1 − T2 exceeds a finite
critical value Tc > 0 can the buoyancy forces overcome the dissipation and
convection currents form.
fluid downwards, in accord with the truism that “hot air rises’’ and “cold air falls.’’
These buoyancy forces constitute the physical mechanism by which the tempera-
ture difference T “drives’’ the air out of equilibrium. As a warm parcel moves
upward, it has to push its way through the surrounding fluid and this motion is
opposed by the dissipative friction force associated with fluid viscosity. Also, as
the parcel rises, it loses heat by thermal conduction to the now cooler surrounding
air, becomes more dense, and the buoyancy force is diminished. Similar dissipative
effects act on a cool descending parcel.
From this microscopic picture, we can understand the experimental fact that
making the temperature difference T positive is a necessary but not sufficient
condition for the air to start moving since the buoyancy forces may not be strong
enough to overcome the dissipative effects of viscosity and conduction. Indeed,
experiment and theory show that only when the temperature difference exceeds a
threshold, a critical value we denote as Tc , will the buoyancy forces be sufficiently
large that the air will spontaneously start to move and a persistent spatiotemporal
structure will appear in the form of convection currents. If the room’s width L is
large compared to its depth d so that the influence of the walls on the bulk fluid
can be ignored, a precise criterion for the onset of convection can be stated in
the form
R > Rc . (1.1)
αgd 3 T
R= , (1.2)
νκ
and the critical value of R has the approximate value
Rc ≈ 1708. (1.3)
The parameters in Eq. (1.2) have the following meaning: g is the gravitational accel-
eration, about 9.8 m/s2 over much of the Earth’s surface; α = −(1/ρ)(∂ρ/∂T )|p
is the fluid’s coefficient of thermal expansion at constant pressure, and measures
the relative change in density ρ as the temperature is varied; d is the uniform depth
of the fluid; T is the uniform temperature difference across the fluid layer; ν is
the fluid’s kinematic viscosity; and κ is the fluid’s thermal diffusivity. Approximate
values of the parameters α, ν, and κ for air, water, and mercury at room temperature
(T = 293 K) and at atmospheric pressure are given in Table 1.1.
The combination of physical parameters in Eq. (1.2) is dimensionless and so has
the same value no matter what physical units are used in any given experiment,
e.g. System Internationale (SI), Centimeter-Gram-Seconds (CGS), or British. This
combination is denoted by the symbol “R’’ and is called the Rayleigh number in
honor of the physicist and applied mathematician Lord Rayleigh who, in 1916, was
the first to identify its significance for determining the onset of convection. The pure
number Rc is called the critical Rayleigh number Rc since it denotes the threshold
that R must exceed for convection to commence. The value Rc can be calculated
directly from the equations that govern the time evolution of a convecting fluid
(the Boussinesq equations) as the criterion when the motionless conducting state of
the fluid first becomes linearly unstable. The general method of this linear stability
analysis is described in Chapter 2.
Despite its dependence on six parameters, you should think of the Rayleigh
number R as simply being proportional to the temperature difference T . The
reason is that all the parameters in Eq. (1.2) except T are approximately constant
in a typical series of convection experiments. Thus the parameters α, ν, and κ in
Eq. (1.1) depend weakly on temperature and are effectively fixed once a particular
fluid is chosen. The acceleration g is fixed once a particular geographical location
is selected for the experiment and the depth of the fluid d is typically fixed once
the convection cell has been designed and is difficult to vary as an experimental
parameter. Only the temperature difference T is easily changed substantially and
so this naturally becomes the experimental control parameter.
You should also note that the numerator αgd 3 T in Eq. (1.2) is related to quanti-
ties that determine the buoyancy force, while the denominator νκ involves quantities
related to the two dissipative mechanisms so Eq. (1.1) indeed states that instability
will not occur until the driving is sufficiently strong compared to the dissipation.
Most nonequilibrium systems have one or more such dimensionless parameters
associated with them and these parameters are key quantities to identify and to
measure when studying a nonequilibrium system.
What kind of dynamics can we expect for the air if the Rayleigh number R is
held constant at some value larger than the critical value Rc ? From Fig. 1.2, we
expect the warm fluid near the floor to rise and the cool fluid near the ceiling to
descend but the entire layer of ascending fluid near the floor cannot pass through the
entire layer of descending fluid near the ceiling because the fluid is approximately
incompressible. What is observed experimentally is pattern formation: the fluid
spontaneously achieves a compromise such that some regions of fluid rise and
neighboring regions descend, leading to the formation of a cellular convection
“pattern’’ in the temperature, velocity, and pressure fields. The distance between
adjacent rising and falling regions turns out to be about the depth of the air. Once
the air begins to convect, the dynamics becomes too complicated to understand by
casual arguments applied to small parcels of air and we need to turn to experiments
to observe what happens and to a deeper mathematical analysis to understand the
experimental results (see Figs. 1.14 and 1.15 below in Section 1.3.2). However, one
last observation can be made. The motion of the fluid parcels inside the experimental
system transport heat and thereby modify the temperature gradient that is felt in a
particular location inside the system. Thus the motion of the medium changes the
balance of driving and dissipation in different parts of the medium, and this is the
reason why the dynamics is nonlinear and often difficult to understand.
The general points we learn from the above discussion about Rayleigh–Bénard
convection are the following. There are mechanisms that can drive a system out of
thermodynamic equilibrium, such as a flux of energy, momentum or matter through
the system. This driving is opposed by one or more dissipative mechanisms such
1.D. The copyright laws of the place where you are located also
govern what you can do with this work. Copyright laws in most
countries are in a constant state of change. If you are outside
the United States, check the laws of your country in addition to
the terms of this agreement before downloading, copying,
displaying, performing, distributing or creating derivative works
based on this work or any other Project Gutenberg™ work. The
Foundation makes no representations concerning the copyright
status of any work in any country other than the United States.
1.E.6. You may convert to and distribute this work in any binary,
compressed, marked up, nonproprietary or proprietary form,
including any word processing or hypertext form. However, if
you provide access to or distribute copies of a Project
Gutenberg™ work in a format other than “Plain Vanilla ASCII” or
other format used in the official version posted on the official
Project Gutenberg™ website (www.gutenberg.org), you must,
at no additional cost, fee or expense to the user, provide a copy,
a means of exporting a copy, or a means of obtaining a copy
upon request, of the work in its original “Plain Vanilla ASCII” or
other form. Any alternate format must include the full Project
Gutenberg™ License as specified in paragraph 1.E.1.
• You pay a royalty fee of 20% of the gross profits you derive
from the use of Project Gutenberg™ works calculated using the
method you already use to calculate your applicable taxes. The
fee is owed to the owner of the Project Gutenberg™ trademark,
but he has agreed to donate royalties under this paragraph to
the Project Gutenberg Literary Archive Foundation. Royalty
payments must be paid within 60 days following each date on
which you prepare (or are legally required to prepare) your
periodic tax returns. Royalty payments should be clearly marked
as such and sent to the Project Gutenberg Literary Archive
Foundation at the address specified in Section 4, “Information
about donations to the Project Gutenberg Literary Archive
Foundation.”
• You comply with all other terms of this agreement for free
distribution of Project Gutenberg™ works.
1.F.
Most people start at our website which has the main PG search
facility: www.gutenberg.org.
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
ebookname.com